JN Information on EB 2010
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 94: 3259-3277, 2005. First published June 8, 2005; doi:10.1152/jn.00481.2005
0022-3077/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
94/5/3259    most recent
00481.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhurov, Y.
Right arrow Articles by Brezina, V.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zhurov, Y.
Right arrow Articles by Brezina, V.

Temperature Compensation of Neuromuscular Modulation in Aplysia

Yuriy Zhurov and Vladimir Brezina

Fishberg Department of Neuroscience, Mount Sinai School of Medicine, New York, New York

Submitted 10 May 2005; accepted in final form 6 June 2005


 ABSTRACT
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Physiological systems that must operate over a range of temperatures often incorporate temperature-compensatory mechanisms to maintain their output within a relatively narrow, functional range of values. We analyze here an example in the accessory radula closer (ARC) neuromuscular system, a representative part of the feeding neuromusculature of the sea slug Aplysia. The ARC muscle's two motor neurons, B15 and B16, release, in addition to ACh that contracts the muscle, modulatory peptide cotransmitters that, through a complex network of effects in the muscle, shape the ACh-induced contractions. It is believed that this modulation is critical in optimizing the performance of the muscle for successful, efficient feeding behavior. However, previous work has shown that the release of the modulatory peptides from the motor neurons decreases dramatically with increasing temperature. From 15 to 25°C, for example, release decreases 20-fold. Yet Aplysia live and feed successfully not only at 15°C, but at 25°C and probably at higher temperatures. Here, working with reduced B15/B16–ARC preparations in vitro as well as a mathematical model of the system, we have found a resolution of this apparent paradox. Although modulator release decreases 20-fold when the temperature is raised from 15 to 25°C, the observed modulation of contraction shape does not decrease at all. Two mechanisms are responsible. First, further downstream within the modulatory network, the modulatory effects themselves—experimentally dissected by exogenous modulator application—have temperature dependencies opposite to that of modulator release, increasing with temperature. Second, the saturating curvature of the dose–response relations within the network diminishes the downstream impact of the decrease of modulator release. Thus two quite distinct mechanisms, one depending on the characteristics of the individual components of the network and the other emerging from the network's structure, combine to compensate for temperature changes to maintain the output of this physiological system.


 INTRODUCTION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The rates of individual biochemical reactions and elementary physiological processes vary intrinsically with temperature. Yet the larger physiological systems that are built from these processes must often maintain their output within a relatively narrow, functional range of values. Animals that are active over a range of environmental temperatures must therefore compensate, if they do not regulate their internal temperature, by building temperature-compensatory mechanisms into the structure of their physiological systems. We analyze here an example in the accessory radula closer (ARC) neuromuscular system of the sea slug Aplysia. This system is tuned by a complex network of neuromodulatory effects that, individually, vary greatly with temperature; yet the overall output of the system does not.

The ARC muscle is one of the muscles of the buccal mass, a complex structure that produces rhythmic, cyclical movements of the animal's food-grasping organ, the radula, in several types of consummatory feeding behavior such as biting, swallowing, and rejection of unsuitable food (Kupfermann 1974Go). The structure of the ARC system is summarized in Fig. 1. The muscle is innervated by two motor neurons, B15 and B16, which, when they fire, release their classical transmitter acetylcholine (ACh) to depolarize and contract the muscle (Cohen et al. 1978Go; Fig. 1, left). In addition, however, the two motor neurons release several families of modulatory peptide cotransmitters that shape the basal ACh-induced contractions (Fig. 1, right; for review and complete references see Brezina et al. 2003aGo, 2005Go; Hooper et al. 1999Go; Kupfermann et al. 1997Go; Weiss et al. 1993Go). Most important are the small cardioactive peptides (SCPs), released from motor neuron B15 (e.g., Cropper et al. 1987aGo; Lloyd et al. 1984Go; Vilim et al. 1996bGo), and the myomodulins (MMs), released from motor neuron B16 (Brezina et al. 1995Go; Cropper et al. 1987bGo, 1991Go; Vilim et al. 2000Go). The SCPs and MMs shape contractions through three main effects on the muscle: 1) they potentiate the contractions by enhancing the muscle's depolarization-activated Ca current that supplies Ca2+ essential for contraction (e.g., Brezina and Weiss 1995Go; Brezina et al. 1994aGo, 1995Go); 2) they depress the contractions by activating in the muscle a K current that opposes the ACh-induced depolarization, activation of the Ca current, and Ca2+ influx (Brezina et al. 1994bGo, 1995Go; Orekhova et al. 2003Go); and 3) they increase the relaxation rate of the contractions probably by modulating the muscle's contractile machinery (Heierhorst et al. 1994Go, 1995Go; Probst et al. 1994Go). The balance of the competing potentiating and depressing effects determines the net modulation of contraction size that, together with the modulation of the relaxation rate, then produces the final shape of the contraction. Similar networks of modulatory effects shape the contractions of the other muscles of the buccal mass (e.g., Church et al. 1993Go; Evans et al. 1999Go; Hurwitz et al. 2000Go; Keating and Lloyd 1999Go).



View larger version (35K):
[in this window]
[in a new window]
 
FIG. 1. Schematic summary of the structure of the B15/B16–ARC neuromuscular system. Variables X(t) are the variables of the dynamic model of the system from Brezina et al. (2003a,b) that we use here in Figs. 2 and 10–12. For details see INTRODUCTION, METHODS, and Reproduction of the observed temperature dependency in a mathematical model of the system in RESULTS.

 
In extensive work over many years, in the real system as well as with a series of increasingly realistic mathematical models, we have analyzed how the modulatory shaping of contractions dynamically tunes the buccal musculature to meet different behavioral demands, to perform efficiently each of the different types of feeding behavior. The modulation is critical for function: without the modulatory tuning, the constraint of the fixed neuromuscular properties would make many sequences of behaviors that the animal might be required to perform inefficient or even completely dysfunctional (Brezina et al. 1996Go, 2000bGo, 2003aGo,bGo, 2005Go).

There is, however, a serious potential difficulty with this scenario. Vilim et al. (1996a)Go found that the release of the modulatory peptide cotransmitters from the motor neurons is extremely sensitive to temperature. Figure 2 replots the data of Vilim et al. (1996a)Go for the release of the SCPs from the ARC motor neuron B15, measured by direct radioimmunoassay of the amount of released peptide (each point is the mean ± SE, smaller than the mean symbol size, of four preparations). Also plotted is the simultaneous release of the buccalins, peptides of another family that are costored in the same dense-core vesicles with the SCPs in B15 (as with the MMs in B16; "Buc" in Fig. 1) and coreleased by the same release process (Vilim et al. 1996bGo, 2000Go). Release of both peptide families declines dramatically with increasing temperature: an increase of 10°C reduces release approximately 20-fold (see Fig. 2 legend).



View larger version (16K):
[in this window]
[in a new window]
 
FIG. 2. Temperature dependency of modulatory peptide cotransmitter release from motor neuron B15. Points replot the experimental data from Fig. 6B of Vilim et al. (1996a). Vilim et al. used radioimmunoassay to directly measure the total amounts of small cardioactive peptide (SCP; filled points) and buccalin (open points) appearing in the perfusate of the accessory radula closer (ARC) muscle when motor neuron B15 was intracellularly stimulated to fire in a physiologically realistic bursting pattern, with intraburst firing frequency fintra = 12 Hz, burst duration dintra = 3.5 s, and interburst interval dinter = 3.5 s, for a total length of firing L = 10 min, at 15, 17, and 19°C. Each point is a mean of n = 4 preparations, with SE smaller than the mean symbol size, however, normalized within each preparation as follows. Because different preparations can release quite different absolute amounts of peptides (Brezina et al. 2000a; Vilim et al. 1996a,b), Vilim et al. expressed all of the data as percentages of the grand total amount of peptide released in the entire series of repetitions of the motor neuron firing pattern at the three temperatures within each preparation. Here we have simply rescaled the plot of Vilim et al. to absolute peptide amounts by multiplying each percentage by the mean grand total amount of peptide released in the 4 preparations (F. S. Vilim, personal communication). Thus the small errors seen here reflect just the variance of the effect of temperature; the much larger variance between preparations has been normalized out. Continuous curves are best fits of our model of the release process (Brezina et al. 2000a, 2003a). To model the temperature dependency, we use the standard concept of the temperature coefficient Q10, relating the value of variable X at temperature T to its value at 15°C by

(1)
With such a coefficient, Q10,r, incorporated into the release equations to scale the instantaneous release r (see Brezina et al. 2003a, Eq. 7 in Fig. 11 legend), analysis of the model (see Brezina et al. 2000a, analysis of Model II) shows that R, the total amount of peptide released by the entire block of firing of length L, as measured in these experiments by Vilim et al., is described by the equation (the equivalent of Eq. 15 of Brezina et al. 2000a)

{J004815p3259eq2}(2)
where S0 is the initial size of the releasable peptide pool, p{infty} = <f>kp+/kp and {tau}p = 1/kp are respectively the steady-state value and the time constant of a slow reaction governing the probability of release, {Phi} = D1–y is the pattern dependency of the release, D = dintra/(dintra + dinter) is the duty cycle of the firing pattern, and <f> = fintraD is the mean firing frequency, with the parameter values [determined from other data of Vilim et al. (1996a,b) by Brezina et al. (2000a, 2003a)] kp+ = 4.0 x 10–10, kp = 3.4 x 10–3 s–1, y = 3, and fintra, dintra, and dinter as above. This leaves just Q10,r, S0,SCP, and S0,Buc as free parameters. Fitting Eq. 2 simultaneously to both the SCP and buccalin data yielded Q10,r = 5.05 x 10–2 ({approx}1/20), S0,SCP = 628 fmol, and S0,Buc = 245 fmol, and the curves shown. [Value of S0,SCP found here is 16% larger, and the value of S0,Buc is 24% larger, than the consensus values of 541 and 198 fmol that were determined by Brezina et al. (2000a), presumably again because the different preparations used released different absolute amounts of peptides.]

 
Figure 2 shows that release of the modulatory peptides is strong at 15°C, a physiological temperature for Aplysia californica, the species studied in all of this work. However, Aplysia do not live only at 15°C. At different times of the year and in different local environments—in deep water as compared with, for example, shallow, sun-heated tide pools—A. californica will encounter water temperatures ranging from well below 15°C to as high as 25°C (Kupfermann and Carew 1974Go). [Other Aplysia species, such as A. fasciata in the Mediterranean Sea or A. dactylomela in the Caribbean, live and feed at up to at least 29°C (Carefoot 1987Go; Gev et al. 1984Go), but the modulatory tuning of the feeding musculature has not been studied in these species.] Aplysia are, of course, poikilotherms that cannot maintain an internal temperature that is appreciably different from that of the surrounding water.

According to Fig. 2, there is little release of the peptide modulators at, say, 25°C.1 Yet Aplysia live and feed at this temperature. Does this mean that the modulatory tuning is, after all, not important for functional, efficient feeding? In this paper we provide a resolution of this difficulty. Even though the release of the modulators declines dramatically with increasing temperature, the overall strength of the modulatory tuning does not decline. This is because the modulatory system incorporates mechanisms, individual effects with the opposite temperature dependency as well as more general mechanisms that are inherent in the network structure of the system, that maintain its output throughout the entire temperature range.


 METHODS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Experiments

