Next Article in Journal
Rapamycin Plus Doxycycline Combination Affects Growth Arrest and Selective Autophagy-Dependent Cell Death in Breast Cancer Cells
Next Article in Special Issue
In Silico Reconstruction of Sperm Chemotaxis
Previous Article in Journal
Isoform-Specific Properties of Orai Homologues in Activation, Downstream Signaling, Physiology and Pathophysiology
Previous Article in Special Issue
The Organization of the Golgi Structures during Drosophila Male Meiosis Requires the Citrate Lyase ATPCL
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 15
10.3390/ijms22158018
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Occurrence of Calcium Oscillations in Human Spermatozoa Is Based on Spatial Signaling Enzymes Distribution

by
Julia Korobkin
1,
Fedor A. Balabin
1,2,
Sergey A. Yakovenko
3,4,
Ekaterina Yu. Simonenko
3 and
Anastasia N. Sveshnikova
1,2,3,5,*
1
Center for Theoretical Problems of Physico-Chemical Pharmacology, Russian Academy of Sciences, 30 Srednyaya Kalitnikovskaya, 109029 Moscow, Russia
2
National Medical Research Center of Pediatric Hematology, Oncology and Immunology Dmitry Rogachev, 1 Samory Mashela St., 117198 Moscow, Russia
3
Faculty of Physics, Lomonosov Moscow State University, 1/2 Leninskie, 119991 Moscow, Russia
4
AltraVita IVF and Gynecology Clinic, 4A Nagornaya, 117186 Moscow, Russia
5
Department of Normal Physiology, Sechenov First Moscow State Medical University, 8/2 Trubetskaya St., 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(15), 8018; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22158018
Submission received: 14 June 2021 / Revised: 15 July 2021 / Accepted: 21 July 2021 / Published: 27 July 2021
(This article belongs to the Special Issue Molecular Mechanisms of Sperm Activation)
Browse Figures
Review Reports Versions Notes

Abstract

:
In human spermatozoa, calcium dynamics control most of fertilization events. Progesterone, present in the female reproductive system, can trigger several types of calcium responses, such as low-frequency oscillations. Here we aimed to identify the mechanisms of progesterone-induced calcium signaling in human spermatozoa. Progesterone-induced activation of fluorophore-loaded spermatozoa was studied by fluorescent microscopy. Two computational models were developed to describe the spermatozoa calcium responses: a homogeneous one based on a system of ordinary differential equations and a three-dimensional one with added space dimensions and diffusion for the cytosolic species. In response to progesterone, three types of calcium responses were observed in human spermatozoa: a single transient rise of calcium concentration in cytosol, a steady elevation, or low-frequency oscillations. The homogenous model provided qualitative description of the oscillatory and the single spike responses, while the three-dimensional model captured the calcium peak shape and the frequency of calcium oscillations. The model analysis demonstrated that an increase in the calcium diffusion coefficient resulted in the disappearance of the calcium oscillations. Additionally, in silico analysis suggested that the spatial distribution of calcium signaling enzymes governs the appearance of calcium oscillations in progesterone-activated human spermatozoa.

1. Introduction

A spermatozoon is a male reproductive cell. Successful fertilization requires triggering of several events in spermatozoa. These include the acrosome reaction, an exocytotic event that results in eversion of the acrosomal vesicle [1], and the hyperactivation, a type of sperm motility characterized by an asymmetrical beating pattern of the sperm tail [2,3].
Calcium is an important second messenger regulating many cellular activities, including fertilization, secretion, neurotransmission, and cell migration [4]. When extracellular signals appear, a cell can increase its cytosolic Ca2+ concentration hundreds of times in milliseconds [5]. Calcium signaling appears to be in the center of sperm activation [6,7,8,9]. In the 1980s, it was demonstrated that cytosolic calcium rise is necessary for the acrosomal reaction [10]. Sperm cytosolic calcium rise was demonstrated in response to zona pellucida recombinant glycoproteins [11], progesterone [10,12], and other types of stimulations [9]. Moreover, spontaneous calcium oscillations occur in up to 50% of the spermatozoa population [13].
In human sperm, intracellular Ca2+ ions play a crucial role in both signal transduction and cell homeostasis. Calcium homeostasis is controlled mainly by the plasma membrane calcium ATPase (PMCA) and the Na+/Ca2+ exchanger (NCX), and calcium plays a crucial role in activating mitochondrial respiration [14]. Mitochondrial abnormalities found in PMCA-deficient sperm [15] suggest a Ca2+ overload due to defective Ca2+ extrusion. However, in bull sperm, mitochondrial respiration is not up-regulated by Ca2+ release [16], and mitochondrial uncoupling does not significantly affect the calcium response occurring in human sperm [17,18].
Progesterone is an endogenous steroid sex hormone involved in the menstrual cycle, pregnancy, and embryogenesis of humans and other species. It is present in the female reproductive tract in nanomolar concentrations and in the human oocyte corona radiata in micromolar concentrations [3].
Different types of responses to progesterone treatment were reported in mammalian spermatozoa—such as single-peak response [6,19,20,21], steady calcium concentration elevation [19,20,21], or calcium oscillations with a period of hundreds of seconds [21,22]. In human sperm, progesterone can trigger acrosome reaction [1] and hyperactivation as well as promote changes in the tyrosine phosphorylation status [22]. Low concentrations of progesterone (10–100 nM) induce sperm motility and activate tyrosine kinases; higher concentrations (1–10 μM) are required to induce acrosome reaction in human spermatozoa [23]. In murine spermatozoa, a transitory increase in intracellular calcium may promote acrosomal exocytosis [19]. In human spermatozoa generating oscillations, the flagella beat mode alternates in synchrony with the oscillation cycle [17]. Low-frequency calcium oscillations observed in a low percentage of human spermatozoa are known to inhibit acrosome reaction [24]. The mechanism of progesterone-induced calcium signaling is not completely clear, although it is known that progesterone modulates the activity of the CatSper (Cation channel of sperm) voltage-gated Ca2+ channel [6], and CatSper is crucial for calcium response [3,25]. In the resting state of the cell, CatSper is inhibited by a plasma membrane lipid, 2-arachidonoylglycerol (2-AG) [26]. Progesterone activates abhydrolase domain-containing protein 2 (ABHD2), an enzyme that cleaves 2-AG into arachidonic acid and glycerol, thus activating CatSper [27].
CatSper activation in spermatozoa induces a calcium signaling pathway, which is similar to those typical for non-excitable somatic cells. It starts with phospholipase C (PLC), which catalyzes inositol-1,4,5-triphosphate (IP3) production from a membrane phospholipid, phosphatidylinositol 4,5-bisphosphate (PIP2). IP3, an important secondary messenger, activates inositol-1,4,5-triphosphate channel-receptors (IP3R) located on the spermatozoon intracellular calcium stores, promoting calcium release into the cytoplasm (Figure 1) [28]. Various phosphatases can promote removal of IP3 and suppression of further signaling (Figure 1). However, there are several unique features in sperm calcium signaling. First, one of the most prominent isoforms of PLC in mammalian spermatozoa is PLCδ, activated solely by the rise of cytosolic calcium concentration [29]. It has the most significant activity level in mature cells amongst all PLC isoforms [30]. There is a lot of evidence that PLCδ participates in the progesterone-induced calcium response in human spermatozoa [31]. It may get activated when calcium concentration increases due to CatSper opening [27]. Second, the origin of the calcium stores in human sperm is not the endoplasmic reticulum, such as in most somatic cells, but there are two separate stores, the acrosome and the redundant nuclear envelope (RNE), a derivative of the nuclear envelope [32]. Interestingly, each of them is spatially distant from CatSper localization sites by approximately 5 μm [33].
To decode the signaling network described above, a computational modeling approach should be used. Several computational models of sperm cell calcium signaling [34,35,36,37], including progesterone-induced [7,38], have been previously developed. These models have successfully described calcium oscillations in sperm cells of various species as a result of positive feedback from calcium on different plasma membrane channels [34,36,37]. Olson et al. have developed a calcium signaling model in response to 8-Br-cAMP and proved the importance of calcium release from the intracellular stores for the calcium response [38], but this model could not describe the calcium oscillations. Li et al. have built a one-dimensional calcium signaling model in response to progesterone [7]. It could not describe the whole range of periods of calcium oscillations in human sperms varying from 40 to 600 s [13,17]. Another homogeneous model by Simons et al. [39] could describe the occurrence of calcium oscillations with a period about 100 s; however, the calcium concentrations in the model were much higher than the experimental ones (about 100–200 nM [40,41]), as were the basal calcium concentrations (50 nM) [42]. There are no studies describing spontaneous calcium oscillations in human sperm with a period of ~500 s. None of those models incorporates the new data on the 2-AG-dependent inhibition of CatSper published by Miller et al. [22].
Here we investigated the mechanisms of progesterone-induced calcium signaling in human spermatozoa using comprehensive computational modeling and single-cell fluorescent microscopy. We showed that the spatial distribution of the signaling enzymes plays an essential role in forming all types of sperm cell calcium responses to progesterone. We theoretically predicted that the distance between the PMCA localization site and calcium stores determines the lag-time of cytosolic calcium response. Additionally, we demonstrated that the rate of calcium diffusion in cytosol influences the period of calcium oscillations, significantly lowering their frequency.

