The Bueno-Orovio–Cherry–Fenton model, also simply called Bueno-Orovio model, is a minimal ionic model for human ventricular cells. It belongs to the category of phenomenological models, because of its characteristic of describing the electrophysiological behaviour of cardiac muscle cells without taking into account in a detailed way the underlying physiology and the specific mechanisms occurring inside the cells.
This mathematical model reproduces both single cell and important tissue-level properties, accounting for physiological action potential development and conduction velocity estimations. It also provides specific parameters choices, derived from parameter-fitting algorithms of the MATLAB Optimization Toolbox, for the modeling of epicardial, endocardial and myd-myocardial tissues. In this way it is possible to match the action potential morphologies, observed from experimental data, in the three different regions of the human ventricles. The Bueno-Orovio–Cherry–Fenton model is also able to describe reentrant and spiral wave dynamics, which occurs for instance during tachycardia or other types of arrhythmias.
From the mathematical perspective, it consists of a system of four differential equations. One PDE, similar to the monodomain model, for an adimensional version of the transmembrane potential, and three ODEs that define the evolution of the so-called gating variables, i.e. probability density functions whose aim is to model the fraction of open ion channels across a cell membrane.
Mathematical modeling
The system of four differential equations reads as follows:
where is the spatial domain and is the final time. The initial conditions are , , , .
refers to the Heaviside function centered in . The adimensional transmembrane potential can be rescaled in mV by means of the affine transformation . , and refer to gating variables, where in particular can be also used as an indication of intracellular calcium concentration (given in the adimensional range [0, 1] instead of molar concentration).
and are the fast inward, slow outward and slow inward currents respectively, given by the following expressions:
All the above-mentioned ionic density currents are partially adimensional and are expressed in .
Different parameters sets, as shown in Table 1, can be used to reproduce the action potential development of epicardial, endocardial and mid-myocardial human ventricular cells. There are some constants of the model, which are not located in Table 1, that can be deduced with the following formulas:
where the temporal constants, i.e. are expressed in seconds, whereas and are adimensional.
The diffusion coefficient results in a value of , which comes from experimental tests on human ventricular tissues.
In order to trigger the action potential development in a certain position of the domain , a forcing term , which represents an externally applied density current, is usually added at the right hand side of the PDE and acts for a short time interval only.
Table 1: values of the parameters for different positions of the human heart Parameter Unity of measure - - - - - - seconds seconds seconds seconds seconds - - seconds seconds seconds seconds seconds seconds - - seconds seconds - - seconds - - Epicardium 0 1.55 0.3 0.13 0.006 0.006 60e-3 1150e-3 1.4506e-3 60e-3 15e-3 65 0.03 200e-3 0.11e-3 400e-3 6e-3 30.0181e-3 0.9957e-3 2.0458 0.65 2.7342e-3 16e-3 2.0994 0.9087 1.8875e-3 0.07 0.94 Endocardium 0 1.56 0.3 0.13 0.2 0.006 75e-3 10e-3 1.4506e-3 6e-3 140e-3 200 0.016 280e-3 0.1e-3 470e-3 6e-3 40e-3 1.2e-3 2 0.65 2.7342e-3 2e-3 2.0994 0.9087 2.9013e-3 0.0273 0.78 Myocardium 0 1.61 0.3 0.13 0.1 0.005 80e-3 1.4506e-3 1.4506e-3 70e-3 8e-3 200 0.016 280e-3 0.078e-3 410e-3 7e-3 91e-3 0.8e-3 2.1 0.6 2.7342e-3 4e-3 2.0994 0.9087 3.3849e-3 0.01 0.5
Weak formulation Assume that refers to the vector containing all the gating variables, i.e. , and contains the corresponding three right hand sides of the ionic model. The Bueno-Orovio–Cherry–Fenton model can be rewritten in the compact form:
Let and be two generic test functions.
To obtain the weak formulation: multiply by the first equation of the model and by the equations for the evolution of the gating variables. Integrating both members of all the equations in the domain :
perform an integration by parts of the diffusive term of the first monodomain-like equation:
Finally the weak formulation reads: Find and , , such that
Numerical discretization There are several methods to discretize in space this system of equations, such as the finite element method (FEM) or isogeometric analysis (IGA).
Time discretization can be performed in several ways as well, such as using a backward differentiation formula (BDF) of order or a Runge–Kutta method (RK).
Space discretization with FEM Let be a tessellation of the computational domain by means of a certain type of elements (such as tetrahedrons or hexahedra), with representing a chosen measure of the size of a single element . Consider the set of polynomials with degree smaller or equal than over an element . Define as the finite dimensional space, whose dimension is . The set of basis functions of is referred to as .
The semidiscretized formulation of the first equation of the model reads: find projection of the solution on , , such that
with , semidiscretized version of the three gating variables, and is the total ionic density current.
The space discretized version of the first equation can be rewritten as a system of non-linear ODEs by setting and :
where , and .
The non-linear term can be treated in different ways, such as using state variable interpolation (SVI) or ionic currents interpolation (ICI). In the framework of SVI, by denoting with and the quadrature nodes and weights of a generic element of the mesh , both and are evaluated at the quadrature nodes:
The equations for the three gating variables, which are ODEs, are directly solved in all the degrees of freedom (DOF) of the tessellation separately, leading to the following semidiscrete form:
Time discretization with BDF (implicit scheme) With reference to the time interval , let be the chosen time step, with number of subintervals. A uniform partition in time is finally obtained.
At this stage, the full discretization of the Bueno-Orovio ionic model can be performed both in a monolithic and segregated fashion. With respect to the first methodology (the monolithic one), at time , the full problem is entirely solved in one step in order to get by means of either Newton method or Fixed-point iterations:
where and are extrapolations of transmembrane potential and gating variables at previous timesteps with respect to , considering as many time instants as the order of the BDF scheme. is a coefficient which depends on the BDF order .
If a segregated method is employed, the equation for the evolution in time of the transmembrane potential and the ones of the gating variables are numerically solved separately:
Firstly, is calculated, using an extrapolation at previous timesteps for the transmembrane potential at the right hand side:
Secondly, is computed, exploiting the value of that has just been calculated:
Another possible segregated scheme would be the one in which is calculated first, and then it is used in the equations for .
See also Monodomain model Bidomain model
References