In Silico Determination of the Effect of Multi-Target Drugs on Calcium Dynamics Signaling Network Underlying Sea Urchin Spermatozoa Motility

The motility of spermatozoa of both Lytechinus pictus and Strongylocentrotus purpuratus sea urchin species is modulated by the egg-derived decapeptide speract via an oscillatory [Ca2+]-dependent signaling pathway. Comprehension of this pathway is hence directly related to the understanding of regulated sperm swimming. Niflumic acid (NFA), a nonsteroidal anti-inflammatory drug alters several ion channels. Though unspecific, NFA profoundly affects how sea urchin sperm respond to speract, increasing the [Ca2+]i oscillation period, amplitude, peak and average level values of the responses in immobilized and swimming cells. A previous logical network model we developed for the [Ca2+] dynamics of speract signaling cascade in sea urchin sperm allows integrated dissection of individual and multiple actions of NFA. Among the channels affected by NFA are: hyperpolarization-activated and cyclic nucleotide gated Na+ channels (HCN), [Ca2+]-dependent Cl− channels (CaCC) and [Ca2+]-dependent K+ channels (CaKC), all present in the sea urchin genome. Here, using our model we investigated the effect of blocking in silico HCN and CaCC channels suggested by experiments. Regarding CaKC channels, arguments can be provided for either their blockage or activation by NFA. Our study yielded two scenarios compliant with experimental observations: i) under CaKC inhibition, this [Ca2+]-dependent K+ channel should be different from the Slo1 channel and ii) under activation of the CaKC channel, another [Ca2+] channel not considered previously in the network is required, such as the pH-dependent CatSper channel. Additionally, our findings predict cause-effect relations resulting from a selective inhibition of those channels. Knowledge of these relations may be of consequence for a variety of electrophysiological studies and have an impact on drug related investigations. Our study contributes to a better grasp of the network dynamics and suggests further experimental work.