Virtual blebbistatin: A robust and rapid software approach to motion artifact removal in optical mapping of cardiomyocytes.

Fluorescent reporters of cardiac electrophysiology provide valuable information on heart cell and tissue function. However, motion artifacts caused by cardiac muscle contraction interfere with accurate measurement of fluorescence signals. Although drugs such as blebbistatin can be applied to stop cardiac tissue from contracting by uncoupling calcium-contraction, their usage prevents the study of excitation-contraction coupling and, as we show, impacts cellular structure. We therefore developed a robust method to remove motion computationally from images of contracting cardiac muscle and to map fluorescent reporters of cardiac electrophysiological activity onto images of undeformed tissue. When validated on cardiomyocytes derived from human induced pluripotent stem cells (iPSCs), in both monolayers and engineered tissues, the method enabled efficient and robust reduction of motion artifact. As with pharmacologic approaches using blebbistatin for motion removal, our algorithm improved the accuracy of optical mapping, as demonstrated by spatial maps of calcium transient decay. However, unlike pharmacologic motion removal, our computational approach allowed direct analysis of calcium-contraction coupling. Results revealed calcium-contraction coupling to be more uniform across cells within engineered tissues than across cells in monolayer culture. The algorithm shows promise as a robust and accurate tool for optical mapping studies of excitation-contraction coupling in heart tissue.

The Balance between Actomyosin Contractility and Microtubule Polymerization Regulates Hierarchical Protrusions That Govern Efficient Fibroblast-Collagen Interactions.

Fibroblasts undergo a critical transformation from an initially inactive state to a morphologically different and contractile state after several hours of being embedded within a physiologically relevant three-dimensional (3D) fibrous collagen-based extracellular matrix (ECM). However, little is known about the critical mechanisms by which fibroblasts adapt themselves and their microenvironment in the earliest stage of cell-matrix interaction. Here, we identified the mechanisms by which fibroblasts interact with their 3D collagen fibrous matrices in the early stages of cell-matrix interaction and showed that fibroblasts use energetically efficient hierarchical micro/nano-scaled protrusions in these stages as the primary means for the transformation and adaptation. We found that actomyosin contractility in these protrusions in the early stages of cell-matrix interaction restricts the growth of microtubules by applying compressive forces on them. Our results show that actomyosin contractility and microtubules work in concert in the early stages of cell-matrix interaction to adapt fibroblasts and their microenvironment to one another. These early stage interactions result in responses to disruption of the microtubule network and/or actomyosin contractility that are opposite to well-known responses to late-stage disruption and reveal insight into the ways that cells adapt themselves and their ECM recursively.

Energy dissipation in quasi-linear viscoelastic tissues, cells, and extracellular matrix.

Characterizing how a tissue's constituents give rise to its viscoelasticity is important for uncovering how hidden timescales underlie multiscale biomechanics. These constituents are viscoelastic in nature, and their mechanics must typically be assessed from the uniaxial behavior of a tissue. Confounding the challenge is that tissue viscoelasticity is typically associated with nonlinear elastic responses. Here, we experimentally assessed how fibroblasts and extracellular matrix (ECM) within engineered tissue constructs give rise to the nonlinear viscoelastic responses of a tissue. We applied a constant strain rate, "triangular-wave" loading and interpreted responses using the Fung quasi-linear viscoelastic (QLV) material model. Although the Fung QLV model has several well-known weaknesses, it was well suited to the behaviors of the tissue constructs, cells, and ECM tested. Cells showed relatively high damping over certain loading frequency ranges. Analysis revealed that, even in cases where the Fung QLV model provided an excellent fit to data, the the time constant derived from the model was not in general a material parameter. Results have implications for design of protocols for the mechanical characterization of biological materials, and for the mechanobiology of cells within viscoelastic tissues.

Investigation of Nanoscopic Phase Separations in Lipid Membranes Using Inverse FCS.

Measurement of the sizes of nanoscopic particles is a difficult challenge, especially in two-dimensional systems such as cell membranes. We have extended inverse fluorescence correlation spectroscopy (iFCS) to endow it with unique advantages for measuring particle size from the nano- to the microscale. We have augmented iFCS with an analysis of moments of fluorescence fluctuations and used it to measure stages of phase separation in model lipid bilayer membranes. We observed two different pathways for the growth of phase domains. In one, nanoscopic gel domains appeared first and then gradually grew to micrometer size. In the other, the domains reached micrometer size quickly, and their number gradually increased. These measurements demonstrate the value of iFCS measurements through their ability, to our knowledge, to provide new information about the mechanism of lipid phase separation and potentially about the physical basis of naturally occurring nanodomains such as lipid rafts.

