
Modelling Neurons & Action Potentials


This page contains various modelling techniques for neurons, action potentials and neuronal stimulation. The techniques range from overly simplified to relatively intricate but still abstract the under workings to simulate the desired behavior. Other writeups exist on my site exploring these details such as modelling of ion channels themselves. The writeups provided here explore:
 Modelling both a UniDirectional and a BiDirectional Action Potential
 Modelling Cardiac Action Potentials
 Modelling Hodgkin and Huxley Neuronal Model using ODE45
 Dynamical Systems Approach to Neuronal Modelling
 Basic Hodgkin Huxley Neuronal Model
Modelling both a UniDirectional and a BiDirectional Action Potential
The following write up investigates the CRRSS equations in modelling a uni and a bidirectional action potential. The first problem develops a Matlab program that will solve for a propagating action potential along a 9cm axon segment. The second problem investigates an extracellular monopole source with an anodic square pulse followed by a cathodic square pulse, varying amplitudes to generate bidirectional action potentials and generating surface plots of membrane voltage as well as a video. The third problem plots various values at the instant the two action potentials have propagated halfway to each end of the axon. The fourth problem experiments with combinations of multiple anodes and cathodes as well as waveshapes and pulsewidths to generate unidirectional propagating action potentials along the axon.
Modelling Cardiac Action Potentials
This writeup studies the movement of ions in a canine ventricular myocyte using the HundRudy dynamic model. The HRd model constructs the canine ventricular epicardial action potential based on the underlying ionic processes. This experiment could provide insight because many arrhythmias are manifested in abnormalities of action potential generation at the cellular level.
Modelling Hodgkin Huxley Neuronal Model using ODE45
This writeup explores the use of ODE45 in modelling the HodgkinHuxley equations.
Overview of each Problem:
Upon closer inspection I believe some of my dynamics may have been off. There should not be an influx of Potassium before its efflux. Also I may have flipped the signs of the ion currents somewhere as. Na rushes in so it should be negative, while K rushes out so it should be positive.
Overview of each Problem:
 Develops a Matlab program that will solve for membrane voltage for a given stimulus. This involves solving a system of four, coupled first order differential equations. ODE45 is used to integrate the model.
 Assuming the resting membrane potential is 65 mV, problem 2 integrates the model until steady state membrane potential is reached (no stimulus is applied). This steadystate value of membrane voltage is then used as the new resting membrane potential and the process is repeated until the initial estimate of membrane potential matches the integrated solution. Then plots are generated of the voltage time course of the first and last model integration.
 Using a depolarizing square pulse with duration of 0.35 ms, the magnitude of the pulse was increased until an action potential was initiated. Various plots for a typical action potential were then generated.
 This problem explores how soon after an action potential can another action potential be initiated by the same depolarizing square pulse, as well as by much stronger pulse of triple the magnitude.
 Use a stimulus in a hyperpolizing way, increasing the duration of the stimulus until an anode break initiation.
Upon closer inspection I believe some of my dynamics may have been off. There should not be an influx of Potassium before its efflux. Also I may have flipped the signs of the ion currents somewhere as. Na rushes in so it should be negative, while K rushes out so it should be positive.