A continuous-binding cross-linker model for passive airway smooth muscle

Although the active properties of airway smooth muscle (ASM) have garnered much modeling attention, the passive mechanical properties are not as well studied. In particular, there are important dynamic effects observed in passive ASM, particularly strain-induced fluidization, which have been observed both experimentally and in models; however, to date these models have left an incomplete picture of the biophysical, mechanistic basis for these behaviors. The well-known Huxley cross-bridge model has for many years successfully described many of the active behaviors of smooth muscle using sliding filament theory; here, we propose to extend this theory to passive biological soft tissue, particularly ASM, using as a basis the attachment and detachment of cross-linker proteins at a continuum of cross-linker binding sites. The resulting mathematical model exhibits strain-induced fluidization, as well as several types of force recovery, at the same time suggesting a new mechanistic basis for the behavior. The model is validated by comparison to new data from experimental preparations of rat tracheal airway smooth muscle. Furthermore, experiments in noncontractile tissue show qualitatively similar behavior, suggesting support for the protein-filament theory as a biomechanical basis for the behavior.     

Could an increase in airway smooth muscle shortening velocity cause airway hyperresponsiveness?

Airway hyperresponsiveness (AHR) is a characteristic feature of asthma. It has been proposed that an increase in the shortening velocity of airway smooth muscle (ASM) could contribute to AHR. To address this possibility, we tested whether an increase in the isotonic shortening velocity of ASM is associated with an increase in the rate and total amount of shortening when ASM is subjected to an oscillating load, as occurs during breathing. Experiments were performed in vitro using 27 rat tracheal ASM strips supramaximally stimulated with methacholine. Isotonic velocity at 20% isometric force (Fiso) was measured, and then the load on the muscle was varied sinusoidally (0.33 ± 0.25 Fiso, 1.2 Hz) for 20 min, while muscle length was measured. A large amplitude oscillation was applied every 4 min to simulate a deep breath. We found that: 1) ASM strips with a higher isotonic velocity shortened more quickly during the force oscillations, both initially (P < 0.001) and after the simulated deep breaths (P = 0.002); 2) ASM strips with a higher isotonic velocity exhibited a greater total shortening during the force oscillation protocol (P < 0.005); and 3) the effect of an increase in isotonic velocity was at least comparable in magnitude to the effect of a proportional increase in ASM force-generating capacity. A cross-bridge model showed that an increase in the total amount of shortening with increased isotonic velocity could be explained by a change in either the cycling rate of phosphorylated cross bridges or the rate of myosin light chain phosphorylation. We conclude that, if asthma involves an increase in ASM velocity, this could be an important factor in the associated AHR.     

Myosin filament polymerization and depolymerization in a model of partial length adaptation in airway smooth muscle

Length adaptation in airway smooth muscle (ASM) is attributed to reorganization of the cytoskeleton, and in particular the contractile elements. However, a constantly changing lung volume with tidal breathing (hence changing ASM length) is likely to restrict full adaptation of ASM for force generation. There is likely to be continuous length adaptation of ASM between states of incomplete or partial length adaption. We propose a new model that assimilates findings on myosin filament polymerization/depolymerization, partial length adaptation, isometric force, and shortening velocity to describe this continuous length adaptation process. In this model, the ASM adapts to an optimal force-generating capacity in a repeating cycle of events. Initially the myosin filament, shortened by prior length changes, associates with two longer actin filaments. The actin filaments are located adjacent to the myosin filaments, such that all myosin heads overlap with actin to permit maximal cross-bridge cycling. Since in this model the actin filaments are usually longer than myosin filaments, the excess length of the actin filament is located randomly with respect to the myosin filament. Once activated, the myosin filament elongates by polymerization along the actin filaments, with the growth limited by the overlap of the actin filaments. During relaxation, the myosin filaments dissociate from the actin filaments, and then the cycle repeats. This process causes a gradual adaptation of force and instantaneous adaptation of shortening velocity. Good agreement is found between model simulations and the experimental data depicting the relationship between force development, myosin filament density, or shortening velocity and length.     

Effect of length oscillations on airway smooth muscle reactivity and cross-bridge cycling