2. Results

2.1. Human Sperm Cells Show Three Different Types of Calcium Responses to Progesterone

To determine the patterns of human sperm cell responses to progesterone, we performed fluorescent microscopy of calcium fluorophore-loaded single human sperm cells (see Methods and Supporting Information Text S1). Calcium concentration was measured using a ratiometric dye Fura-RED, which enabled comprehensive measurements of the signal from moving cells. The sperm heads were immobilized on poly-l-lysine, and the tails were motile and out of the focal plane. For single sperms whose heads were not immobilized properly, the measurement was not performed.
Several types of calcium responses were observed in progesterone-treated human spermatozoa: a single transient calcium increase with the mean width of 160 s and standard deviation of 44 s and mean magnitude of 0.5 μM with standard deviation of 0.2 μM followed by a baseline calcium elevation (Figure 2a–c), a low-frequency oscillatory response (Figure 2d–f), and a steady increase in calcium concentration (Figure 2g,h). In some cases, the transient response had a secondary lower calcium increase starting at 100 ± 20 s (Figure 2c) after progesterone stimulation, and the following calcium concentration was slightly elevated in comparison to non-activated cells (Figure 2b,c). In some cases, the main calcium peak was split (Figure 2b). The oscillatory response was of various types (Figure 2d–f), with a period varying from 80 to 300 s. The mean value was 193 s and the standard deviation was 76 s. Therefore, both the transient response and the calcium oscillations demonstrate a medium level of variance. Approximately half of the cells demonstrated single calcium rise response, while a quarter of cells demonstrated low-frequency oscillations. The total count of responsive cells and count of cells demonstrating each type of response are shown in Figure 2i.

2.2. The Three-Dimensional but Not the Homogeneous Model Provided a Quantitative Description of the Progesterone-Induced Calcium Response

Two computational models incorporating the main biochemical reactions (Figure 1) were constructed to describe the known signaling events (Figure 3a). The models’ numerical parameters (protein copy numbers per cell, turnover numbers, and Michaelis constants of the enzymes) were taken from the data on human proteins (see Supplementary Tables S1–S3 for details). Both homogeneous and three-dimensional models were constructed as a “virtual” sperm cell. To describe the experimentally observed calcium responses (Figure 2), the following parameters of the models were varied: the maximal permeability of CatSper (VCatSper) and the diffusion coefficient for calcium (Dca). The rationale behind the variation of VCatSper was that in the sperm sample incubated in capacitating conditions, there are at least two sperm subpopulations [43], which differ in their membrane potential that affects CatSper (cation channel) activity. Dca is affected by the capacity of the cytosolic buffer [44], particularly by buffering proteins such as 80K-H [45]. In turn, buffering capacity is affected by intracellular pH, which differs between the two subpopulations [43].
The homogenous model (0D) provided a qualitative description of both oscillatory (Figure 3b) and single-peak (Figure 3e) responses. In the 0D model, the maximum width of the calcium peak observed under physiologically reasonable calcium concentrations was 40 s, and the maximum period of calcium oscillations was 30 s. Therefore, the homogeneous model could not describe the experimentally observed peak width (160 ± 44 s) (Figure 3g), oscillations peak width (~100 s), and oscillations period (80–300 s) (Figure 3d) under any parameter value (see Supporting Information Text S11).
The three-dimensional model (3D) produced oscillations with a 250 s period and a 100 s width of each calcium peak (Figure 3c). The single peak response in the 3D model was 130 s wide (Figure 3f). These modeling results corresponded well with the experimental observations (Figure 3d,g).

2.3. Calcium Dynamics in Human Spermatozoa Is Strongly Dependent on the Calcium Diffusion Coefficient

Using the 3D model we developed, a theoretical investigation into the type of calcium response’s dependence on the calcium diffusion coefficient was carried out (Figure 4). It showed that the system might switch its type of response upon variation of calcium diffusion coefficient (Figure 4a). For a fixed diffusion coefficient of IP3 equal to 10 μm2/s, oscillations in the system appeared when the calcium diffusion coefficient was less than 70 μm2/s. In order to investigate the impact of the calcium transition between the principal piece and the head, an auxiliary quasi-heterogeneous model, in which cytosol was divided into two compartments, was built (see Figure S6). This model could describe the experimentally observed low-frequency oscillations (Figure S7), and the period of these oscillations was dependent on the rate coefficient for the calcium transition between the two model compartments (Figure S7). Taken together, these data suggested that some time delay in the movement of calcium from the head to the tail is required for the occurrence of low-frequency oscillations.
Interestingly, oscillations were observed in the 3D model for [IP3] < 0.017 µM (Figure 4e). For fixed [IP3] = 0.013 µM, the period of oscillations and the delay were increasing for lower calcium diffusion coefficients (Figure 4d).
In the main 3D model, the velocity of calcium extrusion by the PMCA appeared to affect the type of response in progesterone-activated human spermatozoa (Figure 4b). The oscillations were observed for PMCA Vmax values between 0.7 and 19.5 μM/s, and the oscillations period depended in a near-linear fashion on the velocity of calcium extrusion by PMCA (Figure 4b). In the model where PMCA was located both in the head and in the principal piece, the period of calcium oscillations was reduced to 60 s (Figure 4d, purple line). We also investigated the impact of the spatial distribution of PLCδ on the response type (Figure 4c, dashed line). In the main model, where the PLCδ activity was restricted to sperm head and MP, an oscillatory response occurred (Figure 4d, red line). For a uniformly distributed PLCδ activity (Figure 4d), a significant delay in calcium response appeared (Figure 4d, dotted line).

2.4. Mechanisms of Spontaneous Oscillations and Calcium Oscillations in Nominally Calcium-Free Medium Revealed by the Model

The heterogeneous model allowed us to take into account the diffusion of IP3 and calcium, thus leading to a description of mechanisms for several experimental phenomena.
First, a low frequency of spontaneous calcium oscillations was described previously [13]. The calcium oscillations with a period of ~500 s could be produced by the heterogeneous model (Figure 5a, dotted curve) with the period (Tsp) and the existence of oscillations depending on the calcium diffusion coefficient (Figure 5b). Interestingly, when the distance (L) between the sperm’s principal piece and the sperm’s store RNE was greater than 7.5 µm, the oscillations disappeared (Figure 5a, solid curve).
The model could also describe the behavior of the sperm cell calcium upon the addition of nominally calcium-free buffer (Figure 5c). In accordance with data from [17], in response to a calcium-free buffer, the frequency of calcium oscillations does not change, while the calcium baseline concentration was lower, and the amplitude of calcium oscillations was higher than in a medium that contained calcium. Details of the calculations of the calcium flux change are given in Supporting Information Text S9.

3. Discussion

