Saucerman, Brunton, Michailova, McCulloch, 2003

Model Status

This CellML version of the Saucerman and McCulloch 2003 model was modified by a post-doc in the Computational Biology lab at The University of Oxford. The model cannot be run in COR due to the presence of "circular arguments" or "DAEs". However this version of the CellML model does run in PCEnv - although it is unstable and will give different simulation results each time it is run (i.e. the model solves, but unreliably). This version of the Saucerman model is the one which is currently being used and discussed on the Physiome tracker. The tracker item is: https://tracker.physiomeproject.org/show_bug.cgi?id=1474.

Model Structure

In cardiac myocytes, the beta-adrenergic signalling network is triggered in response to norepinephrine and epinephrine binding to the Gs-protein coupled beta-adrenergic receptor. Receptor-ligand binding activates the coupled Gs-protein, which in turn activates adenylate cyclase. ATP is converted into the secondary messenger cyclic AMP (cAMP). cAMP promotes dissociation of the protein kinase A (PKA) holoenzyme, whose catalytic subunits go on to phosphorylate a wide range of target proteins, two of which are the L-type calcium channel and phospholamban - which both play essential roles in the regulation of calcium dynamics and transport. In healthy cardiac myocytes, the end result of this signal transduction pathway is to provide coordinated control of contractility, metabolism and gene regulation. However, altered beta-adrenergic signalling may also play a role in the progression of heart failure.

Molecular components of this signalling pathway have been studied in detail as potential therapeutic targets in heart failure. However, the situation is complex, with intracellular compartmentation and functional integration of pathways being suggested as playing an important role in the signalling outcome. Systems level mathematical modelling will help make sense of large quantities of experimental data. Since the Hodgkin-Huxley Squid Axon Model was published in 1952, many ionic models of cell electrophysiology have been developed. Recently systems modelling has been identified as a potential method for improving understanding of signaling networks and the response to genetic and pharmaceutical perturbations.

While to date portions of the mammalian beta-adrenergic signalling network have been described by mathematical models of neuron and HEK-293 cells, no analysis has modelled and validated an entire pathway from ligand to effectors such that neuro-hormonal regulation of cell physiology can be predicted. In the publication described here, Jeffrey Saucerman et al. have developed a systems biology model to investigate the molecular mechanisms underlying the control of the beta-adrenergic signalling network over cardiac myocyte contractility (see figure below). The model was validated using a wide range of experimental data, and model simulations were used to investigate quantitatively the effects of specific molecular perturbations (with potential for pharmaceutical testing). By performing systems analysis the effects of molecular perturbations in the beta-adrenergic signalling network may be understood within the context of integrative physiology.

The model of the signal transduction pathway was embedded within an extension of the Luo-Rudy Ventricular Model II (dynamic), 1994, which was modified for the rabbit ventricular myocyte (Puglisi-Bers Rabbit Ventricular Myocyte Model, 2001). The model was further adapted by inserting the equations for the L-type calcium channel from the Jafri-Rice-Winslow Ventricular Model, 1998, and also the steady-state potassium currents from the Pandit et al. Adult Rat Left Ventricular Myocyte Model, 2001.

The complete original paper reference is cited below:

Modeling beta-adrenergic control of cardiac myocyte contractility in silico , Jeffrey J. Saucerman, Laurence L. Brunton, Anushka P. Michailova, and Andrew D. McCulloch, 2003, Journal of Biological Chemistry , 48, 47997-48003. (Full text (HTML) and PDF versions of the article are available on the Journal of Biological Chemistry website.) PubMed ID: 12972422

Schematic diagram of the integrated model components, including the beta-adrenergic network, calcium handling, and the electrophysiology of the rat ventricular myocyte.

Schematic diagram highlighting the signalling network component of the model.