Search PubMed⌕ Search

Biomedical subjects

G I Zahalak

Publications and source records attributed to G I Zahalak.

At least 19 recordsLinked to original sources

Theoretical model for myocardial trabeculation.

During the morphogenetic process of myocardial trabeculation, most of the cardiac jelly of the initially smooth-walled heart is replaced by sponge-like muscle. The mechanisms that drive and regulate this important process are poorly understood. Using a theoretical model, we examined the possible role that cytoskeletal contraction plays during the initial stages of trabeculation. The myocardium is modeled as a thin viscoelastic membrane consisting of contractile (stress) fibers embedded in an isotropic incompressible matrix, with the interaction of myocardial cells and cardiac jelly fibers providing long-range mechanical effects. The stress fibers are assumed to behave like smooth muscle and to normally operate on the descending limb of their stress-stretch curve. Mechanical instability due to the effectively negative stiffness then leads to the creation of pattern. As a first approximation, computations were carried out for a flat rectangular membrane with stress fibers aligned along a single direction. The computed deformation patterns depend strongly on the magnitude and anisotropy of the long-range effects. Given plausible assumptions about the mechanical properties of the embryonic heart, the model predicts trabecular patterns similar to those observed in the embryo, including the development of circumferential ridges and relatively thin regions ("holes") in the trabecular sheets.

Animals↗

The two-state cross-bridge model of muscle is an asymptotic limit of multi-state models.

The relationship between the two-state model of muscle contraction and multi-state models is examined from the perspective of matched asymptotic expansions, under the assumption that transition rates between attached states are fast compared to those between detached and attached states. A detailed formal analysis of a three-state model reveals that the classic Huxley (1957. Prog. Biophys. Biophys. Chem.7, 225-318) rate equation, as modified for thermodynamic self-consistency by Hill et al. (1975. Biophys. J.15, 335-372), governs the "outer" solution of the three-state equations. Thus, the two-state model remains a valid description of muscle dynamics on physiologically relevant time scales, which are slow compared to millisecond-scale transitions between attached states. But the asymptotic analysis reveals also that the cross-bridge force must be considered to be a nonlinear function of the cross-bridge strain, in contrast to the usual assumption of two-state models. This apparent, or effective, force is determined by both the intrinsic stiffness of the cross-bridge and the equilibrium distribution of cross-bridges among attached states. Further, the asymptotic analysis yields an expression for the energy liberation rate that implies a reduced rate in stretch vs. shortening. Some behaviors of multi-state models that are suggested by the three-state analysis are discussed in qualitative terms.

Biomechanical Phenomena↗

Cell mechanics studied by a reconstituted model tissue.

Tissue models reconstituted from cells and extracellular matrix (ECM) simulate natural tissues. Cytoskeletal and matrix proteins govern the force exerted by a tissue and its stiffness. Cells regulate cytoskeletal structure and remodel ECM to produce mechanical changes during tissue development and wound healing. Characterization and control of mechanical properties of reconstituted tissues are essential for tissue engineering applications. We have quantitatively characterized mechanical properties of connective tissue models, fibroblast-populated matrices (FPMs), via uniaxial stretch measurements. FPMs resemble natural tissues in their exponential dependence of stress on strain and linear dependence of stiffness on force at a given strain. Activating cellular contractile forces by calf serum and disrupting F-actin by cytochalasin D yield "active" and "passive" components, which respectively emphasize cellular and matrix mechanical contributions. The strain-dependent stress and elastic modulus of the active component were independent of cell density above a threshold density. The same quantities for the passive component increased with cell number due to compression and reorganization of the matrix by the cells.

Animals↗

A cell-based constitutive relation for bio-artificial tissues.