Here we studied the mechanisms of progesterone-induced calcium signaling in human spermatozoa using computational modeling and single-cell microscopy. We observed three types of calcium responses: the single transient calcium peak with 160 ± 44 s width, the low-frequency calcium oscillations with a period of 80–300 s, and some of the cells responded to progesterone with steady calcium elevation. A quantitative description of all calcium response types specified above was performed, and two computational models were constructed based on the experimental data (a homogeneous and a three-dimensional one). Only the three-dimensional model could describe the low-frequency calcium oscillations. This type of response in sperm originated from the delay required for calcium to diffuse from the calcium store RNE located in the sperm head to PMCA localized in the sperm flagellum. Not only the shape but also the type of calcium response is affected by the spatial system configuration. It may change from transient to oscillatory upon calcium diffusion coefficient lowering.
The observed calcium mobilization responses of the sperm cell to progesterone correspond well with previously published data. Single peak parameters correspond to peak parameters observed by other authors, i.e., to those of calcium oscillations [19,20] and steady calcium elevation [19,20,21]. The mechanism of calcium oscillations proposed here is essentially the same as in the previous computational works for human [7,38] and other species cells [36,37]; that is, there is a positive feedback activation loop from cytosolic calcium on the plasma membrane or RNE calcium channels.
In this study, we demonstrated that the shape and the type of calcium response are affected by the spatial system configuration. Therefore, we can conclude that the longer the sperm cell, the less frequently it will undergo spontaneous oscillations (Figure 5a). This result correlates with the fact that the sperm length was found to be positively associated with semen characteristics [46]. Interestingly, in murine spermatozoa, which are significantly longer than human ones [47], the oscillatory response was observed more frequently than in human ones (9% [48] vs. 15% [19]). In addition, in our model, the oscillations frequency was higher for higher diffusion coefficients. This result is in accordance with the data from [17], that in the progesterone-activated spermatozoa heated to 37 °C, the oscillations frequency was almost two times higher than in those activated at room temperature, as diffusion rate is increased with the temperature increase. For higher calcium diffusion rates in the model, we observed a disappearance of calcium oscillations (Figure 4a), that have a direct physiological significance, as the temperature (and diffusion rates) in female reproductive tract (up to 37.2 °C) [49] is higher than in testis (around 33 °C) [50] and calcium oscillations prevent the acrosome reaction [13].
The dependence of the calcium response on the spatial system configuration is in accordance with the computational models of Pages et al. [51] and Lechleiter et al. [52]. They state that the shape of calcium response in acinar cells and Xenopus laevis oocytes may change significantly upon coupling with the fluid flow or accounting for calcium diffusion. The system’s geometry also plays an important role in dendritic cells in calcium response [53]. Low-frequency calcium oscillations may also be explained by calcium release delay generated by ER membrane potential [54], and such delay may be generated by diffusion as well. Interestingly, when the quadrilateral distribution of CatSper channels [36] was added to the model, the model response did not change (data not shown).
Approach limitations should also be taken into account. Mass-action law and Hill functions only approximately describe the activity of enzymes; Fick’s diffusion law is also a rough approximation because calcium diffusion in cells is non-linear [55]. Some of the numerical values (for example, a calcium leak coefficient from the RNE) could not be obtained from the literature and had to be estimated. Another possible limitation is that spermatozoa in different functional states may adhere differently to poly-l-lysine. This is discussed in Supporting Information Text S12.

4. Materials and Methods

4.1. Reagents

The sources of materials were as follows: calcium-sensitive cell-permeable fluorescent dyes Fura-2-AM, Fura-RED-AM, and Fluo-4-AM were from Molecular Probes (Eugene, OR, USA); progesterone and DMSO from Sigma-Aldrich (St. Louis, MO, USA); AllGrad® Wash buffer and AllGrad® 90% from Lifeglobal (Guilford, CT, USA). AllGrad® Wash buffer exact electrolyte and supplement composition (1.6 mM calcium, 0.22 mM phosphorus, 2.8 mM potassium, 108 mM chloride, 138 mM sodium mM, 0.24 mM magnesium, 0.18 mM glucose, 4.9 mM lactate, 0.24 mM pyruvate, 5% w/v HSA, gentamicin sulfate, Phenol Red, HEPES, and bicarbonate), was similar to the Global® Total® single-step medium from Lifeglobal (Guilford, CT, USA), as determined in [56].

4.2. Sperm Cells Collection

Sperm samples of 10 adult normozoospermic (according to the World Health Organization) donors were used for the study. The Ethics Committee approved this study for the AltraVita IVF Clinic, decision number #6 from 12 January 2020, and it was conducted according to the Helsinki Declaration’s ethical guidelines. Informed consent was obtained from all subjects. All volunteers were adequately informed about the usage of their clinical and biological data for the experiments before giving their consent.
To evaluate the spermatozoa parameters, the methods routinely used in the AltraVita IVF clinic were performed. Concentration and motility assay was done using a Makler chamber, and evaluation of morphological characteristics was performed using preliminarily stained (with methylene blue and cresyl violet) glass according to Kruger’s criteria (the number of normal spermatozoa is 60–150 million/mL with no less than 70% of them being motile (type a + b) and no less than 13% having a normal morphology).

4.3. Microscopy

Before analysis, spermatozoa were centrifuged using 2-layered AllGrad Wash and incubated for 3 h in bicarbonate and HEPES-buffered culture medium with 5% w/v HSA in tightly closed tubes as described in [57] to induce cells’ capacitation. For microscopy experiments, washed spermatozoa were incubated with 2 μM Fura-RED, AM, or Fluo-8 at 37 °C for 45 min. Cells were resuspended in AllGrad Wash buffer to an approximate density of 108 cells/mL. For epifluorescent microscopy, spermatozoa were immobilized on poly-l-lysine-covered glass slides and investigated in flow chambers [58]. An inverted Nikon Eclipse Ti-E microscope equipped with 100×/1.49 NA TIRF oil objective was used. Cell adhesion was controlled in DIC mode, and fluorescence dynamics were observed in epifluorescence mode. 405 nm and 488 nm lasers were applied sequentially to assess calcium-free Fura-RED and calcium-bound Fura-RED fluorescence in spermatozoa. 488 nm laser was applied in the case of Fluo-8 staining.
Nikon NIS-Elements software was used for microscope image acquisition; Fiji [59] was used for image processing.

4.4. Intracellular Calcium Measurement

To measure calcium concentration in sperm cells, the ratiometric Fura-RED-AM dye was used. The values of Fura-RED fluorescence were converted to cytosolic calcium concentration using the following formula [60]:
C a 2 + = K d F 488 m a x   − F 488 b a c k   F 488 m i n   − F 488 b a c k   R − R m i n   R − R m a x  
where Kd = 140 nM is the dissociation constant of the calcium-Fura-RED complex [61], R is the ratio of fluorescence upon excitation of Fura-RED with λ = 405 nm (calcium-bound dye) to λ = 488 nm (calcium-free dye). F488max and F488min are the maximum and minimum possible fluorescence intensities upon excitation with λ = 488 nm, F488back is the background fluorescence. Rmin and Rmax are the ratios upon the minimum and maximum possible calcium concentrations. Rmax was measured using a 0.025% Triton-X solution in the AllGrad Wash buffer. The maximum ratio detected in Triton-treated sperm cells was 1.9. After that, 5 μM of calcium chelator EGTA was introduced in the solution, and then the cells were incubated for 10 min at room temperature; the minimum ratio detected in spermatozoa was assessed to be 0.2. Sperm cells which exhibited rapid flagella movement were considered viable (Supporting Information Figure S10).

4.5. Computational Model Construction and Validation

Two computational models (a homogeneous and a three-dimensional) were developed to study progesterone-induced activation of human spermatozoa.

4.5.1. The Heterogeneous Model of Progesterone-Induced Activation of Human Spermatozoa

