We developed a mathematical model of the underlying cellular mechanisms responsible for the changes in sinus cycle length (SCL) elicited by vagal stimulation in intact animals. The model incorporated a stimulation-mediated depletion of the releasable pool of acetylcholine (ACh) in the nerve endings, the in vitro reaction kinetics of acetylcholinesterase, and the electrical activity of a pacemaker cell with six membrane ionic currents. SCL increased linearly with the frequency of simulated vagal stimulation, as it does in animal experiments, because the concentration of ACh in the neuroeffector junction ([ACh]) saturated as the frequency of stimulation was increased and because SCL increased geometrically in response to increases in [ACh]. The dependence of SCL on the timing of vagal stimulation in the cardiac cycle resulted, in part, from the dependence of [ACh] on SCL. Simulated vagal stimulation entrained the sinus node because the rate of activation and inactivation of ACh-activated K + channels depended only weakly on membrane potential during diastolic depolarization. SCL increased geometrically with [ACh], because 1) during diastolic depolarization, the amplitude of the ACh-activated K + current was approximately equal to the amplitude of the sum of the other ionic currents, 2) [ACh] was low enough to saturate neither acetylcholinesterase nor the cellular system that activates the ACh-activated K + channels, 3) the pacemaker cell membrane behaved electrotonically like a capacitor, and 4) the sum of all the ionic currents increased linearly with the amplitude of the ACh-activated K + current. (Circulation Research 1989;65:133O-1339) O ur goal is to explain the cellular basis for the changes in sinus cycle length (SCL) elicited by vagal stimulation by using mathematical models of the underlying cellular physiology. In recent papers, we have presented mathematical models of the release of acetylcholine (ACh) from vagal nerve endings, the degradation of ACh in the neuroeffector junctions of the sinus node, and the changes in SCL elicited by ACh.1 . 2 By combining these mathematical models of the underlying cellular physiology, we create a more complete mathematical model that predicts the changes in SCL elicited by vagal stimulation in the intact animal. Our principal goal is to determine how the underlying cellular physiology, which we represented mathematically in our computer model, causes changes in SCL.
Methods
Mathematical ModelsACh release. Figure 1); the letter v was chosen because the releasable pool of ACh may be stored in intracellular vesicles.3 Let r equal the fastest rate of renewal of the releasable pool of ACh. Let p equal the fraction of the releasable pool of ACh, v, that is released by each vagal stimulus. Let u represent the vagal firing pattern. The change in the quantity of ACh that is available for release is represented by the equation: dv/dt=r(l-v)-pvu(t) (1) where u(t)=l if the vagus nerve is stimulated at time t, and u(t)=0 otherwise (t=0 msec at the start of the stimulation). AC...