Fibroblasts Slow Conduction Velocity in a Reconstituted Tissue Model of Fibrotic Cardiomyopathy.

Myocardial function deteriorates over the course of fibrotic cardiomyopathy, due to electrophysiological and mechanical effects of myofibroblasts that are not completely understood. Although a range of experimental model systems and associated theoretical treatments exist at the levels of isolated cardiomyocytes and planar co-cultures of myofibroblasts and cardiomyocytes, interactions between these cell types at the tissue level are less clear. We studied these interactions through an engineered heart tissue (EHT) model of fibrotic myocardium and a mathematical model of the effects of cellular composition on EHT impulse conduction velocity. The EHT model allowed for modulation of cardiomyocyte and myofibroblast volume fractions, and observation of cell behavior in a three-dimensional environment that is more similar to native heart tissue than is planar cell culture. The cardiomyocyte and myofibroblast volume fractions determined the retardation of impulse conduction (spread of the action potential) in EHTs as measured by changes of the fluorescence of the Ca2+ probe, Fluo-2. Interpretation through our model showed retardation far in excess of predictions by homogenization theory, with conduction ceasing far below the fibroblast volume fraction associated with steric percolation. Results point to an important multiscale structural role of myofibroblasts in attenuating impulse conduction in fibrotic cardiomyopathy.

Remodeling by fibroblasts alters the rate-dependent mechanical properties of collagen.

UNLABELLED: The ways that fibroblasts remodel their environment are central to wound healing, development of musculoskeletal tissues, and progression of pathologies such as fibrosis. However, the changes that fibroblasts make to the material around them and the mechanical consequences of these changes have proven difficult to quantify, especially in realistic, viscoelastic three-dimensional culture environments, leaving a critical need for quantitative data. Here, we observed the mechanisms and quantified the mechanical effects of fibroblast remodeling in engineered tissue constructs (ETCs) comprised of reconstituted rat tail (type I) collagen and human fibroblast cells. To study the effects of remodeling on tissue mechanics, stress-relaxation tests were performed on ETCs cultured for 24, 48, and 72h. ETCs were treated with deoxycholate and tested again to assess the ECM response. Viscoelastic relaxation spectra were obtained using the generalized Maxwell model. Cells exhibited viscoelastic damping at two finite time constants over which the ECM showed little damping, approximately 0.2s and 10-30s. Different finite time constants in the range of 1-7000s were attributed to ECM relaxation. Cells remodeled the ECM to produce a relaxation time constant on the order of 7000s, and to merge relaxation finite time constants in the 0.5-2s range into a single time content in the 1s range. Results shed light on hierarchical deformation mechanisms in tissues, and on pathologies related to collagen relaxation such as diastolic dysfunction. STATEMENT OF SIGNIFICANCE: As fibroblasts proliferate within and remodel a tissue, they change the tissue mechanically. Quantifying these changes is critical for understanding wound healing and the development of pathologies such as cardiac fibrosis. Here, we characterize for the first time the spectrum of viscoelastic (rate-dependent) changes arising from the remodeling of reconstituted collagen by fibroblasts. The method also provides estimates of the viscoelastic spectra of fibroblasts within a three-dimensional culture environment. Results are of particular interest because of the ways that fibroblasts alter the mechanical response of collagen at loading frequencies associated with cardiac contraction in humans.

Efficient and optimized identification of generalized Maxwell viscoelastic relaxation spectra.

Viscoelastic relaxation spectra are essential for predicting and interpreting the mechanical responses of materials and structures. For biological tissues, these spectra must usually be estimated from viscoelastic relaxation tests. Interpreting viscoelastic relaxation tests is challenging because the inverse problem is expensive computationally. We present here an efficient algorithm that enables rapid identification of viscoelastic relaxation spectra. The algorithm was tested against trial data to characterize its robustness and identify its limitations and strengths. The algorithm was then applied to identify the viscoelastic response of reconstituted collagen, revealing an extensive distribution of viscoelastic time constants.

Confidence intervals for concentration and brightness from fluorescence fluctuation measurements.