The heterogeneous (three-dimensional, 3D) model was constructed automatically in the VCell software [62] based on the 3D structure of a human sperm cell (Figure S2 depicts the central 2D slice of the cell) and a set of biochemical reactions underlying cell activation (Figure 1). The fluxes for the biochemical reactions were calculated based on laws of chemical kinetics (either mass action law, Henry–Michaelis–Menthen kinetics, or Hill functions). The parameters of the equations were taken from the literature on human proteins (Supplementary Materials Tables S2 and S3). The initial conditions were devised from the copy numbers of sperm proteins, those were estimated using pax-db protein abundance database [63] for human testis or mice sperm cells and data on total cell protein number (Supporting Information Text S8) [64]. The detailed model description, including model reactions and parameters, can be found in Supporting Information Text S3–S6. The full model is a set of 19 differential equations and 52 parameters, with 47 parameters obtained from the literature.
The model sperm cell 3D structure (called “Geometry”) was constructed in the following manner. The Geometry of the model consisted of three membrane compartments: extracellular volume, cytosol, and RNE. The volume and dimensions of the regions were estimated from the sperm microscopy data [32,46,65,66,67]; the details of the process are given in Figure S2, Table S1. The geometry of the three-dimensional model consists of extracellular space (300 µm3); sperm cytosol (18.4 µm3) divided into sperm tail (0.4 μm in diameter [65], 46 μm long [46]), sperm MP (0.8 μm in diameter [65], 7 μm long [66]) and sperm head (4.5 × 3 × 1.5 μm [66]); and sperm calcium store RNE [32] (0.11 µm3) located between sperm head and sperm MP. The diffusion coefficients and initial concentrations are given in Supplementary Materials Tables S1 and S3. Estimation of the five unknown parameters (Supplementary Materials Tables S1–S3) was conducted manually. A detailed description of parameter estimation is given in Supporting Information.
The heterogeneous model was constructed in a modular fashion; namely, two modules (Figure 1) were constructed, and their parameters were estimated independently, and then the modules were merged and validated against sperm cell data. The first module consisted of the CatSper channel, its activators and inhibitors. The module simulated the CatSper opening dynamics upon the treatment of spermatozoa with progesterone. The module included the CatSper channel protein, alpha-beta-hydrolase 2 (ABHD2), 2-arachidonoylglycerol, and arachidonic acid (AA). The localization of CatSper and ABHD2 was restricted to the PM area denoting the sperm cilia membrane, while the 2-arachidonoylglycerol was uniformly distributed in the plasma membrane. ABHD2 gets activated upon binding to a progesterone molecule from the extracellular space (it can be described by the mass-action law, Supplementary Materials Text S3). 2-AG molecules, which inhibit CatSper, get cleaved into glycerol and AA by activated ABHD2, promoting this calcium channel activity (described by the Michaelis–Menthen equation). ABHD2 activity was assumed to be described by Michaelis–Menten kinetics. The CatSper open state probability determined the flux of calcium ions through open CatSper channel into the sperm cytoplasm.
The second module (Calcium module, Figure 1) was focused on sperm cytosolic and store calcium ion concentration and IP3 diffusion and concentration dynamics. IP3 was considered to be synthesized by PLCδ and further dephosphorylated by IP3-5-phosphatase [68,69].
The PLCδ activity was considered to be dependent on calcium [29], while IP3 to activate IPP5P. See Supplementary Materials Text S5 for equations. The equations governing calcium concentration were similar to those in previously published models, except for the consideration of RNE calcium dynamics and use of a different IP3R model [70], and are given in Supplementary Information. The diffusion of IP3, calcium, and calcium buffers was taken into account using Fick’s law. Some of the reactions were restricted to spermatozoa’s specific sites (Figure 3a, details are given in Supplementary Materials Text S2 and Text S5); for example, IP3 generation by PLC occurred only in sperm MP and head. This assumption is based on PIP2 concentration being significantly lower in the principal piece due to voltage-sensitive phosphatase activity [71].

4.5.2. The Homogeneous Model of Progesterone-Induced Activation of Human Spermatozoa

The homogeneous (zero-dimensional, 0D) model of the sperm cell was based on the assumption of well-mixing of the proteins in all compartments of the model. Therefore, the heterogeneous model differed from the homogeneous model because it took into account the diffusion of calcium and IP3. The homogeneous model was constructed automatically in the COPASI software [72] based on the volumes of the compartments of the sperm cell, calculated for the 3D model (Figure S2) and the same set of the biochemical reactions underlying cell activation as the heterogeneous model (Figure 1). The details of the homogeneous model construction are given in Supporting Information, Text S6.

4.5.3. Considerations about the Inclusion of PLCδ into the Models

Human PLCδ catalytic parameters have not been characterized so far. Human sperm cells contain PLCδ [31]. Additional evidence for PLCδ presence and functioning in human sperm cells could be found in Supporting Information, Text S10. To prove the validity of PLCδ (with murine parameters) inclusion in the model, we performed the alignment of the amino acid sequences of murine and human PLCδ by means of BLAST protein alignment (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 15 October 2020), which demonstrated that these proteins are 80–90% similar; therefore, their kinetic parameters should not differ much (Figure S9A,B). Additionally, we performed sensitivity analysis [73] for the parameters of the PLCδ in the homogeneous model and received scaled sensitivities 0.0002 for KPLC and −2 × 10−9 for VPLC (parameter values are given in Supplementary Materials Table S3). This indicates that the possible inaccuracy of parameter estimation would not affect the results of the modeling. Additionally, the variation of the PLCδ parameters in the heterogeneous model (Figure S9) also did not have qualitative effects on the results, and the low-frequency calcium oscillations could be observed for fixed IP3 concentrations (Figure 4E), which correspond to the case when a phospholipase C with low calcium sensitivity (for example, PLCγ) is activated upon addition of progesterone. The acquired results suggest that the model behavior is independent of PLC parameters.

4.5.4. Integration of the Models

The set of partial differential equations was constructed in the VCell software [62] based on the scheme of biochemical reactions (Figure 1). This set with the initial variable values (Table S2) was integrated using the Method of Lines and CVODE in VCell [62]. The corresponding Virtual Cell Model, SpermCalcium, is available in the public domain at http://www.vcell.org/ (accessed on 17 October 2020) under the shared username “Juliajessica”. The homogeneous model was integrated with COPASI [72].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms22158018/s1, Text S1: Calcium concentration measurement in a single cell; Figure S1: Measurement of calcium concentration in a single sperm cell; Text S2: Geometrical region details; Table S1: Three-dimensional model geometry; Figure S2: Three-dimensional model geometry; Text S3: Three-dimensional model design. CatSper module; Table S2: Three-dimensional model design. CatSper module constants; Text S4: CatSper module simulation results; Figure S3: Simulation of CatSper channel flux dependence on progesterone concentration; Figure S4: Simulation of CatSper channel open state probability dependence on time upon activation with 5 µM proges-terone; Text S5: Three-dimensional model design. Calcium module; Table S3: Calcium signaling module parameters; Text S6: Homogeneous model design; Figure S5: Homogeneous model design; Text S7: Quasi-heterogeneous model; Figure S6: Scheme of the quasi-heterogeneous model; Figure S7: Quasi-heterogeneous model computational experiment on progesterone activation of sperm cells; Text S8: Protein number estimation; Text S9: Calculating the change in channel flux upon addition of nominally calcium-free buffer; Figure S8: Proof of the validity of PLCd inclusion in the model; Text S10: Discussion on participation of PLCδ in progesterone activation of human cells; Text S11: Parameter scan for the homogeneous model; Figure S9: Dependence of calcium oscillation maximal concentration on the oscillation frequency for different model parameter value. Figure S10: Vitality check for the spermatozoa adherent to poly-l-lysine. Text S12: Adhesion of different populations of sperm cells to poly-l-lysine

Author Contributions

Conceptualization, E.Y.S. and S.A.Y.; methodology, J.K. and F.A.B.; software, J.K.; validation, J.K.; formal analysis, J.K.; investigation, J.K. and F.A.B.; resources, E.Y.S., S.A.Y. and A.N.S.; data curation, J.K.; writing—original draft preparation, J.K.; writing—review and editing, A.N.S.; visualization, J.K.; supervision, E.Y.S. and A.N.S.; project administration, A.N.S.; funding acquisition, A.N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Russian Foundation for Basic Research, grant number 17-00-00138 and 21-51-10005 and Lomonosov Moscow State University Digital Medicine School.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the IVF Clinic AltraVita Ethics Committee decision number #6 from 12 January 2020.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets generated and/or analyzed during the current study are available in the Google Drive repository, available online https://drive.google.com/drive/folders/1dgNzxcq50eDWuYEwGOzTXQK_bNaDxOUN (accessed on 15 August 2020).

Acknowledgments