Excessive airway narrowing due to airway smooth muscle (ASM) hyperconstriction is a major symptom in many respiratory diseases. In vitro imposition of length oscillations similar to those produced by tidal breathing on contracted ASM have shown to reduce muscle active forces, which is usually attributed to unconfirmed disruption of actomyosin cross-bridges. This research focuses on an in vitro investigation of the effect of mechanical oscillations on ASM reactivity and actomyosin cross-bridges. A computerized organ bath system was used to test maximally precontracted bovine ASM subjected to length oscillations at frequencies in the range of 10-100 Hz superimposed on tidal breathing oscillation. Using an immunofluorescence technique, two specific antibodies against the phospho-serine19 myosin light chain and the α-smooth muscle actin were used to analyze the colocalization between these two filaments. Data were processed using the plug-in "colocalization threshold" of ImageJ 1.43m software. The results demonstrate that both tidal and superimposed length oscillations reduce the active force in contracted ASM for a relatively long term and that the latter enhances the force reduction of the former. This reduction was also found to be frequency and time dependent. Additionally colocalization analysis indicates that length oscillations cause the detachment of the actomyosin connections and that this condition is sustained even after the cessation of the length oscillations.     

The role of a transient potassium current in a bursting neuron model

Presented here is a biophysical cell model which can exhibit low-frequency repetitive activity and bursting behavior. The model is developed from previous models (Av-Ron et al. 1991, 1993) for excitability, oscillations and bursting. A stepwise development of the present model shows the contribution of a transient potassium current (I _A ) to the overall dynamics. By changing a limited set of model parameters one can describe different firing patterns; oscillations with frequencies ranging from 2–200 Hz and a wide range of bursting behaviors in terms of the durations of bursting and quiescence, peak firing frequency and rate of change of the firing frequency.     

Temporal aspects of excitation-contraction coupling in airway smooth muscle.

In airway smooth muscle (ASM), ACh induces propagating intracellular Ca2+ concentration ([Ca2+]i) oscillations (5-30 Hz). We hypothesized that, in ASM, coupling of elevations and reductions in [Ca2+]i to force generation and relaxation (excitation-contraction coupling) is slower than ACh-induced [Ca2+]i oscillations, leading to stable force generation. When we used real-time confocal imaging, the delay between elevated [Ca2+]i and contraction in intact porcine ASM cells was found to be approximately 450 ms. In beta-escin-permeabilized ASM strips, photolytic release of caged Ca2+ resulted in force generation after approximately 800 ms. When calmodulin (CaM) was added, this delay was shortened to approximately 500 ms. In the presence of exogenous CaM and 100 microM Ca2+, photolytic release of caged ATP led to force generation after approximately 80 ms. These results indicated significant delays due to CaM mobilization and Ca2+-CaM activation of myosin light chain kinase but much shorter delays introduced by myosin light chain kinase-induced phosphorylation of the regulatory myosin light chain MLC20 and cross-bridge recruitment. This was confirmed by prior thiophosphorylation of MLC20, in which force generation occurred approximately 50 ms after photolytic release of caged ATP, approximating the delay introduced by cross-bridge recruitment alone. The time required to reach maximum steady-state force was >15 s. Rapid chelation of [Ca2+]i after photolytic release of caged diazo-2 resulted in relaxation after a delay of approximately 1.2 s and 50% reduction in force after approximately 57 s. We conclude that in ASM cells agonist-induced [Ca2+]i oscillations are temporally and spatially integrated during excitation-contraction coupling, resulting in stable force production.     

Cross-bridge model of muscle contraction. Quantitative analysis

We recently presented, in a qualitative manner, a cross-bridge model of muscle contraction which was based on a biochemical kinetic cycle for the actomyosin ATPase activity. This cross-bridge model consisted of two cross-bridge states detached from actin and two cross-bridge states attached to actin. In the present paper, we attempt to fit this model quantitatively to both biochemical and physiological data. We find that the resulting complete cross-bridge model is able to account reasonably well for both the isometric transient data observed when a muscle is subjected to a sudden change in length and for the relationship between the velocity of muscle contraction in vivo and the actomyosin ATPase activity in vitro. This model also illustrates the interrelationship between biochemical and physiological data necessary for the development of a complete cross-bridge model of muscle contraction.     