PREPARATION.  The preparation was the standard preparation used for recording motor neuron-elicited contractions of the ARC muscle (Cohen et al. 1978Go; Weiss et al. 1978Go; for recent use see, e.g., Horn et al. 2004Go; Orekhova et al. 2003Go). Briefly, the preparation consisted of the bilateral buccal ganglia, the two ARC muscles, and the connecting buccal nerves 3 through which the buccal motor neurons B15 and B16 innervate the ARC muscle. The cerebral ganglion, connected to the buccal ganglia by the cerebral–buccal connectives, was also retained. The buccal ganglia (but not the cerebral ganglion) were desheathed. One ARC muscle was pinned out in a separate subchamber and connected to an isotonic transducer (Model 60–3000, Harvard Apparatus, Holliston, MA) to measure the length of the muscle with a light counterbalancing load. The ipsilateral motor neurons B15 and B16 were impaled with standard intracellular microelectrodes and connected to an intracellular amplifier (Axoclamp 2A/B, Axon Instruments, Union City, CA). Under the control of external timers/stimulators (Pulsemaster A300, World Precision Instruments, Sarasota, FL, and Grass S48/S88, Astro-Med, West Warwick, RI), the motor neurons were stimulated with brief current injections to fire spikes in the desired pattern. All signals were sampled and recorded simultaneously by a computer using Digidata 1322A data-acquisition hardware and pCLAMP 8/9 software (Axon Instruments). When the experiment on the first ARC muscle was finished, the ARC muscle on the opposite side of the preparation was also used, whenever possible, for an independent experiment or a matched pair (see RESULTS).

EXPERIMENTAL DESIGN.  Three basic experimental paradigms were used.

1) Single motor neuron self-modulation.  Only one of the two ARC motor neurons, either B15 or B16, was fired in a relatively intense bursting pattern, generally at 10–15 Hz for 2 s every 5 s for 3 min, calculated to release the peptide modulators (see Fig. 3).



View larger version (53K):
[in this window]
[in a new window]
 
FIG. 3. Single motor neuron self-modulation: representative experiment. Motor neuron B16 was fired at 13 Hz for 2 s every 5 s for 3 min, at 15°C (A) and about 22°C (B). Continuous waveform in each case is the instantaneous length of the ARC muscle relative to the baseline length at the beginning of the motor neuron firing (axis at left; decreasing length, thus increasing contraction amplitude, is plotted upward); large black points are measurements of the relaxation rate of the contraction elicited by each motor neuron burst (axis at right; see Fig. 1 and METHODS). Left and right ARC muscles of the same preparation.

 
2) Two motor neuron cross-modulation.  Both motor neurons were fired, one with brief bursts at regular infrequent intervals to elicit "monitoring" contractions of moderate size, the other with a short block of a more intense bursting pattern, interposed between the monitoring contractions once only during the experiment, to release the modulators (see Fig. 4). The monitoring neuron, either B15 or B16, was generally fired for 1 s every 100 s, B15 at 20 or 25 Hz, and B16 at 25 Hz. Such a pattern was likely to have elicited little or no modulator release from the monitoring neuron itself (Brezina et al. 2003bGo; Cropper et al. 1990cGo; Vilim et al. 1996aGo). The modulator-releasing neuron was fired for 4 s every 10 s for 60 s, B15 at 20 Hz, and B16 at 20 or 25 Hz.



View larger version (24K):
[in this window]
[in a new window]
 
FIG. 4. Two motor neuron cross-modulation: representative experiment. Motor neuron B15 was fired at 25 Hz for 1 s every 100 s to elicit monitoring contractions, motor neuron B16 at 25 Hz for 4 s every 10 s for 60 s to release modulators, at 15°C (A) and 25°C (B). Continuous waveform in each of A1 and B1 is the instantaneous length of the ARC muscle (axis at left); large black points are measurements of the relaxation rate of each contraction (axis at right). A2 and B2 expand the contractions 13 from A1 and B1. Left and right ARC muscles of the same preparation.

 
3) Exogenous modulator application.  Motor neuron B15 was fired in the monitoring pattern, generally for 1 s every 100 s at 20 Hz, while exogenous SCPB was applied to the muscle subchamber. Increasing amounts of SCPB, calculated to produce a series of increasing known concentrations in the subchamber, were added cumulatively without wash (Fig. 7). Each concentration was left to work for at least five contractions (>8 min) before the parameters of the contraction were measured. Alternatively, a single concentration of SCPB, 10–7 M, was applied and the time course of the development of its effects was followed for at least ten contractions (>16 min) (Fig. 9, A and B). In some of these experiments the SCPB was then washed out, by perfusing the muscle subchamber (volume {cong} 7 ml) steadily at 1.25 ml/min, while the reversal of the effects was followed with continuing monitoring contractions (Fig. 9C). When the effects had reversed completely (after ≥30 min but sometimes up to 1 h of wash), the application of SCPB was repeated at a different temperature. For additional details see RESULTS and figure legends.



View larger version (23K):
[in this window]
[in a new window]
 
FIG. 7. Modulation by exogenous modulator: dose–response relations for exogenous SCPB. Motor neuron B15 was fired at 20 Hz for 1 s every 100 s to elicit monitoring contractions while increasing concentrations of SCPB were applied to the ARC muscle. Concentrations were applied cumulatively, without intervening wash; each concentration was allowed to work for ≥8 min before the parameters of the final contraction were measured. Five preparations were used, with one of the muscles in each preparation tested at 15°C and the other at 25°C, but for statistical purposes they are here treated as independent experiments. A: representative contractions at 15 and 25°C, in the two muscles of the same preparation. Shown are the final contractions at 10–9, 10–8, 10–7, 10–6, and 10–5 M SCPB ("–9"..."–5"), all superimposed on the same control contraction for comparison. For simplicity, only the motor neuron B15 burst that elicited the control contraction, at 25°C, is shown in each case. B: mean ± SE of the peak contraction amplitude, normalized to the peak amplitude of the control contraction in each muscle, as a function of the SCPB concentration, at 15°C (blue; n = 5) and 25°C (red; n = 5). C: mean ± SE of the relaxation rate as a function of the SCPB concentration, at 15°C (solid blue; n = 5) and 25°C (red; n = 5). Dashed blue curve is the 15°C dose–response relation projected into the range spanned by the 25°C dose–response relation and multiplied by 5. For further explanation see Modulation of the relaxation rate by exogenous modulator is intrinsically stronger at higher temperature in RESULTS.

 


View larger version (23K):
[in this window]
[in a new window]
 
FIG. 9. Modulation by exogenous modulator: time course of the modulation of relaxation rate by 10–7 M SCPB. A: unnormalized time course of the development of the modulation. Means ± SE of the absolute relaxation rates of successive monitoring contractions elicited by motor neuron B15 before and during the application of SCPB (gray block), at 15°C (blue, n = 10) and 25°C (red, n = 10). Labels "Fig. 8, control value" and "Fig. 8, modulated value" indicate where the measurements were taken for Fig. 8. Solid black curves are the best single-exponential fits, constrained in each case to have the time constant {tau} = 179 s obtained from B, to the data at the two temperatures, yielding the magnitudes of the fully developed modulation, 3.24-fold larger at 25 than at 15°C. B: normalized time course of the development of the modulation. Same data as in A but first normalized to between 0 and 100%, the smallest and largest values measured within each experiment, as in Fig. 6. Tested with 2-way ANOVA, the overall difference between the 15 and 25°C conditions was not significant (P > 0.05). Solid black curve is the best single-exponential fit simultaneously to the data at both temperatures, with {tau} = 179 s. C: normalized time course of the reversal of the modulation on washout of the SCPB. Means ± SE of the relaxation rates of monitoring contractions continued after the end of those measured in A and B in a subset of those experiments (n = 3 each at 15 and 25°C), normalized as in B. Statistical significance was tested with 2-way ANOVA followed by pairwise multiple comparisons using the Holm–Sidak test. Overall difference between the 15 and 25°C conditions was highly significant (P < 0.001); ***P < 0.001, **P < 0.01, and *P < 0.05, for the difference between the 15 and 25°C means at a particular time point. Solid black curves are the best single-exponential fits, over the range of times indicated, to the data at the two temperatures, with {tau} = 3,995 s at 15°C and {tau} = 462 s at 25°C.

 
TEMPERATURE.  Temperature was controlled with a Peltier plate under the preparation, and continuously monitored with an IT-23 thermocouple microprobe (Physitemp Instruments, Clifton, NJ) inserted in the muscle subchamber and connected to a digital thermometer (BAT-12, Physitemp Instruments). Because of the relatively large volume of the muscle subchamber, the correspondingly large volume of any perfusing flow, and the long duration of many experiments (about 1 h at each temperature), it was not always possible to maintain the temperature completely steady throughout the entire experiment. However, the temperature generally varied by no more than 1°C. Except in a few experiments done at ambient room temperature (about 22°C; Fig. 3), the temperature was controlled so that the 1°C range included, throughout the experiment, either 15 or 25°C: these are the experiments referred to as having been done at "15°C" or "25°C." When comparing the two temperatures sequentially in the same muscle or in the two muscles of the same preparation, the order in which the temperatures were tested was varied in different preparations; similar results were obtained regardless of the order.

MEASUREMENT OF CONTRACTION PARAMETERS.  The peak amplitude of each contraction was measured relative to the baseline muscle length at the beginning of the motor neuron burst that elicited the contraction, or, in the case of closely spaced bursts that did not allow the muscle to relax fully to the baseline (Fig. 3), the baseline length at the beginning of the first burst. The relaxation rate of each contraction was then evaluated, with a custom-written program in Mathematica (Wolfram Research, Champaign, IL), by fitting a single exponential to the relaxation phase of the contraction from 90 to 33% of the peak amplitude (see Fig. 1, bottom right), or, with the closely spaced contractions, to the beginning of the next motor neuron burst if that came earlier. Over this upper middle portion of the relaxation phase, most contractions were in fact well fitted by single exponentials. The relaxation rate was then the reciprocal of the time constant of the exponential.

Modeling

The dynamic model of the B15/B16–ARC neuromuscular system and its modulation was taken from Brezina et al. (2003aGo,b) and analytically solved or numerically integrated in Mathematica. Specific changes in parameters, dictated by the experimental results, were made as described in the figure legends. The relaxation rate of the simulated contractions was measured in the same way as that of the real contractions.


 RESULTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Single motor neuron self-modulation is not weaker at higher temperature

In the B15/B16–ARC neuromuscular system in vitro, the simplest way to observe the effects of the released modulators is with the single motor neuron self-modulation paradigm, in which just one motor neuron is stimulated both to release the modulators and to elicit the contraction that they modulate. Figure 3 shows a representative experiment of this kind. Motor neuron B16 was stimulated to fire in a repetitive bursting pattern—similar to that which it in fact fires in the intact feeding animal (Cropper et al. 1990aGo; Horn et al. 2004Go; see Fig. 12)—that elicited a sequence of contractions. These contractions progressively changed in shape, in amplitude (continuous waveform; axis at left) and in relaxation rate (large black points; axis at right). Previous studies have strongly suggested that much of the change in amplitude and especially relaxation rate arises from the progressive development of the effects of the modulators that such a pattern of firing also releases (e.g., Brezina et al. 2003bGo; Whim and Lloyd 1990Go).



View larger version (41K):
[in this window]
[in a new window]
 
FIG. 12. Temperature compensation during realistic feeding behavior, simulated with the model. Model was run at 15°C (A) or 25°C (B). Two traces at the top of each of A and B, identical at the two temperatures, are the inputs to the model, the instantaneous firing frequencies of the motor neurons B15 and B16 extracted by Brezina et al. (2005) from the electrical activity recorded with an electrode chronically implanted in the ARC muscle by Horn et al. (2004) during an approximately 2.5-h-long meal, consisting of 749 cycles of the feeding behavior, in an intact, freely feeding animal. Traces below show the resulting trajectories of the most relevant downstream variables of the model (see Fig. 1): the SCP and MM concentrations CSCP and CMM, the modulation of the relaxation rate R, and the modulated contraction amplitude c. C: mean ± SD of the absolute relaxation rate measured from each of the 749 cycles of the meal, at 15°C (blue) and 25°C (red), with the full modulation (filled circles) or with the modulation of the relaxation rate disabled [e.g., with R(t) = 0; empty circles]. Intrinsic temperature dependencies of MM release and of the modulation of the relaxation rate by MM were modeled to be the same as for SCP; all other aspects of the model not already given in legends of Figs. 2, 10, and 11 were as described by Brezina et al. (2003a,b).

 
Because modulator release declines with increasing temperature (Fig. 2), we initially expected that the observed modulation, too, would be weaker at higher temperature. We were therefore surprised to find that this was not the case. Figure 3 shows the same pattern of firing repeated, in the two ARC muscles of the same preparation, at 15°C (A) and at approximately 22°C (B). At the higher temperature, the progressive changes in shape were, if anything, more pronounced. The same result was obtained in two other experiments with motor neuron B16 as well as two experiments with motor neuron B15.

