Experimental measurements from guinea pig ventricular myocytes of dynamic clamp perturbations are used to fit the behavior of the Livshitz Rudy 2009 ventricular model. Contains experimental data from 14 cells (Cell_1 - Cell_14) as well as the Livshitz Rudy 2009 model itself (Cell_15).
References:
- Devenyi, Ryan A., et al. "Differential roles of two delayed rectifier potassium currents in regulation of ventricular action potential duration and arrhythmia susceptibility." The Journal of physiology 595.7 (2017): 2301-2317.
- Livshitz, Leonid, and Yoram Rudy. "Uniqueness and stability of action potential models during rest, pacing, and conduction using problem-solving environment." Biophysical journal 97.5 (2009): 1265-1276.
These files contain the cell membrane capacitance and injected current of the dynamic clamp experiment. The files are organized by row:
- Pacing current injection
- IKr current injection
- IKs current injection
- ICaL current injection
- IK1 current injection
- ICaT current injection
- INaK current injection
- INCX current injection
- Cell membrane capacitance
These files contain cell membrane voltage during each dynamic clamp perturbation. The files are organized by row:
- Static pacing membrane voltage response
- IKr - membrane voltage response
- IKs - membrane voltage response
- ICaL - membrane voltage response
- IK1 - membrane voltage response
- ICaT - membrane voltage response
- INaK - membrane voltage response
- INCX - membrane voltage response
A genetic algorithm and the Livshitz Rudy 2009 model used in the dynamic clamp experiments have been added to this repository as subtrees. If you wish to update the libraries, you must first add the remote repository and merge in any new changes.
Located in the analysis/ directory.
- Nominal_Model_Fitting - The genetic algorithm is used to fit the model to objective data created using the nominal parameters of the model. In short, the ideal fit has an error is ~0 with all parameters set to 1.0.
- Simulation_vs_Objective - A comparison is made between the full dynamic clamp simulation and the simulation run during the GA individual evaluations.
The fitting executable requires 3 arguments:
- Objective file
- Protocol file
- Genetic algorithm settings file
An example command to fit the model to nominal parameters is as follows:
make && ./GA_Fitting objectives/Objective_GA_Protocol_Simulation.dat protocols/Protocol_DynClamp_Simulation.dat settings.ga