Smooth muscle stiffness variation due to external longitudinal oscillations

This research focuses on an in vitro investigation of the stiffness changes of contracted airway smooth muscles (ASM) subjected to external longitudinal oscillations. ASM tissues were dissected from excised pig tracheas and stimulated by a chemical stimulus (acetylcholine, 10(-3) M) to produce maximum contractions. The tissues were then systematically excited with external oscillations. Various frequencies, amplitudes and durations were used in the experiments to determine stiffness changes in response to these variations. Force changes were recorded to reflect the muscle stiffness changes. Two stiffness definitions were used to quantify the results, dynamic stiffness to reflect variations during oscillation and static stiffness to reflect the net effect of oscillation. Under isometric contractions, these two stiffnesses were determined before, during and after oscillations. Incorporating an empirical stiffness equation, a two-dimensional finite element model (FEM) was developed to generalize the tissue responses to oscillation. The main outcomes from this work are: the dynamic stiffness has the tendency to decrease as the frequency and/or amplitude of external oscillation increases; the static stiffness has the tendency of decreasing with an increase in the frequency and/or amplitude of excitation until it reaches almost a constant value for frequencies at and above 25 Hz. The difference in the behavior of the dynamic and static stiffness changes may be attributed to the effect of elasticity and mass inertia that are involved in the dynamic motion. The findings of this research are in agreement with the hypothesis that oscillations exert a direct action on the contractile processes by causing an increased rate of actin-myosin detachments.     

From sarcomere to cell: An efficient algorithm for linking mathematical models of muscle contraction

Two classes of mathematical framework have previously been developed to model active tension generation in contracting muscle. Cross-bridge models of muscle are biophysically based but computationally expensive to solve, and thus unsuitable for embedding in spatially distributed continuum representations. Fading memory models are computationally efficient but provide limited biophysical insight. In this study a novel computational method is proposed for coupling these two frameworks such that biophysical events can be determined and computational tractability maintained. Within the cross-bridge model, the functional forms of the distribution of cross-bridges, as a function of strain in each state, are approximated using the distribution moment approach. Using the variables of area, mean and standard deviation of each distribution, analytic expressions are developed to calculate the temporal dynamics of stiffness, tension and energy. A root finding method is employed to adjust the variables such that the temporal dynamics of the cross-bridge model match those of an equivalent fading memory model. The method is demonstrated for sinusoidal perturbations in length at two frequencies, with an approximate 30-fold increase in computational efficiency over a conventional technique for finding a solution to the cross-bridge model.     

Modeling residual force enhancement with generic cross-bridge models

The interaction of actin and myosin through cross-bridges explains much of muscle behavior. However, some properties of muscle, such as residual force enhancement, cannot be explained by current cross-bridge models. There is ongoing debate whether conceptual cross-bridge models, as pioneered by Huxley (A.F. Huxley, Muscle structure and theories of contraction, Prog. Biophys. Biophys. Chem. 7 (1957) 255) could, if suitably modified, fit experimental data showing residual force enhancement. Here we prove that there are only two ways to explain residual force enhancement with these ‘traditional’ cross-bridge models: the first requires cross-bridges to become stuck on actin (the stuck cross-bridge model) while the second requires that cross-bridges that are pulled off beyond a critical strain enter a ‘new’ unbound state that leads to a new force-producing cycle (the multi-cycle model). Stuck cross-bridge models cannot fit the velocity and stretch amplitude dependence of residual force enhancement, while the multi-cycle models can. The results of this theoretical analysis demonstrate that current kinetic models of cross-bridge action cannot explain the experimentally observed residual force enhancement. Either cross-bridges in the force-enhanced state follow a different kinetic cycle than cross-bridges in a ‘normal’ force state, or the assumptions underlying traditional cross-bridge models must be violated during experiments that show residual force enhancement.     

Modelling concentric contraction of muscle using an improved cross-bridge model.