By using a combination of continuum and statistical mechanics we derive an integral constitutive relation for bio-artificial tissue models consisting of a monodisperse population of cells in a uniform collagenous matrix. This constitutive relation quantitatively models the dependence of tissue stress on deformation history, and makes explicit the separate contribution of cells and matrix to the mechanical behavior of the composite tissue. Thus microscopic cell mechanical properties can be deduced via this theory from measurements of macroscopic tissue properties. A central feature of the constitutive relation is the appearance of "anisotropy tensors" that embody the effects of cell orientation on tissue mechanics. The theory assumes that the tissues are stable over the observation time, and does not in its present form allow for cell migration, reorientation, or internal remodeling. We have compared the predictions of the theory to uniaxial relaxation tests on fibroblast-populated collagen matrices (FPMs) and find that the experimental results generally support the theory and yield values of fibroblast contractile force and stiffness roughly an order of magnitude smaller than, and viscosity comparable to, the corresponding properties of active skeletal muscle. The method used here to derive the tissue constitutive equation permits more sophisticated cell models to be used in developing more accurate representations of tissue properties.

Animals↗

The effects of cross-fiber deformation on axial fiber stress in myocardium.

We incorporated a three-dimensional generalization of the Huxley cross-bridge theory in a finite element model of ventricular mechanics to examine the effect of nonaxial deformations on active stress in myocardium. According to this new theory, which assumes that macroscopic tissue deformations are transmitted to the myofilament lattice, lateral myofilament spacing affects the axial fiber stress. We calculated stresses and deformations at end-systole under the assumption of strictly isometric conditions. Our results suggest that at the end of ejection, nonaxial deformations may significantly reduce active axial fiber stress in the inner half of the wall of the normal left ventricle (18-35 percent at endocardium, depending on location with respect to apex and base). Moreover, this effect is greater in the case of a compliant ischemic region produced by occlusion of the left anterior descending or circumflex coronary artery (26-54 percent at endocardium). On the other hand, stiffening of the remote and ischemic regions (in the case of a two-week-old infarct) lessens the effect of nonaxial deformation on active stress at all locations (9-32 percent endocardial reductions). These calculated effects are sufficiently large to suggest that the influence of nonaxial deformation on active fiber stress may be important, and should be considered in future studies of cardiac mechanics.

Coronary Vessels↗

A distribution-moment model of deactivation in cardiac muscle.

The Distribution Moment (DM) model has simulated experimental data on skeletal muscle, but it has not been used previously to study the mechanics of active contraction in cardiac muscle. In contrast to previous models of striated muscle contraction, all parameters have physical meaning and assumptions concerning biophysical events within the cell are consistent with available data. In order to simulate cardiac muscle deactivation using the DM model it was necessary to make the cross-bridge detachment rates large for large displacements from the neutral equilibrium position of a cross-bridge. To examine the effect of cooperativity on cardiac muscle contraction, we used the DM model's tight coupling scheme with binding of one or two calcium sites regulating contraction. As observed experimentally, our model predicted a reduction of isometric tension development following rapid shortening lengthening transients when contraction is regulated by either one or two calcium binding sites. The predicted deactivating effect increased if the transient was applied late in the twitch when contraction is regulated by two calcium binding sites, but not when it is regulated by one site. This is the first study in which deactivation has been simulated without making any provisions for length-dependent calcium trononin dissociation.

Binding Sites↗

A re-examination of calcium activation in the Huxley cross-bridge model.

This paper investigates mathematical relations between models of calcium activation kinetics and Huxley-type models of cross-bridge dynamics in muscle. It is found that different calcium-activation schemes lead to the same form of generalized Huxley rate equation with calcium activation. [formula: see text] if it is assumed that calcium-troponin interaction rates are fast compared to the rates of transition associated with force-generating cross-bridge states. Calcium affects cross-bridge dynamics by modifying the bonding rate f, but does not affect the number of interacting cross bridges a or the unbonding rate g; this occurs through the appearance in the equation of an activation factor, r, which is a pure function of sarcoplasmic free calcium concentration. In particular, it is shown that both the "tight-coupling" and "loose-coupling" calcium-activation schemes introduced by Zahalak and Ma [1] lead to the same rate equation with the same activation factor; the difference between them appears in the calcium mass-balance equation. While both of these activation models can be made to fit simple twitch and force-velocity data equally well, experimentally observed load-dependent shifts in the free calcium concentration are compatible with the right-coupling scheme, but not with loose coupling.

