Internal viscoelastic loading in cat papillary muscle.

The passive mechanical properties of myocardium were defined by measuring force responses to rapid length ramps applied to unstimulated cat papillary muscles. The immediate force changes following these ramps recovered partially to their initial value, suggesting a series combination of viscous element and spring. Because the stretched muscle can bear force at rest, the viscous element must be in parallel with an additional spring. The instantaneous extension-force curves measured at different lengths were nonlinear, and could be made to superimpose by a simple horizontal shift. This finding suggests that the same spring was being measured at each length, and that this spring was in series with both the viscous element and its parallel spring (Voigt configuration), so that the parallel spring is held nearly rigid by the viscous element during rapid steps. The series spring in the passive muscle could account for most of the series elastic recoil in the active muscle, suggesting that the same spring is in series with both the contractile elements and the viscous element. It is postulated that the viscous element might be coupled to the contractile elements by a compliance, so that the load imposed on the contractile elements by the passive structures is viscoelastic rather than purely viscous. Such a viscoelastic load would give the muscle a length-independent, early diastolic restoring force. The possibility is discussed that the length-independent restoring force would allow some of the energy liberated during active shortening to be stored and released during relaxation.     

An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment

This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola–Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.     

An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

Modeling of summation of individual twitches into unfused tetanus for various types of rat motor units.

Repeated stimulation of motor units (MUs) causes an increase of the force output that cannot be explained by linear summation of equal twitches evoked by the same stimulation pattern. To explain this phenomenon, an algorithm for reconstructing the individual twitches, that summate into an unfused tetanus is described in the paper. The algorithm is based on an analytical function for the twitch course modeling. The input parameters of this twitch model are lead time, contraction and half-relaxation times and maximal force. The measured individual twitches and unfused tetani at 10, 20, 30 and 40 Hz stimulation frequency of three rat motor units (slow, fast resistant to fatigue and fast fatigable) are processed. It is concluded that: (1) the analytical function describes precisely the course of individual twitches; (2) the summation of equal twitches does not follow the results from the experimentally measured unfused tetani, the differences depend on the type of the MU and are bigger for higher values of stimulation frequency and fusion index; (3) the reconstruction of individual twitches from experimental tetanic records can be successful if the tetanus is feebly fused (fusion index up to 0.7); (4) both the maximal forces and time parameters of individual twitches subtracted from unfused tetani change and influence the course of each tetanus. A discrepancy with respect to the relaxation phase was observed between experimental results and model prediction for tetani with fusion index exceeding 0.7. This phase was predicted longer than the experimental one for better fused tetani. Therefore, a separate series of physiological experiments and then, more complex model are necessary for explanation of this distinction.     

A viscoelastic model of mechanically induced and spontaneous contractions of vascular smooth muscle

A new viscoelastic model was developed for the mathematical characterisation of mechanically induced and intrinsic contractional responses of the vascular smooth muscle. To this end, the elastic and viscous analogue elements were supplemented with a new active element generating stress proportional to its momentary elongation. The four-element model consisting of an active element, a parallel viscous element and both series and parallel elastic elements predicted biphasic or damped oscillatory stress relaxation and creep responses which were similar to that found experimentally earlier. Above a certain exciting frequency the model exhibited dissipative and below energy producing behaviour, as indicated by the sign change of the hysteresis loop area. In the case of sinusoidal modulation of the stress generation parameter the model showed parametric resonance, which was regarded as a simulation of intrinsic oscillation of the smooth muscle.     

A model for cytoplasmic rheology consistent with magnetic twisting cytometry.

Magnetic twisting cytometry is gaining wide applicability as a tool for the investigation of the rheological properties of cells and the mechanical properties of receptor-cytoskeletal interactions. Current technology involves the application and release of magnetically induced torques on small magnetic particles bound to or inside cells, with measurements of the resulting angular rotation of the particles. The properties of purely elastic or purely viscous materials can be determined by the angular strain and strain rate, respectively. However, the cytoskeleton and its linkage to cell surface receptors display elastic, viscous, and even plastic deformation, and the simultaneous characterization of these properties using only elastic or viscous models is internally inconsistent. Data interpretation is complicated by the fact that in current technology, the applied torques are not constant in time, but decrease as the particles rotate. This paper describes an internally consistent model consisting of a parallel viscoelastic element in series with a parallel viscoelastic element, and one approach to quantitative parameter evaluation. The unified model reproduces all essential features seen in data obtained from a wide variety of cell populations, and contains the pure elastic, viscoelastic, and viscous cases as subsets.     

Models of skeletal muscle to explain the increase in passive stiffness in desmin knockout muscle