Despite its overwhelming acceptance in muscle research, the cross-bridge theory does not account for all phenomena observed during muscular contractions. A phenomenon which has received much attention in the biomechanics literature, but has evaded convincing explanation and is not accounted for in the formulation of the classic cross-bridge theory, is the persistent aftereffects of muscular length changes on force production. For example, following muscle shortening, the isometric force of a muscle is depressed for a long time period ( > 5 s) compared to the corresponding isometric force following no length change. In the present study, the classic cross-bridge model was modified in two ways in an attempt to account for the force depressions following muscle shortening. First, the steady-state force depressions following shortening were described by a single scalar variable: the work performed by the muscle during shortening; and second, the dynamic, history-dependent cross-bridge properties were described using a fading memory function. The proposed model was developed and tested for shortening of the cat soleus at constant speeds ranging from 4 to 32 mm/s, for shortening at changing speeds, and for shortening of different magnitudes ranging from 2 to 10 mm. The history-dependent forces during shortening and the steady-state force depressions following shortening were well captured with the modified cross-bridge model. The present model contains two mathematically simple adaptations to the classic cross-bridge model, and is the first such model to account for the long-lasting force depressions following muscle shortening using a single scalar variable.     

Stiffness and force in activated frog skeletal muscle fibers.

Single fibers, isolated intact from frog skeletal muscles, were held firmly very near to each end by stiff metal clasps fastened to the tendons. The fibers were then placed horizontally between two steel hooks inserted in eyelets of the tendon clasps. One hook was attached to a capacitance gauge force transducer (resonance frequency up to approximately 50 kHz) and the other was attached to a moving-coil length changer. This allowed us to impose small, rapid releases (complete in less than 0.15 ms) and high frequency oscillations (up to 13 kHz) to one end of a resting or contracting fiber and measure the consequences at the other end with fast time resolution at 4 to 6 degrees C. The stiffness of short fibers (1.8-2.6 mm) was determined directly from the ratio of force to length variations produced by the length changer. The resonance frequency of short fibers was so high (approximately 40 kHz) that intrinsic oscillations were not detectably excited. The stiffness of long fibers, on the other hand, was calculated from measurement of the mechanical resonance frequency of a fiber. Using both short and long fibers, we measured the sinusoids of force at one end of a contracting fiber that were produced by relatively small sinusoidal length changes at the other end. The amplitudes of the sinusoidal length changes were small compared with the size of step changes that produce nonlinear force-extension relations. The sinusoids of force from long fibers changed amplitude and shifted phase with changes in oscillation frequency in a manner expected of a transmission line composed of mass, compliance, and viscosity, similar to that modelled by (Ford, L. E., A. F. Huxley, and R. M. Simmons, 1981, J. Physiol. (Lond.), 311:219-249). A rapid release during the plateau of tetanic tension in short fibers caused a fall in force and stiffness, a relative change in stiffness that putatively was much smaller than that of force. Our results are, for the most part, consistent with the cross-bridge model of force generation proposed by Huxley, A. F., and R. M. Simmons (1971, Nature (Lond.), 213:533-538). However, stiffness in short fibers developed markedly faster than force during the tetanus rise. Thus our findings show the presence of one or more noteworthy cross-bridge states at the onset and during the rise of active tension towards a plateau in that attachment apparently is followed by a relatively long delay before force generation occurs.(ABSTRACT TRUNCATED AT 400 WORDS)     

Information coding by ensembles of resonant neurons

In the present paper, we propose a novel neural procedure for signal processing and coding based on the subthreshold oscillations and resonance of the neural membrane potential that could be used by real neurons to perform frequency spectra analysis and information coding of incoming signals. Taking into account the biophysical properties of the neural membranes, we note that the subthreshold resonant behaviour they exhibit can be used to analyse incoming signals and represent them in the frequency domain. We study the reliability of the representation of signals depending on the biophysical parameters of the neurons, the fault-tolerance of this coding scheme and its robustness against noise and in the presence of spikes. The principal characteristics of our system are the use of the physical phenomenon of neural resonance (rarely considered in the literature for signal coding); it fits well with the biophysical parameters of most neurons that exhibit subthreshold oscillations; it is compatible with experimental data; and it can be easily integrated in a more general model of information processing and coding that includes communication between neurons based on spikes.     

Computer simulation of flagellar movement. III. Models incorporating cross-bridge kinetics.

