Countering reproducibility issues in mathematical models with software engineering techniques: A case study using a one-dimensional mathematical model of the atrioventricular node

DyCeModel: a tool for 1D simulation for distribution of plant hormones controlling tissue patterning

To study the mechanisms of growth and development, it is necessary to analyze the dynamics of the tissue patterning regulators in time and space and to take into account their effect on the cellular dynamics within a tissue. Plant hormones are the main regulators of the cell dynamics in plant tissues; they form gradients and maxima and control molecular processes in a concentration-dependent manner. Here, we present DyCeModel, a software tool implemented in MATLAB for one-dimensional simulation of tissue with a dynamic cellular ensemble, where changes in hormone (or other active substance) concentration in the cells are described by ordinary differential equations (ODEs). We applied DyCeModel to simulate cell dynamics in plant meristems with different cellular structures and demonstrated that DyCeModel helps to identify the relationships between hormone concentration and cellular behaviors. The tool visualizes the simulation progress and presents a video obtained during the calculation. Importantly, the tool is capable of automatically adjusting the parameters by fitting the distribution of the substance concentrations predicted in the model to experimental data taken from the microscopic images. Noteworthy, DyCeModel makes it possible to build models for distinct types of plant meristems with the same ODEs, recruiting specific input characteristics for each meristem. We demonstrate the tool's efficiency by simulation of the effect of auxin and cytokinin distributions on tissue patterning in two types of Arabidopsis thaliana stem cell niches: the root and shoot apical meristems. The resulting models represent a promising framework for further study of the role of hormone-controlled gene regulatory networks in cell dynamics.

A compact multi-functional model of the rabbit atrioventricular node with dual pathways

The atrioventricular node (AVN) is considered a “black box”, and the functioning of its dual pathways remains controversial and not fully understood. In contrast to numerous clinical studies, there are only a few mathematical models of the node. In this paper, we present a compact, computationally lightweight multi-functional rabbit AVN model based on the Aliev-Panfilov two-variable cardiac cell model. The one-dimensional AVN model includes fast (FP) and slow (SP) pathways, primary pacemaking in the sinoatrial node, and subsidiary pacemaking in the SP. To obtain the direction-dependent conduction properties of the AVN, together with gradients of intercellular coupling and cell refractoriness, we implemented the asymmetry of coupling between model cells. We hypothesized that the asymmetry can reflect some effects related to the complexity of the real 3D structure of AVN. In addition, the model is accompanied by a visualization of electrical conduction in the AVN, revealing the interaction between SP and FP in the form of ladder diagrams. The AVN model demonstrates broad functionality, including normal sinus rhythm, AVN automaticity, filtering of high-rate atrial rhythms during atrial fibrillation and atrial flutter with Wenckebach periodicity, direction-dependent properties, and realistic anterograde and retrograde conduction curves in the control case and the cases of FP and SP ablation. To show the validity of the proposed model, we compare the simulation results with the available experimental data. Despite its simplicity, the proposed model can be used both as a stand-alone module and as a part of complex three-dimensional atrial or whole heart simulation systems, and can help to understand some puzzling functions of AVN.