Absence of desmin in skeletal muscle was found to induce an increase in passive stiffness. The present study aimed at developing rheological models of passive muscle to explain this stiffening. Models were elaborated by using experimental data depicting muscle viscoelastic behaviour. The experimental protocol included stepwise extension tests applied on control and desmin knockout soleus muscles from mice. Linear and non-linear models were composed of elastic and viscous elements. They were constructed with the aim at taking the presence or absence of desmin into account by simulating desmin as an elastic element. Furthermore, associated adaptation of connective tissues in absence of desmin was modelled as an additional elastic element. Differences in passive behaviour induced by absence of desmin were predicted by using a linear model and a non-linear one. The non-linear model was selected because: (1) it is able to predict experimental viscoelastic kinetics accounting for the increase in passive stiffness in muscles lacking desmin, (2) its design is consistent with morphological data, and (3) stiffness characteristics of its elements are in accordance with the literature. Finally, this modelling approach demonstrates that both absence of desmin and adaptation of connective tissue are required to explain the increase in passive stiffness in desmin knockout muscles.     

Electrically measuring viscoelastic parameters of adherent cell layers under controlled magnetic forces

Electrical cell-substrate impedance sensing (ECIS) was used to measure the time-dependence and frequency-dependence of impedance for current flowing underneath and between cells. Osteosarcoma cells with a topology similar to a short cylinder (coin-like) surmounted by a dome were used in this study. Application of a small step increase in net vertical stress to the cells (4 and 7 dyn/cm²), via magnetic beads bound to the dorsal (upper) surface, causes an increase in cell body height and an increase in cell-cell separation, as well as stretching of the cell-substrate adhesion bonds. This results in a fast drop in measured resistance (less than 2 s), followed by a slower change with a time constant of 60–150 s. This time constant is about 1.5 times longer at 22 °C than that at 37 °C; it also increases with applied stress. Our frequency scan data, as well as our data for the time course of resistance and capacitance, show that the fast change is associated with both the under-the-cells and between-the-cells resistance. The slower change in resistance mainly reflects the between-the-cells resistance. To obtain viscoelastic parameters from our data we use a simple viscoelastic model comprising viscous and elastic elements (i.e., a dashpot and two springs) for the cell body, and an elastic element (a spring) for the cell-substrate adhesion system. Our results show that the spring constants and the viscosity of the cell body components of this viscoelastic model decrease as the temperature increases, whereas the elastic modulus of cell-substrate adhesion increases with temperature. At 37 °C, for the cell body we obtain a value of about 10⁵ P for the viscous element of the viscoelastic model, and a spring constant expressed in units of an elastic modulus of about 10⁴ dyn/cm² for the spring in series with the viscous element, with another spring with a modulus of about 2×10³ dyn/cm² in parallel with these. In comparable units, we have a modulus for the cell-substrate adhesion system of about 3×10³ dyn/cm². Received: 23 March 1998 / Revised version: 23 June 1998 / Accepted: 1 July 1998     

Mechanical properties of the passive myocardium: experiment and a mathematical model

A three-dimensional mathematical model of the rheological properties of a morphofunctional unit of the myocardium has been constructed, which consists of transversal and longitudinal elastic elements and tilted viscoelastic elements hinged without friction. The parameters of the viscosity and elasticity of structural elements of the model do not depend on the deformation value. The model makes it possible to adequately describe the features of the viscoelastic behavior of isolated samples of the passive myocardium and of the myocardium under longitudinal strain. A good agreement between the calculated and experimental data for an intact preparation of the rat myocardium and a preparation with removed cardiomyocytes has been shown.     

Hybrid duplex: a novel method to study the contractile function of heterogeneous myocardium

In an earlier study, we experimentally mimicked the effects of mechanical interaction between different regions of the ventricular wall by allowing pairs of independently maintained cardiac muscle fibers to interact mechanically in series or in parallel. This simple physiological model of heterogeneous myocardium, which has been termed "duplex," has provided new insight into basic effects of cardiac electromechanical heterogeneity. Here, we present a novel "hybrid duplex," where one of the elements is an isolated cardiac muscle and the other a "virtual cardiac muscle." The virtual muscle is represented by a computational model of cardiomyocyte electromechanical activity. We present in detail the computer-based digital control system that governs the mechanical interaction between virtual and biological muscle, the software used for data analysis, and working implementations of the model. Advantages of the hybrid duplex method are discussed, and experimental recordings are presented for illustration and as proof of the principle.     

A branching process model of gene amplification following chromosome breakage.

