Computing conduction velocity restitution

Author: Jason Bayer jason.bayer@ihu-liryc.fr

This example demonstrates how to compute conduction velocity restitution in cardiac tissue. To run the experiments of this tutorial do

Introduction

Conduction velocity restitution is an important property of cardiac tissue. As pacing frequency is increased, conduction velocity will become slower. For this tutorial, the user will shown how to construct conduction velocity restitution curves to describe the conduction properties of cardiac tissue in response to various pacing protocols.

Problem setup

1D cable model

A 1.0 cm cable of epicardial ventricular myocytes is used to generate a CV restitution curve for a user defined pacing protocol. The model domain was discretized with linear finite elements with an average edge length of 0.01 cm.

Ionic model

This tutorial uses the most recent version of the ten Tusscher ionic model for human ventricular myocytes [Tusscher2006] that was modified for the study [Bayer2016] . This ionic model is labeled GTT2_fast in openCARP's LIMPET library.

Pacing protocol

The left side of the 1D cable model is paced with 5-ms-long stimuli at twice capture amplitude for an S1S2 restitution pacing protocol defined by the user inputs.The user sets the cycle length and number of beats for S1 pacing, and the range of cycle lengths to apply for the S2.

Conduction velocity

Activation times are computed for each S2 beat of the pacing protocol using the openCARP option LATs (see tutorial X). CV is then computed along the cable by taking the difference in activation times at the locations 0.25 cm and 0.75 cm divided by the distance between the two points.

To run the experiments of this tutorial do

Usage

The following optional arguments are available (default values are indicated):

After running run.py, the CI and CV are output into the ASCII file CVrestitution_GTT2_fast_bcl_.....

If the program is run with the --visualize option, the CV restitution curve in the file above will be plotted using pythons plotting functions.

Tasks

  1. Determine the minimum S2 CI for the default parameters
  2. Determine the effect of decreasing Gil by 75% on the minimum S2 CI.
  3. Determine the effect of increasing Gil by 50% on the minimum S2 CI.

Solutions to the tasks

  1. The minimum CI for the default parameters is 325 ms.
  1. Reducing Gil by 75% slows down CV and shifts the minimum S2 CI to 300 ms.
  1. Increasing Gil by 50% speeds up CV and shifts the minimum S2 CI to 375 ms.

References

Bayer2016

Bayer JD, Lalani GG, EJ Vigmond, SM Narayan, NA Trayanova. Mechanisms linking electrical alternans and clinical ventricular arrhythmia in human heart failure. Heart Rhythm, 13(9):1922-1931, 2016. [Pubmed]

© Copyright 2020 openCARP project    Supported by DFG    Contact    Imprint and data protection