We did not analyze these results in any detail because, although the single motor neuron paradigm is simple, its results are difficult to interpret. As the motor neuron firing continues, all of the system components shown in Fig. 1—modulator release, modulator concentrations, the various modulator effects, as well as processes that shape the basal contraction such as presynaptic facilitation and potentiation of ACh release (Cohen et al. 1978Go) and postsynaptic summation of the contraction waveform itself—continue to evolve at different rates (Brezina et al. 2003aGo), in continuing interaction with each other. Each feature of the contraction shape thus has multiple interpretations.

Two motor neuron cross-modulation is not weaker at higher temperature

More readily interpretable are the results of the two motor neuron cross-modulation paradigm, first introduced in this system by Whim and Lloyd (1990)Go. One motor neuron is stimulated to fire brief bursts at regular infrequent intervals, which release little or no modulator themselves (Brezina et al. 2003bGo; Cropper et al. 1990cGo; Vilim et al. 1996aGo) but elicit regular "monitoring" contractions. The other motor neuron is then fired in a relatively short, intense bursting pattern to release the modulators. This pattern inevitably elicits a contraction, too, but afterward the contraction rapidly relaxes while the effects of the released modulators, some of which are very long lasting (Brezina et al. 2003aGo; see following text), persist to modulate the continuing monitoring contractions.2 Figure 4A1 shows a representative experiment. After the modulator-releasing firing of motor neuron B16, the monitoring contractions elicited by motor neuron B15 were increased in both amplitude (continuous waveform; axis at left) and relaxation rate (large black points; axis at right), and remained so for many minutes. From what is known about the system (see INTRODUCTION), both effects were presumably caused by the MMs released from motor neuron B16. When, in other experiments, the roles of the monitoring and modulator-releasing motor neurons were reversed, similar increases in contraction amplitude and relaxation rate were observed, in that case presumably caused by the SCPs released from motor neuron B15. Previous work has shown that the SCPs and MMs increase the relaxation rate in an identical fashion and have identical potentiating effects on contraction amplitude (Fig. 1; Brezina et al. 1996Go, 2003aGo).

Figure 4, A and B shows the same experiment repeated, in the two ARC muscles of the same preparation, at 15 and 25°C. Figure 4, A2 and B2 compares in detail three contractions expanded from A1 and B1, in particular a control contraction before the modulator-releasing firing (contraction 1) with the modulated contraction about 2 min after the end of the modulator-releasing firing (contraction 2), at each temperature. Clearly, the observed modulation was not weaker at the higher temperature. The modulation of the relaxation rate, especially, was significantly stronger.

Quantification of the modulation of relaxation rate

In our analysis of the results of the two motor neuron cross-modulation experiments, we focused primarily on the modulation of the relaxation rate and only secondarily on that of contraction amplitude. Even in the cross-modulation paradigm, the modulation of contraction amplitude is still a composite process with multiple interpretations. On the potentiating effect of the SCPs and MMs, mediated by enhancement of Ca current in the muscle, is superimposed a simultaneous depressing effect, mediated by activation of a K current (INTRODUCTION, Fig. 1). Furthermore, together with the SCPs and MMs both motor neurons B15 and B16 corelease the buccalins, which depress contractions (Cropper et al. 1988Go, 1990bGo; Orekhova et al. 2003Go; Vilim et al. 1994Go). Thus no net change in contraction amplitude, for example, could mean no modulator release, or conversely robust modulator release with the potentiating effect of the released SCPs, or MMs, balanced by an equally strong depressing effect of the SCPs or MMs, or of the coreleased buccalins.3 To fully understand how the modulation of contraction amplitude is affected by temperature, these various effects would have to be studied in isolation, each with its own paradigm in a different type of reduced preparation (see, e.g., Brezina et al. 1994aGo,bGo, 1995Go, 2003aGo; Orekhova et al. 2003Go).

In contrast, the modulation of the relaxation rate is, as far as is known, a unitary process, which furthermore is activated, in an identical fashion, just by the SCPs and MMs: the buccalins do not modulate the relaxation rate (Cropper et al. 1988Go; Vilim et al. 1994Go). Analysis of the relaxation rate, however, has its own difficulty. In different preparations, the initial, "unmodulated" relaxation rate of the ARC muscle—the control relaxation rate against which the modulation is then evaluated—can be very different (Brezina et al. 1995Go; see following text and Figs. 5 and 8). We attempted to manipulate the control relaxation rate by altering the load on the muscle, with limited success. Where possible we also made matched pairs of measurements, at 15 and 25°C, in the two muscles of the same preparation or successively in the same muscle. This reduced the variability somewhat, but considerable variability still remained.



View larger version (23K):
[in this window]
[in a new window]
 
FIG. 5. Two motor neuron cross-modulation: summary of measurements of the modulation of relaxation rate. A: modulated relaxation rate, of the monitoring contraction about 2 min after the end of the modulator-releasing firing (i.e., contraction 2 in Fig. 4), is plotted against the control relaxation rate, of the monitoring contraction before the modulator-releasing firing (contraction 1 in Fig. 4). Motor neuron B16 was fired to elicit the monitoring contractions and motor neuron B15 to release the modulators, at 15°C (filled blue circles, n = 15) or 25°C (red circles, n = 5), or, conversely, B15 was fired to elicit the monitoring contractions and B16 to release the modulators, at 15°C (blue squares, n = 8) or 25°C (red squares, n = 6). Green lines join matched pairs of measurements, at 15 and 25°C, in the two ARC muscles of the same preparation. Open blue circles are measurements from experiments in which motor neuron B16 was fired to elicit the monitoring contractions and motor neuron B15 to release the modulators, at 15°C, but in the presence of 10–5 M exogenous SCPB (n = 6). B: same measurements of control and modulated relaxation rate as in A plotted along one dimension for statistical comparison. Same symbols are used as in A; the motor neurons B15 and B16 are pooled. Box in each column indicates the mean ± SE. Thin black lines join the corresponding control and modulated values that are plotted against each other in A; the thick black lines join the means. Statistical significance was tested with 2-way ANOVA followed by pairwise multiple comparisons using the Holm–Sidak test. Comparing the control with the modulated mean under each of the three conditions, there were highly significant differences (P < 0.001) at 15 and at 25°C, but no significant difference (P > 0.05) at 15°C in the presence of exogenous SCPB. Comparing between conditions, there was no significant difference between the control means (P > 0.05, "n.s."), but a highly significant difference between the modulated means (P < 0.001, "***"), at 15 and 25°C. For further explanation see Quantification of the modulation of relaxation rate in RESULTS.

 


View larger version (26K):
[in this window]
[in a new window]
 
FIG. 8. Modulation by exogenous modulator: magnitude of the modulation of relaxation rate by 10–7 M SCPB. Figure constructed along the lines of Fig. 5. A: modulated relaxation rate, of the final monitoring contraction elicited by motor neuron B15 after the SCPB had been applied to the muscle and the modulation had fully developed (>16 min after the SCPB application: see Fig. 9A), is plotted against the control relaxation rate. Green lines join matched pairs of measurements, made in the same muscle at 15°C (blue circle; n = 10) and 25°C (red circle; n = 10) (tested in either order in different muscles) with two applications of SCPB separated by wash. B: same measurements of control and modulated relaxation rate as in A plotted along one dimension for statistical comparison. Statistical significance was tested with 2-way ANOVA followed by pairwise multiple comparisons using the Holm–Sidak test. There were highly significant differences (P < 0.001) between the control and modulated means at each temperature as well as between temperatures ("***").

 
For this reason, we have first of all plotted in Fig. 5A, from each of the 40 cross-modulation experiments that we performed, simply the modulated relaxation rate measured about 2 min after the end of the modulator-releasing firing (i.e., that of contraction 2 in Fig. 4) against the control relaxation rate (that of contraction 1 in Fig. 4). There is considerable scatter of the control relaxation rate values, along the horizontal dimension of Fig. 5A, as of the modulated values, along the vertical dimension. No modulation is represented by a diagonal line through the points of equal control and modulated value (gray line marked "0% increase in relaxation rate"). The blue symbols indicate the experiments performed at 15°C and the red symbols at 25°C, with modulator release from motor neuron B15 (circles) or B16 (squares). (The two motor neurons gave similar results and are lumped together in Fig. 5B and in the following discussion.) Figure 5B then shows the same measurements of control and modulated relaxation rate as in Fig. 5A plotted separately along one dimension for statistical comparison. The box in each column indicates the mean ± SE. The thin black lines join the corresponding control and modulated values that are plotted against each other in Fig. 5A; the thick black lines join the means. Together, the plots in Fig. 5, A and B, show a number of noteworthy features.

First, it can be seen that the modulator-releasing firing always brought about an increase—never any significant decrease—in the relaxation rate. Essentially all points in Fig. 5A lie above the 0% diagonal line; the lines joining the corresponding control and modulated values in Fig. 5B all slant upward. In individual experiments the increase could be small (e.g., experiment a in Fig. 5A, lying close to the 0% diagonal line) or large (experiment b, with >3,000% increase in the relaxation rate). The mean increase, however, was statistically highly significant at both 15 and 25°C (see Fig. 5 legend).

Previous work has suggested that a major factor that determines the size of the increase in relaxation rate in this kind of experiment is the control relaxation rate that the modulation acts on (Brezina et al. 1995Go). If the control relaxation rate is small, the absolute increase (the change in units of s–1), and still more the relative increase (the change in %), is large. If the control relaxation rate is already large, the absolute and the relative increases are small. In either case, the modulation tends to bring the relaxation rate to about the same absolute level. A natural explanation is that there is some upper limit, some maximal rate, beyond which the muscle cannot be made to relax any faster by the modulation under those conditions. To demonstrate this again here, we performed six experiments with the two motor neuron cross-modulation paradigm, firing motor neuron B15 to release the SCPs, at 15°C, but in the presence of a saturating concentration of one of the SCPs, 10–5 M SCPB (see Brezina et al. 1996Go; Lloyd et al. 1984Go). These experiments are shown by the open blue circles in Fig. 5, A and B. The exogenous SCP increased—essentially tripled—the relaxation rate to >1.5 s–1, whereupon the modulator-releasing firing had little further effect. All of these experiments lie close to the 0% diagonal line in Fig. 5A; there was no statistically significant difference between their control and modulated means in Fig. 5B (see legend). We can reasonably take the mean modulated relaxation rate in these experiments, about 1.8 s–1, to be the presumed maximal relaxation rate at 15°C. This upper limit is indicated by the horizontal dashed line across Fig. 5, A and B. Note that all of the other experiments at 15°C fall below the line.

Consider now the experiments at 25°C (red symbols) in Fig. 5, A and B. Clearly, these often had modulated relaxation rates considerably larger than any experiment at 15°C, exceeding by a substantial margin the limit just discussed. The green lines in Fig. 5A indicate matched pairs of experiments, at 15 and 25°C, in the two muscles of the same preparation. In some cases at least, the two muscles had similar control relaxation rates—for example, the pair cd—and in each of these cases the modulated relaxation rate was larger at 25 than at 15°C. Finally, although in Fig. 5B the mean control relaxation rate was somewhat larger at 25 than at 15°C (see also Figs. 79), statistical analysis (see legend) showed that the difference was not statistically significant. The mean modulated relaxation rate, however, was substantially larger at 25 than at 15°C, and this difference was highly significant.

