1
|
Garrett AS, Means SA, Roesler MW, Miller KJW, Cheng LK, Clark AR. Modeling and experimental approaches for elucidating multi-scale uterine smooth muscle electro- and mechano-physiology: A review. Front Physiol 2022; 13:1017649. [PMID: 36277190 PMCID: PMC9585314 DOI: 10.3389/fphys.2022.1017649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Accepted: 09/21/2022] [Indexed: 11/13/2022] Open
Abstract
The uterus provides protection and nourishment (via its blood supply) to a developing fetus, and contracts to deliver the baby at an appropriate time, thereby having a critical contribution to the life of every human. However, despite this vital role, it is an under-investigated organ, and gaps remain in our understanding of how contractions are initiated or coordinated. The uterus is a smooth muscle organ that undergoes variations in its contractile function in response to hormonal fluctuations, the extreme instance of this being during pregnancy and labor. Researchers typically use various approaches to studying this organ, such as experiments on uterine muscle cells, tissue samples, or the intact organ, or the employment of mathematical models to simulate the electrical, mechanical and ionic activity. The complexity exhibited in the coordinated contractions of the uterus remains a challenge to understand, requiring coordinated solutions from different research fields. This review investigates differences in the underlying physiology between human and common animal models utilized in experiments, and the experimental interventions and computational models used to assess uterine function. We look to a future of hybrid experimental interventions and modeling techniques that could be employed to improve the understanding of the mechanisms enabling the healthy function of the uterus.
Collapse
|
2
|
López Garza G, Castellanos NP, Godínez R. Cell-to-cell modeling of the interface between atrial and sinoatrial anisotropic heterogeneous nets. Comput Biol Chem 2017; 68:245-259. [PMID: 28460307 DOI: 10.1016/j.compbiolchem.2017.04.008] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Accepted: 04/17/2017] [Indexed: 10/19/2022]
Abstract
The transition between sinoatrial cells and atrial cells in the heart is not fully understood. Here we focus on cell-to-cell mathematical models involving typical sinoatrial cells and atrial cells connected with experimentally observed conductance values. We are interested mainly in the geometry of the microstructure of the conduction paths in the sinoatrial node. We show with some models that appropriate source-sink relationships between atrial and sinoatrial cells may occur according to certain geometric arrangements.
Collapse
Affiliation(s)
- Gabriel López Garza
- Mathematics Department, Universidad Autónoma Metropolitana, M. City, Mexico.
| | - Norma P Castellanos
- Electric Engineering Department, Universidad Autónoma Metropolitana, M. City, Mexico
| | - Rafael Godínez
- Electric Engineering Department, Universidad Autónoma Metropolitana, M. City, Mexico
| |
Collapse
|
3
|
Sheldon RE, Mashayamombe C, Shi SQ, Garfield RE, Shmygol A, Blanks AM, van den Berg HA. Alterations in gap junction connexin43/connexin45 ratio mediate a transition from quiescence to excitation in a mathematical model of the myometrium. J R Soc Interface 2015; 11:20140726. [PMID: 25401181 DOI: 10.1098/rsif.2014.0726] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The smooth muscle cells of the uterus contract in unison during delivery. These cells achieve coordinated activity via electrical connections called gap junctions which consist of aggregated connexin proteins such as connexin43 and connexin45. The density of gap junctions governs the excitability of the myometrium (among other factors). An increase in gap junction density occurs immediately prior to parturition. We extend a mathematical model of the myometrium by incorporating the voltage-dependence of gap junctions that has been demonstrated in the experimental literature. Two functional subtypes exist, corresponding to systems with predominantly connexin43 and predominantly connexin45, respectively. Our simulation results indicate that the gap junction protein connexin45 acts as a negative modulator of uterine excitability, and hence, activity. A network with a higher proportion of connexin45 relative to connexin43 is unable to excite every cell. Connexin45 has much more rapid gating kinetics than connexin43 which we show limits the maximum duration of a local burst of activity. We propose that this effect regulates the degree of synchronous excitation attained during a contraction. Our results support the hypothesis that as labour approaches, connexin45 is downregulated to allow action potentials to spread more readily through the myometrium.
Collapse
|
4
|
Lammers WJEP, Stephen B, Al-Sultan MA, Subramanya SB, Blanks AM. The location of pacemakers in the uteri of pregnant guinea pigs and rats. Am J Physiol Regul Integr Comp Physiol 2015; 309:R1439-46. [PMID: 26377559 DOI: 10.1152/ajpregu.00187.2015] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2015] [Accepted: 09/01/2015] [Indexed: 01/31/2023]
Abstract
The pregnant uterus is a smooth muscle organ whose pattern of contraction is dictated by the propagation of electrical impulses. Such electrical activity may originate from one or more pacemakers, but the location of these sites has not yet been determined. To detect the location of the pacemaker in the gravid uterus, two approaches were used: 1) determine the site from where the contraction started using isolated uteri from the pregnant guinea pig, and videotape their contractions; and 2) record, in isolated uteri from pregnant term rats, with 240 extracellular electrodes simultaneously, and determine where the electrical bursts started. In both the contractile and electrophysiological experiments, there was not a single, specific pacemaker area. However, most contractions (guinea pig 87%) and bursts (rat 76%) started close to the mesometrial border (mean 2.7 ± 4.0 mm SD in guinea pigs and 1.3 ± 1.4 mm in rats). In addition, in the rat, most sites of initiations were located closer to the ovarial end of the horn (mean distance from the ovarial end 6.0 ± 6.2 mm SD), whereas such an orientation was not seen in the guinea pig. In both guinea pig and rat uteri at term, there is not one specific pacemaker area. Rather, contractile and electrical activity may arise from any site, with the majority starting close to the mesometrial border. Furthermore, in the rat, most activities started at the ovarial end of the horn. This may suggest a slightly different pattern of contraction in both species.
Collapse
Affiliation(s)
- Wim J E P Lammers
- Department of Physiology, College of Medicine and Health Sciences, United Arab Emirates University, Al Ain, United Arab Emirates; and
| | - Betty Stephen
- Department of Physiology, College of Medicine and Health Sciences, United Arab Emirates University, Al Ain, United Arab Emirates; and
| | - Mahmood Ahmed Al-Sultan
- Department of Physiology, College of Medicine and Health Sciences, United Arab Emirates University, Al Ain, United Arab Emirates; and
| | - Sandeep B Subramanya
- Department of Physiology, College of Medicine and Health Sciences, United Arab Emirates University, Al Ain, United Arab Emirates; and
| | - Andrew M Blanks
- Division of Translational and Systems Medicine, Warwick Medical School, Clinical Sciences Research Laboratory, Coventry, United Kingdom
| |
Collapse
|
5
|
Simulating uterine contraction by using an electro-chemo-mechanical model. Biomech Model Mechanobiol 2015; 15:497-510. [PMID: 26162461 DOI: 10.1007/s10237-015-0703-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2015] [Accepted: 07/02/2015] [Indexed: 10/23/2022]
Abstract
Contractions of uterine smooth muscle cells consist of a chain of physiological processes. These contractions provide the required force to expel the fetus from the uterus. The inclusion of these physiological processes is, therefore, imperative when studying uterine contractions. In this study, an electro-chemo-mechanical model to replicate the excitation, activation, and contraction of uterine smooth muscle cells is developed. The presented modeling strategy enables efficient integration of knowledge about physiological processes at the cellular level to the organ level. The model is implemented in a three-dimensional finite element setting to simulate uterus contraction during labor in response to electrical discharges generated by pacemaker cells and propagated within the myometrium via gap junctions. Important clinical factors, such as uterine electrical activity and intrauterine pressure, are predicted using this simulation. The predictions are in agreement with clinically measured data reported in the literature. A parameter study is also carried out to investigate the impact of physiologically related parameters on the uterine contractility.
Collapse
|
6
|
Sheldon RE, Baghdadi M, McCloskey C, Blanks AM, Shmygol A, van den Berg HA. Spatial heterogeneity enhances and modulates excitability in a mathematical model of the myometrium. J R Soc Interface 2013; 10:20130458. [PMID: 23843249 DOI: 10.1098/rsif.2013.0458] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The muscular layer of the uterus (myometrium) undergoes profound changes in global excitability prior to parturition. Here, a mathematical model of the myocyte network is developed to investigate the hypothesis that spatial heterogeneity is essential to the transition from local to global excitation which the myometrium undergoes just prior to birth. Each myometrial smooth muscle cell is represented by an element with FitzHugh-Nagumo dynamics. The cells are coupled through resistors that represent gap junctions. Spatial heterogeneity is introduced by means of stochastic variation in coupling strengths, with parameters derived from physiological data. Numerical simulations indicate that even modest increases in the heterogeneity of the system can amplify the ability of locally applied stimuli to elicit global excitation. Moreover, in networks driven by a pacemaker cell, global oscillations of excitation are impeded in fully connected and strongly coupled networks. The ability of a locally stimulated cell or pacemaker cell to excite the network is shown to be strongly dependent on the local spatial correlation structure of the couplings. In summary, spatial heterogeneity is a key factor in enhancing and modulating global excitability.
Collapse
Affiliation(s)
- Rachel E Sheldon
- MOAC Doctoral Training Centre, University of Warwick, Coventry CV4 7AL, UK.
| | | | | | | | | | | |
Collapse
|
7
|
Benson AP, Aslanidi OV, Zhang H, Holden AV. The canine virtual ventricular wall: A platform for dissecting pharmacological effects on propagation and arrhythmogenesis. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2008; 96:187-208. [PMID: 17915298 DOI: 10.1016/j.pbiomolbio.2007.08.002] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
We have constructed computational models of canine ventricular cells and tissues, ultimately combining detailed tissue architecture and heterogeneous transmural electrophysiology. The heterogeneity is introduced by modifying the Hund-Rudy canine cell model in order to reproduce experimentally reported electrophysiological properties of endocardial, midmyocardial (M) and epicardial cells. These models are validated against experimental data for individual ionic current and action potential characteristics, and their rate dependencies. 1D and 3D heterogeneous virtual tissues are constructed, with detailed tissue architecture (anisotropy and orthotropy, due to fibre orientation and sheet structure) of the left ventricular wall wedge extracted from a diffusion tensor imaging data set. The models are used to study the effects of tissue heterogeneity and class III drugs on transmural propagation and tissue vulnerability to re-entry. We have determined relationships between the transmural dispersion of action potential duration (APD) and the vulnerable window in the 1D virtual ventricular wall, and demonstrated how changes in the transmural heterogeneity, and hence tissue vulnerability, can lead to generation of re-entry in the 3D ventricular wedge. Two class III drugs with opposite qualitative effects on transmural APD heterogeneity are considered: d-sotalol that increases transmural APD dispersion, and amiodarone that decreases it. Simulations with the 1D virtual ventricular wall show that under d-sotalol conditions the vulnerable window is substantially wider compared to amiodarone conditions, primarily in the epicardial region where unidirectional conduction block persists until the adjacent M cells are fully repolarised. Further simulations with the 3D ventricular wedge have shown that ectopic stimulation of the epicardial region results in generation of sustained re-entry under d-sotalol conditions, but not under amiodarone conditions or in control. Again, APD increase in M cells was identified as the major contributor to tissue vulnerability--re-entry was initiated primarily due to ectopic excitation propagating around the unidirectional conduction block in the M cell region. This suggests an electrophysiological mechanism for the anti- and proarrhythmic effects of the class III drugs: the relative safety of amiodarone in comparison to d-sotalol can be explained by relatively low transmural APD dispersion, and hence, a narrow vulnerable window and low probability of re-entry in the tissue.
Collapse
Affiliation(s)
- Alan P Benson
- Computational Biology Laboratory, Institute of Membrane and Systems Biology, Worsley Building, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK.
| | | | | | | |
Collapse
|
8
|
Garfield RE, Maner WL. Physiology and electrical activity of uterine contractions. Semin Cell Dev Biol 2007; 18:289-95. [PMID: 17659954 PMCID: PMC2048588 DOI: 10.1016/j.semcdb.2007.05.004] [Citation(s) in RCA: 145] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2007] [Accepted: 05/03/2007] [Indexed: 11/22/2022]
Abstract
Presently, there is no effective treatment for preterm labor. The most obvious reason for this anomaly is that there is no objective manner to evaluate the progression of pregnancy through steps leading to labor, either at term or preterm. Several techniques have been adopted to monitor labor, and/or to diagnose labor, but they are either subjective or indirect, and they do not provide an accurate prediction of when labor will occur. With no method to determine preterm labor, treatment might never improve. Uterine electromyography (EMG) methods may provide such needed diagnostics.
Collapse
Affiliation(s)
- Robert E Garfield
- University of Texas Medical Branch, Department of Obstetrics and Gynecology, Division of Reproductive Sciences, 301 University, Route 1062, Galveston, TX 77555, United States.
| | | |
Collapse
|
9
|
Benson AP, Halley G, Li P, Tong WC, Holden AV. Virtual cell and tissue dynamics of ectopic activation of the ventricles. CHAOS (WOODBURY, N.Y.) 2007; 17:015105. [PMID: 17411262 DOI: 10.1063/1.2404634] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Cardiac ventricular cells and tissues are normally excitable, and are activated by propagating waves of excitation that are initiated in the specialized pacemaking region of the heart. However, isolated or repetitive activity can be initiated at abnormal (ectopic) sites in the ventricles. To trigger an endogenous ectopic beat, there must be a compact focus of cells with changed membrane excitation parameters and kinetics, which initiate activity by after-depolarizations triggered by propagating activity, or that have bifurcated into autorhythmicity. This ectopic focus needs to be large enough, and adequately coupled, to drive the surrounding tissue. We investigate the initiation of ectopic excitation in computational models of human ventricular cells triggered by after-depolarizations and by up/down regulation of specific membrane conductance systems, and the propagation and evolution of ectopic activity in homogeneous or heterogeneous and isotropic, anisotropic, or orthotropic tissues.
Collapse
Affiliation(s)
- Alan P Benson
- Computational Biology Laboratory, Institute of Membrane and Systems Biology, and Cardiovascular Research Institute, Worsley Building, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| | | | | | | | | |
Collapse
|
10
|
Bernard S, Cajavec B, Pujo-Menjouet L, Mackey MC, Herzel H. Modelling transcriptional feedback loops: the role of Gro/TLE1 in Hes1 oscillations. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1155-70. [PMID: 16608701 DOI: 10.1098/rsta.2006.1761] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
The transcriptional repressor Hes1, a basic helix-loop-helix family protein, periodically changes its expression in the presomitic mesoderm. Its periodic pattern of expression is retained in a number of cultured murine cell lines. In this paper, we introduce an extended mathematical model for Hes1 oscillatory expression that includes regulation of Hes1 transcription by Drosophila Groucho (Gro) or its vertebrate counterpart, the transducine-like enhancer of split/Groucho-related gene product 1 (TLE1). Gro/TLE1 is a necessary corepressor required by a number of DNA-binding transcriptional repressors, including Hes1. Models of direct repression via Hes1 typically display an expression overshoot after transcription initiation which is not seen in the experimental data. However, numerical simulation and theoretical predictions of our model show that the cofactor Gro/TLE1 reduces the overshoot and is thus necessary for a rapid and finely tuned response of Hes1 to activation signals. Further, from detailed linear stability and numerical bifurcation analysis and simulations, we conclude that the cooperativity coefficient (h) for Hes1 self-repression should be large (i.e. h>or=4). Finally, we introduce the characteristic turnaround duration, and show that for our model the duration of the repression loop is between 40 and 60 min.
Collapse
Affiliation(s)
- Samuel Bernard
- Institute for Theoretical Biology, Humboldt University, Invalidenstrasse 43, 10115 Berlin, Germany
| | | | | | | | | |
Collapse
|
11
|
Biktasheva IV, Simitev RD, Suckley R, Biktashev VN. Asymptotic properties of mathematical models of excitability. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1283-98. [PMID: 16608708 DOI: 10.1098/rsta.2006.1770] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
We analyse small parameters in selected models of biological excitability, including Hodgkin-Huxley (Hodgkin & Huxley 1952 J. Physiol.117, 500-544) model of nerve axon, Noble (Noble 1962 J. Physiol.160, 317-352) model of heart Purkinje fibres and Courtemanche et al. (Courtemanche et al. 1998 Am. J. Physiol.275, H301-H321) model of human atrial cells. Some of the small parameters are responsible for differences in the characteristic time-scales of dynamic variables, as in the traditional singular perturbation approaches. Others appear in a way which makes the standard approaches inapplicable. We apply this analysis to study the behaviour of fronts of excitation waves in spatially extended cardiac models. Suppressing the excitability of the tissue leads to a decrease in the propagation speed, but only to a certain limit; further suppression blocks active propagation and leads to a passive diffusive spread of voltage. Such a dissipation may happen if a front propagates into a tissue recovering after a previous wave, e.g. re-entry. A dissipated front does not recover even when the excitability restores. This has no analogy in FitzHugh-Nagumo model and its variants, where fronts can stop and then start again. In two spatial dimensions, dissipation accounts for breakups and self-termination of re-entrant waves in excitable media with Courtemanche et al. kinetics.
Collapse
Affiliation(s)
- I V Biktasheva
- Department of Computer Science, University of Liverpool Liverpool, L69 3BX, UK.
| | | | | | | |
Collapse
|
12
|
Gavaghan D, Garny A, Maini PK, Kohl P. Mathematical models in physiology. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1099-106. [PMID: 16608698 DOI: 10.1098/rsta.2006.1757] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Computational modelling of biological processes and systems has witnessed a remarkable development in recent years. The search-term (modelling OR modeling) yields over 58000 entries in PubMed, with more than 34000 since the year 2000: thus, almost two-thirds of papers appeared in the last 5-6 years, compared to only about one-third in the preceding 5-6 decades. The development is fuelled both by the continuously improving tools and techniques available for bio-mathematical modelling and by the increasing demand in quantitative assessment of element inter-relations in complex biological systems. This has given rise to a worldwide public domain effort to build a computational framework that provides a comprehensive theoretical representation of integrated biological function-the Physiome. The current and next issues of this journal are devoted to a small sub-set of this initiative and address biocomputation and modelling in physiology, illustrating the breadth and depth of experimental data-based model development in biological research from sub-cellular events to whole organ simulations.
Collapse
Affiliation(s)
- David Gavaghan
- Oxford University Computing Laboratory, Wolfson Building, Parks Road, Oxford OX1 3QD, UK.
| | | | | | | |
Collapse
|
13
|
Calcium Oscillations and Ectopic Beats in Virtual Ventricular Myocytes and Tissues: Bifurcations, Autorhythmicity and Propagation. ACTA ACUST UNITED AC 2005. [DOI: 10.1007/11494621_31] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|