Excitation-contraction coupling/mitochondrial energetics (ECME) model (Cortassa et al. 2006)

 Download zip file 
Help downloading and running models
Accession:105383
"An intricate network of reactions is involved in matching energy supply with demand in the heart. This complexity arises because energy production both modulates and is modulated by the electrophysiological and contractile activity of the cardiac myocyte. Here, we present an integrated mathematical model of the cardiac cell that links excitation-contraction coupling with mitochondrial energy generation. The dynamics of the model are described by a system of 50 ordinary differential equations. The formulation explicitly incorporates cytoplasmic ATP-consuming processes associated with force generation and ion transport, as well as the creatine kinase reaction. Changes in the electrical and contractile activity of the myocyte are coupled to mitochondrial energetics through the ATP, Ca21, and Na1 concentrations in the myoplasmic and mitochondrial matrix compartments. ..."
Reference:
1 . Cortassa S, Aon MA, Marbán E, Winslow RL, O'Rourke B (2003) An integrated model of cardiac mitochondrial energy metabolism and calcium dynamics. Biophys J 84:2734-55 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell; Electrogenic pump;
Brain Region(s)/Organism:
Cell Type(s): Heart cell;
Channel(s): I L high threshold; I Sodium; I Potassium; Na/Ca exchanger; I_SERCA;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: C or C++ program;
Model Concept(s): Activity Patterns; Temporal Pattern Generation; Signaling pathways; Calcium dynamics;
Implementer(s):
Search NeuronDB for information about:  I L high threshold; I Sodium; I Potassium; Na/Ca exchanger; I_SERCA;
#pragma once

/*     ----------------------------------------------------

         NOTICE OF COPYRIGHT AND OWNERSHIP OF SOFTWARE

         Copyright 2004, The Johns Hopkins University
            School of Medicine. All rights reserved.
			For research use only; commercial use prohibited.
			Distribution without permission of Raimond L. Winslow
			not permitted. rwinslow@bme.jhu.edu

         Name of Program: Guinea Pig C++: Coupled, Algbraic, BB, MCA
         Version: Documented Version, version 1.0.1
         Date: August 2004

       -----------------------------------------------------  
*/

#include "integrator.h"
#include "model.h"
#include "fileIO.h"

#include <sundialstypes.h>	/* definitions of types realtype and             */
							/* integertype, and the constant FALSE           */
#include <cvode.h>			/* prototypes for CVodeMalloc, CVode, and CVodeFree, */
							/* constants OPT_SIZE, BDF, NEWTON, SV, SUCCESS,     */
							/* NST, NFE, NSETUPS, NNI, NCFN, NETF                */
#include <cvdense.h>		/* prototype for CVDense, constant DENSE_NJE         */
#include <nvector.h>		/* definitions of type N_Vector and macro N_VIth,    */
#include <nvector_serial.h>	/* definitions of type N_Vector and macro N_VIth,    */
							/* prototypes for N_VNew, N_VFree                    */
#include <dense.h>			/* definitions of type DenseMat, macro DENSE_ELEM    */

#include <stdlib.h>  
#include <iostream>
using namespace std;

class IntegratorCVode :
	public Integrator
{
public:
	IntegratorCVode(void);
	virtual ~IntegratorCVode(void);

	void setParameters(fileIO& data);

	void setupCVode(Model *model);

	void integrateModel(Model *model, fileIO *out, int runNumber);

protected:
	void iterateToTime(double time);
	void refreshCVode(Model *model);

	//Internal variables
	double cvode_t;
	bool usingErrorWeights;
	N_Vector ew;
	N_Vector cvode_y;
	N_Vector cvode_ic;
	void* cvode_mem;
	double* ropt;
	long* iopt;
	M_Env machEnv;

	//Parameters
	double reltol;
	double abstol;
	double step_max;
	double step_min;
	double start_time;
	double stepsize;
};

Loading data, please wait...