Synaptic release was simulated using a Simulink sequential storage model with three vesicular pools. Modeling was modular and easily extendable to the systems with greater number of vesicular pools, parallel input, or time-varying parameters. Given an input (short or long tetanic trains, patterned or random stimulation) and the storage model, the vesicular release, the replenishment of various vesicular pools, and the vesicular content of all pools could be simulated for the time-invariant and time-varying storage systems. From the input stimuli and either a noiseless or a noisy output, the parameters of such storage systems could also be estimated using the optimization technique that minimizes in the least square sense the error between the observed release and the predicted release. All parameters of the storage model could be evaluated with sufficiently long input-output data pairs. Not surprisingly, the parameters characterizing the processes near the release locus, such as the fractional release and the size of the immediately available pool and its coupling to the small store, as well as the state variables associated with the immediately available pool, such as its vesicular content and replenishment, could be determined with fewer stimuli. The possibility of estimating parameters with random inputs extends the applicability of the method to in vivo synapses with the physiological inputs. The parameter estimation was also possible under the time-variant, but slowly changing, conditions as well as for open systems that are part of larger vesicular storage systems but whose parameters can either not be reliably determined or are of no interest. The quality of parameter estimation was monitored continuously by comparing the observed and predicted output and/or estimated parameters with the true values. Finally, the method was tested experimentally using the rat phrenic-diaphragm neuromuscular junction.

No abstract

No abstract

No abstract

Tracking spectral changes of rapidly varying signals is a demanding task. In this study, we explore on Monte Carlo-simulated glutamate-activated AMPA patch and synaptic currents whether a wavelet analysis offers such a possibility. Unlike Fourier methods that determine only the frequency content of a signal, the wavelet analysis determines both the frequency and the time. This is owing to the nature of the basis functions, which are infinite for Fourier transforms (sines and cosines are infinite), but are finite for wavelet analysis (wavelets are localized waves). In agreement with previous reports, the frequency of the stationary patch current fluctuations is higher for larger currents, whereas the mean-variance plots are parabolic. The spectra of the current fluctuations and mean-variance plots are close to the theoretically predicted values. The median frequency of the synaptic and nonstationary patch currents is, however, time dependent, though at the peak of synaptic currents, the median frequency is insensitive to the number of glutamate molecules released. Such time dependence demonstrates that the "composite spectra" of the current fluctuations gathered over the whole duration of synaptic currents cannot be used to assess the mean open time or effective mean open time of AMPA channels. The current (patch or synaptic) versus median frequency plots show hysteresis. The median frequency is thus not a simple reflection of the overall receptor saturation levels and is greater during the rise phase for the same saturation level. The hysteresis is due to the higher occupancy of the doubly bound state during the rise phase and not due to the spatial spread of the saturation disk, which remains remarkably constant. Albeit time dependent, the variance of the synaptic and nonstationary patch currents can be accurately determined. Nevertheless the evaluation of the number of AMPA channels and their single current from the mean-variance plots of patch or synaptic currents is not highly accurate owing to the varying number of the activatable AMPA channels caused by desensitization. The spatial nonuniformity of open, bound, and desensitized AMPA channels, and the time dependence and spatial nonuniformity of the glutamate concentration in the synaptic cleft, further reduce the accuracy of estimates of the number of AMPA channels from synaptic currents. In conclusion, wavelet analysis of nonstationary fluctuations of patch and synaptic currents expands our ability to determine accurately the variance and frequency of current fluctuations, demonstrates the limits of applicability of techniques currently used to evaluate the single channel current and number of AMPA channels, and offers new insights into the mechanisms involved in the generation of unitary quantal events at excitatory central synapses.