Original Research ARTICLE
Computational modeling of inhibition of voltage-gated Ca channels: identification of different effects on uterine and cardiac action potentials
- Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, UK
The uterus and heart share the important physiological feature whereby contractile activation of the muscle tissue is regulated by the generation of periodic, spontaneous electrical action potentials (APs). Preterm birth arising from premature uterine contractions is a major complication of pregnancy and there remains a need to pursue avenues of research that facilitate the use of drugs, tocolytics, to limit these inappropriate contractions without deleterious actions on cardiac electrical excitation. A novel approach is to make use of mathematical models of uterine and cardiac APs, which incorporate many ionic currents contributing to the AP forms, and test the cell-specific responses to interventions. We have used three such models—of uterine smooth muscle cells (USMC), cardiac sinoatrial node cells (SAN), and ventricular cells—to investigate the relative effects of reducing two important voltage-gated Ca currents—the L-type (ICaL) and T-type (ICaT) Ca currents. Reduction of ICaL (10%) alone, or ICaT (40%) alone, blunted USMC APs with little effect on ventricular APs and only mild effects on SAN activity. Larger reductions in either current further attenuated the USMC APs but with also greater effects on SAN APs. Encouragingly, a combination of ICaL and ICaT reduction did blunt USMC APs as intended with little detriment to APs of either cardiac cell type. Subsequent overlapping maps of ICaL and ICaT inhibition profiles from each model revealed a range of combined reductions of ICaL and ICaT over which an appreciable diminution of USMC APs could be achieved with no deleterious action on cardiac SAN or ventricular APs. This novel approach illustrates the potential for computational biology to inform us of possible uterine and cardiac cell-specific mechanisms. Incorporating such computational approaches in future studies directed at designing new, or repurposing existing, tocolytics will be beneficial for establishing a desired uterine specificity of action.
Computational modeling of an action potential (AP) of an electrically excitable cell was first developed in 1952 with the landmark study of neurons (Hodgkin and Huxley, 1952). Its success led to the development of other models such as the tonic AP in cardiac cells (Noble, 1960) and the bursting AP in β-pancreatic cells (Chay and Keizer, 1983). In the intervening years there has been an enormous advance in our understanding of the cardiac physiome (Noble, 2007; Schmitz et al., 2011; Noble et al., 2012) and computational analysis of electrical excitability now runs hand-in-hand with physiological experimentation in heart research (Crampin et al., 2004; Bassingthwaighte et al., 2009; Masumiya et al., 2009; Nikolaidou et al., 2012; Zhang et al., 2012). Many computational models exist to describe in considerable detail cardiac cell-specific excitation-contraction properties, including the biophysical details of the constituent ion currents and calcium fluxes. These include multicellular tissue and organ maps of spatiotemporal electrical and calcium wave propagation (Rudy, 2000; Zhang et al., 2000; Kleber and Rudy, 2004; Severi et al., 2009; Aslanidi et al., 2011b; Atkinson et al., 2011). Mathematical models are continuously being developed and applied to predicting the risks of pathophysiological phenomena (e.g., the likelihood of dyssynchronous activation and fibrillation) (Benson et al., 2008; Bishop and Plank, 2012; Cherry et al., 2012; Kharche et al., 2012; Behradfar et al., 2014) as well as the potential beneficial effects of drugs and treatments (Levin et al., 2002; Muzikant and Penland, 2002; Davies et al., 2012; Mirams et al., 2012; di Veroli et al., 2013, 2014).
The uterus is also an electrically excitable tissue whose contractile function is determined by episodic spontaneous APs and calcium fluxes. Although our comprehension of the electrophysiological basis of uterine AP formation lags behind that of cardiac muscle, there is an increasing awareness that computational approaches, such as those which have been applied so extensively to cardiac muscle, may foster advances in this matter (Taggart et al., 2007; Aslanidi et al., 2011a; Tong et al., 2011; Sharp et al., 2013). For example, we have developed a biophysically-detailed uterine smooth muscle cell (USMC) model validated against experimental data and it can describe many different uterine AP forms and the corresponding intracellular calcium changes (Tong et al., 2011).
Utilizing computational models for the examination of uterine APs offers the additional possibility of predicting the likely actions of drugs that target ion channels or exchangers with the intention of attenuating premature uterine contractions. Hitherto, such tocolytics have been used clinically without prior extensive in silico assessment of effectiveness. Uterine APs, and the resultant contraction of smooth muscle cells, are markedly dependent upon the activation of a prominent voltage-gated inward (depolarizing) current, the long-lasting L-type calcium current (ICaL). Nifedipine, an L-type voltage-gated calcium channel blocker, is currently used as a tocolytic treatment (RCOG, 2011). This treatment is useful for delaying labor in the short term (Conde-Agudelo et al., 2011). However, it is not without adverse side effects, in particular on maternal cardiovascular performance (van Veen et al., 2005; Guclu et al., 2006; Gaspar and Hajagos-Toth, 2013), and it is not presently recommended for longer-term use, nor for women with cardiac disease (RCOG, 2011). This alerts one to another necessary consideration of tocolytic drugs intended to limit uterine APs namely, what possible actions may there be on cardiac electrical excitability?
Another voltage-gated Ca current, in addition to ICaL, that, theoretically, may be a suitable target for developing tocolytic drugs is the short-lasting, transient T-type Ca channel current (ICaT). T-type calcium currents are believed to play a role in pacemaking in many cell types (Mangoni et al., 2006; Perez-Reyes et al., 2009). ICaT has been observed in uterine tissues (Inoue et al., 1990; Young et al., 1991, 1993; Blanks et al., 2007) and putative blockers of ICaT reduce in vitro uterine contractions (Lee et al., 2009). ICaT has a different current-voltage (I-V) profile from ICaL and, unlike the latter, there is considered to be little ICaT present in ventricular cardiomyocytes although it has been suggested to contribute to sinoatrial (SA) node APs (Ono and Iijima, 2010; Mesirca et al., 2014).
Our over-arching objective was to utilize computational models of uterine and cardiac cells to theoretically test the possible cell-specific effects of inhibiting voltage-gated Ca entry (ICaL and ICaT) in a manner that may be anticipated to occur with drugs targeting these pathways. Using publicly-available computational models of uterine and cardiac APs we have performed a series of simulation experiments to investigate the following questions:
(i) Are there comparable effects on uterine and cardiac APs of reducing ICaL?
(ii) Are there similar effects on uterine and cardiac APs of reducing ICaT?
(iii) Does combined reduction of ICaL and ICaT have differential actions on uterine and cardiac APs?
Uterine Smooth Muscle Cell Model
For these simulation studies we used our previously published USMC model (Tong et al., 2011). A schematic of the ionic currents contributing to the model is shown in Figure 1A. The model source code, parameter values and the full description are provided in Tong et al. (2011). All the USMC AP simulations were started at the same resting state. The numerical values of all the dynamical variables at this state—the initial conditions—are provided in the Supplementary Materials.
Figure 1. Schematic diagrams of the ionic currents in each of the three cell models. (A) Uterine smooth muscle cell model (USMC) (Tong et al., 2011), (B) rabbit sinoatrial nodal cell model (SAN) (Garny et al., 2003), (C) modified guinea pig ventricular cell model (LRd00) (Faber and Rudy, 2000). The lollipop sign indicates channels that are voltage-gated. The L-type Ca channel and T-type Ca channel, carrying ICaL and ICaT, respectively, are indicated in red.
Cardiac Cell Models
There are numerous computational models of cardiac cells and many are listed in the CellML model repository (http://models.cellml.org/). They represent different parts of the heart from the central and peripheral sinoatrial nodal cells (Zhang et al., 2000; Garny et al., 2003), the atrial cardiomyocytes (Aslanidi et al., 2009), the atrial-ventricular node and His-bundle (Inada et al., 2009), Purkinje fiber cells (Stewart et al., 2009; Li and Rudy, 2011) and ventricular cardiomyocytes (Faber and Rudy, 2000). These cell types can be grouped, rather roughly, into those of a contractile or conducting function. In our preliminary simulations using these cardiac models, we selected the ventricular and sinoatrial nodal cell models that showed the greatest propensity to alter AP form following a reduction of Ca current thus forming the most sensitive situation for comparing the same maneuvre to that on the USMC model.
For the sinoatrial nodal cell, we used the Garny et al. (2003) rabbit sinoatrial nodal cell model (SAN) and its 0D-capable version configurations for the SAN central cell. A schematic of the ionic currents contributing to the model is shown in Figure 1B. The model equations and parameter values and are listed in the CellML model repository. The initial conditions for the SAN model are provided in the Supplementary Materials.
For the ventricular cardiomyocyte cell, we used as a platform the Faber and Rudy (2000) guinea pig cell model (LRd00) and its M cell configurations. A schematic of the ionic current components of the model is shown in Figure 1C. The model equations, parameter values and initial conditions are listed in the CellML model repository. A copy of the model source code can also be downloaded from http://rudylab.wustl.edu/research/cell/methodology/cellmodels/LRd/code.htm. However, as detailed below in the Results, we had to modify this model in order to incorporate a more robust description of ICaT. The LRd00 modified model was paced for 20 min at 2 Hz to allow it to reach a stable state with its M cell configurations. Then, the values of all the dynamical variables in this stable state were saved: these stable state conditions of the modified LRd00 model are provided in the Supplementary Materials.
Simulations of Uterine and Cardiac APs with Reduced ICaL and ICaT
The effects of reduced ICaL and ICaT on the uterine and cardiac APs were assessed and quantified. The ICaL and ICaT were reduced by multiplying their maximal conductances with a scaling factor between 0 and 1 and we assumed the same proportional reduction would occur in all three cell types. As the APs between the three unique cell types have different forms and characteristics, it is difficult to use a common measure to assess their AP behaviors. Therefore, we chose to assess a characteristic that best reflects the main function of each of the cell types: the bursting activities in the uterine cell, the pacemaking ability of the SAN cell and the action potential duration (APD) of the ventricular cell. For the USMC model, after adjusting the levels of ICaL and ICaT, an AP was induced by a 5 s current clamp at −0.5 pA pF−1 and the maximal peak membrane voltage (Vpeak) reached after the first initial spike was measured. This assessment gives an indication of the presence or absence of bursting spikes in a USMC AP. For the SAN model, after adjusting the levels of ICaL and ICaT, the model was simulated for 10 s and the pacemaking frequency was measured from the APs in the last 2 s of simulations. For the LRd00 model, an AP was induced by a 0.5 ms current pulse at −90 pA pF−1 and the APD at 90% repolarization level (APD90) was assessed.
All simulations were computed using XPPAUT Version 6 (Ermentrout, 2002) in a Dell Optiplex PC with an Intel(R) Core(TM) i7-3770 CPU @ 3.40 GHz. For the USMC model, all simulations were computed with a fixed time step of 0.02 ms and the Euler method. For the SAN model, all simulations were computed with a fixed time step of 0.1 ms and the fourth-order Runge–Kutta method. For the LRd00 model, all simulations were computed with a fixed time step of 0.002 ms and the fourth-order Runge–Kutta method.
Validation of the Cell Models
We first verified each source code of the three individual cell models with their respective published results. Both the USMC and the SAN models were validated in this regard. Surprisingly, we found that the LRd00 model could not be validated so in respect of ICaT: the simulated current tracings (Figure 2A), and the corresponding current-voltage (I-V) relationship (Figure 2B), of the LRd00 model did not match the experimental results that they were rated against (Balke et al., 1992). The simulations produced current responses at membrane voltages that were much farther to the right (positive) than observed in the experimental data or, indeed, to that anticipated from other ICaT models. When we traced back the original formulation of this ICaT to Zeng et al. (1995), and the experimental values from which its kinetics details were based to Droogmans and Nilius (1989), we did not find any typographical errors in the model equations or the experimental values to explain the discrepancy. As an aim of our study was to investigate the effects of altering ICaT, it was essential that the models incorporated the correct formulations. We were left, therefore, with no option but to attempt to modify the ICaT in the LRd00 model such that it did reflect more closely the experimental data. For those interested in the details of this process, our attempts to modify the model were undertaken as follows: we had first tried, but failed, to produce a reasonable ICaT model by either adjusting the LRd00 ICaT equations or reformulating new equations using data from Droogmans and Nilius (1989). We next tried substituting the LRd00 ICaT with a validated ICaT model from Li and Rudy (2011), which was developed for canine Purkinje fiber cells. The simulated current tracings with this canine ICaT resemble the ICaT current tracings in Balke et al. (1992) but the kinetics were too fast and the peak of the I-V relationship was still too positive compared to the experimental data from Balke et al. (1992). Using the ICaT from Li and Rudy (2011) as our template, we corrected the differences between the simulations and the experimental data with the modifications that resulted in Equations (1–8) for the new ventricular ICaT that are listed in the Supplementary Materials. As discussed by Droogmans and Nilius (1989), we also found that the kinetics of the ventricular ICaT currents were best described by a gating product of b2g where b and g are gating variables for activation and inactivation (Equation 1 in the Supplementary Materials). The I-V data in Balke et al. (1992) was matched with a half-activation at −50 mV in the activation steady-state function (Equation 3 in the Supplementary Materials). Although this half-activation value deviated from the reported value for ICaT (Droogmans and Nilius, 1989), it is similar to those reported values from clonal Cav3.1 expression data (Serrano et al., 1999; Hering et al., 2003). As the kinetics of the canine ICaT model were too fast compared to the ventricular cell experimental data within the LRd00 model (i.e., the data from Balke et al., 1992), we slowed down the activation and inactivation by scaling the time constants (Equations 5, 6 in the Supplementary Materials). The resultant new ICaT model satisfactorily described the ventricular ICaT data (Balke et al., 1992) and we replaced the original LRd00 ICaT with this new ICaT model. The ratio between peak ICaT and ICaL in guinea pig ventricular cells are reported to be between 0.13 and 0.19 (Balke et al., 1992; Masumiya et al., 1998; Zorn-Pauly et al., 2004). Therefore, we chose a value for the maximal conductance of our new ICaT model (gCaT) so that the ratio of the maximum ICaT: ICaL in the modified LRd00 cell model was 0.15.
Figure 2. Incorporation of a modified T-type calcium current model for LRd00. The original ICaT contained in the LRd00 cardiac ventricular cell model did not behave like a T-type calcium current. The simulated raw current tracings (A) and the I-V current plots (B) appeared at holding membrane potentials farther to the right than of the experimental data of Balke et al. (1992) that was used as a source experimental dataset. Simulated current tracings at five different voltage steps from a holding potential of −90 mV (A) and simulated peak current-voltage (I-V) relationships (B) of the original and modified ICaT currents are shown in comparison to the experimental data from Balke et al. (1992). All current tracings were normalized to their maximal peak current from their I-V relationships. The modified ICaT now closely resembled the experimental data. (Experimental tracings and data, adapted with permission from Figure 4 from Balke et al. (1992); copyright 1992 The Physiological Society.)
Normal USMC, SAN, and Ventricular APs
USMC APs can take a variety of forms of variable durations. Bursting type USMC AP, wherein an initial spike of an AP followed by rapid and repetitive fluctuations in membrane potential that persist for the whole AP duration, are often observed. These may assist in the maintenance of long-duration contractions of many tens of seconds that are a feature of the uterus during labor. Therefore, a promising feature of a tocolytic would be one that could dampen, or delay, the bursting in the USMC APs with minimum effects on the APs of the cardiac SAN and the ventricular cells.
Figure 3A shows the simulated APs from the three cell models under control conditions before any maneuvre. The USMC model showed a resting membrane potential (RMP) of −55.5 mV before stimulation. During the evoked USMC AP, burstings occurred throughout the AP with a frequency between 2 and 2.54 Hz and the amplitudes of the repetitive spikes were around 40 mV. The Vpeak after the initial spike was −7 mV. The SAN cell is autorhythmic where periodic APs occur without external stimulation. The frequency and amplitude of the pacemaking activities were 3.15 Hz and 75 mV (between −56 and 19 mV), respectively. These values were the same throughout the whole 10 s simulation. The LRd00 ventricular cell showed a comparatively hyperpolarized (negative) RMP at −87 mV. A single AP was evoked in respond to a stimulus and the ventricular APD90 was 165 ms. These quantified characteristics are the reference values before any maneuvre.
Figure 3. Effects of reducing ICaL in uterine and cardiac cells. Effects on the action potential forms and characteristics with reduced ICaL in different cell types: the USMC cell model (left), the SAN cell model (center), and the LRd00 ventricular cell model (right). (A) Control configuration at full ICaL and the effects of reducing ICaL by (B) 10% and (C) 20%. The negative sign indicated reduction.
Effects of Reducing ICaL in Uterine and Cardiac Cells
ICaL is a major depolarization current in USMCs and thus a logical target for tocolytics. We examined the effect of reducing ICaL by 10% (Figure 3B). In the USMC model, this change dampened the initial repetitive spikes of the USMC AP and reduced the duration of the burstings by around 60%. The amplitudes of the remaining, late-onset, bursting spikes were reduced by 30% to about 28 mV and the Vpeak of the USMC AP was lowered to −14.8 mV. A less severe effect of 10% less ICaL was evident in the APs from the SAN and the LRd00 models. The amplitude of the SAN pacemaking APs was modestly reduced by 12% and the frequency remained unchanged. The LRd00 APD90 was slightly reduced by 5%. When ICaL was reduced by 20% (Figure 3C), the burstings of the USMC AP were suppressed completely. However, the SAN cell also stopped pacemaking. The LRd00 APD90 was slightly reduced by ~10% at 150 ms. These data indicated that although reducing some ICaL suppressed uterine bursting APs, this maneuvre also affected both the cardiac pacemaker APs and the ventricular AP.
Effects of Reducing ICaT in Uterine and Cardiac Cells
ICaT may be one of the SAN pacemaking currents in the heart and may also be involved in the generation of APs in USMCs. The effects of our maneuvres to reduce ICaT in each of the three cell types are shown in Figure 4. The control cases of the three models are displayed in Figure 4A. Reduction of ICaT by 40% (Figure 4B) delayed the onset of the burstings in the USMC AP. However, both the Vpeak and the amplitude of the later-onset burstings remained the same as the control. With 40% less ICaT, the pacemaking by the SAN cell slowed down to about 2.94 Hz with the peak potential at a slightly lower level of 15 mV. Reduction of ICaT by 40% did not affect the LRd00 APD90. When more ICaT was reduced, by 80% (Figure 4C), the onset of the USMC burstings was further delayed with only one spike appearing at the end of the AP. However, the amplitude of this spike was similar in size to those of the control. Also, the RMP level became more hyperpolarized at ~−57 mV. With this large reduction of ICaT (by 80%), the SAN pacemaking APs slowed down to 2.63–2.7 Hz. The pacemaking potential in the SAN model remained at the same level at −56 mV but the peak potential was lowered to 12 mV. At this level, the SAN model also showed a longer (~0.5 s) transient during simulation. With 80% less ICaT, the LRd00 APD90 remained unchanged. Without ICaT (Figure 4D), the RMP level of the USMC became more hyperpolarized at ~−58 mV and the bursting of the USMC AP stopped. The SAN pacemaking also ceased without ICaT but the LRd00 AP was unaffected. These data indicated that reduction of ICaT suppressed the uterine bursting AP with no impact on the ventricular AP. However, depending on the magnitude of ICaT reduction, there was also an effect on the SAN APs.
Figure 4. Effects of reducing ICaT in uterine and cardiac cells. Effects on the action potential forms and characteristics with reduced ICaT in different cell types: the USMC cell model (left), the SAN cell model (center), and the LRd00 ventricular cell model (right). (A) Control configuration at full ICaT and the effects of reducing ICaT by: (B) 40%, (C) 80%, and (D) 100%. The negative sign indicated reduction.
Effects of Reducing ICaL and ICaT in Uterine and Cardiac Cells
The above data suggested that there was some promise in reducing either ICaL or ICaT and effecting an inhibition of USMC APs with little deleterious action on ventricular APs. The manouevres that reduced ICaL or ICaT and had mild effects on SAN AP may be tolerable. However, reductions in either current that effected a larger inhibition of uterine APs, as one may desire the action of a tocolytic in our experimental setting to be, also had dramatic actions on SAN AP form. The outcomes of reducing ICaL or ICaT in the USMC models were slightly different: reducing ICaT altered the RMP and the onset of spike bursting whereas changes in ICaL had a predominant action on spike amplitude. This begat the question, could a combined reduction of ICaL and ICaT exert the desired action on USMC APs with little impact on cardiac cell APs? An example is illustrated in Figure 5. The control cases of the three models are shown in Figure 5A and the effects of reducing only ICaL by 10%, or only ICaT by 40%, are illustrated again in Figures 5B,C, respectively. This enables comparison to the effect arising when both ICaL and ICaT were reduced, by 10 and 40%, respectively, (Figure 5D). In this case, the burstings of the USMC AP were completely suppressed (Figure 5D). With this combination, the pacemaking frequency of the SAN APs remained the same at 3.125 Hz and the amplitude was reduced by 20% (from −52 to 7.7 mV). The LRd00 APD90 was slightly reduced by ~4% to 158 ms. Compared to the cases when only either ICaL or ICaT were reduced alone to completely suppress the burstings of the USMC AP (Figures 3C, 4D), this combined approach clearly performed better while preserving many of the properties of the cardiac APs.
Figure 5. Effects of reducing ICaL and ICaT in uterine and cardiac cells. Effects on the action potential forms and characteristics when both ICaL and ICaT were proportionally reduced in different cell types: the USMC cell model (left), the SAN cell model (center), and the LRd00 ventricular cell model (right). (A) Control configuration before reducing either of the currents. (B) Reducing ICaL by 10% only. (C) Reducing ICaT by 40% only. (D) Reducing ICaL by 10% and ICaT by 40%. The negative sign indicated reduction.
We plotted two-parameter maps of percentage current inhibitions—of ICaL vs. ICaT—in each cell model to examine the effective ranges of different combinations of ICaL and ICaT reductions on the properties of uterine and cardiac APs (Figure 6). After the initial spike of a USMC AP, the peak of the bursting spikes typically reach values between −20 and 0 mV whereas without bursting, the plateau AP voltage usually stayed below −30 mV. Therefore, we assessed the presence or absence of burstings in a USMC AP from the level of Vpeak when the proportions of ICaL and ICaT were changed and the results of the parameter plots are color-coded in Figure 6. In Figure 6A, the white region indicated the combinations of ICaL and ICaT that would produce a USMC AP with some burstings and the shaded region indicated the combinations of currents with which the USMC cell would not produce such activity. The USMC AP example shown previously in Figure 5D, with 10% reduction of ICaL and 40% reduction of ICaT, lay just within the non-bursting domain (as indicated by the green star). For the SAN cell, the pacemaker frequency was monitored as one changed the proportions of ICaL and ICaT and the results are shown in Figure 6B. The colored domain indicated pacemaking occurrence and the white region suggested the parameter combinations would result in no pacemaking. The frequencies of the pacemaking APs were similar throughout the pacemaking zone and the transition between the two regions was steep (not shown). The SAN example shown in Figure 5D with 10% reduction of ICaL and 40% reduction of ICaT was within the pacemaking region. For the LRd00 ventricular cell, the APD90 was monitored as the proportions of ICaL and ICaT were changed and the results are shown in Figure 6C. Changes in the APD90 occurred only with changes in ICaL, indicating that ICaT did not influence the APD of the ventricular cell model. From overlapping these parameter maps, we derived the boundaries of combinations of ICaL and ICaT for attaining specific sets of AP forms in the uterine and cardiac cells as shown in Figure 6D. The colored lines separate different properties of the uterine and cardiac cells obtained from the simulations in Figures 6A–C. The blue line traces the level of Vpeak at −30 mV in the parameter space of ICaL and ICaT of the USMC cell model and it separates the bursting and quiescent conditions. The red line separates the pacemaking conditions with ICaL and ICaT in the SAN cell model from the non-pacemaking. We allowed for 10% variations in the APD90 of the LRd00 ventricular cells and this threshold at 150 ms was indicated by the gray line. These three criteria were superimposed on the same parameter map of ICaL and ICaT and revealed an overlapping area that satisfied all three conditions. In other words, any combination of reduced ICaL and ICaT within this area can produce our desired outcome for tocolytic treatment, namely, effective suppression of USMC APs with minimal effects on cardiac APs.
Figure 6. Summarized effects of different combinations of reduced ICaL and ICaT on the APs of the three different cell types. The shaded areas in the 2-parameter maps indicate the type of AP that would be produced with the different combinations of proportionally reduced ICaL and ICaT: (A) suppressed bursting in the USMC cell model, (B) pacemaking activities in the SAN cell model, and (C) less than 10% change in the action potential duration at 90% repolarization level (APD90) in the LRd00 ventricular cell model. The green star illustrates the corresponding position on each 2-parameter map of the example shown in Figure 5D, i.e., with reductions of 10% ICaL and 40% ICaT. (D) Combinations of ICaL and ICaT (green area) that would suppress bursting in the USMC without affecting the pacemaking in the SAN cell or the APD90 in the LRd00 ventricular cell.
The above experiments have considered the outcomes on AP form when imposing the same proportional reduction of ICaL and ICaT in all three cell types. However, another interesting possibility would be the same absolute quantities of currents being reduced in all three cell types. This consideration arises from examining the range of current densities of ICaL and ICaT reported from the different cell types as noted in Table 1. Cells of the cardiac conduction system tend to have much more ICaT than contractile cardiomyocytes and, when comparing between uterine and cardiac cell types, the total ICaT found in USMCs amounts to around 42% of that reported from guinea pig SAN. We re-performed our simulations to test the effects of reducing the same quantities of currents in each of the three cell type models. This resulted in a transformation of the XY axes of the parameter maps in Figure 6 to current densities and the result is shown in Figure 7. Now, the simulations reveal a larger domain within which one may reduce ICaL and ICaT to produce effective inhibition of USMC model APs with little detriment to cardiac cell model APs. In particular, it is possible to completely block the USMC bursting by large reduction of the USMC ICaT without affecting the SAN pacemaking. Therefore, our first experimental protocol of examining the effects of proportional current reductions likely represents the most stringent case.
Figure 7. Effects of the same absolute amounts of ICaL and ICaT reductions on the APs of the three different cell types. A transformation of the simulation results in Figure 6D using the maximum values from published experimental current densities of ICaL and ICaT (Table 1) in human uterine smooth muscle cells (ICaL = 4.5 pA pF−1 and ICaT = 2.32 pA pF−1), guinea pig SAN cells (ICaL = 4.6 pA pF−1 and ICaT = 5.4 pA pF−1), and guinea pig ventricular cells (ICaL = 12.5 pA pF−1 and ICaT = 1.97 pA pF−1). The dashed lines are the experimental limits of ICaL and ICaT that can be reduced in the different cell types. The solid color lines and the green area represented the same conditions as in Figure 6.
Muscle contractions are determined by the form and duration of the preceding APs and the corresponding calcium fluxes. ICaL is involved in the depolarization of APs and is a major source of calcium influx in many excitable cell types including uterine smooth muscle and cardiac muscle cells. As it has been shown in experimental studies (Kodama et al., 1997; Terrar et al., 2007; Lee et al., 2009; Young and Bemis, 2009) and in this study, reducing ICaL affects the amplitude and the duration of the AP and/or contractile force in these cell types. These properties made ICaL an obvious choice as a target for tocolytic treatment. Indeed, nifedipine is used clinically as a second line tocolytic for short periods and within tight dosage recommendations; caution in its use is warranted because of the possibility of side effects including palpitations and hypotension (van Geijn et al., 2005; RCOG, 2011). This alerts one to a necessary consideration of the possible actions on cardiac electrical excitability by tocolytic drugs intended to limit uterine APs.
From our simulations, the USMC model seems to be more susceptible to reduced ICaL than the two cardiac (SAN and ventricular) cell models. However, the influence on the cardiac SAN model of reducing ICaL appears rather steep and so, as one attempts to increase the USMC effects by greater ICaL reductions, there too is a potentially damaging influence on SAN activity and, by extension, heart pacing. This may relate to the side effects on cardiac functions and the tight dosage guidelines recommended for nifedipine use as a tocolytic (RCOG, 2011).
With its potential role in facilitating the onset of USMC APs, ICaT is an alternative tocolytic target. The RMP of USMCs of around −55 mV is close to the voltage level for the window current from ICaT (Blanks et al., 2007). However, experimental studies in cardiac SAN cells illustrate that maneuvres designed to reduce ICaT result in alterations of SAN APs (Hagiwara et al., 1988; Masumiya et al., 1998; Tanaka et al., 2008). In our study, modest reduction of ICaT hyperpolarizes the USMC model RMP and diminishes the bursting spikes with a mild effect on the SAN model cell AP. However, as above for ICaL, further reductions in ICaT with the purpose of abrogating USMC model APs resulted in deleterious actions on the SAN cell model APs.
In contrast, our data suggests that reducing both ICaL and ICaT can confer advantages for achieving our objective. This strategy may facilitate delaying/slowing the onset of the USMC APs. In our simulations, this combined approach clearly works more effectively than reducing either of the currents alone. Furthermore, we have identified a range of desirable combinations of ICaL and ICaT inhibition for maximizing the desirable action on USMC model APs without attracting an adverse effect on cardiac SAN or ventricular cell model APs.
If this approach is to work, then there will have to be an agent that blocks ICaT with at least as good selectivity as nifedipine for ICaL. Mibefradil has been suggested to be a selective T-type Ca channel blocker but there are several studies indicating that it can affect both L-type and T-type calcium channels (Masumiya et al., 1999; Protas and Robinson, 2000; de Paoli et al., 2002). From such a seeming disadvantage, could one utilize the cross-channel inhibition of a particular drug for the purposes described above? i.e., might the use of one drug eventuate the situation of beneficially inhibiting both ICaL and ICaT? Mibefradil, at 1 μM, was shown to block 55% of ICaT in rabbit sinonatrial nodal cells (Protas and Robinson, 2000) but also and 64% of ICaL; at 3 μM, it blocked 28% of ICaT and 15% of ICaL in guinea pig ventricular cells (de Paoli et al., 2002); at 10 μM, it inhibited 90% ICaT in guinea pig ventricular cells (Masumiya et al., 1999) and also 40% ICaL. If we compare these combinations of proportional inhibitions of ICaL and ICaT against our results in Figure 6, all of these combinations fall outside our desired proportional mix of these two currents for a useful tocolytic treatment. Instead, based on our simulation results, if mibefradil was used as tocolytic treatments at these concentrations, it might reduce uterine bursting APs but it would likely affect cardiac pacemaking functions as well. Indeed, mibefradil, at 1 μM, was shown to attenuate the contractile forces of uterine muscle strips from late pregnant rats (IC50 ~1 μM, Lee et al., 2009). However, at the same concentration, mibefradil also reversibly stopped the pacemaking activities from sinoatrial nodal cells (Protas and Robinson, 2000). Thus, it would seem that the dual actions, and limited channel-specific discretion, of mibefradil on ICaL and ICaT would not be advantageous to our purpose in both in silico and in vitro experimental settings. In addition, it has been reported to have inhibitory actions on a number of non-Ca channels including INa(1.5), a major Na current in cardiomyocytes (McNulty and Hanck, 2004). This highlights another issue, therefore, in the prosecution of these studies in the longer-term, i.e., that the development of inhibitors with greater selectivity for T-type Ca channels would be a major advance.
In this in silico assessment of potential tocolycis, not only did we consider the effectiveness of action on the uterus, but we also assessed the possible side effects to different cells of the heart. This approach is only possible as the computational models of both organs are developed with sufficient biophysical details. The most obvious limitation for this kind of assessment is that it depends on the depth and accuracy of biophysical details in the computational models. We can only include quantitative details from experimental data but often we are forced to make assumptions and educated guesses for the unknowns. Indeed, we have seen in our study that close scrutiny of existing models, however well established, can reveal anomalies that need addressing—in our case the ICaT characteristics in the ventricular cell model—before experimental simulations with particular purposes can be tested. This is the benefit of the iterative process. As computational models continue to be tested, and evolved to incorporate new and relevant biophysical details, the accuracy of the quantitative predictions will improve.
Our work brings to the fore the notion that, instead of searching for a single tocolytic compound with high specificity, there may be merit in considering a cocktail of more than one compound. The present recommendations in the UK are to avoid combinations of tocolytic drugs for fear of increasing the risk of side effects (RCOG, 2011). However, this caution arises from a lack of depth to the background research data and the paucity of currently available tocolytic options. The progression to labor has been described as a modular accumulation of (patho)physiological systems—or MAPS—and it is such modularity that is a major challenge to overcome in seeking to prevent or inhibit labor (Mitchell and Taggart, 2009). The suspicion is that once module A of the labor process is inhibited, there is module B, or module C or more that eventually facilitates a similar function and outcome of labor and delivery. If this is so, then the possibilities offered by combination drug strategies must continue to be explored.
Our assessment only considered effects at the cell level. However, to fully evaluate the actions of any tocolytic compounds, we need to also consider their effects at the tissue and organ levels. Currently there is no validated computational organ model for the uterus that would serve well this purpose although it is recognized by many that efforts toward this are required (Aslanidi et al., 2011a; Sharp et al., 2013). We are under no illusion that there is a need for much more “wet” experimentation—particularly in USMCs—to furnish computational model improvement, and to test with increasing rigor and clarity, important hypotheses of relevance to tocolysis.
Our approach herein can be regarded as a useful platform to be built upon for assessing the potential of tocolytics that act upon ion channels or electrogenic ion exchangers. Using the in silico approach described above enables future research to assess in parallel the potential benefits of attenuating premature uterine contractions vs. the risks of deleterious actions on cardiac electrical excitation.
Encouraging such in silico assessments at the early stages of laboratory or clinical studies will foster an integrated and iterative process between mathematical modeling and experimentation such that each informs the other. Related to this, it can alert one to the possibility of otherwise unforeseen actions of drugs and inform the subsequent protocol development for experimentation. Not only will this be beneficial in the examination of new physiological mechanisms and the actions of novel drugs but it will be useful, as indicated in this study, for the investigation of repurposing the application of existing drugs. It is to be hoped that similar approaches may form part of forthcoming scientific and clinical scientific strategies.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by an MRC Bioinformatics Training Fellowship G0902091.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fphys.2014.00399/abstract
Aslanidi, O., Atia, J., Benson, A. P., van den Berg, H. A., Blanks, A. M., Choi, C., et al. (2011a). Towards a computational reconstruction of the electrodynamics of premature and full term human labour. Prog. Biophys. Mol. Biol. 107, 183–192. doi: 10.1016/j.pbiomolbio.2011.07.004
Aslanidi, O. V., Boyett, M. R., Dobrzynski, H., Li, J., and Zhang, H. (2009). Mechanisms of transition from normal to reentrant electrical activity in a model of rabbit atrial tissue: interaction of tissue heterogeneity and anisotropy. Biophys. J. 96, 798–817. doi: 10.1016/j.bpj.2008.09.057
Aslanidi, O. V., Colman, M. A., Stott, J., Dobrzynski, H., Boyett, M. R., Holden, A. V., et al. (2011b). 3D virtual human atria: a computational platform for studying clinical atrial fibrillation. Prog. Biophys. Mol. Biol. 107, 156–168. doi: 10.1016/j.pbiomolbio.2011.06.011
Atkinson, A., Inada, S., Li, J., Tellez, J. O., Yanni, J., Sleiman, R., et al. (2011). Anatomical and molecular mapping of the left and right ventricular His-Purkinje conduction networks. J. Mol. Cell. Cardiol. 51, 689–701. doi: 10.1016/j.yjmcc.2011.05.020
Balke, C. W., Rose, W. C., Marban, E., and Wier, W. G. (1992). Macroscopic and unitary properties of physiological ion flux through T-type Ca2+channels in guinea-pig heart cells. J. Physiol. 456, 247–265.
Benson, A. P., Aslanidi, O. V., Zhang, H., and Holden, A. V. (2008). The canine virtual ventricular wall: a platform for dissecting pharmacological effects on propagation and arrhythmogenesis. Prog. Biophys. Mol. Biol. 96, 187–208. doi: 10.1016/j.pbiomolbio.2007.08.002
Bishop, M. J., and Plank, G. (2012). The role of fine-scale anatomical structure in the dynamics of reentry in computational models of the rabbit ventricles. J. Physiol. 590, 4515–4535. doi: 10.1113/jphysiol.2012.229062
Blanks, A. M., Zhao, Z.-H., Shmygol, A., Bru-Mercier, G., Astle, S., and Thornton, S. (2007). Characterization of the molecular and electrophysiological properties of the T-type calcium channel in human myometrium. J. Physiol. 581, 915–926. doi: 10.1113/jphysiol.2007.132126
Camara, A. K., Begic, Z., Kwok, W. M., and Bosnjak, Z. J. (2001). Differential modulation of the cardiac L- and T-type calcium channel currents by isoflurane. Anesthesiology 95, 515–524. doi: 10.1097/00000542-200108000-00038
Cherry, E. M., Fenton, F. H., and Gilmour, R. F. Jr. (2012). Mechanisms of ventricular arrhythmias: a dynamical systems-based perspective. Am. J. Physiol. Heart Circ. Physiol. 302, H2451–H2463. doi: 10.1152/ajpheart.00770.2011
Conde-Agudelo, A., Romero, R., and Kusanovic, J. P. (2011). Nifedipine in the management of preterm labor: a systematic review and metaanalysis. Am. J. Obstet. Gynecol. 204, 134.e1–134.e20. doi: 10.1016/j.ajog.2010.11.038
Crampin, E. J., Halstead, M., Hunter, P., Nielsen, P., Noble, D., Smith, N., et al. (2004). Computational physiology and the physiome project. Exp. Physiol. 89, 1–26. doi: 10.1113/expphysiol.2003.026740
Davies, M. R., Mistry, H. B., Hussein, L., Pollard, C. E., Valentin, J.-P., Swinton, J., et al. (2012). An in silico canine cardiac midmyocardial action potential duration model as a tool for early drug safety assessment. Am. J. Physiol. Heart Circ. Physiol. 302, H1466–H1480. doi: 10.1152/ajpheart.00808.2011
de Paoli, P., Cerbai, E., Koidl, B., Kirchengast, M., Sartiani, L., and Mugelli, A. (2002). Selectivity of different calcium antagonists on T- and L-type calcium currents in guinea-pig ventricular myocytes. Pharmacol. Res. 46, 491–497. doi: 10.1016/S1043661802002360
di Veroli, G. Y., Davies, M. R., Zhang, H., Abi-Gerges, N., and Boyett, M. R. (2013). High-throughput screening of drug-binding dynamics to HERG improves early drug safety assessment. Am. J. Physiol. Heart Circ. Physiol. 304, H104–H117. doi: 10.1152/ajpheart.00511.2012
di Veroli, G. Y., Davies, M. R., Zhang, H., Abi-Gerges, N., and Boyett, M. R. (2014). hERG inhibitors with similar potency but different binding kinetics do not pose the same proarrhythmic risk: Implications for drug safety assessment. J. Cardiovasc. Electrophysiol. 25, 197–207. doi: 10.1111/jce.12289
Faber, G. M., and Rudy, Y. (2000). Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study. Biophys. J. 78, 2392–2404. doi: 10.1016/S0006-3495(00)76783-X
Garny, A., Kohl, P., Hunter, P. J., Boyett, M. R., and Noble, D. (2003). One-dimensional rabbit sinoatrial node models: benefits and limitations. J. Cardiovasc. Electrophysiol. 14, S121–S132. doi: 10.1046/j.1540.8167.90301.x
Gaspar, R., and Hajagos-Toth, J. (2013). Calcium channel blockers as tocolytics: principles of their actions, adverse effects and therapeutic combinations. Pharmaceuticals (Basel) 6, 689–699. doi: 10.3390/ph6060689
Guclu, S., Gol, M., Saygili, U., Demir, N., Sezer, O., and Baschat, A. A. (2006). Nifedipine therapy for preterm labor: effects on placental, fetal cerebral and atrioventricular doppler parameters in the first 48 hours. Ultrasound Obstet. Gynecol. 27, 403–408. doi: 10.1002/uog.2728
Inada, S., Hancox, J. C., Zhang, H., and Boyett, M. R. (2009). One-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells. Biophys. J. 97, 2117–2127. doi: 10.1016/j.bpj.2009.06.056
Inoue, Y., Nakao, K., Okabe, K., Izumi, H., Kanda, S., Kitamura, K., et al. (1990). Some electrical properties of human pregnant myometrium. Am. J. Obstet. Gynecol. 162, 1090–1098. doi: 10.1016/0002-9378(90)91322-4
James, A. F., Tominaga, T., Okada, Y., and Tominaga, M. (1996). Distribution of cAMP-activated chloride current and CFTR mRNA in the guinea pig heart. Circ. Res. 79, 201–207. doi: 10.1161/01.RES.79.2.201
Jones, K., Shmygol, A., Kupittayanant, S., and Wray, S. (2004). Electrophysiological characterization and functional importance of calcium-activated chloride channel in rat uterine myocytes. Pflugers Arch. 448, 36–43. doi: 10.1007/s00424-003-1224-7
Kharche, S., Adeniran, I., Stott, J., Law, P., Boyett, M. R., Hancox, J. C., et al. (2012). Pro-arrhythmogenic effects of the S140G KCNQ1 mutation in human atrial fibrillation—insights from modelling. J. Physiol. 590, 4501–4514. doi: 10.1113/jphysiol.2012.229146
Knock, G. A., and Aaronson, P. I. (1999). Calcium antagonistic properties of the cyclooxygenase-2 inhibitor nimesulide in human myometrial myocytes. Br. J. Pharmacol. 127, 1470–1478. doi: 10.1038/sj.bjp.0702685
Kodama, I., Nikmaram, M. R., Boyett, M. R., Suzuki, R., Honjo, H., and Owen, J. M. (1997). Regional differences in the role of the Ca2+ and Na+ currents in pacemaker activity in the sinoatrial node. Am. J. Physiol. 272, H2793–H2806.
Kojima, A., Ito, Y., Kitagawa, H., Matsuura, H., and Nosaka, S. (2014). Direct negative chronotropic action of desflurane on sinoatrial node pacemaker activity in the guinea pig heart. Anesthesiology 120, 1400–1413. doi: 10.1097/ALN.0000000000000165
Kojima, A., Kitagawa, H., Omatsu-Kanbe, M., Matsuura, H., and Nosaka, S. (2012). Inhibitory effects of sevoflurane on pacemaking activity of sinoatrial node cells in guinea-pig heart. Br. J. Pharmacol. 166, 2117–2135. doi: 10.1111/j.1476-5381.2012.01914.x
Kusaka, M., and Sperelakis, N. (1995). Inhibition of L-type calcium current by genistein, a tyrosine kinase inhibitor, in pregnant rat myometrial cells. Biochim. Biophys. Acta 1240, 196–200. doi: 10.1016/0005-2736(95)00191-3
Lee, S.-E., Ahn, D.-S., and Lee, Y.-H. (2009). Role of T-type Ca channels in the spontaneous phasic contraction of pregnant rat uterine smooth muscle. Korean J. Physiol. Pharmacol. 13, 241–249. doi: 10.4196/kjpp.2009.13.3.241
Levin, J. M., Penland, R. C., Stamps, A. T., and Cho, C. R. (2002). Using in silico biology to facilitate drug development, in 'In Silico' Simulation of Biological Processes: Novartis Foundation Symposium, Vol. 247, eds G. Bock and J. A. Goode (Chichester: John Wiley & Sons, Ltd). doi: 10.1002/0470857897.ch18
Li, P., and Rudy, Y. (2011). A model of canine purkinje cell electrophysiology and Ca(2+) cycling: rate dependence, triggered activity, and comparison to ventricular myocytes. Circ. Res. 109, 71–79. doi: 10.1161/CIRCRESAHA.111.246512
Mangoni, M. E., Couette, B., Marger, L., Bourinet, E., Striessnig, J., and Nargeot, J. (2006). Voltage-dependent calcium channels and cardiac pacemaker activity: from ionic currents to genes. Prog. Biophys. Mol. Biol. 90, 38–63. doi: 10.1016/j.pbiomolbio.2005.05.003
Masumiya, H., Kase, J., Tanaka, Y., Tanaka, H., and Shigenobu, K. (1999). Effects of mibefradil, a selective T-type Ca2+ channel antagonist, on sino-atrial node and ventricular myocardia. Res. Commun. Mol. Pathol. Pharmacol. 104, 321–329.
Masumiya, H., Oku, Y., and Okada, Y. (2009). Inhomogeneous distribution of action potential characteristics in the rabbit sino-atrial node revealed by voltage imaging. J. Physiol. Sci. 59, 227–241. doi: 10.1007/s12576-009-0032-z
Masumiya, H., Shijuku, T., Tanaka, H., and Shigenobu, K. (1998). Inhibition of myocardial L- and T-type Ca2+ currents by efonidipine: possible mechanism for its chronotropic effect. Eur. J. Pharmacol. 349, 351–357. doi: 10.1016/S0014-2999(98)00204-0
Mirams, G. R., Davies, M. R., Cui, Y., Kohl, P., and Noble, D. (2012). Application of cardiac electrophysiology simulations to pro-arrhythmic safety testing. Br. J. Pharmacol. 167, 932–945. doi: 10.1111/j.1476-5381.2012.02020.x
Mitchell, B. F., and Taggart, M. J. (2009). Are animal models relevant to key aspects of human parturition? Am. J. Physiol. Regul. Integr. Comp. Physiol. 297, R525–R545. doi: 10.1152/ajpregu.00153.2009
Miyoshi, H., Urabe, T., and Fujiwara, A. (1991). Electrophysiological properties of membrane currents in single myometrial, cells isolated from pregnant rats. Pflugers Arch. 419, 386–393. doi: 10.1007/BF00371121
Nikolaidou, T., Aslanidi, O. V., Zhang, H., and Efimov, I. R. (2012). Structure-function relationship in the sinus and atrioventricular nodes. Pediatr. Cardiol. 33, 890–899. doi: 10.1007/s00246-012-0249-0
Okabe, K., Inoue, Y., and Soeda, H. (1999). Estradiol inhibits Ca2+ and K+ channels in smooth muscle cells from pregnant rat myometrium. Eur. J. Pharmacol. 376, 101–108. doi: 10.1016/S0014-2999(99)00353-2
Ono, K., and Iijima, T. (2005). Pathophysiological significance of T-type Ca2+ channels: properties and functional roles of T-type Ca2+ channels in cardiac pacemaking. J. Pharmacol. Sci. 99, 197–204. doi: 10.1254/jphs.FMJ05002X2
Perez-Reyes, E., Van Deusen, A. L., and Vitko, I. (2009). Molecular pharmacology of human Cav3.2 T-type Ca2+ channels: block by antihypertensives, antiarrhythmics, and their analogs. J. Pharmacol. Exp. Ther. 328, 621–627. doi: 10.1124/jpet.108.145672
Rosati, B., Dun, W., Hirose, M., Boyden, P. A., and McKinnon, D. (2007). Molecular basis of the T- and L-type Ca2+ currents in canine Purkinje fibres. J. Physiol. 579, 465–471. doi: 10.1113/jphysiol.2006.127480
Royal College of Obstetricians and Gynaecologists (RCOG). (2011). Tocolysis for Women in Preterm Labour. RCOG Green-Top Guideline No. 1b. Available online at: http://www.rcog.org.uk/womens-health/clinical-guidance/tocolytic-drugs-women-preterm-labour-green-top-1b
Schmitz, J. P. J., Vanlier, J., van Riel, N. A. W., and Jeneson, J. A. L. (2011). Computational modeling of mitochondrial energy transduction. Crit. Rev. Biomed. Eng. 39, 363–377. doi: 10.1615/CritRevBiomedEng.v39.i5.20
Severi, S., Corsi, C., Rocchetti, M., and Zaza, A. (2009). Mechanisms of beta-adrenergic modulation of I(Ks) in the guinea-pig ventricle: insights from experimental and model-based analysis. Biophys. J. 96, 3862–3872. doi: 10.1016/j.bpj.2009.02.017
Shmigol, A. V., Eisner, D. A., and Wray, S. (1998). Properties of voltage-activated [Ca2+]i transients in, single smooth muscle cells isolated from pregnant rat uterus. J. Physiol. 511, 803–811. doi: 10.1111/j.1469-7793.1998.803bg.x
Stewart, P., Aslanidi, O. V., Noble, D., Noble, P. J., Boyett, M. R., and Zhang, H. (2009). Mathematical models of the electrical action potential of Purkinje fibre cells. Philos. Trans. A Math. Phys. Eng. Sci. 367, 2225–2255. doi: 10.1098/rsta.2008.0283
Taggart, M. J., Blanks, A., Kharche, S., Holden, A., Wang, B., and Zhang, H. (2007). Towards understanding the myometrial physiome: approaches for the construction of a virtual physiological uterus. BMC Pregnancy Childbirth, 7 Suppl. 1:S3. doi: 10.1186/1471-2393-7-S1-S3
Takeda, K., Yamagishi, R., Masumiya, H., Tanaka, H., and Shigenobu, K. (2004). Effect of cilnidipine on L- and T-type calcium currents in guinea pig ventricle and action potential in rabbit sinoatrial node. J. Pharmacol. Sci. 95, 398–401. doi: 10.1254/jphs.SCJ04001X
Tanaka, H., Komikado, C., Namekata, I., Nakamura, H., Suzuki, M., Tsuneoka, Y., et al. (2008). Species difference in the contribution of T-type calcium current to cardiac pacemaking as revealed by r(-)-efonidipine. J. Pharmacol. Sci. 107, 99–102. doi: 10.1254/jphs.SC0070405
Terrar, D. A., Wilson, C. M., Graham, S. G., Bryant, S. M., and Heath, B. M. (2007). Comparison of guinea-pig ventricular myocytes and dog Purkinje fibres for in vitro assessment of drug-induced delayed repolarization. J. Pharmacol. Toxicol. Methods 56, 171–185. doi: 10.1016/j.vascn.2007.04.005
Tong, W.-C., Choi, C. Y., Kharche, S., Holden, A. V., Zhang, H., and Taggart, M. J. (2011). A computational model of the ionic currents, Ca2+ dynamics and action potentials underlying contraction of isolated uterine smooth muscle. PLoS ONE 6:e18685. doi: 10.1371/journal.pone.0018685
van Veen, A. J., Pelinck, M. J., van Pampus, M. G., and Erwich, J. J. H. M. (2005). Severe hypotension and fetal death due to tocolysis with nifedipine. BJOG 112, 509–510. doi: 10.1111/j.1471-0528.2004. 00480.x
Young, R. C., and Bemis, A. (2009). Calcium-activated chloride currents prolongs the duration of contractions in pregnant rat myometrial tissue. Reprod. Sci. 16, 734–739. doi: 10.1177/1933719109334965
Young, R. C., Smith, L. H., and Anderson, N. Jr. (1991). Passive membrane properties and inward calcium current of human uterine smooth muscle cells. Am. J. Obstet. Gynecol. 164, 1132–1139. doi: 10.1016/0002-9378(91)90601-M
Young, R. C., Smith, L. H., and McLaren, M. D. (1993). T-type and L-type calcium currents in freshly dispersed human uterine smooth muscle cells. Am. J. Obstet. Gynecol. 169, 785–792. doi: 10.1016/0002-9378(93)90006-5
Zeng, J., Laurita, K. R., Rosenbaum, D. S., and Rudy, Y. (1995). Two components of the delayed rectifier K+ current in ventricular myocytes of the guinea pig type. Theoretical formulation and their role in repolarization. Circ. Res. 77, 140–152. doi: 10.1161/01.RES.77.1.140
Zhang, H., Butters, T., Adeniran, I., Higham, J., Holden, A. V., Boyett, M. R., et al. (2012). Modeling the chronotropic effect of isoprenaline on rabbit sinoatrial node. Front. Physiol. 3:241. doi: 10.3389/fphys.2012.00241
Zhang, H., Holden, A. V., Kodama, I., Honjo, H., Lei, M., Varghese, T., et al. (2000). Mathematical models of action potentials in the periphery and center of the rabbit sinoatrial node. Am. J. Physiol. Heart Circ. Physiol. 279, H397–H421.
Keywords: computational modeling, uterus, cardiac, calcium channel, tocolytics
Citation: Tong W-C, Ghouri I and Taggart MJ (2014) Computational modeling of inhibition of voltage-gated Ca channels: identification of different effects on uterine and cardiac action potentials. Front. Physiol. 5:399. doi: 10.3389/fphys.2014.00399
Received: 13 June 2014; Accepted: 26 September 2014;
Published online: 16 October 2014.
Edited by:Raheela N. Khan, University of Nottingham, UK
Reviewed by:Reka Albert, Pennsylvania State University, USA
Andreas Bergdahl, Concordia University, Canada
Copyright © 2014 Tong, Ghouri and Taggart. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wing-Chiu Tong and Michael J. Taggart, Institute of Cellular Medicine, Newcastle University, 3rd Floor, William Leech Building, Medical School Framlington Place Newcastle upon Tyne, NE2 4HH, UK e-mail: email@example.com; firstname.lastname@example.org