Actins↗

Non-axial muscle stress and stiffness.

A generalization is developed of the classic two-state Huxley cross-bridge model to account for non-axial active stress and stiffness. The main ingredients of the model are: (i) a relation between the general three-dimensional deformation of an element of muscle and the deformations of the cross-bridges, that assumes macroscopic deformation is transmitted to the myofibrils, (ii) radial as well as axial cross-bridge stiffness, and (iii) variations of the attachment and detachment rates with lateral filament spacing. The theory leads to a generalized Huxley rate equation on the bond-distribution function, n(zeta, theta, t), of the form [equation: see text] where the Dij are the components of the relative velocity gradient and rho and ñ are functions of the polar angle, theta, and time that describe, respectively, the deformation of the myofilament lattice and the distribution of accessible actin sites (both of these functions can be calculated from the macroscopic deformation). Explicit expressions, in terms of n, are derived for the nine components of the active stress tensor, and the 21 non-vanishing components of the active stiffness tensor; the active stress tensor is found to be unsymmetric. The theory predicts that in general non-axial deformations will modify active axial stress and stiffness, and also give rise to non-axial (e.g., shearing) components. Under most circumstances the magnitudes of the non-axial stress and stiffness components will be small compared with the axial and, further, the effects of non-axial deformation rates will be small compared with those of the axial rate. Large transverse deformations may, however, greatly reduce the axial force and stiffness. The theory suggests a significant mechanical role for the non-contractile proteins in muscle, namely that of equilibrating the unsymmetric active stresses. Some simple applications of the theory are provided to illustrate its physical content.

Animals↗

The effect of hyperosmolality on the rate of heat production of quiescent trabeculae isolated from the rat heart.

We have measured the rate of heat production of isolated, quiescent, right ventricular trabeculae of the rat under isosmotic and hyperosmotic conditions, using a microcalorimetric technique. In parallel experiments, we measured force production and intracellular calcium concentration ([Ca2+]i). The rate of resting heat production under isosmotic conditions (mean +/- SEM, n = 32) was 100 +/- 7 mW (g dry wt)-1; it increased sigmoidally with osmolality, reaching a peak that was about four times the isosmotic value at about twice normal osmotic pressure. The hyperosmotic thermal response was: (a) abolished by anoxia, (b) attenuated by procaine, (c) insensitive to verapamil, ouabain, and external calcium concentration, and (d) absent in chemically skinned trabeculae bathed in low-Ca2+ "relaxing solution." Active force production was inhibited at all osmolalities above isosmotic. Passive (tonic) force increased to, at most, 15% of the peak active force developed under isosmotic conditions while [Ca2+]i increased, at most, 30% above its isosmotic value. We infer that hyperosmotic stimulation of resting cardiac heat production reflects, in large part, greatly increased activity of the sarcoplasmic reticular Ca2+ ATPase in the face of increased efflux via a procaine-inhibitable Ca(2+)-release channel.

Animals↗

The variation of isometric energy rates with muscle length: a distribution-moment model analysis.

The Distribution-Moment Model of skeletal muscle, which has been enhanced recently to make possible the calculation of chemical energy release (E) and heat production (H) rates [1], is applied to isometric muscle. Under steady-state isometric conditions the model predicts a simple relation between the energy rates and the muscle length, namely (E/Emax) = (H/Hmax) = [1 + B alpha(symbol see text)]/[1 + B], where (symbol see text) is the ratio of muscle length to the "optimal" length at which maximal isometric tension is produced, and alpha (symbol see text) is a function numerically equal to the ratio of the tetanic isometric force to its maximum value. The single dimensionless constant in this relation, B, can be calculated from model parameters characterizing muscle dynamics at the optimum length, and has a value near unity for frog sartorius at 0 degrees C. The predicted behavior is shown to agree reasonably well with experimental measurements of heat production and phosphocreatine (PCr) hydrolysis. The model relates the isometric energy rates to PCr hydrolysis in (1) cross-bridge interactions, and (2) calcium pumping into the sarcoplasmic reticulum.