The theory of photon count histogram (PCH) analysis describes the distribution of fluorescence fluctuation amplitudes due to populations of fluorophores diffusing through a focused laser beam and provides a rigorous framework through which the brightnesses and concentrations of the fluorophores can be determined. In practice, however, the brightnesses and concentrations of only a few components can be identified. Brightnesses and concentrations are determined by a nonlinear least-squares fit of a theoretical model to the experimental PCH derived from a record of fluorescence intensity fluctuations. The χ(2) hypersurface in the neighborhood of the optimum parameter set can have varying degrees of curvature, due to the intrinsic curvature of the model, the specific parameter values of the system under study, and the relative noise in the data. Because of this varying curvature, parameters estimated from the least-squares analysis have varying degrees of uncertainty associated with them. There are several methods for assigning confidence intervals to the parameters, but these methods have different efficacies for PCH data. Here, we evaluate several approaches to confidence interval estimation for PCH data, including asymptotic standard error, likelihood joint-confidence region, likelihood confidence intervals, skew-corrected and accelerated bootstrap (BCa), and Monte Carlo residual resampling methods. We study these with a model two-dimensional membrane system for simplicity, but the principles are applicable as well to fluorophores diffusing in three-dimensional solution. Using simulated fluorescence fluctuation data, we find the BCa method to be particularly well-suited for estimating confidence intervals in PCH analysis, and several other methods to be less so. Using the BCa method and additional simulated fluctuation data, we find that confidence intervals can be reduced dramatically for a specific non-Gaussian beam profile.

Physically-induced cytoskeleton remodeling of cells in three-dimensional culture.

Characterizing how cells in three-dimensional (3D) environments or natural tissues respond to biophysical stimuli is a longstanding challenge in biology and tissue engineering. We demonstrate a strategy to monitor morphological and mechanical responses of contractile fibroblasts in a 3D environment. Cells responded to stretch through specific, cell-wide mechanisms involving staged retraction and reinforcement. Retraction responses occurred for all orientations of stress fibers and cellular protrusions relative to the stretch direction, while reinforcement responses, including extension of cellular processes and stress fiber formation, occurred predominantly in the stretch direction. A previously unreported role of F-actin clumps was observed, with clumps possibly acting as F-actin reservoirs for retraction and reinforcement responses during stretch. Responses were consistent with a model of cellular sensitivity to local physical cues. These findings suggest mechanisms for global actin cytoskeleton remodeling in non-muscle cells and provide insight into cellular responses important in pathologies such as fibrosis and hypertension.

Fibrocartilage tissue engineering: the role of the stress environment on cell morphology and matrix expression.

Although much is known about the effects of uniaxial mechanical loading on fibrocartilage development, the stress fields to which fibrocartilaginous regions are subjected to during development are mutiaxial. That fibrocartilage develops at tendon-to-bone attachments and in compressive regions of tendons is well established. However, the three-dimensional (3D) nature of the stresses needed for the development of fibrocartilage is not known. Here, we developed and applied an in vitro system to determine whether fibrocartilage can develop under a state of periodic hydrostatic tension in which only a single principal component of stress is compressive. This question is vital to efforts to mechanically guide morphogenesis and matrix expression in engineered tissue replacements. Mesenchymal stromal cells in a 3D culture were exposed to compressive and tensile stresses as a result of an external tensile hydrostatic stress field. The stress field was characterized through mechanical modeling. Tensile cyclic stresses promoted spindle-shaped cells, upregulation of scleraxis and type one collagen, and cell alignment with the direction of tension. Cells experiencing a single compressive stress component exhibited rounded cell morphology and random cell orientation. No difference in mRNA expression of the genes Sox9 and aggrecan was observed when comparing tensile and compressive regions unless the medium was supplemented with the chondrogenic factor transforming growth factor beta3. In that case, Sox9 was upregulated under static loading conditions and aggrecan was upregulated under cyclic loading conditions. In conclusion, the fibrous component of fibrocartilage could be generated using only mechanical cues, but generation of the cartilaginous component of fibrocartilage required biologic factors in addition to mechanical cues. These studies support the hypothesis that the 3D stress environment influences cell activity and gene expression in fibrocartilage development.

Fatty acid 2-hydroxylase mediates diffusional mobility of Raft-associated lipids, GLUT4 level, and lipogenesis in 3T3-L1 adipocytes.