A computer simulation procedure is used to analyze the generation of propagated bending waves by flagellar models in which active sliding is generated by a cycle of cross-bridge activity. Two types of cross-bridge cycle have been examined in detail. In both cycles, cross-bridge attachment is followed immediately by a configurational change in the cross-bridge, which transfers energy to a stretched elastic element and generates a shearing force between the filaments. In the first model, which has cross-bridge behavior close to current ideas about cross-bridge behavior in muscle, cross-bridge attachment is proportional to curvature of the flagellum and detachment is an exponential decay process. The configurational change is equivalent to an angular deviation of pi/5 radians. In the second type of cross-bridge cycle, cross-bridge attachment occurs rapidly when a critical curvature is reached, and detachment occurs when a critical curvature in the opposite direction is reached. With this cycle, an unrealistically large angular deviation of the cross-bridges, equivalent to 3.0 radians, is required to obtain bending waves of normal amplitude. Both models generate bending wave patterns similar to those obtained in earlier work. However, the behavior of the second type of cross-bridge model more closely matches the actual behavior of flagella under experimental conditions: the chemical turnover rate per beat cycle remains constant as the viscosity is increased, and reduction in the number of active cross-bridges can cause a reduction in beat frequency, with little change in amplitude or wavelength.     

Effects of length oscillation on the subsequent force development in swine tracheal smooth muscle.

It has been shown that deep inspiration (DI) taken before application of bronchoconstricting stimuli causes a reduction in the subsequent bronchoconstriction; a fast DI has a greater inhibitory effect than a slow DI. We hypothesize that periodic length changes imposed on a relaxed airway smooth muscle (ASM) would attenuate subsequent bronchoconstriction by disrupting the organization of the contractile apparatus, and this could be an important mechanism for the observed bronchoprotective effect of DI and tidal breathing. Length oscillations of different amplitude, frequency, and duration were applied to a relaxed muscle. The effects of such perturbations on force development were then assessed. Results show that oscillations reduce the subsequent force generation and that the magnitude of force reduction is proportional to amplitude and duration of the length oscillation. After the oscillation, isometric force recovered to the preoscillation level in a series of isometric contractions, and the rate of recovery was facilitated by frequent stimulation. The in vitro behavior of ASM found in this study could account for the observed temporary reduction in bronchoconstriction subsequent to a DI.     

描述肝细胞中两类不同特性钙离子浓度振荡

细胞内第二信使钙离子通常以浓度振荡的方式转导多种生理学信息,影响细胞分化、成熟和凋亡等各种生理过程。肝细胞实验中看到在一定浓度范围的激动剂刺激下,细胞质钙离子浓度的变化可呈现出很不相同的图像。例如顺脱羟肾上腺素刺激下,可出现简单的周期振荡,振荡频率随激动剂浓度的不同有变化;在腺苷三磷酸刺激下,随激动剂浓度从低到高,胞质钙离子浓度的变化可以从开始出现简单振荡,到形成复杂的多峰间歇振荡。肝细胞中钙离子浓度振荡的峰值多在500nmol/L到800nmol/L范围。给出一个四为量数学模型的改进形式,可以模拟肝细胞中钙离子浓度从简单振荡向复杂振荡的变化。数值计算给出与实验结果比较一致的振荡波形和振幅。     

Role of MgATP and MgADP in the cross-bridge kinetics in chemically skinned rabbit psoas fibers. Study of a fast exponential process (C)

The role of the substrate (MgATP) and product (MgADP) molecules in cross-bridge kinetics is investigated by small amplitude length oscillations (peak to peak: 3 nm/cross-bridge) and by following amplitude change and phase shift in tension time courses. The range of discrete frequencies used for this investigation is 0.25-250 Hz, which corresponds to 0.6-600 ms in time domain. This report investigates the identity of the high frequency exponential advance (process C), which is equivalent to "phase 2" of step analysis. The experiments are performed in maximally activated (pCa 4.5-5.0) single fibers from chemically skinned rabbit psoas fibers at 20 degrees C and at the ionic strength 195 mM. The rate constant 2 pi c deduced from process (C) increases and saturates hyperbolically with an increase in MgATP concentration, whereas the same rate constant decreases monotonically with an increase in MgADP concentration. The effects of MgATP and MgADP are opposite in all respects we have studied. These observations are consistent with a cross-bridge scheme in which MgATP and MgADP are in rapid equilibria with rigorlike cross-bridges, and they compete for the substrate site on myosin heads. From our measurements, the association constants are found to be 1.4 mM-1 for MgATP and 2.8 mM-1 for MgADP. We further deduced that the composite second order rate constant of MgATP binding to cross-bridges and subsequent isomerization/dissociation reaction to be 0.57 x 10(6)M-1s-1.     