It thus appears that in the two motor neuron cross-modulation experiments the modulation of the relaxation rate was not only not weaker at the higher temperature, but on the contrary was significantly stronger.

The modulation of relaxation rate reverses more rapidly at higher temperature

In Fig. 4 it appears that, although the modulation of the relaxation rate was stronger at 25 than at 15°C, it reversed more rapidly after the end of the modulator-releasing firing (as, apparently, did the modulation of contraction amplitude). The group data in Fig. 6, averaging together the time courses of the modulation of the relaxation rate in all of the experiments (without exogenous SCP) in Fig. 5, bear this out. The faster reversal of the modulation at 25 than at 15°C was statistically highly significant (see Fig. 6 legend).



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 6. Two motor neuron cross-modulation: reversal of the modulation of relaxation rate. In each of the experiments in Fig. 5 (without exogenous SCP), the time series of the relaxation rates measured from each of the monitoring contractions from before the modulator-releasing firing to about 16 min after it (i.e., the span shown in Fig. 4, A1 and B1) was normalized to between 0 and 100%, the smallest and largest values measured in that particular experiment. Plotted here are the means ± SE of the normalized time series at 15°C (blue; n = 23) and 25°C (red; n = 11). Statistical significance was tested with 2-way ANOVA followed by pairwise multiple comparisons using the Holm–Sidak test. Overall difference between the 15 and 25°C conditions was highly significant (P < 0.001); ***P < 0.001 and *P < 0.05, for the difference between the 15 and 25°C means at a particular time point.

 
Modulation of the relaxation rate by exogenous modulator is intrinsically stronger at higher temperature

So far, with both of our motor neuron stimulation paradigms, we have studied the modulation, that of the relaxation rate in particular, as one overall input–output process, Motor neuron firing -> Relaxation rate modulation, and determined what we will call its cumulative temperature dependency. The overall process is composed of a sequence of more elementary steps: Motor neuron firing -> Modulator release -> Modulator concentration -> Relaxation rate modulation (Fig. 1). Each of these steps will have its own intrinsic temperature dependency, which will then combine with that of the other steps to give the cumulative temperature dependency of the overall process. To understand the origin of the cumulative temperature dependency, we must therefore know the intrinsic temperature dependency of the individual steps. We already know that of the first step, Motor neuron firing -> Modulator release (Fig. 2). What about the intrinsic temperature dependency of the subsequent steps, in particular the last step, Modulator concentration -> Relaxation rate modulation, the actual modulatory effect itself?4



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 10. Two motor neuron cross-modulation reproduced with a mathematical model incorporating temperature dependency. Here an experiment like that in Fig. 4 was simulated, with motor neuron B16 fired to elicit the monitoring contractions and motor neuron B15 to release the modulator (SCP). f, firing frequency; c, contraction (dimensionless). Model was taken from Brezina et al. (2003a,b) and used without modification except for the incorporation of temperature dependency for two of its steps. Intrinsic temperature dependency of modulator release was modeled from the experimental data in Fig. 2 as described in Fig. 2 legend. Intrinsic temperature dependency of the modulation of the relaxation rate was modeled from the data in Figs. 79 as follows. In the model the modulator-induced increase in relaxation rate is described by the variable R (Fig. 1), which ranges from 0 to a maximal value, Rmax (Brezina et al. 2003a). R is then converted to an increase in the absolute relaxation rate of the muscle, in units of s–1, through a coupling factor {delta}R (see Eq. 1h of Brezina et al. 2003b; here {delta}R = 0.05). Dynamics of R, when controlled by a single modulator such as SCP, are governed by the first-order differential equation

(3)
where is the concentration of SCP, CSCP, modified by the Hill coefficient hR = 0.7, kR+ and kR are rate constants, and t is time. After a step at t = 0 to a new, maintained CSCP(t ≥0) [with an intrinsic time constant of <30 s (Brezina et al. 2003a), the step in CSCP can be considered quasi-instantaneous on the timescale of the experiments in Fig. 9], Eq. 3 predicts a single-exponential rise (or fall) in R

(4)
with steady-state value

(5)
where KR {equiv} kR/kR+, and time constant

(6)
(See Eq. 5 of Brezina et al. 2003a.) To preserve the shape of the steady-state dose–response relation (Eq. 5) as modeled by Brezina et al. (2003a), furthermore apparently at all temperatures (Fig. 7C), the value KR = 3.64 x 10–6 of Brezina et al. (2003a) must be preserved at all temperatures. This restricts the choice of kR+ and kR and requires both to have the same temperature coefficient, Q10,kR. When CSCP(t ≥0) = 0, as in Fig. 9C, Eq. 6 yields {tau}R = 1/kR. Thus the values of kR at 15 and 25°C, and so as well as, through the relationship KR = kR/kR+ = 3.64 x 10–6 , the corresponding values of kR+, can in principle be immediately obtained from the time constants of the fitted exponentials in Fig. 9C. Ratio of these time constants implies {approx} 8.6. However, the absolute values of the time constants are not entirely consistent with some of the previous values compiled by Brezina et al. (2003a), and in general not all of the available data can be completely accommodated within the simplifying framework of Eqs. 46 as discussed by Brezina et al. (2003a). We therefore selected a consensus set of values of kR+, kR, and that appeared to represent the best compromise between the previous data of Brezina et al. (2003a) and the time constants in Fig. 9C, as well as Fig. 9, A and B, here. We set kR+(15°C) = 152.5 s–1 , kR(15°C) = 5.55 x 10–4 s–1 (so that, still, KR = 3.64 x 10–6 ), and, conservatively, = 5. Finally, the temperature-dependent scaling of the magnitude of the modulation seen in Figs. 7C and 9A requires an additional temperature coefficient of Rmax, . Figure 7C suggests {approx} 5 and Fig. 9A {approx} 3.2; we set, conservatively, = 3, with Rmax(15°C) = 100%.

 


View larger version (22K):
[in this window]
[in a new window]
 
FIG. 11. Second, structural mechanism of temperature compensation: analysis with the model. A: cumulative steady-state temperature dependency at successive levels in the sequence Motor neuron firing -> Modulator release -> Modulator concentration -> Relaxation rate modulation, illustrated with motor neuron B15 and SCP. (Motor neuron B16 and MM would present a qualitatively identical picture, although quantitatively the curves would be slightly different.) With input firing of motor neuron B15 at a steady frequency fB15 = 10 Hz, the model was solved analytically for the steady states of rSCP, CSCP, and R, the relevant successive variables in the sequence (Fig. 1). These steady states are given by the equations [taken from Brezina et al. (2003a) or derived from the master dynamic equations therein]

(7)
[SSCP, the size of the releasable pool of SCP, was fixed at the initial size S0,SCP to eliminate the progressive depletion of the pool that would otherwise prevent the system from reaching any nonzero steady state (Brezina et al. 2003a)]

(8)
where v = 10 µl is the effective volume into which the SCP is released and kC,SCP = 0.1 s–1 is the rate constant of SCP removal from v, and either

(9)
without any intrinsic temperature dependency, or


(10)

incorporating the intrinsic temperature dependency of R that was modeled in Fig. 10 legend. (Variables and parameters not defined here have already been introduced, and the specific parameter values used here have been given, in legends of Figs. 2 and 10.) Equations 710 are plotted here as a function of the temperature T, normalized in each case to the magnitude at 15°C; Eqs. 7 and 8, when normalized both yielding the same curve, are shown by the light gray curve; Eq. 9, by the dark gray curve; and Eq. 10, by the black curve. Inset extends the temperature range of the main plot to 40°C. B: intrinsic steady-state dose–response relation of the step Modulator concentration -> Relaxation rate modulation, i.e., plot of R{infty}, as a percentage of Rmax or the temperature-scaled Rmax (either Eq. 9 or Eq. 10), as a function of CSCP [or CMM: the intrinsic dose–response relations of the two modulators are identical (Brezina et al. 1996, 2003a)]. For further explanation see Second, structural mechanism of temperature compensation: saturating curvature of dose–response relations in RESULTS.

 
To examine this, we performed experiments in which we fired motor neuron B15 to elicit regular monitoring contractions as in the two motor neuron cross-modulation paradigm, but then, instead of stimulating motor neuron B16 to release the modulators, we exogenously applied different steady concentrations of a modulator, SCPB, directly onto the ARC muscle. Figure 7 shows the results of one series of these experiments in which we applied cumulatively increasing concentrations of SCPB to determine its dose–response relations at 15 and at 25°C. Figure 7A shows representative contractions at the two temperatures in the presence of 10–9, 10–8, 10–7, 10–6, and 10–5 ("–9"..."–5") M SCPB, superimposed in each case on a control contraction for comparison. Figure 7B then plots the dose–response relations for the net modulation of peak contraction amplitude, and Fig. 7C the dose–response relations for the modulation of the relaxation rate, averaged from five experiments each at 15°C (blue) and 25°C (red).

On contraction amplitude, previous work has shown that the potentiating (Ca-current) and depressing (K-current) effects of the modulators combine asymmetrically so as to give net potentiation at low modulator concentrations and net depression at high concentrations (see Brezina et al. 1995Go, 1996Go, 2003aGo). It appears from Fig. 7, A and B, that both the net potentiation at low SCPB concentrations and the net depression at high SCPB concentrations—and therefore presumably, both of the underlying potentiating and depressing effects—were larger at 25 than at 15°C. Indeed, high SCPB concentrations often produced at 25°C a very substantial net depression (Fig. 7A), such as is only rarely seen at lower temperatures with the SCPs—although it is seen with the MMs, which activate a much larger amplitude of the K current in the muscle (Brezina et al. 1995Go, 1996Go).5

The modulation of the relaxation rate, too, was larger at 25 than at 15°C (Fig. 7C). Indeed, it was much larger: as Fig. 7C shows, it was necessary to scale the 15°C dose–response relation up approximately fivefold before it reproduced, quite well at all concentrations, the magnitude of the dose–response relation obtained at 25°C.6

In another series of these experiments, we applied just a single steady concentration of the modulator, 10–7 M SCPB, and followed systematically the time course of the development of the modulation, and in some experiments also that of its reversal on washout of the SCPB, at both 15 and 25°C. From these experiments, Fig. 8 first presents, in the same way as did Fig. 5, the magnitude of the fully developed modulation of the relaxation rate. Figure 8A plots the modulated relaxation rate against the control relaxation rate in each experiment at 15°C (blue circles) and 25°C (red circles). The green lines join matched pairs of measurements. (In this series of experiments the matched measurements were made not just in the same preparation, but in the very same muscle with two applications of SCPB separated by wash.) Figure 8B then plots the same measurements, and their means, for statistical comparison. Again, the modulated relaxation rate was always much larger at 25 than at 15°C. At 15°C the modulated relaxation rate was smaller, but at 25°C considerably larger, than the presumed maximal relaxation rate attainable at 15°C (horizontal dashed line, taken from Fig. 5). In Fig. 8B, the difference between the mean modulated relaxation rates at the two temperatures was statistically highly significant (see Fig. 8 legend). (In this series of experiments, unlike in Fig. 5, the difference between the mean control relaxation rates—larger at 25 than at 15°C—was also significant.)