Straight chain fatty acid alpha-oxidation increases during differentiation of 3T3-L1 adipocytes, leading to a marked accumulation of odd chain length fatty acyl moieties. Potential roles of this pathway in adipocyte differentiation and lipogenesis are unknown. Mammalian fatty acid 2-hydroxylase (FA2H) was recently identified and suggested to catalyze the initial step of straight chain fatty acid alpha-oxidation. Accordingly, we examined whether FA2H modulates adipocyte differentiation and lipogenesis in mature adipocytes. FA2H level markedly increases during differentiation of 3T3-L1 adipocytes, and small interfering RNAs against FA2H inhibit the differentiation process. In mature adipocytes, depletion of FA2H inhibits basal and insulin-stimulated glucose uptake and lipogenesis, which are partially rescued by the enzymatic product of FA2H, 2-hydroxy palmitic acid. Expression of fatty-acid synthase and SCD1 was decreased in FA2H-depleted cells, and levels of GLUT4 and insulin receptor proteins were reduced. 2-Hydroxy fatty acids are enriched in cellular sphingolipids, which are components of membrane rafts. Accelerated diffusional mobility of raft-associated lipids was shown to enhance degradation of GLUT4 and insulin receptor in adipocytes. Consistent with this, depletion of FA2H appeared to increase raft lipid mobility as it significantly accelerated the rates of fluorescence recovery after photobleaching measurements of lipid rafts labeled with Alexa 488-conjugated cholera toxin subunit B. Moreover, the enhanced recovery rates were partially reversed by treatment with 2-hydroxy palmitic acid. In conclusion, our findings document the novel role of FA2H in adipocyte lipogenesis possibly by modulation of raft fluidity and level of GLUT4.

Stretch-activated force shedding, force recovery, and cytoskeletal remodeling in contractile fibroblasts.

The stress fiber network within contractile fibroblasts structurally reinforces and provides tension, or "tone", to tissues such as those found in healing wounds. Stress fibers have previously been observed to polymerize in response to mechanical forces. We observed that, when stretched sufficiently, contractile fibroblasts diminished the mechanical tractions they exert on their environment through depolymerization of actin filaments then restored tissue tension and rebuilt actin stress fibers through staged Ca(++)-dependent processes. These staged Ca(++)-modulated contractions consisted of a rapid phase that ended less than a minute after stretching, a plateau of inactivity, and a final gradual phase that required several minutes to complete. Active contractile forces during recovery scaled with the degree of rebuilding of the actin cytoskeleton. This complementary action demonstrates a programmed regulatory mechanism that protects cells from excessive stretch through choreographed active mechanical and biochemical healing responses.

A simplified approach to quasi-linear viscoelastic modeling.

The fitting of quasi-linear viscoelastic (QLV) constitutive models to material data often involves somewhat cumbersome numerical convolution. A new approach to treating quasi-linearity in 1-D is described and applied to characterize the behavior of reconstituted collagen. This approach is based on a new principle for including nonlinearity and requires considerably less computation than other comparable models for both model calibration and response prediction, especially for smoothly applied stretching. Additionally, the approach allows relaxation to adapt with the strain history. The modeling approach is demonstrated through tests on pure reconstituted collagen. Sequences of "ramp-and-hold" stretching tests were applied to rectangular collagen specimens. The relaxation force data from the "hold" was used to calibrate a new "adaptive QLV model" and several models from literature, and the force data from the "ramp" was used to check the accuracy of model predictions. Additionally, the ability of the models to predict the force response on a reloading of the specimen was assessed. The "adaptive QLV model" based on this new approach predicts collagen behavior comparably to or better than existing models, with much less computation.

Cellular and matrix contributions to tissue construct stiffness increase with cellular concentration.

The mechanics of bio-artificial tissue constructs result from active and passive contributions of cells and extracellular matrix (ECM). We delineated these for a fibroblast-populated matrix (FPM) consisting of chick embryo fibroblast cells in a type I collagen ECM through mechanical testing, mechanical modeling, and selective biochemical elimination of tissue components. From a series of relaxation tests, we found that contributions to overall tissue mechanics from both cells and ECM increase exponentially with the cell concentration. The force responses in these relaxation tests exhibited a logarithmic decay over the 3600 second test duration. The amplitudes of these responses were nearly linear with the amplitude of the applied stretch. The active component of cellular forces rose dramatically for FPMs containing higher cell concentrations.

Incremental mechanics of collagen gels: new experiments and a new viscoelastic model.

Paired incremental uniaxial step (i.e., relaxation) and ramp tests were conducted simultaneously on four (nominally) identical samples of type I collagen gel, over a direct strain range 0 < epsilon < 0.2. The paired step and ramp responses could not both be predicted by a simple viscoelastic constitutive relation (either linear or Fung-type), but could be predicted reasonably accurately by a general nonlinear viscoelastic relation with a strain-dependent relaxation spectrum, of the form sigma(t) = f(t)-infinity g(t-tau,epsilon)[d(epsilon)(tau)/d(tau)]d(tau). Based on a four-term exponential-series approximation, we measured the stiffness moduli and time constants of the relaxation function, g(t,epsilon), for the four gel samples that we tested, and found that the time constants were independent of strain but the moduli increased strongly with strain. Further, we found that the time constants did not vary across the four gels, but the moduli varied by a factor of about 2 across the gels. Some additional tests show features of the response of collagen gels to cycles of application and removal of loading.