Biomechanical Phenomena↗

A distribution-moment model of energetics in skeletal muscle.

In this paper we develop a theory for calculating the chemical energy liberation and heat production of a skeletal muscle subjected to an arbitrary history of stimulation, loading, and length variation. This theory is based on and complements the distribution-moment (DM) model of muscle [Zahalak and Ma, J. biomech. Engng 112, 52-62 (1990)]. The DM model is a mathematical approximation of the A. F. Huxley cross-bridge theory and represents a muscle in terms of five (normalized) state variables: A, the muscle length, c, the sarcoplasmic free calcium concentration, and Q0, Q1, Q2, the first three moments of the actin-myosin bond-distribution function (which, respectively, have macroscopic interpretations as the muscle stiffness, force, and elastic energy stored in the contractile tissue). From this model are derived two equations which predict the chemical energy liberation and heat production rates in terms of the five DM state variables, and which take account of the following factors: (1) phosphocreatine hydrolysis associated with cross-bridge cycling; (2) phosphocreatine hydrolysis associated with sarcoplasmic-reticulum pumping of calcium; (3) passive calcium flux across the sarcoplasmic-reticulum membrane; (4) calcium-troponin bonding; (5) cross-bridge bonding at zero strain; (6) cross-bridge strain energy; (7) tendon strain energy; and (8) external work. Using estimated parameters appropriate for a frog sartorius at 0 degree C, the energy rates are calculated for several experiments reported in the literature, and reasonable agreement is found between our model and the measurements. (The selected experiments are confined to the plateau of the isometric length-tension curve, although our theory admits arbitrary length variations.) The two most important contributions to the energy rates are phosphocreatine hydrolysis associated with cross-bridge cycling and with sarcoplasmic-reticulum calcium pumping, and these two contributions are approximately equal under tetanic, isometric, steady-state conditions. The contribution of the calcium flux across the electrochemical potential gradient at the sarcoplasmic-reticulum membrane was found to be small under all conditions examined, and can be neglected. Long-term fatigue and oxidative recovery effects are not included in this theory. Also not included is the so-called 'unexplained energy' presumably associated with reactions which have not yet been identified. Within these limitations our model defines clear quantitative interrelations between the activation, mechanics, and energetics in muscle, and permits rational estimates of the energy production to be calculated for arbitrary programs of muscular work.

Actins↗

Muscle activation and contraction: constitutive relations based directly on cross-bridge kinetics.

This paper reviews recent work aimed at deriving tractable constitutive relations for skeletal muscle from biophysical cross-bridge theories. Discussion is focused on a model proposed previously by the first author (the Distribution-Moment Model), which emphasizes the important role of the moments of the actin-myosin bond-distribution function. The theory leads to a relatively simple third order state variable model for contraction dynamics in which the state variables are the three lowest order moments of the bond-distribution function; further, these three moments have simple macroscopic interpretations as muscle stiffness, force, and elastic energy. New results are presented on the formulation of a compatible model for excitation-contraction coupling, and this model requires the introduction of only one more state variable--the free calcium concentration.

Actins↗

Extensional recovery of an intact erythrocyte from a tank-treading motion.

Normal human erythrocytes suspended in shear flow are stretched into quasi ellipsoidal forms while their membranes rotate smoothly (tank-treading). Following abrupt cessation of shear the cells recover their discoidal shapes approximately exponentially, in the manner of a Kelvin-Voigt (K-V) solid. To test the hypothesis that the recovery process is membrane-controlled, the effects of initial deformation, cytoplasmic viscosity and membrane surface-to-volume ratio were studied. It was concluded that the membrane dynamics dominates the transient shape recovery, and that the characteristic recovery time is dependent on the initial deformation. Hence, the usual simplified analysis based on retraction of a plane sheet of K-V material with constant moduli appears to be an inadequate treatment of transient whole cell recovery.

