In silico determination of the effect of multi-target drugs on calcium dynamics signaling network underlying sea urchin spermatozoa motility

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

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.

Strategies to Enhance Logic Modeling-Based Cell Line-Specific Drug Synergy Prediction

Discrete dynamical modeling shows promise in prioritizing drug combinations for screening efforts by reducing the experimental workload inherent to the vast numbers of possible drug combinations. We have investigated approaches to predict combination responses across different cancer cell lines using logic models generated from one generic prior-knowledge network representing 144 nodes covering major cancer signaling pathways. Cell-line specific models were configured to agree with baseline activity data from each unperturbed cell line. Testing against experimental data demonstrated a high number of true positive and true negative predictions, including also cell-specific responses. We demonstrate the possible enhancement of predictive capability of models by curation of literature knowledge further detailing subtle biologically founded signaling mechanisms in the model topology. In silico model analysis pinpointed a subset of network nodes highly influencing model predictions. Our results indicate that the performance of logic models can be improved by focusing on high-influence node protein activity data for model configuration and that these nodes accommodate high information flow in the regulatory network.

Modular analysis of the control of flagellar Ca2+-spike trains produced by CatSper and CaV channels in sea urchin sperm

Intracellular calcium ([Ca²⁺]_i) is a basic and ubiquitous cellular signal controlling a wide variety of biological processes. A remarkable example is the steering of sea urchin spermatozoa towards the conspecific egg by a spatially and temporally orchestrated series of [Ca²⁺]_i spikes. Although this process has been an experimental paradigm for reproduction and sperm chemotaxis studies, the composition and regulation of the signalling network underlying the cytosolic calcium fluctuations are hitherto not fully understood. Here, we used a differential equations model of the signalling network to assess which set of channels can explain the characteristic envelope and temporal organisation of the [Ca²⁺]_i-spike trains. The signalling network comprises an initial membrane hyperpolarisation produced by an Upstream module triggered by the egg-released chemoattractant peptide, via receptor activation, cGMP synthesis and decay. Followed by downstream modules leading to intraflagellar pH (pH_i), voltage and [Ca²⁺]_i fluctuations. The Upstream module outputs were fitted to kinetic data on cGMP activity and early membrane potential changes measured in bulk cell populations. Two candidate modules featuring voltage-dependent Ca²⁺-channels link these outputs to the downstream dynamics and can independently explain the typical decaying envelope and the progressive spacing of the spikes. In the first module, [Ca²⁺]_i-spike trains require the concerted action of a classical Ca_V-like channel and a potassium channel, BK (Slo1), whereas the second module relies on pH_i-dependent CatSper dynamics articulated with voltage-dependent neutral sodium-proton exchanger (NHE). We analysed the dynamics of these two modules alone and in mixed scenarios. We show that the [Ca²⁺]_i dynamics observed experimentally after sustained alkalinisation can be reproduced by a model featuring the CatSper and NHE module but not by those including the pH-independent Ca_V and BK module or proportionate mixed scenarios. We conclude in favour of the module containing CatSper and NHE and highlight experimentally testable predictions that would corroborate this conclusion.

Fertilisation in marine invertebrates, such as the sea urchin, occurs during broadcast spawning events in which males and females of co-localised species ejaculate sperm and spawn eggs synchronously. During these events, spermatozoa have to find and fertilise conspecific eggs in the midst of all the other ones, which is a remarkable navigation and mating choice achievement. Sperm cells do this by navigating towards the source of species-specific peptides released by the egg, steered by spatial and temporally orchestrated peaks in intracellular calcium concentration that trigger sudden reorientations. How these calcium spikes are regulated and timed remains elusive. Different calcium channels have been implicated by indirect experimental evidence giving rise to a complex network of putative interacting components. We gained insight into the structure and function of this network by modelling it as a set of candidate modules that could be studied separately. By using this ‘divide and conquer’ approach to the complexity of the network, we could characterise the potential dynamics of each module and confront these dynamics with specific quantitative data. Our results indicate that the channel mediating calcium signals in sea urchin sperm is likely CatSper, a calcium channel necessary for human male fertility.

Network model predicts that CatSper is the main Ca2+ channel in the regulation of sea urchin sperm motility

Spermatozoa sea urchin swimming behaviour is regulated by small peptides from the egg outer envelope. Speract, such a peptide, after binding to its receptor in Strongylocentrotus purpuratus sperm flagella, triggers a signaling pathway that culminates with a train of intracellular calcium oscillations, correlated with changes in sperm swimming pattern. This pathway has been widely studied but not fully characterized. Recent work on Arbacia punctulata sea urchin spermatozoa has documented the presence of the Ca²⁺ CatSper channel in their flagella and its involvement in chemotaxis. However, if other calcium channels participate in chemotaxis remains unclear. Here, based on an experimentally-backed logical network model, we conclude that CatSper is fundamental in the S. purpuratus speract-activated sea urchin sperm signaling cascade, although other Ca²⁺ channels could still be relevant. We also present for the first time experimental corroboration of its active presence in S. purpuratus sperm flagella. We argue, prompted by in silico knock-out calculations, that CatSper is the main generator of calcium oscillations in the signaling pathway and that other calcium channels, if present, have a complementary role. The approach adopted here allows us to unveil processes, which are hard to detect exclusively by experimental procedures.