Figure 9 summarizes the time course of the modulation of the relaxation rate in these experiments. Figure 9A shows simply the averaged absolute relaxation rates at each time point during the development of the modulation in all of the experiments at 15°C (blue) and at 25°C (red). Figure 9B shows the same data first normalized within each experiment, as in Fig. 6, to discard magnitude differences and allow comparison just of the time course shapes at the two temperatures. Figure 9C shows, similarly normalized, the time courses of the reversal of the modulation on washout of the SCPB. To characterize the time courses quantitatively, we fitted the experimental data with best single-exponential fits (solid black curves) as follows. First, statistical analysis (see Fig. 9 legend) showed that the normalized time courses of the development of the modulation in Fig. 9B were not significantly different at 15 and 25°C. We therefore fitted all of the data in Fig. 9B simultaneously with the same exponential, with time constant {tau} = 179 s. With this time constant, fitting the unnormalized data in Fig. 9A then yielded the magnitudes of the fully developed modulation at 15 and 25°C. In this series of experiments, the modulation was about 3.2-fold larger at 25 than at 15°C. (This result was already implicit in Fig. 8, which used the same data.) Finally, statistical analysis—and just simple inspection—showed that the time courses of the reversal of the modulation in Fig. 9C were very significantly different at 15 and 25°C. We therefore fitted the data in Fig. 9C, over the range of times shown, with two different exponentials at the two temperatures, yielding {tau} = 3,995 s at 15°C and {tau} = 462 s at 25°C. Thus the modulation of the relaxation rate by the exogenous modulator was not only intrinsically larger at the higher temperature, but, just like the modulation by the endogenously released modulator in Figs. 4 and 6, it reversed more rapidly when the modulator was removed.

Reproduction of the observed temperature dependency in a mathematical model of the system

A dynamic model of the B15/B16–ARC neuromuscular system, including the basal ACh-induced contraction and the major components of the modulation—all of the variables X(t) in Fig. 1—was previously constructed from a variety of experimental data by Brezina et al. (2003aGo,b). We now modified the model to incorporate into it temperature dependency for the two steps that we have focused on in this paper. From the experimental data in Fig. 2, we modeled the intrinsic temperature dependency of the release of the modulators (for details see Fig. 2 legend). The continuous curves in Fig. 2 show the modeled release, decreasing, like the experimentally measured release, about 20-fold for a 10°C increase in temperature. From the experimental data in Figs. 7 9, we then modeled the intrinsic temperature dependency of the modulation of the relaxation rate (for details see Fig. 10 legend). Figure 7C suggests that the magnitude of the modulation might increase about fivefold, whereas Fig. 9A suggests about 3.2-fold, for a 10°C increase in temperature. Conservatively, we modeled it as increasing threefold. We also modeled a fivefold increase in the rate constants of the modulation, so that the reversal of the modulation, in particular, would proceed faster at higher tem-perature as observed in Fig. 9C. Again, the fivefold increase was conservative: Fig. 9C suggests a greater than eightfold increase.

With the incorporation of the intrinsic temperature dependencies of these two steps, the model was then well able to reproduce the cumulative temperature dependency that we observed, in particular, in our two motor neuron cross-modulation experiments. Figure 10 shows a simulation of an experiment like that in Fig. 4. After the firing of the modulator-releasing motor neuron, the modeled relaxation rate in Fig. 10, just like the experimentally measured rate in Fig. 4, increases considerably more at 25 than at 15°C, but then falls back to the basal rate much faster.

The modulation of contraction amplitude in the simulation in Fig. 10—not only the net modulation but each of its component effects, which can be inspected in the model—is smaller at 25 than at 15°C, whereas in the real experiment in Fig. 4 it is not smaller, and may even be somewhat larger. Because we did not characterize separately the intrinsic temperature dependencies of each of the component effects, we could not incorporate them into the model as we did that of the modulation of the relaxation rate, but presumably the modeling of such increases of the effects with temperature as were observed in Fig. 7, A and B, would largely correct the deficiency seen in Fig. 10. Even without this correction, nevertheless, it is remarkable in Fig. 10 that the modulation of contraction amplitude is not grossly incorrect. Even though the modulator release decreases 20-fold as the temperature increases from 15 to 25°C, the modulation of contraction amplitude decreasesonly moderately. This is because the system incorporates a second type of temperature-compensating mechanism, quite distinct from the intrinsic temperature dependencies of the various steps, as the following analysis of the model reveals.

Second, structural mechanism of temperature compensation: saturating curvature of dose–response relations

For consistency we will discuss the second mechanism as it operates, once again, in the particular case of the modulation of the relaxation rate, but from its principles it will be seen that the mechanism will operate similarly for all of the modulatory effects, all those that contribute to the modulation of contraction amplitude as well. In the case of the relaxation rate, that there must be an additional mechanism becomes clear when we consider the magnitudes of the intrinsic temperature dependencies that we modeled above. Modulator release decreases 20-fold for a 10°C increase in temperature, whereas the modulation of the relaxation rate increases, intrinsically, only threefold. Cumulatively, when the modulation is caused by the released modulator, we should therefore expect the modulation still to decrease, to a first approximation, approximately 20/threefold—more than sixfold—when the temperature is raised from 15 to 25°C. Instead, we see in Fig. 10 that the modulation does not decrease, and even increases. An additional temperature-compensating mechanism must therefore be operating in the model and, presumably, in the real system. Furthermore, it is not a mechanism that we ourselves have explicitly incorporated into the model. It must be implicit in the structure of the model network—present, for example, even in the previous model of Brezina et al. (2003aGo,b), which had no explicit temperature dependency at all.

To demonstrate the second mechanism more formally, in Fig. 11A we have analytically solved the model equations to plot the cumulative steady-state temperature dependency from the level of motor neuron firing to each of the subsequent levels in the sequence leading to the modulation of the relaxation rate. In other words, for the same representative pattern of motor neuron firing throughout (see Fig. 11 legend), Fig. 11A plots the temperature dependency of modulator release (i.e., of the step Motor neuron firing -> Modulator release; light gray curve), modulator concentration (Motor neuron firing -> Modulator concentration; light gray curve again, identical to that for modulator release), and finally the modulation of the relaxation rate (Motor neuron firing -> Relaxation rate modulation), either without (dark gray curve) or with (black curve) the incorporation of the intrinsic temperature dependency of the step Modulator concentration -> Relaxation rate modulation that we modeled above. The light gray curve of the temperature dependency of modulator release is simply the curve already seen in Fig. 2, along which modulator release decreases 20-fold for a 10°C increase in temperature. Because modulator concentration is, in this model, simply a linear reflection of modulator release (see Eq. 8 in Fig. 11 legend), its cumulative temperature dependency is identical. However, the next level of the model, the modulation of the relaxation rate, does not at all follow this dramatic decrease with temperature. Even without any intrinsic temperature dependency of its own (dark gray curve), the modulation of the relaxation rate decreases only very moderately. The incorporation of the intrinsic temperaturedependency of the modulatory effect then converts the moderate decrease into a positive increase (black curve). Foldwise, however, it is the second mechanism, revealed in the large difference between the cumulative temperature dependency of the modulator concentration and that of the modulation of the relaxation rate, even without its own intrinsic temperature dependency, that is most significant under these circumstances.

What precisely is the nature of this second mechanism? Consider Fig. 11B, which shows simply the intrinsic steady-state dose–response relation of the modulatory effect, the step Modulator concentration -> Relaxation rate modulation. This quasi-hyperbolic relation is everywhere curved concave-down.7 Consequently any change in the input modulator concentration will result in a smaller relative change in the output relaxation rate modulation. For example, the point labeled "Actual at 15°C" is the steady-state point reached with the motor neuron firing used in Fig. 11A at 15°C. When the temperature is now increased to 25°C, there is a 20-fold decrease in modulator concentration and, if the modulation of the relaxation rate simply reflected the modulator concentration linearly, we would also expect a 20-fold decrease in the modulation of the relaxation rate, to the point labeled "Expected at 25°C with linear dose–response relation." Instead, along the curved dose–response relation, the modulation of the relaxation rate decreases much less, to the point labeled "Actual at 25°C."

The inset of Fig. 11A extends the main plot to higher temperatures. Even at much higher temperatures, together the two temperature-compensating mechanisms operating in the system—the "first" mechanism of compensating intrinsic temperature dependency and the second, structural mechanism—do a remarkably good job of maintaining the overall output, the modulation of the relaxation rate. From 15 to about 27°C, as already mentioned, the modulation (black curve) actually increases. Above 27°C it begins to decrease, but even at 40°C—where modulator release and concentration have decreased 1,745-fold relative to 15°C—it still has not decreased below that at 15°C. Only above 40°C does this happen.

Temperature compensation during realistic feeding behavior

The power of the second, structural mechanism will decrease as the motor neuron firing frequency and thus modulator release and modulator concentration decrease, because the system will be operating less along the saturated, horizontal limb and more along the vertical—although still concave-down—limb of the dose–response relation in Fig. 11B. The firing frequencies of the motor neurons B15 and B16 recorded with chronically implanted electrodes in intact, feeding animals (Cropper et al. 1990aGo,cGo; Horn et al. 2004Go) are lower than the frequency used in Fig. 11A, and the released modulators probably reach concentrations only of the order of 100 nM (Brezina et al. 2003aGo). Furthermore, the analysis in Fig. 11 was done with a constant, regular motor neuron firing pattern in the steady state, whereas the firing patterns recorded in vivo are extremely variable and irregular (Horn et al. 2004Go; Zhurov et al. 2004Go) and the modulated neuromuscular system generally operates far from the steady state (Brezina et al. 2003aGo,bGo, 2005Go). Is the temperature compensation still effective under these circumstances?

To answer this question, we ran our model with the motor neuron firing patterns recorded in vivo. Specifically, we used the firing patterns of the motor neurons B15 and B16 from an entire approximately 2.5-h-long meal, consisting of 749 cycles of the feeding behavior, recorded in an intact, freely feeding animal by Horn et al. (2004)Go and analyzed by Brezina et al. (2005)Go. Figure 12 shows the results. The two traces at the top of each of A and B are the 2.5-h-long firing patterns of the two motor neurons, the inputs to the model. The traces below then show the resulting trajectories of the most relevant downstream variables of the model, that is, the concentrations of SCP and MM, the modulation of the relaxation rate, and the final modulated contraction amplitude, at 15°C (A) and at 25°C (B). Clearly, increasing the temperature from 15 to 25°C decreased the modulator concentrations so as to make them essentially invisible on the scale used in Fig. 12. Yet the modulation of the relaxation rate decreased hardly at all, and neither did the modulated contraction amplitude. (Note again that the model incorporated no intrinsic temperature dependency for any component of the modulation of contraction amplitude; with such incorporation, the modulated contraction amplitude would probably have been as large at 25 as at 15°C or even larger.) Figure 12C summarizes measurements, comparable to those made in Figs. 3, 4, and 10, of the absolute relaxation rate of the simulated contraction in each of the 749 cycles of the meal. This final, most realistic measure of the relaxation rate was likewise well maintained at 25°C at a high level (red) not appreciably different from that at 15°C (blue). Furthermore, this maintenance was the result of the modulation: when the modulation of the relaxation rate was disabled in the model, the absolute relaxation rate decreased severalfold, at both 15 and 25°C (empty circles).8 It thus appeared that the temperature-compensating mechanisms operating in the model, and presumably in the real system, are indeed very effective even with the frequencies and patterns in which the motor neurons actually fire in real behavior.


 DISCUSSION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Temperature compensation in the B15/B16–ARC neuromuscular system of Aplysia

The basic experimental result of this report is that, although modulator release from the ARC motor neurons decreases dramatically with increasing temperature (Fig. 2), the modulation of ARC muscle contraction amplitude and relaxation rate brought about by the released modulators does not decrease at all, and may even increase (Figs. 35).

The apparent paradox of this result requires us to reexamine our assumptions. Two potential straightforward explanations must be considered. First, it could be that, after all, modulator release does not decrease with increasing temperature. The data in Fig. 2 are sparse. However, they are part of an extensive, solid body of work by Vilim et al. (1996aGo,bGo, 2000Go) in which modulator release was directly measured by radioimmunoassay rather than merely inferred from the subsequent effects of the released modulators. Furthermore, one of the early effects of the released modulators, elevation of cyclic AMP (cAMP) in the ARC muscle, has likewise been reported to decrease with increasing temperature (Whim and Lloyd 1990Go; see further below).