The authors are grateful to M. A. Panteleev for his kind advice during the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Okabe, M. The Acrosome Reaction: A Historical Perspective. In Sperm Acrosome Biogenesis and Function during Fertilization; Springer: Cham, Switzerland, 2016; pp. 1–13. [Google Scholar] [CrossRef]
  2. Armon, L.; Eisenbach, M. Behavioral Mechanism during Human Sperm Chemotaxis: Involvement of Hyperactivation. PLoS ONE 2011, 6, e28359. [Google Scholar] [CrossRef]
  3. Lishko, P.V.; Botchkina, I.L.; Kirichok, Y. Progesterone activates the principal Ca2+ channel of human sperm. Nature 2011, 471, 387–391. [Google Scholar] [CrossRef] [PubMed]
  4. Clapham, D.E. Calcium Signaling. Cell 2007, 131, 1047–1058. [Google Scholar] [CrossRef] [Green Version]
  5. Merritt, J.E.; Rink, T.J. Rapid increases in cytosolic free calcium in response to muscarinic stimulation of rat parotid acinar cells. J. Biol. Chem. 1987, 262, 4958–4960. [Google Scholar] [CrossRef]
  6. Strünker, T.; Goodwin, N.; Brenker, C.; Kashikar, N.; Weyand, I.; Seifert, R.; Kaupp, U.B. The CatSper channel mediates progesterone-induced Ca2+ influx in human sperm. Nature 2011, 471, 382–386. [Google Scholar] [CrossRef]
  7. Li, L.-F.; Xiang, C.; Zhu, Y.-B.; Qin, K.-R. Modeling of progesterone-induced intracellular calcium signaling in human spermatozoa. J. Theor. Biol. 2014, 351, 58–66. [Google Scholar] [CrossRef] [PubMed]
  8. Correia, J.; Michelangeli, F.; Publicover, S. Regulation and roles of Ca2+ stores in human sperm. Reproduction 2015, 150, R65–R76. [Google Scholar] [CrossRef] [Green Version]
  9. Sánchez-Tusie, A.A.; Vasudevan, S.R.; Churchill, G.C.; Nishigaki, T.; Treviño, C.L. Characterization of NAADP-Mediated Calcium Signaling in Human Spermatozoa. Biochem. Biophys. Res. Commun. 2014, 443, 531–536. [Google Scholar] [CrossRef]
  10. Blackmore, P.F.; Beebe, S.J.; Danforth, D.R.; Alexander, N. Progesterone and 17 Alpha-Hydroxyprogesterone. Novel Stimulators of Calcium Influx in Human Sperm. J. Biol. Chem. 1990, 265, 1376–1380. [Google Scholar] [CrossRef]
  11. José, O.; Hernández-Hernández, O.; Chirinos, M.; González-González, M.E.; Larrea, F.; Almanza, A.; Felix, R.; Darszon, A.; Treviño, C.L. Recombinant human ZP3-induced sperm acrosome reaction: Evidence for the involvement of T- and L-type voltage-gated calcium channels. Biochem. Biophys. Res. Commun. 2010, 395, 530–534. [Google Scholar] [CrossRef]
  12. Yeste, M.; Fernández-Novell, J.M.; Ramió-Lluch, L.; Estrada, E.; Rocha, L.G.; Pérez, J.; Álvaro, C.; Blanco, M.T.M.; Concha, I.I.; Ramírez, A.; et al. Intracellular calcium movements of boar spermatozoa during ‘in vitro’ capacitation and subsequent acrosome exocytosis follow a multiple-storage place, extracellular calcium-dependent model. Andrology 2015, 3, 729–747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mata-Martínez, E.; Darszon, A.; Treviño, C.L. pH-dependent Ca+2 oscillations prevent untimely acrosome reaction in human sperm. Biochem. Biophys. Res. Commun. 2018, 497, 146–152. [Google Scholar] [CrossRef]
  14. Magnus, G.; Keizer, J. Minimal Model of P-Cell Mitochondrial Ca2+ Handling. Model. Physiol. 1997, 273, C717–C733. [Google Scholar] [CrossRef]
  15. Okunade, G.W.; Miller, M.L.; Pyne, G.J.; Sutliff, R.L.; O’Connor, K.T.; Neumann, J.C.; Andringa, A.; Miller, D.A.; Prasad, V.; Doetschman, T.; et al. Targeted Ablation of Plasma Membrane Ca2+-ATPase (PMCA) 1 and 4 Indicates a Major Housekeeping Function for PMCA1 and a Critical Role in Hyperactivated Sperm Motility and Male Fertility for PMCA. J. Biol. Chem. 2004, 279, 33742–33750. [Google Scholar] [CrossRef] [Green Version]
  16. Gualtieri, R.; Talevi, R. Selection of highly fertilization-competent bovine spermatozoa through adhesion to the Fallopian tube epithelium in vitro. Reproduction 2003, 125, 251–258. [Google Scholar] [CrossRef]
  17. Harper, C.; Barratt, C.; Publicover, S.J. Stimulation of Human Spermatozoa with Progesterone Gradients to Simulate Approach to the Oocyte. J. Biol. Chem. 2004, 279, 46315–46325. [Google Scholar] [CrossRef] [Green Version]
  18. De Oliveira, G.M.; Lefièvre, L.; Ford, C.; Herrero, M.B.; Barratt, C.; Connolly, T.J.; Nash, K.; Morales-Garcia, A.; Kirkman-Brown, J.C.; Publicover, S. Mobilisation of Ca2+ stores and flagellar regulation in human sperm by S-nitrosylation: A role for NO synthesised in the female reproductive tract. Development 2008, 135, 3677–3686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Romarowski, A.; Sánchez-Cárdenas, C.; Gomez, H.V.R.; Molina, L.D.C.P.; Trevino, C.L.; Hernández-Cruz, A.; Darszon, A.; Buffone, M.G. A Specific Transitory Increase in Intracellular Calcium Induced by Progesterone Promotes Acrosomal Exocytosis in Mouse Sperm1. Biol. Reprod. 2016, 94, 63. [Google Scholar] [CrossRef]
  20. Achikanu, C.; Pendekanti, V.; Teague, R.; Publicover, S. Effects of pH manipulation, CatSper stimulation and Ca2+-store mobilization on [Ca2+]i and behaviour of human sperm. Hum. Reprod. 2018, 33, 1802–1811. [Google Scholar] [CrossRef]
  21. Servin-Vences, M.R.; Tatsu, Y.; Ando, H.; Guerrero, A.; Yumoto, N.; Darszon, A.; Nishigaki, T. A caged progesterone analog alters intracellular Ca2+ and flagellar bending in human sperm. Reproduction 2012, 144, 101–109. [Google Scholar] [CrossRef] [Green Version]
  22. Tesarik, J.; Mendoza, C. Insights into the Function of a Sperm-Surface Progesterone Receptor: Evidence of Ligand-Induced Receptor Aggregation and the Implication of Proteolysis. Exp. Cell Res. 1993, 205, 111–117. [Google Scholar] [CrossRef] [PubMed]
  23. Sagare-Patil, V.; Galvankar, M.; Satiya, M.; Bhandari, B.; Gupta, S.; Modi, D. Differential concentration and time dependent effects of progesterone on kinase activity, hyperactivation and acrosome reaction in human spermatozoa. Int. J. Androl. 2012, 35, 633–644. [Google Scholar] [CrossRef] [PubMed]
  24. Sánchez-Cárdenas, C.; Servin-Vences, R.; José, O.; Trevino, C.L.; Hernández-Cruz, A.; Darszon, A. Acrosome Reaction and Ca2+ Imaging in Single Human Spermatozoa: New Regulatory Roles of [Ca2+]i. Biol. Reprod. 2014, 91, 67. [Google Scholar] [CrossRef]
  25. Xia, J.; Reigada, D.; Mitchell, C.H.; Ren, D. CATSPER Channel-Mediated Ca2+ Entry into Mouse Sperm Triggers a Tail-to-Head Propagation. Biol. Reprod. 2007, 77, 551–559. [Google Scholar] [CrossRef] [Green Version]
  26. Miller, M.R.; Mannowetz, N.; Iavarone, A.T.; Safavi, R.; Gracheva, E.O.; Smith, J.F.; Hill, R.; Bautista, D.M.; Kirichok, Y.; Lishko, P.V. Unconventional endocannabinoid signaling governs sperm activation via the sex hormone progesterone. Science 2016, 352, 555–559. [Google Scholar] [CrossRef] [Green Version]
  27. Vyklicka, L.; Lishko, P.V. Dissecting the signaling pathways involved in the function of sperm flagellum. Curr. Opin. Cell Biol. 2020, 63, 154–161. [Google Scholar] [CrossRef]
  28. Tomes, C.N. The proteins of exocytosis: Lessons from the sperm model. Biochem. J. 2015, 465, 359–370. [Google Scholar] [CrossRef]
  29. Allen, V.; Swigart, P.; Cheung, R.; Cockcroft, S.; Katan, M. Regulation of inositol lipid-specific phospholipase Cδ by changes in Ca2+ ion concentrations. Biochem. J. 1997, 327, 545–552. [Google Scholar] [CrossRef]
  30. Fukami, K.; Yoshida, M.; Inoue, T.; Kurokawa, M.; Fissore, R.A.; Yoshida, N.; Mikoshiba, K.; Takenawa, T. Phospholipase Cδ4 is required for Ca2+ mobilization essential for acrosome reaction in sperm. J. Cell Biol. 2003, 161, 79–88. [Google Scholar] [CrossRef] [PubMed]
  31. Arenaza, I.U.; Osinalde, N.; Akimov, V.; Puglia, M.; Candenas, L.; Pinto, F.M.; Muñoa-Hoyos, I.; Gianzo, M.; Matorras, R.; Irazusta, J.; et al. Phosphoproteomic and Functional Analyses Reveal Sperm-specific Protein Changes Downstream of Kappa Opioid Receptor in Human Spermatozoa. Mol. Cell. Proteom. 2019, 18, S118–S131. [Google Scholar] [CrossRef] [Green Version]
  32. Ho, H.-C.; Suarez, S.S. Characterization of the Intracellular Calcium Store at the Base of the Sperm Flagellum That Regulates Hyperactivated Motility. Biol. Reprod. 2003, 68, 1590–1596. [Google Scholar] [CrossRef] [Green Version]
  33. Van Der Horst, G.; Maree, L.; Kotzé, S.H.; O’Riain, M.J. Sperm structure and motility in the eusocial naked mole-rat, Heterocephalus glaber: A case of degenerative orthogenesis in the absence of sperm competition? BMC Evol. Biol. 2011, 11, 351. [Google Scholar] [CrossRef] [Green Version]
  34. Priego-Espinosa, D.A.; Darszon, A.; Guerrero, A.; González-Cota, A.L.; Nishigaki, T.; Martínez-Mekler, G.; Carneiro, J. Modular analysis of the control of flagellar Ca2+-spike trains produced by CatSper and CaV channels in sea urchin sperm. PLoS Comput. Biol. 2020, 16, e1007605. [Google Scholar] [CrossRef] [PubMed]
  35. Espinal-Enríquez, J.; Priego-Espinosa, D.A.; Darszon, A.; Beltran, C.; Martínez-Mekler, G. Network model predicts that CatSper is the main Ca2+ channel in the regulation of sea urchin sperm motility. Sci. Rep. 2017, 7, 4236. [Google Scholar] [CrossRef] [Green Version]
  36. Espinal-Enríquez, J.; Darszon, A.; Guerrero, A.; Martínez-Mekler, G. In Silico Determination of the Effect of Multi-Target Drugs on Calcium Dynamics Signaling Network Underlying Sea Urchin Spermatozoa Motility. PLoS ONE 2014, 9, e104451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Aguilera, L.U.; Galindo, B.E.; Sánchez, D.; Santillan, M. What Is the Core Oscillator in the Speract-Activated Pathway of the Strongylocentrotus purpuratus Sperm Flagellum? Biophys. J. 2012, 102, 2481–2488. [Google Scholar] [CrossRef] [Green Version]
  38. Olson, S.D.; Suarez, S.S.; Fauci, L.J. A Model of CatSper Channel Mediated Calcium Dynamics in Mammalian Spermatozoa. Bull. Math. Biol. 2010, 72, 1925–1946. [Google Scholar] [CrossRef]
  39. Simons, J.; Fauci, L. A Model for the Acrosome Reaction in Mammalian Sperm. Bull. Math. Biol. 2018, 80, 2481–2501. [Google Scholar] [CrossRef]
  40. Bedu-Addo, K.; Barratt, C.; Kirkman-Brown, J.C.; Publicover, S. Patterns of [Ca2+]i mobilization and cell response in human spermatozoa exposed to progesterone. Dev. Biol. 2007, 302, 324–332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Espino, J.; Mediero, M.; Lozano, G.M.; Bejarano, I.; Ortiz, Á.; García, J.F.; Pariente, J.A.; Rodríguez, A.B. Reduced levels of intracellular calcium releasing in spermatozoa from asthenozoospermic patients. Reprod. Biol. Endocrinol. 2009, 7, 11. [Google Scholar] [CrossRef] [Green Version]
  42. Brewis, I.A.; Morton, I.E.; Mohammad, S.N.; Browes, C.E.; Moore, H.D. Measurement of intracellular calcium concentration and plasma membrane potential in human spermatozoa using flow cytometry. J. Androl. 2000, 21, 238–249. [Google Scholar] [PubMed]
  43. López-González, I.; Torres-Rodríguez, P.; Sánchez-Carranza, O.; Solís-López, A.; Santi, C.; Darszon, A.; Treviño, C. Membrane hyperpolarization during human sperm capacitation. Mol. Hum. Reprod. 2014, 20, 619–629. [Google Scholar] [CrossRef] [Green Version]
  44. Matthews, E.A.; Edietrich, D. Buffer mobility and the regulation of neuronal calcium domains. Front. Cell. Neurosci. 2015, 9, 48. [Google Scholar] [CrossRef]
  45. Naaby-Hansen, S.; Diekman, A.; Shetty, J.; Flickinger, C.J.; Westbrook, A.; Herr, J.C. Identification of calcium-binding proteins associated with the human sperm plasma membrane. Reprod. Biol. Endocrinol. 2010, 8, 6. [Google Scholar] [CrossRef] [Green Version]
  46. Mossman, J.A.; Pearson, J.T.; Moore, H.D.; Pacey, A. Variation in mean human sperm length is linked with semen characteristics. Hum. Reprod. 2013, 28, 22–32. [Google Scholar] [CrossRef]
  47. Cummins, J.M.; Woodall, P.F. On mammalian sperm dimensions. Reproduction 1985, 75, 153–175. [Google Scholar] [CrossRef]
  48. Kirkman-Brown, J.C.; Barratt, C.L.R.; Publicover, S.J. Slow calcium oscillations in human spermatozoa. Biochem. J. 2004, 378, 827–832. [Google Scholar] [CrossRef] [PubMed]
  49. Samples, J.T.; Abrams, R.M. Reliability of Urine Temperature as a Measurement of Basal Body Temperature. JOGN Nurs. 1984, 13, 319–323. [Google Scholar] [CrossRef] [PubMed]
  50. Goldstein, M.; Eid, J.-F. Elevation of Intratesticular and Scrotal Skin Surface Temperature in Men with Varicocele. J. Urol. 1989, 142, 743–745. [Google Scholar] [CrossRef]
  51. Vera-Sigüenza, E.; Pages, N.; Rugis, J.; Yule, D.I.; Sneyd, J. A Mathematical Model of Fluid Transport in an Accurate Reconstruction of Parotid Acinar Cells. Bull. Math. Biol. 2018, 81, 699–721. [Google Scholar] [CrossRef] [PubMed]
  52. Lechleiter, J.; Girard, S.; Clapham, D.; Peralta, E. Subcellular patterns of calcium release determined by G protein-specific residues of muscarinic receptors. Nature 1991, 350, 505–508. [Google Scholar] [CrossRef]
  53. Bell, M.; Bartol, T.; Sejnowski, T.; Rangamani, P. Dendritic spine geometry and spine apparatus organization govern the spatiotemporal dynamics of calcium. J. Gen. Physiol. 2019, 151, 1017–1034. [Google Scholar] [CrossRef] [Green Version]
  54. Marhl, M.; Schuster, S.; Brumen, M.; Heinrich, R. Modelling the interrelations between calcium oscillations and ER membrane potential oscillations. Biophys. Chem. 1997, 63, 221–239. [Google Scholar] [CrossRef]
  55. Zador, A.; Koch, C. Linearized models of calcium dynamics: Formal equivalence to the cable equation. J. Neurosci. 1994, 14, 4705–4715. [Google Scholar] [CrossRef] [Green Version]
  56. Morbeck, D.E.; Baumann, N.A.; Oglesbee, D. Composition of single-step media used for human embryo culture. Fertil. Steril. 2017, 107, 1055–1060.e1. [Google Scholar] [CrossRef] [Green Version]
  57. De Jonge, C.J.; Barratt, C. Methods for the Assessment of Sperm Capacitation and Acrosome Reaction Excluding the Sperm Penetration Assay. Adv. Struct. Saf. Stud. 2012, 927, 113–118. [Google Scholar] [CrossRef]
  58. Balabin, F.A.; Morozova, D.; Mayorov, A.S.; Martyanov, A.A.; Panteleev, M.A.; Sveshnikova, A. Clusterization of Inositol Trisphosphate Receptors Determines the Shape of the Calcium Oscillation Peak in Platelet Cytosol. Mosc. Univ. Phys. Bull. 2018, 73, 526–533. [Google Scholar] [CrossRef]
  59. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Grynkiewicz, G.; Poenie, M.; Tsien, R.Y. A new generation of Ca2+ indicators with greatly improved fluorescence properties. J. Biol. Chem. 1985, 260, 3440–3450. [Google Scholar] [CrossRef]
  61. Naraghi, M. T-jump study of calcium binding kinetics of calcium chelators. Cell Calcium 1997, 22, 255–268. [Google Scholar] [CrossRef] [Green Version]
  62. Moraru, I.; Morgan, F.; Li, Y.; Loew, L.; Schaff, J.C.; Lakshminarayana, A.; Slepchenko, B.M.; Gao, F.; Blinov, M. Virtual Cell modelling and simulation software environment. IET Syst. Biol. 2008, 2, 352–362. [Google Scholar] [CrossRef] [Green Version]
  63. Wang, M.; Weiss, M.; Simonovic, M.; Haertinger, G.; Schrimpf, S.P.; Hengartner, M.; von Mering, C. PaxDb, a Database of Protein Abundance Averages Across All Three Domains of Life. Mol. Cell. Proteom. 2012, 11, 492–500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Milo, R. What is the total number of protein molecules per cell volume? A call to rethink some published values. BioEssays 2013, 35, 1050–1055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Serres, C.; Escalier, D.; David, G. Ultrastructural morphometry of the human sperm flagellum with a stereological analysis of the lengths of the dense fibres. Biol. Cell 1984, 49, 153–161. [Google Scholar] [CrossRef]
  66. Sunanda, P.; Panda, B.; Dash, C.; Padhy, R.N.; Routray, P. An illustration of human sperm morphology and their functional ability among different group of subfertile males. Andrology 2018, 6, 680–689. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Laufer, N.; Segal, S.; Yaffe, H.; Svartz, H.; Grover, N. Volume and Shape of Normal Human Spermatozoa. Fertil. Steril. 1977, 28, 456–458. [Google Scholar] [CrossRef]
  68. Thomas, P.; Meizel, S. Phosphatidylinositol 4,5-bisphosphate hydrolysis in human sperm stimulated with follicular fluid or progesterone is dependent upon Ca2+ influx. Biochem. J. 1989, 264, 539–546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Hellsten, E.; Evans, J.P.; Bernard, D.J.; Jänne, P.A.; Nussbaum, R.L. Disrupted Sperm Function and Fertilin β Processing in Mice Deficient in the Inositol Polyphosphate 5-Phosphatase Inpp5b. Dev. Biol. 2001, 240, 641–653. [Google Scholar] [CrossRef] [Green Version]
  70. Cai, X.; Li, X.; Qi, H.; Wei, F.; Chen, J.; Shuai, J. Comparison of gating dynamics of different IP3R channels with immune algorithm searching for channel parameter distributions. Phys. Biol. 2016, 13, 056005. [Google Scholar] [CrossRef]
  71. Kawai, T.; Miyata, H.; Nakanishi, H.; Sakata, S.; Morioka, S.; Sasaki, J.; Watanabe, M.; Sakimura, K.; Fujimoto, T.; Sasaki, T.; et al. Polarized PtdIns(4,5)P2 distribution mediated by a voltage-sensing phosphatase (VSP) regulates sperm motility. Proc. Natl. Acad. Sci. USA 2019, 116, 26020–26028. [Google Scholar] [CrossRef] [PubMed]
  72. Hoops, S.; Sahle, S.; Gauges, R.; Lee, C.; Pahle, J.; Simus, N.; Singhal, M.; Xu, L.; Mendes, P.; Kummer, U. COPASI—A COmplex PAthway SImulator. Bioinformatics 2006, 22, 3067–3074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Saltelli, A.; Ratto, M.; Tarantola, A.S.; Campolongo, F. Sensitivity Analysis for Chemical Models. Chem. Rev. 2005, 105, 2811–2828. [Google Scholar] [CrossRef] [PubMed]
