No Arabic abstract
Ventricular tachycardia (VT) and ventricular fibrillation (VF) are lethal rhythm disorders, which are associated with the occurrence of abnormal electrical scroll waves in the heart. Given the technical limitations of imaging and probing, the in situ visualization of these waves inside cardiac tissue remains a challenge. Therefore, we must, perforce, rely on in-silico simulations of scroll waves in mathematical models for cardiac tissue to develop an understanding of the dynamics of these waves in mammalian hearts. We use direct numerical simulations of the Hund-Rudy-Dynamic (HRD) model, for canine ventricular tissue, to examine the interplay between electrical scroll-waves and conduction and ionic inhomogeneities, in anatomically realistic canine ventricular geometries with muscle-fiber architecture. We find that millimeter-sized, distributed, conduction inhomogeneities cause a substantial decrease in the scroll wavelength, thereby increasing the probability for wave breaks; by contrast, single, localized, medium-sized ($simeq $ cm) conduction inhomogeneities, exhibit the potential to suppress wave breaks or enable the self-organization of wave fragments into stable, intact scrolls. We show that ionic inhomogeneities, both distributed or localised, suppress scroll-wave break up. The dynamics of a stable rotating wave is not affected significantly by such inhomogeneities, except at high concentrations of distributed inhomogeneities, which can cause a partial break up of scroll waves. Our results indicate that inhomogeneities in the canine ventricular tissue are less arrhythmogenic than inhomogeneities in porcine ventricular tissue, for which an earlier in silico study has shown that the inhomogeneity-induced suppression of scroll waves is a rare occurrence.
We conduct a systematic,direct-numerical-simulation study,in mathematical models for ventricular tissue,of the dependence of spiral-and scroll-wave dynamics on $G_{Kr}$, the maximal conductance of the delayed rectifier Potassium current($I_{Kr}$) channel,and the parameter $gamma_{Cao}$,which determines the magnitude and shape of the current $I_{CaL}$ for the L-type calcium-current channel,in both square and anatomically realistic,whole-ventricle simulation domains using canine and human models. We use ventricular geometry with fiber-orientation details and employ a physiologically realistic model for a canine ventricular myocyte. We restrict ourselves to an HRD-model parameter regime, which does not produce spiral- and scroll-wave instabilities because of other,well-studied causes like a very sharp action-potential-duration-restitution (APDR) curve or early after depolarizations(EADs) at the single-cell level. We find that spiral- or scroll-wave dynamics are affected predominantly by a simultaneous change in $I_{CaL}$ and $I_{Kr}$,rather than by a change in any one of these currents;other currents do not have such a large effect on these wave dynamics in this parameter regime of the HRD model.We obtain stability diagrams in the $G_{Kr} -gamma_{Cao}$ plane.In the 3D domain,the geometry of the domain supports the confinement of the scroll waves and makes them more stable compared to their spiral-wave counterparts in 2D domain. We have also carried out a comparison of our HRD results with their counterparts for the human-ventricular TP06 model and have found important differences. In both these models,to make a transition,(from broken-wave to stable-scroll states or vice versa) we must simultaneously increase $I_{Kr}$ and decrease $I_{CaL}$;a modification of only one of these currents is not enough to effect this transition.
Cilia are elastic hairlike protuberances of the cell membrane found in various unicellular organisms and in several tissues of most living organisms. In some tissues such as the airway tissues of the lung, the coordinated beating of cilia induce a fluid flow of crucial importance as it allows the continuous cleaning of our bronchia, known as mucociliary clearance. While most of the models addressing the question of collective dynamics and metachronal wave consider homogeneous carpets of cilia, experimental observations rather show that cilia clusters are heterogeneously distributed over the tissue surface. The purpose of this paper is to investigate the role of spatial heterogeneity on the coherent beating of cilia using a very simple one dimensional model for cilia known as the rower model. We systematically study systems consisting of a few rowers to hundreds of rowers and we investigate the conditions for the emergence of collective beating. When considering a small number of rowers, a phase drift occurs, hence a bifurcation in beating frequency is observed as the distance between rowers clusters is changed. In the case of many rowers, a distribution of frequencies is observed. We found in particular the pattern of the patchy structure that shows the best robustness in collective beating behavior, as the density of cilia is varied over a wide range.
Spiral waves of excitation in cardiac tissue are associated with life-threatening cardiac arrhythmias. It is, therefore, important to study the electrophysiological factors that affect the dynamics of these spiral waves. By using an electrophysiologically detailed mathematical model of a myocyte (cardiac cell), we study the effects of cellular parameters, such as membrane-ion-channel conductances, on the properties of the action-potential (AP) of a myocyte. We then investigate how changes in these properties, specifically the upstroke velocity and the AP duration (APD), affect the frequency $omega$ of a spiral wave in the mathematical model that we use for human-ventricular tissue. We find that an increase (decrease) in this upstroke-velocity or a decrease (increase) in the AP duration increases (decreases) $omega$. We also study how other intercellular factors, such as the fibroblast-myocyte coupling, diffusive coupling strength, and the effective number of neighboring myocytes, modulate $omega$. Finally, we demonstrate how a spiral wave can drift to a region with a high density of fibroblasts. Our results provide a natural explanation for the anchoring of spiral waves in highly fibrotic regions in fibrotic hearts.
In this work, we introduce a novel computational framework that we developed to use numerical simulations to investigate the complexity of brain tissue at a microscopic level with a detail never realised before. Directly inspired by the advances in computational neuroscience for modelling brain cells, we propose a generative model that enables us to simulate molecular diffusion within realistic digitalised brain cells, such as neurons and glia, in a completely controlled and flexible fashion. We validate our new approach by showing an excellent match between the morphology and simulated DW-MR signal of the generated digital model of brain cells and those of digital reconstruction of real brain cells from available open-access databases. We demonstrate the versatility and potentiality of the framework by showing a select set of examples of relevance for the DW-MR community. Further development is ongoing, which will support even more realistic conditions like dense packing of numerous 3D complex cell structures and varying cell surface permeability.
Marine species reproduce and compete while being advected by turbulent flows. It is largely unknown, both theoretically and experimentally, how population dynamics and genetics are changed by the presence of fluid flows. Discrete agent-based simulations in continuous space allow for accurate treatment of advection and number fluctuations, but can be computationally expensive for even modest organism densities. In this report, we propose an algorithm to overcome some of these challenges. We first provide a thorough validation of the algorithm in one and two dimensions without flow. Next, we focus on the case of weakly compressible flows in two dimensions. This models organisms such as phytoplankton living at a specific depth in the three-dimensional, incompressible ocean experiencing upwelling and/or downwelling events. We show that organisms born at sources in a two-dimensional time-independent flow experience an increase in fixation probability.