Second, the possibility has to be considered that the observed changes in contraction amplitude and relaxation rate are not really a result of the released modulators. A very extensive body of evidence (reviewed by Brezina et al. 2003aGo; Weiss et al. 1993Go; see INTRODUCTION) argues, however, that they are. The modulatory peptide genes are expressed by the motor neurons; the peptides are processed and are contained in dense-core vesicles in the motor neuron terminals within the ARC muscle; the modulators are released from these terminals onto the muscle fibers with physiological frequencies and patterns of firing of the motor neurons; exogenous application of the modulators, at the concentrations likely be achieved by the motor neuron firing, mimics the changes in contraction shape observed with the motor neuron firing. It has not yet been possible to carry out in a completely satisfactory manner the final demonstration, that blocking one of the earlier steps in this sequence—release of the modulators, ideally—also blocks the observed changes in contraction shape, because of the lack of specific pharmacological or other tools. (Indeed, our initial motivation for studying temperature dependency was to develop such a tool, in the expectation that, if increased temperature so dramatically decreases the release of the modulators, it should of course decrease all of their effects too, indeed for practical purposes eliminate them to reveal just the basal, unmodulated contractions of the muscle. Evidently, this simple logic can fail.) However, the gap has to a large extent been bridged by modeling of the sort that we have continued in this paper, which has previously shown that, when the intrinsic characteristics of the individual steps in the modulatory sequence are modeled individually, joined together, and activated by the patterns of motor neuron firing that are used in experiments or recorded in vivo, the cumulative changes in contraction shape predicted by the model largely reproduce those that are actually observed (Brezina et al. 2003aGo,bGo, 2005Go). As a last point, we note that it would be difficult to explain how the two motor neuron cross-paradigm in Fig. 4 works if it is not through the release from each motor neuron of some modulatory substance that has long-lasting, general effects on contractions, even those elicited by the other motor neuron.

We believe therefore that we can rule out these two essentially artifactual explanations for our results. We can rule them out all the more easily because we have in this paper developed what we believe is the correct resolution of the apparent paradox. Experimentally, we have found that the final modulatory effects have very considerable intrinsic temperature dependencies of their own, opposite to that of the modulator release (Figs. 7 9), so that as the modulator release decreases with increasing temperature, the modulatory effects intrinsically increase. At the same time, our modeling has revealed another, very different mechanism of temperature compensation that is built into the structure of the system. Because of the curvature of the intrinsic dose–response relations of the modulatory effects, a given decrease in modulator release and concentration decreases the modulatory effects much less (Fig. 11). The model incorporating both of these intrinsic mechanisms of temperature compensation is well able to reproduce the cumulative temperature dependency of the modulation that is brought about, through the entire sequence of modulatory steps, by the motor neuron firing (Fig. 10). That is, even though one of the steps in the sequence—the release of the modulators—undergoes a dramatic decrease as the temperature increases, the final modulation that is observed does not decrease at all. Even though there are large intrinsic temperature dependencies within the sequence, the cumulative final temperature dependency is small.

We have focused here principally on the modulation of the relaxation rate and neglected that of contraction amplitude, largely for technical reasons: the modulation of the relaxation rate is, to all appearances, a single effect whereas that of contraction amplitude is distinctly composite. However, it is likely that the modulation of contraction amplitude is temperature compensated by the same two mechanisms. There are indications in our data that the opposing effects that combine to give the net modulation of contraction amplitude—in particular, the enhancement of Ca current that potentiates contraction amplitude and the activation of K current that depresses it—each have their own intrinsic temperature dependency, increasing, like the modulation of the relaxation rate, with increasing temperature (Fig. 7B). The Ca- and K-current effects' intrinsic dose–response relations, too, are shaped like that of the modulation of the relaxation rate (Brezina et al. 1996Go, 2003aGo).

Clearly, each of the "elementary" steps into which we have analytically decomposed the modulatory sequence here is itself a black box containing lower-level processes, each with its own intrinsic temperature dependency. This is implied already by the fact that the rate of one of our "elementary" processes, that is, the modulator release, decreases with increasing temperature. On thermodynamic grounds, we should expect all truly elementary reaction rates to increase with temperature. The temperature dependency of the modulator release is therefore presumably the resultant of the temperature dependencies of opposing lower-level processes, each increasing with temperature but with the "backward" rate—itself probably the resultant of still lower-level processes—increasing more steeply than the "forward" rate. A similar explanation can be given for our attribution of temperature dependency to the magnitude of the modulation (Figs. 7C and 9A; Fig. 10 legend), which is not a rate at all but presumably is the resultant of the rates of lower-level processes.

For the modulation of the relaxation rate, some of the lower-level detail can be sketched in. The step Modulator concentration -> Relaxation rate modulation probably proceeds, in its early stages, through a conventional sequence of binding of the modulators, the SCPs and MMs, to their respective receptors in the membrane of the ARC muscle fibers, activation of a G protein such as Gs, activation of adenylyl cyclase, and elevation of cAMP in the muscle fibers (Brezina et al. 1994aGo; Hooper et al. 1994bGo; Lloyd et al. 1984Go; Whim and Lloyd 1989Go, 1990Go). Whim and Lloyd (1990)Go reported that, when a fixed concentration of SCP is exogenously applied, the resulting elevation of cAMP is identical at 16 and 22°C. If so, then this entire early sequence of events, up to the elevation of cAMP, has no intrinsic temperature dependency and cannot be the locus of the temperature dependency that we have discovered here in Figs. 79. cAMP then activates cAMP-dependent protein kinase (PKA; Hooper et al. 1994aGo,bGo), which then phosphorylates the giant muscle protein twitchin (Heierhorst et al. 1994Go, 1995Go, 1996Go; Probst et al. 1994Go; see also Siegman et al. 1997Go, 1998Go; Yamada et al. 2001Go). The phosphorylation/dephosphorylation state of twitchin in the ARC muscle has been implicated in the control of the relaxation rate (Probst et al. 1994Go), although the precise mechanism remains unclear. Twitchin is a structural protein in the muscle, but it is also a kinase that, under the control of signals such as the cAMP-dependent phosphorylation as well as Ca2+ acting through Ca2+-binding proteins such as calmodulin and S100 (Heierhorst et al. 1996Go), will phosphorylate, directly or indirectly by phosphorylating other kinases, other muscle proteins. With regard to the temperature dependency, the activity of some molluscan PKAs has been reported to be temperature dependent (e.g., Bardales et al. 2004Go), but probably much of the temperature dependency that we have discovered will be found to originate somewhere within the late, complex steps involving twitchin. Certainly the very slow reversal of the modulation of the relaxation rate, which is strongly temperature dependent (Fig. 9C), is not likely to be rate limited by any such early process as a decay in free Ca2+ concentration, cAMP level, or even PKA activity, all of which decay much faster (Probst et al. 1994Go; Whim and Lloyd 1990Go; see Brezina et al. 2003aGo), but may well reflect a subsequent slow process such as dephosphorylation.

Temperature affects all biological processes, and a full appreciation of how the B15/B16–ARC neuromuscular system functions throughout the temperature range that the animal may encounter will require knowledge of the temperature dependency also of all of the other important components of the system. For example, how does temperature affect the shape of the basal, unmodulated contractions? [There are indications that increasing temperature decreases their amplitude (Vilim et al. 1996aGo) and perhaps increases their relaxation rate (Figs. 7C, 8, and 9A here).] How does temperature affect the firing frequencies and patterns of the motor neurons? Nothing is known furthermore about how the animal may acclimatize to different temperatures in the wild. Nevertheless, we believe that our results here have demonstrated, in basic outline at least, how the combination of the temperature-compensating mechanisms that are built into the structure of the system can maintain its output—the strength of the modulatory tuning and, as a consequence, the modulated contraction shape—relatively constant despite large changes in temperature. In our simulation with realistic motor neuron firing patterns (Fig. 12), the overall strength of the modulatory tuning varies hardly at all from 15 to 25°C and perhaps even higher (Fig. 11A, inset). Because the different dimensions of the contraction shape—amplitude and relaxation rate, for instance—are the result of different processes with somewhat different cumulative temperature dependency, the shape will not remain exactly identical at all temperatures, and the extent of its change will furthermore depend on the motor neuron firing frequency and pattern. However, these small residual changes in contraction shape are unlikely to matter functionally. In this neuromuscular system, the real motor neuron firing patterns and resulting contraction shapes do not appear to be precisely tailored to the specific behavioral demands in any particular cycle of the feeding behavior. Rather, they exhibit large quasi-random variability—much larger than the residual changes with temperature—from one cycle to the next (as can be seen in Fig. 12) so as to diversify, it is conjectured, the feeding movements in the face of environmental uncertainty (Brezina et al. 2005Go; Horn et al. 2004Go).

Temperature compensation in other systems

Because temperature affects all biological reactions, all animals that are active over a range of temperatures, unless they regulate their internal temperature, face this same problem with respect to each of their physiological systems, large and small: how to make the overall output of the system built from individual temperature-sensitive reactions sufficiently temperature insensitive to maintain function throughout the temperature range. A well-known example is that of circadian oscillators, which cannot allow the period of the clock to change if the temperature changes (see, e.g., Roenneberg and Merrow 1999Go). The fact that temperature affects all reactions, and furthermore has—unlike, say, a neuromodulator—universal access to them all, at the same time provides plenty of material from which temperature compensation can be built. In the systems that have been investigated, it is interesting to note the extreme heterogeneity of the types of processes whose intrinsic temperature dependencies complement each other to maintain the overall output. On different time scales, the temperature dependencies of gene expression, biochemical and signal-transduction pathways, biophysical processes in neurons, synaptic physiology, muscle contractility, neuronal firing patterns and higher-order control strategies of the nervous system, and behavior can all play complementary roles (e.g., Feder and Hofmann 1999Go; Heitler and Edwards 1998Go; Montgomery and Macdonald 1990Go; Neumeister et al. 2000Go; Podrabsky and Somero 2004Go; Pörtner 2002Go; Rall and Woledge 1990Go; Rome 1990Go). In addition to this relatively self-evident mechanism of complementation of intrinsic temperature dependencies, it is recognized (e.g., Montgomery and Macdonald 1990Go; see following text) that the way in which the processes are coupled together—the structure of the network—also plays a major role. Temperature-compensatory mechanisms are thus of two types, an instance of each of which we have analyzed in this paper.

General lessons

Changing temperature is just one kind of stress that biological networks must resist (Roenneberg and Merrow 1999Go), and the mechanisms of temperature compensation have much in common with the general mechanisms that have been identified in many networks that make their output robust against dynamic noise, phenotypic parameter variation, and environmental perturbations of all kinds. Both the intrinsic characteristics of the processes within the network and the dynamic nonlinear properties that emerge from the way in which the processes are connected together are important. For example, recent work has particularly emphasized the role of negative feedback in the stabilization of circadian oscillators (e.g., Cheng et al. 2001Go; Gonze et al. 2002Go; Smolen et al. 2001Go; Vilar et al. 2002Go). However, as we have shown here, even in the absence of feedback merely the typical saturating shapes of the dose–response relations in the network exert a stabilizing influence, and this influence will grow as more such steps are connected in sequential cascades (see, e.g., Kenakin 1997Go).