Ijms 22 08018 g001 550
Figure 1. Scheme of the main biochemical reactions during sperm cell activation by progesterone. We suggest a division of the reactions into two modules—the CatSper module and the Calcium module. (a) CatSper module. Without activation, the CatSper channel is inhibited by the plasma membrane lipid 2-arachidonoylglycerol (2-AG). Progesterone activates an abhydrolase domain-containing protein 2 (ABHD2), an enzyme that cleaves 2-AG into arachidonic acid and glycerol, and thus activates CatSper. (b) Calcium module. Calcium ions enter cytosol through open CatSper channels and then diffuse into spermatozoon midpiece (MP), where a specific isoform of phospholipase C, PLCδ, is located. This enzyme activity is upregulated by calcium; PLCδ catalyzes IP3 production from the membrane phospholipid phosphatidylinositol 4,5-bisphosphate. IP3 activates IP3 channel-receptors (IP3R) located on the spermatozoon calcium store, RNE (redundant nuclear envelope), and thus induces calcium release into the cell cytosol. Plasma Membrane Calcium ATPase (PMCA) present in the sperm cilia membrane extrudes calcium into extracellular media. Signaling Pathway Calcium ATPase (SPCA) pumps calcium into the calcium stores RNE. Moreover, calcium buffers such as calmodulin (CaM) located in the cytosol and calreticulin (CalRet) located in the endoplasmic reticulum can reversibly bind calcium. IP3 is metabolized by inositol polyphosphate 5-phosphatase (IPP5P).
Figure 1. Scheme of the main biochemical reactions during sperm cell activation by progesterone. We suggest a division of the reactions into two modules—the CatSper module and the Calcium module. (a) CatSper module. Without activation, the CatSper channel is inhibited by the plasma membrane lipid 2-arachidonoylglycerol (2-AG). Progesterone activates an abhydrolase domain-containing protein 2 (ABHD2), an enzyme that cleaves 2-AG into arachidonic acid and glycerol, and thus activates CatSper. (b) Calcium module. Calcium ions enter cytosol through open CatSper channels and then diffuse into spermatozoon midpiece (MP), where a specific isoform of phospholipase C, PLCδ, is located. This enzyme activity is upregulated by calcium; PLCδ catalyzes IP3 production from the membrane phospholipid phosphatidylinositol 4,5-bisphosphate. IP3 activates IP3 channel-receptors (IP3R) located on the spermatozoon calcium store, RNE (redundant nuclear envelope), and thus induces calcium release into the cell cytosol. Plasma Membrane Calcium ATPase (PMCA) present in the sperm cilia membrane extrudes calcium into extracellular media. Signaling Pathway Calcium ATPase (SPCA) pumps calcium into the calcium stores RNE. Moreover, calcium buffers such as calmodulin (CaM) located in the cytosol and calreticulin (CalRet) located in the endoplasmic reticulum can reversibly bind calcium. IP3 is metabolized by inositol polyphosphate 5-phosphatase (IPP5P).
Ijms 22 08018 g001
Ijms 22 08018 g002 550
Figure 2. Responses of calcium in human spermatozoa to progesterone at 5 μM. (a–c) A single transient peak of calcium concentration, 160 ± 44 s wide and 0.5 ± 0.2 μM high, was observed in most progesterone-activated spermatozoa. Sometimes the main peak was followed by a minor calcium elevation (c) starting in ~100 s after progesterone activation. The peak width and amplitude were typical for 36 out of 64 cells. (d–f) Different examples of low-frequency (T~100 s) calcium oscillations detected in progesterone-activated spermatozoa. The oscillation amplitude and frequency were typical for 18 out of 64 cells. (g,h) A steady calcium increase was observed in some progesterone-treated cells. This type of response was observed in 12 out of 64 cells. (i) The histogram of different types of calcium responses in human spermatozoa shows the diversity of calcium responses to activation, with no particular dominating group.
Figure 2. Responses of calcium in human spermatozoa to progesterone at 5 μM. (a–c) A single transient peak of calcium concentration, 160 ± 44 s wide and 0.5 ± 0.2 μM high, was observed in most progesterone-activated spermatozoa. Sometimes the main peak was followed by a minor calcium elevation (c) starting in ~100 s after progesterone activation. The peak width and amplitude were typical for 36 out of 64 cells. (d–f) Different examples of low-frequency (T~100 s) calcium oscillations detected in progesterone-activated spermatozoa. The oscillation amplitude and frequency were typical for 18 out of 64 cells. (g,h) A steady calcium increase was observed in some progesterone-treated cells. This type of response was observed in 12 out of 64 cells. (i) The histogram of different types of calcium responses in human spermatozoa shows the diversity of calcium responses to activation, with no particular dominating group.
Ijms 22 08018 g002
Ijms 22 08018 g003 550
Figure 3. Computational models design and typical simulation results in comparison with experimental data. Progesterone concentration for both the model and the experiment was set at 5 µM. (a) Three-dimensional model design. PMCA, as well as CatSper protein location, is restricted to sperm cilia membrane; the corresponding reactions occur only in the specified location. IP3 production by PLCδ is also assumed to be happening only in the sperm MP and head. (b) Oscillatory response in the homogenous model. The oscillation period (30 s) described by the homogeneous model is significantly lower than the experimentally observed one (100 s). CatSper maximum conductance VCatSper was set to 0.20 µM/s. Progesterone concentration was set to 5 µM at t = 25 s. (c) Oscillatory response in the 3D model. CatSper maximum conductance VCatSper was set to 0.7 µM/s. Dca was set to 18 µm2/s. Progesterone concentration was set to 5 µM at t = 25 s. The model oscillatory response period (200 s) was close to the experimentally observed one. (d) Experimentally observed calcium oscillations. (e) Homogeneous model single-spike response. The peak width in the homogeneous model (40 s) is significantly lower than that observed in the experiment (160 s). CatSper maximum conductance VCatSper was set equal to 0.5 µM/s. Progesterone concentration was set to 5 µM at t = 25 s. (f) Three-dimensional model of a single peak response. The model single peak width (130 s) was close to the experimentally observed one (160 s). CatSper maximum conductance VCatSper was set to 1.6 µM/s. Dca was set to 23.5 µm2/s. Progesterone concentration was set from 0 to 5 µM at t = 25 s. (g) An experimentally observed single peak (see details of the experiments in Figure 2 legend and Methods).
Figure 3. Computational models design and typical simulation results in comparison with experimental data. Progesterone concentration for both the model and the experiment was set at 5 µM. (a) Three-dimensional model design. PMCA, as well as CatSper protein location, is restricted to sperm cilia membrane; the corresponding reactions occur only in the specified location. IP3 production by PLCδ is also assumed to be happening only in the sperm MP and head. (b) Oscillatory response in the homogenous model. The oscillation period (30 s) described by the homogeneous model is significantly lower than the experimentally observed one (100 s). CatSper maximum conductance VCatSper was set to 0.20 µM/s. Progesterone concentration was set to 5 µM at t = 25 s. (c) Oscillatory response in the 3D model. CatSper maximum conductance VCatSper was set to 0.7 µM/s. Dca was set to 18 µm2/s. Progesterone concentration was set to 5 µM at t = 25 s. The model oscillatory response period (200 s) was close to the experimentally observed one. (d) Experimentally observed calcium oscillations. (e) Homogeneous model single-spike response. The peak width in the homogeneous model (40 s) is significantly lower than that observed in the experiment (160 s). CatSper maximum conductance VCatSper was set equal to 0.5 µM/s. Progesterone concentration was set to 5 µM at t = 25 s. (f) Three-dimensional model of a single peak response. The model single peak width (130 s) was close to the experimentally observed one (160 s). CatSper maximum conductance VCatSper was set to 1.6 µM/s. Dca was set to 23.5 µm2/s. Progesterone concentration was set from 0 to 5 µM at t = 25 s. (g) An experimentally observed single peak (see details of the experiments in Figure 2 legend and Methods).
Ijms 22 08018 g003
Ijms 22 08018 g004 550
Figure 4. Influence of model parameters on the calcium oscillation response in the 3D model. (a) Variation in the calcium ion diffusion coefficient led to the disappearance of the oscillation at higher Dca values. The solid line represents steady-state calcium concentration. The dashed and dotted line represents the lowest concentration upon calcium oscillation). The dashed line represents the highest concentration upon calcium oscillation. (b) Variation of the Vmax for PMCA. The solid black line denotes the steady-state or maximum calcium concentration upon calcium oscillation. The dashed and dotted black line denotes the minimum calcium concentration during oscillations. (c) Dependence of the response type on PLCδ distribution. For a PLCδ activity restricted to sperm head and MP, an oscillatory response occurred (red line). For a PLCδ activity distributed uniformly, a significant delay in calcium response appeared (dashed line). For PMCA activity distributed uniformly, the oscillations period was decreased to 60 s (purple line). (d) For fixed IP3 concentration of 0.013 µm, the oscillations period and the delay were increasing for lower calcium diffusion coefficients. (e) Dependence of cytosolic calcium concentration on (fixed) IP3 concentration. Dca = 10 µm2/s for model plots from (c,d). VCatSper = 0.7 µM/s for model plots from (a–e).
Figure 4. Influence of model parameters on the calcium oscillation response in the 3D model. (a) Variation in the calcium ion diffusion coefficient led to the disappearance of the oscillation at higher Dca values. The solid line represents steady-state calcium concentration. The dashed and dotted line represents the lowest concentration upon calcium oscillation). The dashed line represents the highest concentration upon calcium oscillation. (b) Variation of the Vmax for PMCA. The solid black line denotes the steady-state or maximum calcium concentration upon calcium oscillation. The dashed and dotted black line denotes the minimum calcium concentration during oscillations. (c) Dependence of the response type on PLCδ distribution. For a PLCδ activity restricted to sperm head and MP, an oscillatory response occurred (red line). For a PLCδ activity distributed uniformly, a significant delay in calcium response appeared (dashed line). For PMCA activity distributed uniformly, the oscillations period was decreased to 60 s (purple line). (d) For fixed IP3 concentration of 0.013 µm, the oscillations period and the delay were increasing for lower calcium diffusion coefficients. (e) Dependence of cytosolic calcium concentration on (fixed) IP3 concentration. Dca = 10 µm2/s for model plots from (c,d). VCatSper = 0.7 µM/s for model plots from (a–e).
Ijms 22 08018 g004
Ijms 22 08018 g005 550
Figure 5. Miscellaneous experimental data described by the heterogeneous model. (a) Spontaneous calcium oscillations were demonstrated by the 3D model with distance L between the sperm’s principal piece and sperm’s store RNE of 5 μm (black line), but not 7.5 μm (dashed line). (b) Dependence of the model spontaneous oscillations period on the calcium diffusion coefficient Dca. (c) 3D model calcium oscillations in progesterone-activated spermatozoa upon addition of a nominally calcium-free medium. The moment of setting the conditions to one corresponding to nominally-free calcium medium is shown by the red arrow.
Figure 5. Miscellaneous experimental data described by the heterogeneous model. (a) Spontaneous calcium oscillations were demonstrated by the 3D model with distance L between the sperm’s principal piece and sperm’s store RNE of 5 μm (black line), but not 7.5 μm (dashed line). (b) Dependence of the model spontaneous oscillations period on the calcium diffusion coefficient Dca. (c) 3D model calcium oscillations in progesterone-activated spermatozoa upon addition of a nominally calcium-free medium. The moment of setting the conditions to one corresponding to nominally-free calcium medium is shown by the red arrow.
Ijms 22 08018 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Korobkin, J.; Balabin, F.A.; Yakovenko, S.A.; Simonenko, E.Y.; Sveshnikova, A.N. Occurrence of Calcium Oscillations in Human Spermatozoa Is Based on Spatial Signaling Enzymes Distribution. Int. J. Mol. Sci. 2021, 22, 8018. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22158018

AMA Style

Korobkin J, Balabin FA, Yakovenko SA, Simonenko EY, Sveshnikova AN. Occurrence of Calcium Oscillations in Human Spermatozoa Is Based on Spatial Signaling Enzymes Distribution. International Journal of Molecular Sciences. 2021; 22(15):8018. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22158018

Chicago/Turabian Style

Korobkin, Julia, Fedor A. Balabin, Sergey A. Yakovenko, Ekaterina Yu. Simonenko, and Anastasia N. Sveshnikova. 2021. "Occurrence of Calcium Oscillations in Human Spermatozoa Is Based on Spatial Signaling Enzymes Distribution" International Journal of Molecular Sciences 22, no. 15: 8018. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22158018

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop