Mechanical stress analysis of microfluidic environments designed for isolated biological cell investigations

Computational simulation of applying mechanical vibration to mesenchymal stem cell for mechanical modulation toward bone tissue engineering

Evaluation of cell response to mechanical stimuli at in vitro conditions is known as one of the important issues for modulating cell behavior. Mechanical stimuli, including mechanical vibration and oscillatory fluid flow, act as important biophysical signals for the mechanical modulation of stem cells. In the present study, mesenchymal stem cell (MSC) consists of cytoplasm, nucleus, actin, and microtubule. Also, integrin and primary cilium were considered as mechanoreceptors. In this study, the combined effect of vibration and oscillatory fluid flow on the cell and its components were investigated using numerical modeling. The results of the FEM and FSI model showed that the cell response (stress and strain values) at the frequency of 30 H z mechanical vibration has the highest value. The achieved results on shear stress caused by the fluid flow on the cell showed that the cell experiences shear stress in the range of 0 . 1 - 10 Pa . Mechanoreceptors that bind separately to the cell surface, can be highly stimulated by hydrodynamic pressure and, therefore, can play a role in the mechanical modulation of MSCs at in vitro conditions. The results of this research can be effective in future studies to optimize the conditions of mechanical stimuli applied to the cell culture medium and to determine the mechanisms involved in mechanotransduction.

Microfluidic Surgery in Single Cells and Multicellular Systems

Microscale surgery on single cells and small organisms has enabled major advances in fundamental biology and in engineering biological systems. Examples of applications range from wound healing and regeneration studies to the generation of hybridoma to produce monoclonal antibodies. Even today, these surgical operations are often performed manually, but they are labor intensive and lack reproducibility. Microfluidics has emerged as a powerful technology to control and manipulate cells and multicellular systems at the micro- and nanoscale with high precision. Here, we review the physical and chemical mechanisms of microscale surgery and the corresponding design principles, applications, and implementations in microfluidic systems. We consider four types of surgical operations: (1) sectioning, which splits a biological entity into multiple parts, (2) ablation, which destroys part of an entity, (3) biopsy, which extracts materials from within a living cell, and (4) fusion, which joins multiple entities into one. For each type of surgery, we summarize the motivating applications and the microfluidic devices developed. Throughout this review, we highlight existing challenges and opportunities. We hope that this review will inspire scientists and engineers to continue to explore and improve microfluidic surgical methods.

Oxygen mass transfer and shear stress effects on Pseudomonas putida BCRC 14365 growth to improve bioreactor design and performance

In this work, the experimental evidence is presented for two basic issues including oxygen mass transfer and shear analysis on the microorganism containing medium on the most prominent sections of the bioreactor. Computational fluid dynamics (CFD) methodology reproduces shear rate values for specific impeller designs using the commercial code (Fluent 6.2). CFD calculates volumetric mass transfer coefficient based on the Higbie's penetration theory. Four types of impeller are used. The spherical probe is used to measure flow hydrodynamic parameters to obtain shear rate by electro-diffusion (ED) method. The obtained results are validated experimentally and it is shown that a fully axial pattern impeller represents more enhanced results than partially axial and radial. In this regard, experimental results for volumetric oxygen mass transfer coefficient (k _l a) confirm CFD predictions by acceptable deviations of 2.65, 8.90, and 9.20 for 0.15, 0.2, and 0.3 VVM, respectively. These results collaboratively indicate that LIGHTNIN-C 200 type operates more efficiently by reflecting the flow to the bottom corner stagnation areas with the minimum tolerable shear and the most velocity distribution uniformity. Furthermore, the values of k _l a improve by aeration rate. Conversely, increasing the rotational speed of impeller creates difficulties for cell growth due to the generated harsh shear condition. CFD provide a better understanding of how operational and geometrical variables may be manipulated to achieve a moderate shear rate and acceptable level of mass transfer.

Engineering of a synthetic quadrastable gene network to approach Waddington landscape and cell fate determination

The process of cell fate determination has been depicted intuitively as cells travelling and resting on a rugged landscape, which has been probed by various theoretical studies. However, few studies have experimentally demonstrated how underlying gene regulatory networks shape the landscape and hence orchestrate cellular decision-making in the presence of both signal and noise. Here we tested different topologies and verified a synthetic gene circuit with mutual inhibition and auto-activations to be quadrastable, which enables direct study of quadruple cell fate determination on an engineered landscape. We show that cells indeed gravitate towards local minima and signal inductions dictate cell fates through modulating the shape of the multistable landscape. Experiments, guided by model predictions, reveal that sequential inductions generate distinct cell fates by changing landscape in sequence and hence navigating cells to different final states. This work provides a synthetic biology framework to approach cell fate determination and suggests a landscape-based explanation of fixed induction sequences for targeted differentiation.

DOI:http://dx.doi.org/10.7554/eLife.23702.001

Cells in animals use a process called differentiation to specialize into specific cell types such as skin cells and liver cells. Proteins called transcription factors drive particular steps in differentiation by controlling the activity of specific genes. Many transcription factors interact with each other to form complex networks that regulate gene activity to determine the fate of a cell and control the whole differentiation process. Some individual gene networks can program cells to become any one of several different cell fates, a feature known as multistability.

In the 1950s, a scientist called Conrad Waddington proposed the concept of an “epigenetic landscape” to describe how the fate of a cell is decided as an animal develops. The cell, depicted as a ball, rolls down a rugged landscape and has the option of taking several different routes. Each route will eventually lead to a distinct cell fate. As the ball moves down the hill, the choice of routes and final destinations becomes more limited. Theoretical approaches have been used to understand how gene regulatory networks shape the epigenetic landscape of an animal. However, few studies have experimentally tested the findings of the theoretical approaches and it is not clear how environmental inputs help to determine which path a cell will take.

Although bacteria cells do not generally specialize into particular cell types, bacteria cells can use multistability in transcription factor networks to switch between different behaviors or “states” in response to cues from the environment. Wu et al. used a bacterium called E. coli as a model to investigate how a gene network called MINPA from mammals, which is involved in differentiation and is believed to show multistability, can guide cells to adopt different states. The work combined experimental and mathematical approaches to design, construct and test an artificial version of the MINPA gene network in E. coli.

The experiments showed that MINPA could direct the cells to adopt four different stable states in which the cells produced fluorescent proteins of different colors. With the help of mathematical modeling, Wu et al. charted how the landscape of cell states changed when external chemical cues were applied. Exposing the cells to several cues in particular orders guided the cells to different final states.

The findings of Wu et al. shed new light on how the fate of a cell is determined and provide a theoretical framework for understanding the complex networks that control cell differentiation. This could help develop new ways of directing cell fate that could ultimately be used to generate cells to treat human diseases.

DOI:http://dx.doi.org/10.7554/eLife.23702.002

Guiding pancreatic beta cells to target electrodes in a whole-cell biosensor for diabetes

We are developing a cell-based bioelectronic glucose sensor that exploits the multi-parametric sensing ability of pancreatic islet cells for the treatment of diabetes. These cells sense changes in the concentration of glucose and physiological hormones and immediately react by generating electrical signals. In our sensor, signals from multiple cells are recorded as field potentials by a micro-electrode array (MEA). Thus, cell response to various factors can be assessed rapidly and with high throughput. However, signal quality and consequently overall sensor performance rely critically on close cell-electrode proximity. Therefore, we present here a non-invasive method of further exploiting the electrical properties of these cells to guide them towards multiple micro-electrodes via electrophoresis. Parameters were optimized by measuring the cell's zeta potential and modeling the electric field distribution. Clonal and primary mouse or human β-cells migrated directly to target electrodes during the application of a 1 V potential between MEA electrodes for 3 minutes. The morphology, insulin secretion, and electrophysiological characteristics were not altered compared to controls. Thus, cell manipulation on standard MEAs was achieved without introducing any external components and while maintaining the performance of the biosensor. Since the analysis of the cells' electrical activity was performed in real time via on-chip recording and processing, this work demonstrates that our biosensor is operational from the first step of electrically guiding cells to the final step of automatic recognition. Our favorable results with pancreatic islets, which are highly sensitive and fragile cells, are encouraging for the extension of this technique to other cell types and microarray devices.

Fluid-Structure Interactions Analysis of Shear-Induced Modulation of a Mesenchymal Stem Cell: An Image-Based Study

Although effects of biochemical modulation of stem cells have been widely investigated, only recent advances have been made in the identification of mechanical conditioning on cell signaling pathways. Experimental investigations quantifying the micromechanical environment of mesenchymal stem cells (MSCs) are challenging while computational approaches can predict their behavior due to in vitro stimulations. This study introduces a 3D cell-specific finite element model simulating large deformations of MSCs. Here emphasizing cell mechanical modulation which represents the most challenging multiphysics phenomena in sub-cellular level, we focused on an approach attempting to elicit unique responses of a cell under fluid flow. Fluorescent staining of MSCs was performed in order to visualize the MSC morphology and develop a geometrically accurate model of it based on a confocal 3D image. We developed a 3D model of a cell fixed in a microchannel under fluid flow and then solved the numerical model by fluid-structure interactions method. By imposing flow characteristics representative of vigorous in vitro conditions, the model predicts that the employed external flow induces significant localized effective stress in the nucleo-cytoplasmic interface and average cell deformation of about 40%. Moreover, it can be concluded that a lower strain level is made in the cell by the oscillatory flow as compared with steady flow, while same ranges of effective stress are recorded inside the cell in both conditions. The deeper understanding provided by this study is beneficial for better design of single cell in vitro studies.

Cytoskeletal strains in modeled optohydrodynamically stressed healthy and diseased biological cells

Controlled external chemomechanical stimuli have been shown to influence cellular and tissue regeneration/degeneration, especially with regards to distinct disease sequelae or health maintenance. Recently, a unique three-dimensional stress state was mathematically derived to describe the experimental stresses applied to isolated living cells suspended in an optohydrodynamic trap (optical tweezers combined with microfluidics). These formulae were previously developed in two and three dimensions from the fundamental equations describing creeping flows past a suspended sphere. The objective of the current study is to determine the full-field cellular strain response due to the applied three-dimensional stress environment through a multiphysics computational simulation. In this investigation, the multiscale cytoskeletal structures are modeled as homogeneous, isotropic, and linearly elastic. The resulting computational biophysics can be directly compared with experimental strain measurements, other modeling interpretations of cellular mechanics including the liquid drop theory, and biokinetic models of biomolecule dynamics. The described multiphysics computational framework will facilitate more realistic cytoskeletal model interpretations, whose intracellular structures can be distinctly defined, including the cellular membrane substructures, nucleus, and organelles.

Geometrical effects in microfluidic-based microarrays for rapid, efficient single-cell capture of mammalian stem cells and plant cells

In this paper, a detailed numerical and experimental investigation into the optimisation of hydrodynamic micro-trapping arrays for high-throughput capture of single polystyrene (PS) microparticles and three different types of live cells at trapping times of 30 min or less is described. Four different trap geometries (triangular, square, conical, and elliptical) were investigated within three different device generations, in which device architecture, channel geometry, inter-trap spacing, trap size, and trap density were varied. Numerical simulation confirmed that (1) the calculated device dimensions permitted partitioned flow between the main channel and the trap channel, and further, preferential flow through the trap channel in the absence of any obstruction; (2) different trap shapes, all having the same dimensional parameters in terms of depth, trapping channel lengths and widths, main channel lengths and widths, produce contrasting streamline plots and that the interaction of the fluid with the different geometries can produce areas of stagnated flow or distorted field lines; and (3) that once trapped, any motion of the trapped particle or cell or a shift in its configuration within the trap can result in significant increases in pressures on the cell surface and variations in the shear stress distribution across the cell's surface. Numerical outcomes were then validated experimentally in terms of the impact of these variations in device design elements on the percent occupancy of the trapping array (with one or more particles or cells) within these targeted short timeframes. Limitations on obtaining high trap occupancies in the devices were shown to be primarily a result of particle aggregation, channel clogging and the trap aperture size. These limitations could be overcome somewhat by optimisation of these device design elements and other operational variables, such as the average carrier fluid velocity. For example, for the 20 μm polystyrene microparticles, the number of filled traps increased from 32% to 42% during 5-10 min experiments in devices with smaller apertures. Similarly, a 40%-60% reduction in trapping channel size resulted in an increase in the amount of filled traps, from 0% to almost 90% in 10 min, for the human bone marrow derived mesenchymal stem cells, and 15%-85% in 15 min for the human embryonic stem cells. Last, a reduction of the average carrier fluid velocity by 50% resulted in an increase from 80% to 92% occupancy of single algae cells in traps. Interestingly, changes in the physical properties of the species being trapped also had a substantial impact, as regardless of the trap shape, higher percent occupancies were observed with cells compared to single PS microparticles in the same device, even though they are of approximately the same size. This investigation showed that in microfluidic single cell capture arrays, the trap shape that maximizes cell viability is not necessarily the most efficient for high-speed single cell capture. However, high-speed trapping configurations for delicate mammalian cells are possible but must be optimised for each cell type and designed principally in accordance with the trap size to cell size ratio.

Cell-to-cell variability in deformations across compressed myoblasts

Many biological consequences of external mechanical loads applied to cells depend on localized cell deformations rather than on average whole-cell-body deformations. Such localized intracellular deformations are likely to depend, in turn, on the individual geometrical features of each cell, e.g., the local surface curvatures or the size of the nucleus, which always vary from one cell to another, even within the same culture. Our goal here was to characterize cell-to-cell variabilities in magnitudes and distribution patterns of localized tensile strains that develop in the plasma membrane (PM) and nuclear surface area (NSA) of compressed myoblasts, in order to identify resemblance or differences in mechanical performances across the cells. For this purpose, we utilized our previously developed confocal microscopy-based three-dimensional cell-specific finite element modeling methodology. Five different C2C12 undifferentiated cells belonging to the same culture were scanned confocally and modeled, and were then subjected to compression in the simulation setting. We calculated the average and peak tensile strains in the PM and NSA, the percentage of PM area subjected to tensile strains above certain thresholds and the coefficient of variation (COV) in average and peak strains. We found considerable COV values in tensile strains developing at the PM and NSA (up to ~35%) but small external compressive deformations induced greater variabilities in intracellular strains across cells compared to large deformations. Interestingly, the external deformations needed to cause localized PM or NSA strains exceeding each threshold were very close across the different cells. Better understanding of variabilities in mechanical performances of cells-either of the same type or of different types-is important for interpreting experimental data in any experiments involving delivery of mechanical loads to cells.

Finite-element modeling of viscoelastic cells during high-frequency cyclic strain

Mechanotransduction refers to the mechanisms by which cells sense and respond to local loads and forces. The process of mechanotransduction plays an important role both in maintaining tissue viability and in remodeling to repair damage; moreover, it may be involved in the initiation and progression of diseases such as osteoarthritis and osteoporosis. An understanding of the mechanisms by which cells respond to surrounding tissue matrices or artificial biomaterials is crucial in regenerative medicine and in influencing cellular differentiation. Recent studies have shown that some cells may be most sensitive to low-amplitude, high-frequency (i.e., 1-100 Hz) mechanical stimulation. Advances in finite-element modeling have made it possible to simulate high-frequency mechanical loading of cells. We have developed a viscoelastic finite-element model of an osteoblastic cell (including cytoskeletal actin stress fibers), attached to an elastomeric membrane undergoing cyclic isotropic radial strain with a peak value of 1,000 µstrain. The results indicate that cells experience significant stress and strain amplification when undergoing high-frequency strain, with peak values of cytoplasmic strain five times higher at 45 Hz than at 1 Hz, and peak Von Mises stress in the nucleus increased by a factor of two. Focal stress and strain amplification in cells undergoing high-frequency mechanical stimulation may play an important role in mechanotransduction.

Volumetric stress-strain analysis of optohydrodynamically suspended biological cells

Ongoing investigations are exploring the biomechanical properties of isolated and suspended biological cells in pursuit of understanding single-cell mechanobiology. An optical tweezer with minimal applied laser power has positioned biologic cells at the geometric center of a microfluidic cross-junction, creating a novel optohydrodynamic trap. The resulting fluid flow environment facilitates unique multiaxial loading of single cells with site-specific normal and shear stresses resulting in a physical albeit extensional state. A recent two-dimensional analysis has explored the cytoskeletal strain response due to these fluid-induced stresses [Wilson and Kohles, 2010, "Two-Dimensional Modeling of Nanomechanical Stresses-Strains in Healthy and Diseased Single-Cells During Microfluidic Manipulation," J Nanotechnol Eng Med, 1(2), p. 021005]. Results described a microfluidic environment having controlled nanometer and piconewton resolution. In this present study, computational fluid dynamics combined with multiphysics modeling has further characterized the applied fluid stress environment and the solid cellular strain response in three dimensions to accompany experimental cell stimulation. A volumetric stress-strain analysis was applied to representative living cell biomechanical data. The presented normal and shear stress surface maps will guide future microfluidic experiments as well as provide a framework for characterizing cytoskeletal structure influencing the stress to strain response.

Manipulation of Suspended Single Cells by Microfluidics and Optical Tweezers

Chondrocytes and osteoblasts experience multiple stresses in vivo. The optimum mechanical conditions for cell health are not fully understood. This paper describes the optical and microfluidic mechanical manipulation of single suspended cells enabled by the μPIVOT, an integrated micron resolution particle image velocimeter (μPIV) and dual optical tweezers instrument (OT). In this study, we examine the viability and trap stiffness of cartilage cells, identify the maximum fluid-induced stresses possible in uniform and extensional flows, and compare the deformation characteristics of bone and muscle cells. These results indicate cell photodamage of chondrocytes is negligible for at least 20 min for laser powers below 30 mW, a dead cell presents less resistance to internal organelle rearrangement and deforms globally more than a viable cell, the maximum fluid-induced shear stresses are limited to ~15 mPa for uniform flows but may exceed 1 Pa for extensional flows, and osteoblasts show no deformation for shear stresses up to 250 mPa while myoblasts are more easily deformed and exhibit a modulated response to increasing stress. This suggests that global and/or local stresses can be applied to single cells without physical contact. Coupled with microfluidic sensors, these manipulations may provide unique methods to explore single cell biomechanics.

A Distinct Catabolic to Anabolic Threshold Due to Single-Cell Static Nanomechanical Stimulation in a Cartilage Biokinetics Model

Understanding physicochemical interactions during biokinetic regulation will be critical for the creation of relevant nanotechnology supporting cellular and molecular engineering. The impact of nanoscale influences in medicine and biology can be explored in detail through mathematical models as an in silico testbed. In a recent single-cell biomechanical analysis, the cytoskeletal strain response due to fluid-induced stresses was characterized (Wilson, Z. D., and Kohles, S. S., 2010, "Two-Dimensional Modeling of Nanomechanical Strains in Healthy and Diseased Single-Cells During Microfluidic Stress Applications," J. Nanotech. Eng. Med., 1(2), p. 021005). Results described a microfluidic environment having controlled nanometer and piconewton resolution for explorations of multiscale mechanobiology. In the present study, we constructed a mathematical model exploring the nanoscale biomolecular response to that controlled microenvironment. We introduce mechanical stimuli and scaling factor terms as specific input values for regulating a cartilage molecule synthesis. Iterative model results for this initial multiscale static load application have identified a transition threshold load level from which the mechanical input causes a shift from a catabolic state to an anabolic state. Modeled molecule homeostatic levels appear to be dependent upon the mechanical stimulus as reflected experimentally. This work provides a specific mathematical framework from which to explore biokinetic regulation. Further incorporation of nanomechanical stresses and strains into biokinetic models will ultimately lead to refined mechanotransduction relationships at the cellular and molecular levels.

Two-Dimensional Modeling of Nanomechanical Strains in Healthy and Diseased Single-Cells During Microfluidic Stress Applications

Investigations in cellular and molecular engineering have explored the impact of nanotechnology and the potential for monitoring and control of human diseases. In a recent analysis, the dynamic fluid-induced stresses were characterized during microfluidic applications of an instrument with nanometer and picoNewton resolution as developed for single-cell biomechanics (Kohles, S. S., Nève, N., Zimmerman, J. D., and Tretheway, D. C., 2009, "Stress Analysis of Microfluidic Environments Designed for Isolated Biological Cell Investigations," ASME J. Biomech. Eng., 131(12), p. 121006). The results described the limited stress levels available in laminar, creeping-flow environments, as well as the qualitative cellular strain response to such stress applications. In this study, we present a two-dimensional computational model exploring the physical application of normal and shear stress profiles (with 0.1, 1.0, and 10.0 Pa peak amplitudes) potentially available within uniform and extensional flow states. The corresponding cellular strains and strain patterns were determined within cells modeled with healthy and diseased mechanical properties (5.0-0.1 kPa moduli, respectively). Strain energy density results integrated over the volume of the planar section indicated a strong mechanical sensitivity involving cells with disease-like properties. In addition, ex vivo microfluidic environments creating in vivo stress states would require freestream flow velocities of 2-7 mm/s. Knowledge of the nanomechanical stresses-strains necessary to illicit a biologic response in the cytoskeleton and cellular membrane will ultimately lead to refined mechanotransduction relationships.