At the same time, the dynamic interconnectedness of such networks makes it difficult to use simple linear logic to predict the outcome of their perturbation. Fox and Lloyd (2001)Go and Koh and Weiss (2005)Go, citing the fact that modulator release decreases dramatically with increasing temperature in the ARC system, used temperature to try to separate, in another buccal neuromuscular system and in the CNS of Aplysia, the subsequent effects of modulators released there from basal, modulator-independent effects. Effects that likewise decreased with increasing temperature were identified as potential effects of the released modulators; effects that did not decrease were not. We, too, initially sought to use this logic in the ARC system itself. There, however, the logic clearly fails: by it, we would have had to conclude that the effects on contraction amplitude and relaxation rate, very well established as caused by the modulators by numerous other lines of evidence, were not modulator dependent. Instead, our studies have led us to discover the temperature-compensating mechanisms that we have described in this paper. In view of the existence of such mechanisms, the use of temperature to separate modulator-dependent from modulator-independent processes requires reconsideration.


 GRANTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was funded by National Institute of Neurological Disorders and Stroke Grant R01 NS-41497 to V. Brezina.


 ACKNOWLEDGMENTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank A. J. Susswein, T. H. Carefoot, and T. Capo for information on the physiological temperatures for Aplysia, and we thank the reviewers of the paper for helpful suggestions.


 FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1 As Fig. 2 implies, release of the peptide modulators has not been directly measured at temperatures >19°C. However, measurements of downstream effects of the released modulators such as the SCP-induced elevation of cAMP in the ARC muscle (Fig. 1) strongly suggest that the release continues to decrease up to at least 22°C (e.g., Whim and Lloyd 1990Go), although one lesson of our work here is that such indirect evidence must be treated with caution because the downstream effects are by no means a straightforward, linear readout of the release (see DISCUSSION). Back

2 The persisting modulation is not simply a long-lasting consequence of the (often quite large) contraction elicited by the modulator-releasing firing. This was shown by Whim and Lloyd (1990)Go, who reversibly applied hexamethonium, a blocker of the ACh receptors in the muscle that mediate the contraction, during the period of the modulator-releasing firing. This completely blocked the contraction elicited by the modulator-releasing firing, but when the hexamethonium was then washed out to bring back the monitoring contractions elicited by the other motor neuron, these were found to be modulated as usual. Back

3 Indeed, although usually the modulator-releasing firing increased the amplitude of the monitoring contractions as in Fig. 4, occasionally, when the modulator-releasing motor neuron was B16, the amplitude was at first decreased, for the first one or two monitoring contractions after the modulator-releasing firing, before increasing in the usual way. Based on what is known quantitatively about the modulatory effects and their interaction (Brezina et al. 1995Go, 1996Go, 2003aGo), this was almost certainly a sign, not of weak modulation, but of strong modulation. The potentiating (Ca-current) and depressing (K-current) effects combine asymmetrically so as to give net potentiation at low modulator concentrations and net depression at high concentrations (see, e.g., Fig. 7, A and B). The depressing effect has fast intrinsic dynamics, and so will reverse rapidly after the modulator-releasing firing as the released modulator concentration falls, whereas the potentiating effect, with slow dynamics, will persist. Finally, the depressing effect is activated particularly strongly by the MMs, released from motor neuron B16. Thus, very likely, the transient net depression occasionally observed after the modulator-releasing firing of motor neuron B16 reflected exceptionally strong release of the MMs and activation of the depressing (K-current) effect. By extension, the depressing effect probably contributed to the net modulation of contraction amplitude much more generally, even when it was hidden in the usual net potentiation. Back

4 Clearly, each of the "elementary" steps mentioned here can be broken down further and, in particular, the last step Modulator concentration -> Relaxation rate modulation can be broken down to Modulator concentration -> Muscle cAMP level -> Relaxation rate modulation (Fig. 1). Nevertheless, we study here the Modulator concentration -> Relaxation rate modulation step as a unit, primarily because we can easily manipulate the input to this step, modulator concentration, but not so easily the muscle cAMP level. For the same reason, also, the muscle cAMP level is not an explicit part of our model of the system (see Brezina et al. 2003aGo) with which we seek to reproduce our experimental results in Figs. 10 and 11. Back

5 Potentially, the changes in contraction amplitude were a problem for the interpretation of the modulation of the relaxation rate common because the relaxation rate depends to some degree on the amplitude: small contractions tend to relax faster than large contractions (Brezina et al. 2003bGo). However, in these and indeed in most experiments in this paper, the changes in amplitude were relatively small (the magnitudes in Figs. 4 and 7, A and B, are typical), considerably smaller than the severalfold potentiation, or complete suppression of contractions, that can occur (see Brezina et al. 1996Go). In some of the experiments in Fig. 9, we explicitly tested to what extent the changes in contraction amplitude could account for the changes in relaxation rate, by adjusting, once the modulation had developed, the firing frequency of the contraction-eliciting motor neuron so as to return the amplitude of the modulated contraction back to the control amplitude. The change in relaxation rate nevertheless remained very substantial. Although there was indeed an inverse relationship between contraction amplitude and relaxation rate, its slope was relatively shallow, so that the changes in contraction amplitude could account for only a small part of the observed changes in relaxation rate. We therefore made no special attempt to allow for the changing contraction amplitude in our analysis of the modulation of the relaxation rate in this paper. Back

6 Our treatment of the relaxation-rate data in Figs. 79 is based on the simplest conceptual model—then formalized in our actual model in Figs. 1012—of the effect of temperature: that it is only the modulation that is temperature dependent. Effects of temperature on the basal relaxation rate are neglected. With the data in Fig. 5 this treatment appears to be adequate, but the data in Figs. 79 suggest that the basal relaxation rate, too, has a significant temperature dependency. A completely realistic model would need to postulate a separate temperature dependency for each process; our model here lumps the temperature dependency of the basal relaxation rate into that of the modulation. Functionally (e.g., in Fig. 12), the result will be similar in either case: because the modulation is coupled, through the firing of the motor neurons, to the basal contraction, the final modulated relaxation rate will have a similar cumulative temperature dependency no matter how the cumulative temperature dependency has been decomposed into the intrinsic temperature dependencies of the processes within the system. Back

7 The dose–response relation (Eqs. 9 and 10 in Fig. 11 legend) is strictly hyperbolic only if the Hill coefficient h, here hR, is equal to 1; here hR = 0.7. The concave-down curvature of the dose–response relation, and thus the consequences discussed here, holds for h ≤ 1; for h > 1, it holds only over part of the domain of the dose–response relation. Back

8 Even without the modulation of the relaxation rate, the relaxation rates were slightly different at 15 and 25°C in Fig. 12C because the modulation of contraction amplitude, still present, also altered to some extent the relaxation rate, in the model as in reality (see footnote 5). Back

Address for reprint requests and other correspondence: V. Brezina, Department of Neuroscience, Box 1218, Mt. Sinai School of Medicine, 1 Gustave L. Levy Place, New York, NY 10029 (E-mail Vladimir.Brezina{at}mssm.edu)


 REFERENCES
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Bardales JR, Díaz-Enrich MJ, Ibarguren I, and Villamarín JA. Isoforms of cAMP-dependent protein kinase in the bivalve mollusk Mytilus galloprovincialis: activation by cyclic nucleotides and effect of temperature. Arch Biochem Biophys 432: 71–78, 2004.[Web of Science][Medline]

Brezina V, Bank B, Cropper EC, Rosen S, Vilim FS, Kupfermann I, and Weiss KR. Nine members of the myomodulin family of peptide cotransmitters at the B16–ARC neuromuscular junction of Aplysia. J Neurophysiol 74: 54–72, 1995.[Abstract/Free Full Text]

Brezina V, Church PJ, and Weiss KR. Temporal pattern dependence of neuronal peptide transmitter release: models and experiments. J Neurosci 20: 6760–6772, 2000a.[Abstract/Free Full Text]

Brezina V, Evans CG, and Weiss KR. Enhancement of Ca current in the accessory radula closer muscle of Aplysia californica by neuromodulators that potentiate its contractions. J Neurosci 14: 4393–4411, 1994a.[Abstract]

Brezina V, Evans CG, and Weiss KR. Activation of K current in the accessory radula closer muscle of Aplysia californica by neuromodulators that depress its contractions. J Neurosci 14: 4412–4432, 1994b.[Abstract]

Brezina V, Horn CC, and Weiss KR. Modeling neuromuscular modulation in Aplysia. III. Interaction of central motor commands and peripheral modulatory state for optimal behavior. J Neurophysiol 93: 1523–1556, 2005.[Abstract/Free Full Text]

Brezina V, Orekhova IV, and Weiss KR. Functional uncoupling of linked neurotransmitter effects by combinatorial convergence. Science 273: 806–810, 1996.[Abstract]

Brezina V, Orekhova IV, and Weiss KR. Optimization of rhythmic behaviors by modulation of the neuromuscular transform. J Neurophysiol 83: 260–279, 2000b.[Abstract/Free Full Text]

Brezina V, Orekhova IV, and Weiss KR. Neuromuscular modulation in Aplysia. I. Dynamic model. J Neurophysiol 90: 2592–2612, 2003a.[Abstract/Free Full Text]

Brezina V, Orekhova IV, and Weiss KR. Neuromuscular modulation in Aplysia. II. Modulation of the neuromuscular transform in behavior. J Neurophysiol 90: 2613–2628, 2003b.[Abstract/Free Full Text]

Brezina V and Weiss KR. Ca2+-activated K current in the ARC muscle of Aplysia. J Neurophysiol 73: 993–1004, 1995.[Abstract/Free Full Text]

Carefoot TH. Aplysia: its biology and ecology. Oceanogr Mar Biol Ann Rev 25: 167–284, 1987.

Cheng P, Yang Y, and Liu Y. Interlocked feedback loops contribute to the robustness of the Neurospora circadian clock. Proc Natl Acad Sci USA 98: 7408–7413, 2001.[Abstract/Free Full Text]

Church PJ, Whim MD, and Lloyd PE. Modulation of neuromuscular transmission by conventional and peptide transmitters released from excitatory and inhibitory motor neurons in Aplysia. J Neurosci 13: 2790–2800, 1993.[Abstract]

Cohen JL, Weiss KR, and Kupfermann I. Motor control of buccal muscles in Aplysia. J Neurophysiol 41: 157–180, 1978.[Free Full Text]

Cropper EC, Kupfermann I, and Weiss KR. Differential firing patterns of the peptide-containing cholinergic motor neurons B15 and B16 during feeding behavior in Aplysia. Brain Res 522: 176–179, 1990a.[CrossRef][Web of Science][Medline]

Cropper EC, Lloyd PE, Reed W, Tenenbaum R, Kupfermann I, and Weiss KR. Multiple neuropeptides in cholinergic motor neurons of Aplysia: evidence for modulation intrinsic to the motor circuit. Proc Natl Acad Sci USA 84: 3486–3490, 1987a.[Abstract/Free Full Text]

Cropper EC, Miller MW, Tenenbaum R, Kolks MAG, Kupfermann I, and Weiss KR. Structure and action of buccalin: a modulatory neuropeptide localized to an identified small cardioactive peptide-containing cholinergic motor neuron of Aplysia californica. Proc Natl Acad Sci USA 85: 6177–6181, 1988.[Abstract/Free Full Text]

Cropper EC, Miller MW, Vilim FS, Tenenbaum R, Kupfermann I, and Weiss KR. Buccalin is present in the cholinergic motor neuron B16 of Aplysia and it depresses accessory radula closer muscle contractions evoked by stimulation of B16. Brain Res 512: 175–179, 1990b.[CrossRef][Web of Science][Medline]

Cropper EC, Price D, Tenenbaum R, Kupfermann I, and Weiss KR. Release of peptide cotransmitters from a cholinergic motor neuron under physiological conditions. Proc Natl Acad Sci USA 87: 933–937, 1990c.[Abstract/Free Full Text]

Cropper EC, Tenenbaum R, Kolks MAG, Kupfermann I, and Weiss KR. Myomodulin: a bioactive neuropeptide present in an identified cholinergic buccal motor neuron of Aplysia. Proc Natl Acad Sci USA 84: 5483–5486, 1987b.[Abstract/Free Full Text]