A Cross-Bridge Cycle with Two Tension-Generating Steps Simulates Skeletal Muscle Mechanics

We examined whether cross-bridge cycle models with one or two tension-generating steps can account for the force-velocity relation of and tension response to length steps of frog skeletal muscle. Transition-state theory defined the strain dependence of the rate constants. The filament stiffness was non-Hookean. Models were refined against experimental data by simulated annealing and downhill simplex runs. Models with one tension-generating step were rejected, as they had a low efficiency and fitted the experimental data relatively poorly. The best model with two tension-generating steps (stroke distances 5.6 and 4.6 nm) and a cross-bridge stiffness of 1.7 pN/nm gave a good account of the experimental data. The two tensing steps allow an efficiency of up to 38% during shortening. In an isometric contraction, 54.7% of the attached heads were in a pre-tension-generating state, 44.5% of the attached heads had undergone the first tension-generating step, and only 0.8% had undergone both tension-generating steps; they bore 34%, 64%, and 2%, respectively, of the isometric tension. During slow shortening, the second tensing step made a greater contribution. During lengthening, up to 93% of the attached heads were in a pre-tension-generating state yet bore elevated tension by being dragged to high strains before detaching.     

Thixotropy and Rheopexy of Muscle Fibers Probed Using Sinusoidal Oscillations

Length changes of muscle fibers have previously been shown to result in a temporary reduction in fiber stiffness that is referred to as thixotropy. Understanding the mechanism of this thixotropy is important to our understanding of muscle function since there are many instances in which muscle is subjected to repeated patterns of lengthening and shortening. By applying sinusoidal length changes to one end of single permeabilized muscle fibers and measuring the force response at the opposite end, we studied the history-dependent stiffness of both relaxed and activated muscle fibers. For length change oscillations greater than 1 Hz, we observed thixotropic behavior of activated fibers. Treatment of these fibers with EDTA and blebbistatin, which inhibits myosin-actin interactions, quashed this effect, suggesting that the mechanism of muscle fiber thixotropy is cross-bridge dependent. We modeled a half-sarcomere experiencing sinusoidal length changes, and our simulations suggest that thixotropy could arise from force-dependent cross-bridge kinetics. Surprisingly, we also observed that, for length change oscillations less than 1 Hz, the muscle fiber exhibited rheopexy. In other words, the stiffness of the fiber increased in response to the length changes. Blebbistatin and EDTA did not disrupt the rheopectic behavior, suggesting that a non-cross-bridge mechanism contributes to this phenomenon.     

A Cross-Bridge Cycle with Two Tension-Generating Steps Simulates Skeletal Muscle Mechanics

We examined whether cross-bridge cycle models with one or two tension-generating steps can account for the force-velocity relation of and tension response to length steps of frog skeletal muscle. Transition-state theory defined the strain dependence of the rate constants. The filament stiffness was non-Hookean. Models were refined against experimental data by simulated annealing and downhill simplex runs. Models with one tension-generating step were rejected, as they had a low efficiency and fitted the experimental data relatively poorly. The best model with two tension-generating steps (stroke distances 5.6 and 4.6 nm) and a cross-bridge stiffness of 1.7 pN/nm gave a good account of the experimental data. The two tensing steps allow an efficiency of up to 38% during shortening. In an isometric contraction, 54.7% of the attached heads were in a pre-tension-generating state, 44.5% of the attached heads had undergone the first tension-generating step, and only 0.8% had undergone both tension-generating steps; they bore 34%, 64%, and 2%, respectively, of the isometric tension. During slow shortening, the second tensing step made a greater contribution. During lengthening, up to 93% of the attached heads were in a pre-tension-generating state yet bore elevated tension by being dragged to high strains before detaching.