We have devised a mathematical model of gene amplification utilizing recent experimental observations concerning dihydrofolate reductase (DHFR) gene amplification in CHO cells. The mathematical model, based on a biological model which proposes that acentric elements are the initial intermediates in gene amplification, includes the following features: (1) initiation of amplification by chromosomal breakage to produce an acentric structure; (2) replication of acentric DNA, once per cell cycle; (3) dissociation of replicated acentric DNA; (4) unequal segregation of acentric DNA fragments to daughter cells at mitosis; (5) subsequent reintegration of acentric fragments into chromosomes. These processes are assumed to be independent for each element present in a cell at a given time. Thus, processes of unequal segregation and integration may occur in parallel, not necessarily in a unique sequence, and may be reiterated in one or multiple cell cycles. These events are described mathematically as a Galton-Watson branching process with denumerable infinity of object types. This mathematical model qualitatively and quantitatively reproduces the major elements of the dynamical behavior of DHFR genes observed experimentally. The agreement between the mathematical model and the experimental data lends credence to the biological model proposed by Windle et al. (1991), including the importance of chromosome breakage and subsequent gene deletion resulting from resection of the broken chromosome ends as initial events in gene amplification.     

Simulation of post-tetanic potentiation and fatigue in muscle using a visco-elastic model

Post-tetanic potentiation (PTP) in single motor units was simulated using a simple visco-elastic model. Single isometric twitches and unfused tetani were obtained using a wide range of physiological input rates. Values of model parameters were chosen to simulate contraction times close to those of fast and slow muscle fibers. PTP has been attributed either to i) an augmented plateau level of active state or ii) an increase in time constant of active state decay. Our results show that a prolonged decay time of active state can account for most of the experimental data obtained in amphibian and mammalian preparations. In particular, potentiation is more marked in unfused tetani than in single twitches. Moreover the model accounts for PTP even in the case of a reduction of active state plateau due to fatigue.     

A rheological motor model for vertebrate skeletal muscle in due consideration of non-linearity

A rheological motor model that satisfies the major mechanical properties of the skeletal muscle is proposed. The model consists of two Maxwell elements and a Voigt element connected in parallel with each other and has a force generator in it. The model well explains the mechanical behavior in quick and slow recovery phases in the isometric contraction of the muscle and achieves a sufficient isotonic shortening speed. The energy liberation of the motor in isotonic contraction is calculated and a mechanism of control is proposed, which operates so as to decrease the dissipated energy by altering the weights of the elastic and viscous constants in Maxwell elements. And thereby it becomes possible for the motor to possess non-linearity in energy liberation and load-velocity relation alike in muscle. The model would be a base model to be utilized for analyzing the kinetics of human macrosystems and/or for modeling the human neuromuscular system of motion control.     

Viscoelastic nature of isolated tomato (Solanum lycopersicum) fruit cuticles: a mathematical model

Viscoelastic behaviour of isolated tomato fruit cuticle (CM) is well known and extensively described. Temperature and hydration conditions modify the mechanical properties of CM. Mechanical data from previous transient‐creep analysis developed in tomato fruit cuticle under different temperature and hydration conditions have been used to propose a rheological model that describes the viscoelastic nature of CM. As a composite material, the biomechanical behaviour of the plant cuticle will depend not only on the mechanical characteristics of the individual components by themselves but also on the sum of them. Based on this previous information, we proposed a two‐element model to describe the experimental behaviour: an elastic hookean element connected in parallel to a viscous element or Voigt element that will describe the mechanical behaviour of the isolated CM and cutin under the studied conditions. The main parameters of the model, E₁ and E₂ will reflect the elastic and viscoelastic behaviour of the cuticle. Relationship between these physical parameters and the change in CM properties were discussed in order to elucidate the rheological processes taking place in CM. This model describes both the influence of temperature and hydration and the behaviour of the isolated cutin and the inferred contribution of the cuticle fraction of polysaccharides when the whole cuticle is tested.     

Parallel and serial processes in the human oculomotor system: bimodal integration and express saccades

Saccadic reaction times (SRTs) were analyzed in the context of stochastic models of information processing (e.g., Townsend and Ashby 1983) to reveal the processing architecture(s) underlying integrative interactions between visual and auditory inputs and the mechanisms of express saccades. The results support the following conclusions. Bimodal (visual + auditory) targets are processed in parallel, and facilitate SRT to an extent that exceeds levels attainable by probability summation. This strongly implies neural summation between elements responding to spatially aligned visual and auditory inputs in the human oculomotor system. Second, express saccades are produced within a separable processing stage that is organized in series with that responsible for intersensory integration. A model is developed that implements this combination of parallel and serial processing. The activity in parallel input channels is summed within a sensory stage which is organized in series with a pre-motor and motor stage. The time course of each subprocess is considered a random variable, and different experimental manipulations can selectively influence different stages. Parallels between the model and physiological data are explored.     

Dynamic stiffness profiles in the left ventricle.

