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.

Fluorescence Correlation Spectroscopy and Photon Counting Histograms in Finite, Bounded Domains.

Analysis of fluctuations arising as fluorescent particles pass through a focused laser beam has enabled quantitative characterization of a broad range of molecular kinetic processes. Two key mathematical frameworks that have enabled these quantifications are fluorescence correlation spectroscopy (FCS) and photon counting histogram (PCH) analysis. Although these frameworks are effective and accurate when the focused laser beam is well approximated by an infinite Gaussian beam with a waist that is small compared to the size of the region over which the fluorescent particles can diffuse, they cannot be applied to situations in which this region is bounded at the nanoscale. We therefore derived general forms of the FCS and PCH frameworks for bounded systems. The finite-domain form of FCS differs from the classical form in its boundary and initial conditions and requires development of a new Fourier space solution for fitting data. Our finite-domain FCS predicts simulated data accurately and reduces to a previous model for the special case when the system is much larger than the Gaussian beam and can be considered to be infinite. We also derived the PCH form for the bounded systems. Our approach enables estimation of the concentration of diffusing fluorophores within a finite domain for the first time, to our knowledge. The method opens the possibility of quantification of kinetics in several systems for which this has never been possible.

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.

A review on qualifications and cost effectiveness of induced pluripotent stem cells (IPSCs)-induced cardiomyocytes in drug screening tests.

Human induced pluripotent stem cells (hIPSCs) have initiated a higher degree of successes in disease modelling, preclinical evaluation of drug therapy and pharmaco-toxicological testing. Since the discovery of iPSCs in 2006, many advanced techniques have been introduced to differentiate iPSCs to cardiomyocytes, which have been progressively improved. The disease models from iPSC-induced cardiomyocytes (iPSC-CM) have been successfully helping to study a variety of cardiac diseases such as long QT syndrome, drug-induced long QT, different cardiomyopathies related to mutations in mitochondria or desmosomal proteins and other rare genetic diseases. IPSC-CMs have also been used to screen the role of chemicals in cardiovascular drug discovery and individualisation of drug dosages. In this review, the quality of current procedures for characterisation and maturation of iPSC-CM lines will be discussed. Also, we will focus on time efficiency and cost of standard differentiation methods after reprogramming.

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.

Short communication: vascular smooth muscle cell stiffness as a mechanism for increased aortic stiffness with aging.

RATIONALE: Increased aortic stiffness, an important feature of many vascular diseases, eg, aging, hypertension, atherosclerosis, and aortic aneurysms, is assumed because of changes in extracellular matrix (ECM). OBJECTIVE: We tested the hypothesis that the mechanisms also involve intrinsic stiffening of vascular smooth muscle cells (VSMCs). METHODS AND RESULTS: Stiffness was measured in vitro both by atomic force microscopy (AFM) and in a reconstituted tissue model, using VSMCs from aorta of young versus old male monkeys (Macaca fascicularis) (n=7/group), where aortic stiffness increases by 200% in vivo. The apparent elastic modulus was increased (P<0.05) in old (41.7+/-0.5 kPa) versus young (12.8+/-0.3 kPa) VSMCs but not after disassembly of the actin cytoskeleton with cytochalasin D. Stiffness of the VSMCs in the reconstituted tissue model was also higher (P<0.05) in old (23.3+/-3.0 kPa) than in young (13.7+/-2.4 kPa). CONCLUSIONS: These data support the novel concept, not appreciated previously, that increased vascular stiffness with aging is attributable not only to changes in ECM but also to intrinsic changes in VSMCs.

Fluorescence correlation spectroscopy: past, present, future.

In recent years fluorescence correlation spectroscopy (FCS) has become a routine method for determining diffusion coefficients, chemical rate constants, molecular concentrations, fluorescence brightness, triplet state lifetimes, and other molecular parameters. FCS measures the spatial and temporal correlation of individual molecules with themselves and so provides a bridge between classical ensemble and contemporary single-molecule measurements. It also provides information on concentration and molecular number fluctuations for nonlinear reaction systems that complement single-molecule measurements. Typically implemented on a fluorescence microscope, FCS samples femtoliter volumes and so is especially useful for characterizing small dynamic systems such as biological cells. In addition to its practical utility, however, FCS provides a window on mesoscopic systems in which fluctuations from steady states not only provide the basis for the measurement but also can have important consequences for the behavior and evolution of the system. For example, a new and potentially interesting field for FCS studies could be the study of nonequilibrium steady states, especially in living cells.

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.

The biochemical response of the heart to hypertension and exercise.

Mechanical stress on the heart can lead to crucially different outcomes. Exercise is beneficial because it causes heart muscle cells to enlarge (hypertrophy). Chronic hypertension also causes hypertrophy, but in addition it causes an excessive increase in fibroblasts and extracellular matrix (fibrosis), death of cardiomyocytes and ultimately heart failure. Recent research shows that stimulation of physiological (beneficial) hypertrophy involves several signaling pathways, including those mediated by protein kinase B (also known as Akt) and the extracellular-signal-regulated kinases 1 and 2 (ERK1/2). Hypertension, beta-adrenergic stimulation and agonists such as angiotensin II (Ang II) activate not only ERK1/2 but also p38 and the Jun N-terminal kinase (JNK), leading to pathological heart remodeling. Despite this progress, the mechanisms that activate fibroblasts to cause fibrosis and those that differentiate between exercise and hypertension to produce physiological and pathological responses, respectively, remain to be established.

A model for positive feedback control of the transformation of fibroblasts to myofibroblasts.

The phenotypic conversion of normal fibroblasts to myofibroblasts is central to normal wound healing and to pathological fibrosis that can occur in the heart and many other tissues. The transformation occurs in two stages. The first stage is driven mainly by mechanical changes such as increased stiffness of the heart due to hypertension and cellular contractility. The second stage requires both increasing stiffness and biochemical factors such as the growth factor, TGFß. As more and more cells convert from weakly contractile fibroblasts to strongly contractile myofibroblasts, the stiffness of the ventricular muscle increases. We propose a simple model for the establishment of non-equilibrium steady states with different compositions of fibroblasts and myofibroblasts. Under some conditions a positive feedback loop resulting from the increasing stiffness caused by increasing numbers of myofibroblasts can produce a bifurcation between steady states with low and high myofibroblast content. We illustrate the large mechanical differences between normal fibroblasts and myofibroblasts with measurements in engineered tissue constructs.

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.

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.

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.

Discrete quasi-linear viscoelastic damping analysis of connective tissues, and the biomechanics of stretching.

The time- and frequency-dependent properties of connective tissue define their physiological function, but are notoriously difficult to characterize. Well-established tools such as linear viscoelasticity and the Fung quasi-linear viscoelastic (QLV) model impose forms on responses that can mask true tissue behavior. Here, we applied a more general discrete quasi-linear viscoelastic (DQLV) model to identify the static and dynamic time- and frequency-dependent behavior of rabbit medial collateral ligaments. Unlike the Fung QLV approach, the DQLV approach revealed that energy dissipation is elevated at a loading period of ∼10s. The fitting algorithm was applied to the entire loading history on each specimen, enabling accurate estimation of the material's viscoelastic relaxation spectrum from data gathered from transient rather than only steady states. The application of the DQLV method to cyclically loading regimens has broad applicability for the characterization of biological tissues, and the results suggest a mechanistic basis for the stretching regimens most favored by athletic trainers.

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.