Biomechanical Phenomena↗

Determination of cellular mechanical properties by cell poking, with an application to leukocytes.

In this paper we review the cell-poking technique as an approach for investigating the mechanical properties of living cells. We first summarize the rationale for the technique and the mainly qualitative results obtained so far. Then we provide a technical description of the instrument as it is configured at present. This is followed by a discussion of the current status of analytical results available for interpreting cell-poking measurements. In the final section we apply these results to an analysis of unmodulated and modulated lymphocytes and neutrophils, and conclude that the mechanical response of these leukocytes to indentation is not consistent with simple models developed by previous investigators on the basis of micropipette-aspiration experiments.

Animals↗

Cell poking: quantitative analysis of indentation of thick viscoelastic layers.

A recently introduced device, the cell poker, measures the force required to indent the exposed surface of a cell adherent to a rigid substratum. The cell poker has provided phenomenological information about the viscoelastic properties of several different types of cells, about mechanical changes triggered by external stimuli, and about the role of the cytoskeleton in these mechanical functions. Except in special cases, however, it has not been possible to extract quantitative estimates of viscosity and elasticity moduli from cell poker measurements. This paper presents cell poker measurements of well characterized viscoelastic polymeric materials, polydimethylsiloxanes of different degrees of polymerization, in a simple shape, a flat, thick layer, which for our purposes can be treated as a half space. Analysis of the measurements in terms of a linear viscoelasticity theory yields viscosity values for three polymer samples in agreement with those determined by measurements on a macroscopic scale. Theoretical analysis further indicates that the measured limiting static elasticity of the layers may result from the tension generated at the interface between the polymer and water. This work demonstrates the possibility of obtaining quantitative viscoelastic material properties from cell poker measurements and represents the first step in extending these quantitative studies to more complicated structures including cells.

Cell Adhesion↗

A practical electrode-array myoprocessor for surface electromyography.

A new myoprocessor is described which produces a relatively smooth and accurate measure of a muscle's electrical activity during both slow and rapid movements. This is achieved partly by spatial averaging of weakly correlated EMG signals from electrodes distributed over the surface of the muscle.

Electromyography↗

Deduction of intrinsic mechanical properties of the erythrocyte membrane from observations of tank-treading in the rheoscope.

Measurements of the dimensions and membrane rotational frequency of individual erythrocytes steadily tank-treading in a Rheoscope are used to deduce the surface shear viscosity (eta m) and the shear elastic modulus (mu m) of the membrane. Previously published algorithms (Trans-Son-Tay et al., Biophys. J. 46: 65, 1984, and 51: 915, 1987) plus an assumed area-conserving membrane velocity field (Secomb and Skalak, Q. J. Mech. Appl. Math. XXXV 2: 233, 1982) are applied to calculate eta m as a function of the second invariant of the surface strain rate and mu m as a function of the second invariant of membrane strain. The results indicate density-related increases in membrane stiffness and viscosity, shear-thinning viscous behavior, and strain-stiffening elastic behavior.

Algorithms↗

Membrane stress and internal pressure in a red blood cell freely suspended in a shear flow.

Presented is an algorithm for the approximate calculation of the membrane stress distribution and the internal pressure of a steadily tank-treading red cell. The algorithm is based on an idealized ellipsoidal model of the tank-treading cell (Keller, S.R., and R. Skalak, 1982, J. Fluid Mech., 120:27-47) joined with experimental observations of projected length, width, and tank-treading frequency. The results are inexact because the membrane shape and velocity are assumed a priori, rather than being determined via appropriate material constitutive relations for the membrane; these results are, nevertheless, believed to be approximately correct, and show that internal pressure builds up slowly as cell elongation increases, rising more rapidly as the deformed cell approaches the limiting geometry of a prolate ellipsoid. The maximum shear stress resultant in the membrane was found to be below but approaching the yield point range at the highest shear rate applied.

Algorithms↗