No Arabic abstract
The analytical theory of our earlier study (Mortensen et al. (2021), Mathematical Medicine and Biology, 38(1), pp. 106-131) is extended to address the outstanding cases of fibroblast barrier distribution and myocyte strait distribution. In particular, closed-form approximations to the resting membrane potential and to the critical parameter values for propagation are derived for these two non-uniform fibroblast distributions and are in good agreement with numerical estimates.
The electrical coupling between myocytes and fibroblasts and the spacial distribution of fibroblasts within myocardial tissues are significant factors in triggering and sustaining cardiac arrhythmias but their roles are poorly understood. This article describes both direct numerical simulations and an asymptotic theory of propagation and block of electrical excitation in a model of atrial tissue with myocyte-fibroblast coupling. In particular, three idealised fibroblast distributions are introduced: uniform distribution, fibroblast barrier and myocyte strait, all believed to be constituent blocks of realistic fibroblast distributions. Primary action potential biomarkers including conduction velocity, peak potential and triangulation index are estimated from direct simulations in all cases. Propagation block is found to occur at certain critical values of the parameters defining each idealised fibroblast distribution and these critical values are accurately determined. An asymptotic theory proposed earlier is extended and applied to the case of a uniform fibroblast distribution. Biomarker values are obtained from hybrid analytical-numerical solutions of coupled fast-time and slow-time periodic boundary value problems and compare well to direct numerical simulations. The boundary of absolute refractoriness is determined solely by the fast-time problem and is found to depend on the values of the myocyte potential and on the slow inactivation variable of the sodium current ahead of the propagating pulse. In turn, these quantities are estimated from the slow-time problem using a regular perturbation expansion to find the steady state of the coupled myocyte-fibroblast kinetics. The asymptotic theory gives a simple analytical expression that captures with remarkable accuracy the block of propagation in the presence of fibroblasts.
In this paper, we derive an effective macroscale description suitable to describe the growth of biological tissue within a porous tissue-engineering scaffold. As in our recent work (Holden textit{et al.} A multiphase multiscale model for nutrient limited tissue growth, The ANZIAM Journal, 2018, doi:10.1017/S1446181118000044) the underlying tissue dynamics is described as a multiphase mixture, thereby naturally accommodating features such as interstitial growth and active cell motion. Via a linearisation of the underlying multiphase model (whose nonlinearity poses significant challenge for such analyses), we obtain, by means of multiple-scales homogenisation, a simplified macroscale model that nevertheless retains explicit dependence on both the microscale scaffold structure and the tissue dynamics. The model we obtain comprises Darcy flow, and differential equations for the volume fraction of cells within the scaffold and the concentration of nutrient, required for growth. These are coupled to underlying Stokes-type cell problems that provide permeability tensors to parameterise the macroscale description. In Holden textit{et al.}, the cell problems retain macroscale dependence, posing significant computational challenges; here, we obtain a decoupled system whereby the quasi-steady cell-problems may be solved separately from the macroscale description, thereby greatly reducing the complexity associated with fully-coupled multiscale descriptions. Moreover, we indicate how the formulation is influenced by a set of alternative microscale boundary conditions.S
Mathematical models of cardiac electrical excitation are increasingly complex, with multiscale models seeking to represent and bridge physiological behaviours across temporal and spatial scales. The increasing complexity of these models makes it computationally expensive to both evaluate long term (>60 seconds) behaviour and determine sensitivity of model outputs to inputs. This is particularly relevant in models of atrial fibrillation (AF), where individual episodes last from seconds to days, and inter-episode waiting times can be minutes to months. Potential mechanisms of transition between sinus rhythm and AF have been identified but are not well understood, and it is difficult to simulate AF for long periods of time using state-of-the-art models. In this study, we implemented a Moe-type cellular automaton on a novel, topologically correct surface geometry of the left atrium. We used the model to simulate stochastic initiation and spontaneous termination of AF, arising from bursts of spontaneous activation near pulmonary veins. The simplified representation of atrial electrical activity reduced computational cost, and so permitted us to investigate AF mechanisms in a probabilistic setting. We computed large numbers (~10^5) of sample paths of the model, to infer stochastic initiation and termination rates of AF episodes using different model parameters. By generating statistical distributions of model outputs, we demonstrated how to propagate uncertainties of inputs within our microscopic level model up to a macroscopic level. Lastly, we investigated spontaneous termination in the model and found a complex dependence on its past AF trajectory, the mechanism of which merits future investigation.
A continuum model of epithelial tissue mechanics was formulated using cellular-level mechanical ingredients and cell morphogenetic processes, including cellular shape changes and cellular rearrangements. This model can include finite deformation, and incorporates stress and deformation tensors, which can be compared with experimental data. Using this model, we elucidated dynamical behavior underlying passive relaxation, active contraction-elongation, and tissue shear flow. This study provides an integrated scheme for the understanding of the mechanisms that are involved in orchestrating the morphogenetic processes in individual cells, in order to achieve epithelial tissue morphogenesis.
Within developing embryos, tissues flow and reorganize dramatically on timescales as short as minutes. This includes epithelial tissues, which often narrow and elongate in convergent extension movements due to anisotropies in external forces or in internal cell-generated forces. However, the mechanisms that allow or prevent tissue reorganization, especially in the presence of strongly anisotropic forces, remain unclear. We study this question in the converging and extending Drosophila germband epithelium, which displays planar polarized myosin II and experiences anisotropic forces from neighboring tissues, and we show that in contrast to isotropic tissues, cell shape alone is not sufficient to predict the onset of rapid cell rearrangement. From theoretical considerations and vertex model simulations, we predict that in anisotropic tissues two experimentally accessible metrics of cell patterns, the cell shape index and a cell alignment index, are required to determine whether an anisotropic tissue is in a solid-like or fluid-like state. We show that changes in cell shape and alignment over time in the Drosophila germband predict the onset of rapid cell rearrangement in both wild-type and snail twist mutant embryos, where our theoretical prediction is further improved when we also account for cell packing disorder. These findings suggest that convergent extension is associated with a transition to more fluid-like tissue behavior, which may help accommodate tissue shape changes during rapid developmental events.