Cropper EC, Vilim FS, Alevizos A, Tenenbaum R, Kolks MAG, Rosen S, Kupfermann I, and Weiss KR. Structure, bioactivity, and cellular localization of myomodulin B: a novel Aplysia peptide. Peptides 12: 683–690, 1991.[CrossRef][Web of Science][Medline]

Evans CG, Vilim FS, Harish O, Kupfermann I, Weiss KR, and Cropper EC. Modulation of radula opener muscles in Aplysia. J Neurophysiol 82: 1339–1351, 1999.[Abstract/Free Full Text]

Feder ME and Hofmann GE. Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol 61: 243–282, 1999.[CrossRef][Web of Science][Medline]

Fox LE and Lloyd PE. Evidence that post-tetanic potentiation is mediated by neuropeptide release in Aplysia. J Neurophysiol 86: 2845–2855, 2001.[Abstract/Free Full Text]

Gev S, Achituv Y, and Susswein AJ. Seasonal determinants of the life cycle in two species of Aplysia found in shallow waters along the Mediterranean coasts of Israel. J Exp Mar Biol Ecol 74: 67–73, 1984.[CrossRef]

Gomez D, Halloy J, and Goldbeter A. Robustness of circadian rhythms with respect to molecular noise. Proc Natl Acad Sci USA 99: 673–678, 2002.[Abstract/Free Full Text]

Heierhorst J, Kobe B, Feil SC, Parker MW, Benian BM, Weiss KR, and Kemp BE. Ca2+/S100 regulation of giant protein kinases. Nature 380: 636–639, 1996.[CrossRef][Medline]

Heierhorst J, Probst WC, Kohanski RA, Buku A, and Weiss KR. Phosphorylation of myosin regulatory light chains by the molluscan twitchin kinase. Eur J Biochem 233: 426–431, 1995.[Web of Science][Medline]

Heierhorst J, Probst WC, Vilim FS, Buku A, and Weiss KR. Autophosphorylation of molluscan twitchin and interaction of its kinase domain with calcium/calmodulin. J Biol Chem 269: 21086–21093, 1994.[Abstract/Free Full Text]

Heitler WJ and Edwards DH. Effect of temperature on a voltage-sensitive electrical synapse in crayfish. J Exp Biol 201: 503–513, 1998.[Abstract/Free Full Text]

Hooper SL, Brezina V, Cropper EC, and Weiss KR. Flexibility of muscle control by modulation of muscle properties. In: Beyond Neurotransmission: Neuromodulation and Its Importance for Information Processing, edited by Katz PS. Oxford, UK: Oxford Univ. Press, 1999, p. 241–274.

Hooper SL, Probst WC, Cropper EC, Kupfermann I, and Weiss KR. SCP application or B15 stimulation activates cAPK in the ARC muscle of Aplysia. Brain Res 657: 337–341, 1994a.[CrossRef][Web of Science][Medline]

Hooper SL, Probst WC, Cropper EC, Kupfermann I, and Weiss KR. Myomodulin application increases cAMP and activates cAMP-dependent protein kinase in the accessory radula closer muscle of Aplysia. Neurosci Lett 179: 167–170, 1994b.[CrossRef][Web of Science][Medline]

Horn CC, Zhurov Y, Orekhova IV, Proekt A, Kupfermann I, Weiss KR, and Brezina V. Cycle-to-cycle variability of neuromuscular activity in Aplysia feeding behavior. J Neurophysiol 92: 157–180, 2004.[Abstract/Free Full Text]

Hurwitz I, Cropper EC, Vilim FS, Alexeeva V, Susswein AJ, Kupfermann I, and Weiss KR. Serotonergic and peptidergic modulation of the buccal mass protractor muscle (I2) in Aplysia. J Neurophysiol 84: 2810–2820, 2000.[Abstract/Free Full Text]

Keating C and Lloyd PE. Differential modulation of motor neurons that innervate the same muscle but use different excitatory transmitters in Aplysia. J Neurophysiol 82: 1759–1767, 1999.[Abstract/Free Full Text]

Kenakin T. Pharmacologic Analysis of Drug-Receptor Interaction (3rd ed). Philadelphia, PA: Lippincott–Raven, 1997.

Koh H-Y and Weiss KR. Peptidergic contribution to posttetanic potentiation at a central synapse of Aplysia. J Neurophysiol 94: 1281–1286, 2005.[Abstract/Free Full Text]

Kupfermann I. Feeding behavior in Aplysia: a simple system for the study of motivation. Behav Biol 10: 1–26, 1974.[CrossRef][Web of Science][Medline]

Kupfermann I, Brezina V, Cropper EC, Deodhar D, Probst WC, Rosen SC, Vilim FS, and Weiss KR. Reconfiguration of the peripheral plant during various forms of feeding behaviors in the mollusc Aplysia. In: Neurons, Networks, and Motor Behavior, edited by Stein PSG, Grillner S, Selverston AI, and Stuart DG. Cambridge, MA: MIT Press, 1997, p. 217–222.

Kupfermann I and Carew TJ. Behavior patterns of Aplysia californica in its natural environment. Behav Biol 12: 317–337, 1974.[CrossRef][Web of Science][Medline]

Lloyd PE, Kupfermann I, and Weiss KR. Evidence for parallel actions of a molluscan neuropeptide and serotonin in mediating arousal in Aplysia. Proc Natl Acad Sci USA 81: 2934–2937, 1984.[Abstract/Free Full Text]

Montgomery JC and Macdonald JA. Effects of temperature on nervous system: implications for behavioral performance. Am J Physiol Regul Integr Comp Physiol 259: R191–R196, 1990.[Abstract/Free Full Text]

Neumeister H, Ripley B, Preuss T, and Gilly WF. Effects of temperature on escape jetting in the squid Loligo opalescens. J Exp Biol 203: 547–557, 2000.[Abstract]

Orekhova IV, Alexeeva V, Church PJ, Weiss KR, and Brezina V. Multiple presynaptic and postsynaptic sites of inhibitory modulation by myomodulin at ARC neuromuscular junctions of Aplysia. J Neurophysiol 89: 1488–1502, 2003.[Abstract/Free Full Text]

Podrabsky JE and Somero GN. Changes in gene expression associated with acclimation to constant temperatures and fluctuating daily temperatures in an annual killifish Austrofundulus limnaeus. J Exp Biol 207: 2237–2254, 2004.[Abstract/Free Full Text]

Pörtner HO. Physiological basis of temperature-dependent biogeography: trade-offs in muscle design and performance in polar ectotherms. J Exp Biol 205: 2217–2230, 2002.[Abstract/Free Full Text]

Probst WC, Cropper EC, Heierhorst J, Hooper SL, Jaffe H, Vilim F, Beushausen S, Kupfermann I, and Weiss KR. cAMP-dependent phosphorylation of Aplysia twitchin may mediate modulation of muscle contractions by neuropeptide cotransmitters. Proc Natl Acad Sci USA 91: 8487–8491, 1994.[Abstract/Free Full Text]

Rall JA and Woledge RC. Influence of temperature on mechanics and energetics of muscle contraction. Am J Physiol Regul Integr Comp Physiol 259: R197–R203, 1990.[Abstract/Free Full Text]

Roenneberg T and Merrow M. Circadian systems and metabolism. J Biol Rhythms 14: 449–459, 1999.[Abstract/Free Full Text]

Rome LC. Influence of temperature on muscle recruitment and muscle function in vivo. Am J Physiol Regul Integr Comp Physiol 259: R210–R222, 1990.[Abstract/Free Full Text]

Siegman MJ, Funabara D, Kinoshita S, Watabe S, Hartshorne DJ, and Butler TM. Phosphorylation of a twitchin-related protein controls catch and calcium sensitivity of force production in invertebrate smooth muscle. Proc Natl Acad Sci USA 95: 5383–5388, 1998.[Abstract/Free Full Text]

Siegman MJ, Mooers SU, Li C, Narayan S, Trinkle-Mulcahy L, Watabe S, Hartshorne DJ, and Butler TM. Phosphorylation of a high molecular weight (~600 kDa) protein regulates catch in invertebrate smooth muscle. J Muscle Res Cell Motil 18: 655–670, 1997.[CrossRef][Web of Science][Medline]

Smolen P, Baxter DA, and Byrne JH. Modeling circadian oscillations with interlocking positive and negative feedback loops. J Neurosci 21: 6644–6656, 2001.[Abstract/Free Full Text]

Vilar JMG, Kueh HY, Barkai N, and Leibler S. Mechanisms of noise-resistance in genetic oscillators. Proc Natl Acad Sci USA 99: 5988–5992, 2002.[Abstract/Free Full Text]

Vilim FS, Cropper EC, Price DA, Kupfermann I, and Weiss KR. Release of peptide cotransmitters in Aplysia: regulation and functional implications. J Neurosci 16: 8105–8114, 1996a.[Abstract/Free Full Text]

Vilim FS, Cropper EC, Price DA, Kupfermann I, and Weiss KR. Peptide cotransmitter release from motorneuron B16 in Aplysia californica: costorage, corelease, and functional implications. J Neurosci 20: 2036–2042, 2000.[Abstract/Free Full Text]

Vilim FS, Cropper EC, Rosen SC, Tenenbaum R, Kupfermann I, and Weiss KR. Structure, localization, and action of buccalin B: a bioactive peptide from Aplysia. Peptides 15: 959–969, 1994.[CrossRef][Web of Science][Medline]

Vilim FS, Price DA, Lesser W, Kupfermann I, and Weiss KR. Costorage and corelease of modulatory peptide cotransmitters with partially antagonistic actions on the accessory radula closer muscle of Aplysia californica. J Neurosci 16: 8092–8104, 1996b.[Abstract/Free Full Text]

Weiss KR, Brezina V, Cropper EC, Heierhorst J, Hooper SL, Probst WC, Rosen SC, Vilim FS, and Kupfermann I. Physiology and biochemistry of peptidergic cotransmission in Aplysia. J Physiol (Paris) 87: 141–151, 1993.

Weiss KR, Cohen JL, and Kupfermann I. Modulatory control of buccal musculature by a serotonergic neuron (metacerebral cell) in Aplysia. J Neurophysiol 41: 181–203, 1978.[Free Full Text]

Whim MD and Lloyd PE. Frequency-dependent release of peptide cotransmitters from identified cholinergic motor neurons in Aplysia. Proc Natl Acad Sci USA 86: 9034–9038, 1989.[Abstract/Free Full Text]

Whim MD and Lloyd PE. Neuropeptide cotransmitters released from an identified cholinergic motor neuron modulate neuromuscular efficacy in Aplysia. J Neurosci 10: 3313–3322, 1990.[Abstract]

Yamada A, Yoshio M, Kojima H, and Oiwa K. An in vitro assay reveals essential protein components for the "catch" state of invertebrate smooth muscle. Proc Natl Acad Sci USA 98: 6635–6640, 2001.[Abstract/Free Full Text]

Zhurov Y, Weiss KR, and Brezina V. Intraburst variability of motor neuron spike timing shapes contractions of the ARC muscle of Aplysia. Abstract 537.5, 2004 Abstract Viewer. Washington, DC: Society for Neuroscience, online, 2004.




This article has been cited by other articles:


Home page
J. Neurophysiol.Home page
I. Hurwitz, A. Harel, S. Markowitz, O. Noy, and A. J. Susswein
Control of Feeding in Aplysia With Ad Libitum Access to Food: Presence of Food Increases the Intervals Between Feeding Bouts
J Neurophysiol, January 1, 2006; 95(1): 106 - 118.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
J. T. Birmingham
Simple Mechanism for Stabilizing Motor Output. Focus on "Temperature Compensation of Neuromuscular Modulation in Aplysia"
J Neurophysiol, November 1, 2005; 94(5): 2997 - 2998.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
94/5/3259    most recent
00481.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhurov, Y.
Right arrow Articles by Brezina, V.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zhurov, Y.
Right arrow Articles by Brezina, V.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2005 by the The American Physiological Society.