Diastolic pressure-volume (P-V) curves were calculated on a beat-to-beat basis in the open-chest, pentobarbital-anesthetized dog, using the technique of direct transmitral flow measurement previously described. P-V curves were constructed and the slope (dP/dV) was plotted vs. pressure and time. dP/dV was used as an index of stiffness in each heart and its instantaneous changes with time were followed throughout the diastolic period. The end-diastolic P-V relation based on points from successive cycles during volume loading was found to be exponential. In contrast, the instantaneous P-V relation during any one diastolic period was not exponential. That is, the dynamic dP/dV vs. pressure plot was nonlinear. In the normal heart, stiffness was characterized in early diastole by a negative dP/dV as the ventricle continued to relax, and then frequently decreased prior to a second stiffness rise with atrial augmentation. These findings can be explained by a model containing an element whose deformation is rate dependent, i.e., a parallel viscous element. Stiffness profiles in mitral stenosis where dynamic effects are minimized substantiate this conclusion.     

Effect of colchicine on viscoelastic properties of neutrophils

The effect of colchicine (15-60 micrograms/ml) on the viscoelastic properties of human neutrophils was studied by the micropipette technique. The small deformation of the neutrophil in response to a step aspiration pressure was analyzed by using a three-element model in which an elastic element, K1, is in parallel with a Maxwell element composed of another elastic element, K2, in series with a viscous element, mu. Colchicine treatment of neutrophils caused decreases in K2 and mu without affecting K1. The results indicate that the integrity of the microtubules plays a significant role in providing the viscoelastic resistance (as represented by the Maxwell element in the model) of neutrophils to deforming stress.     

A macroscopic ansatz to deduce the Hill relation

In this study, we derive the hyperbolic force-velocity relation of concentric muscular contraction, first formulated empirically by A.V. Hill in 1938, from three essential model assumptions: (1) the structural assembly of three well-known elements - i.e. active, parallel damping, and serial - fulfilling a force equilibrium, (2) the parallel damping coefficient explicitly depending on muscle force output and three parameters, and (3) the kinematic gearing ratio between active and serial element being assigned to a parameter. The energy source within the muscle represented by the force of the active element is an additional fifth parameter. As a result we find the Hill “constants” A and B as functions of our five model parameters. Using A and B values from literature on experimental data, we predict heat power release of our model. By calculating enthalpy rate and mechanical efficiency, we compare the model heat power to predictions from another Hill-type model, to Hill's original findings, and to findings from modern muscle heat measurements. We reconsider why the biggest share of heat rate during isometric contractions (maintenance heat) and the velocity-dependent heat rate during concentric contractions in addition to maintenance heat rate (shortening heat rate) may be traced back to the same mechanism represented by the kinematic gearing ratio. Namely, we suggest that the serial element transfers attachment-detachment fluctuations of actin-myosin crossbridges within one sarcomere to others in the same sarcomere and to those in parallel and in series. Numerically, in case of negligible passive muscular damping, we find the ratio between A and isometric force (relative A) to depend exclusively on the kinematic gearing ratio, whereas the maintenance heat rate scales with the square of relative A. Moreover, this mechanical coupling internal to the muscle fibres may also be behind the macroscopic force dependency of the overall parallel damping coefficient.     

Effects of transmural electrical heterogeneities and electrotonic interactions on the dispersion of cardiac repolarization and action potential duration: A simulation study

It has been shown in the literature that myocytes isolated from the ventricular walls at various intramural depths have different action potential durations (APDs). When these myocytes are embedded in the ventricular wall, their inhomogeneous properties affect the sequence of repolarization and the actual distribution of the APDs in the entire wall. In this article, we implement a mathematical model to simulate the combined effect of (a) the non-homogeneous intrinsic membrane properties (in particular the non-homogeneous APDs) and (b) the electrotonic currents that modulate the APDs when the myocytes are embedded in the ventricular myocardium. In particular, we study the effect of (a) and (b) on the excitation and repolarization sequences and on the distribution of APDs in the ventricles. We implement a Monodomain tissue representation that includes orthotropic anisotropy, transmural fiber rotation and homogeneous or heterogeneous transmural intrinsic membrane properties, modeled according to the phase I Luo-Rudy membrane ionic model. Three-dimensional simulations are performed in a cartesian slab with a parallel finite element solver employing structured isoparametric trilinear finite elements in space and a semi-implicit adaptive method in time. Simulations of excitation and repolarization sequences elicited by epicardial or endocardial pacing show that in a homogeneous slab the repolarization pathways approximately follow the activation sequence. Conversely, in the heterogeneous cases considered in this study, we observed two repolarization wavefronts that started from the epi and the endocardial faces respectively and collided in the thickness of the wall and in one case an additional repolarization wave starting from an intramural site. Introducing the heterogeneities along the transmural epi-endocardial direction affected both the repolarization sequence and the APD dispersion, but these effects were clearly discernible only in transmural planes. By contrast, in planes parallel to epi- and endocardium the APD distribution remained remarkably similar to that observed in the homogeneous model. Therefore, the patterns of the repolarization sequence and APD dispersion on the epicardial surface (or any other intramural surface parallel to it) do not reveal the uniform transmural heterogeneity.