# Model of the Mammalian ET cell Computational model and files to recreate the results from the paper "A Computational Model of the Mammalian External Tufted Cell" https://doi.org/10.1016/j.jtbi.2018.10.003 Author: Ryan Viertel ## Usage: data = ET(input,sampling_rate); input: input vector, if no input then just a vector of zeros sampling_rate: rate at which the input vector should be sampled. 1000 for milisecond returns the following struct: data.T - time vector data.X - ODE variables at each time step * data.X(:,1) - Membrane Potential * data.X(:,2) - nK * data.X(:,3) - hNaP * data.X(:,4) - hH * data.X(:,5) - mLVA * data.X(:,6) - hLVA * data.X(:,7) - mBK * data.X(:,8) - Calcium * data.X(:,9) - nHVK data.events - list of spike events data.which - event type * 1 - spike * 2 - burst start * 3 - burst end data.current - system currents * data.current(:,1) = transient sodium * data.current(:,2) = fast potassium * data.current(:,3) = leak * data.current(:,4) = persistent sodium * data.current(:,5) = hyperpolarization activated * data.current(:,6) = LVA calcium * data.current(:,7) = HVA calcium * data.current(:,8) = large conductance potassium * data.current(:,9) = HVK current ## example ### create the input vector input = zeros(1,5000); ### run the model data = ET(input,1000); ### plot the voltage trace plot(data.T,data.X(:,1)) ## ME-PCM The code used to sample the model throughout parameter space to determine stability and investigate the effect of model parameters on model output is found in the ME-PCM directory ## xpp The ODE file to recreate the bifurcation diagram is found in the xpp directory