Analyzing the functional divergence of Slo1 and Slo3 channel subfamilies

Physiological role of potassium channels in mammalian germ cell differentiation, maturation, and capacitation

BACKGROUND

Ion channels are essential for differentiation and maturation of germ cells, and even for fertilization in mammals. Different types of potassium channels have been identified, which are grouped into voltage-gated channels (Kv), ligand-gated channels (Kligand ), inwardly rectifying channels (Kir ), and tandem pore domain channels (K2P ).

MATERIAL-METHODS

The present review includes recent findings on the role of potassium channels in sperm physiology of mammals.

RESULTS-DISCUSSION

While most studies conducted thus far have been focused on the physiological role of voltage- (Kv1, Kv3, and Kv7) and calcium-gated channels (SLO1 and SLO3) during sperm capacitation, especially in humans and rodents, little data about the types of potassium channels present in the plasma membrane of differentiating germ cells exist. In spite of this, recent evidence suggests that the content and regulation mechanisms of these channels vary throughout spermatogenesis. Potassium channels are also essential for the regulation of sperm cell volume during epididymal maturation and for preventing premature membrane hyperpolarization. It is important to highlight that the nature, biochemical properties, localization, and regulation mechanisms of potassium channels are species-specific. In effect, while SLO3 is the main potassium channel involved in the K+ current during sperm capacitation in rodents, different potassium channels are implicated in the K+ outflow and, thus, plasma membrane hyperpolarization during sperm capacitation in other mammalian species, such as humans and pigs.

CONCLUSIONS

Potassium conductance is essential for male fertility, not only during sperm capacitation but throughout the spermiogenesis and epididymal maturation.

Pharmacological Evidence Suggests That Slo3 Channel Is the Principal K+ Channel in Boar Spermatozoa

Sperm ion channels are associated with the quality and type of flagellar movement, and their differential regulation is crucial for sperm function during specific phases. The principal potassium ion channel is responsible for the majority of K⁺ ion flux, resulting in membrane hyperpolarization, and is essential for sperm capacitation-related signaling pathways. The molecular identity of the principal K⁺ channel varies greatly between different species, and there is a lack of information about boar K⁺ channels. We aimed to determine the channel identity of boar sperm contributing to the primary K⁺ current using pharmacological dissection. A series of Slo1 and Slo3 channel modulators were used for treatment. Sperm motility and related kinematic parameters were monitored using a computer-assisted sperm analysis system under non-capacitated conditions. Time-lapse flow cytometry with fluorochromes was used to measure changes in different intracellular ionic concentrations, and conventional flow cytometry was used to determine the acrosome reaction. Membrane depolarization, reduction in acrosome reaction, and motility parameters were observed upon the inhibition of the Slo3 channel, suggesting that the Slo3 gene encodes the main K⁺ channel in boar spermatozoa. The Slo3 channel was localized on the sperm flagellum, and the inhibition of Slo3 did not reduce sperm viability. These results may aid potential animal-model-based extrapolations and help to ameliorate motility and related parameters, leading to improved assisted reproductive methods in industrial livestock production.

Homozygous mutation in SLO3 leads to severe asthenoteratozoospermia due to acrosome hypoplasia and mitochondrial sheath malformations

BACKGROUND

Potassium channels are important for the structure and function of the spermatozoa. As a potassium transporter, the mSlo3 is essential for male fertility as Slo3 knockout male mice were infertile with the series of functional defects in sperm cells. However, no pathogenic variant has been detected in human SLO3 to date. Here we reported a human case with homozygous SLO3 mutation. The function of SLO3 in human sperm and the corresponding assisted reproductive strategy are also investigated.

METHODS

We performed whole-exome sequencing analysis from a large cohort of 105 patients with asthenoteratozoospermia. The effects of the variant were investigated by quantitative RT-PCR, western blotting, and immunofluorescence assays using the patient spermatozoa. Sperm morphological and ultrastructural studies were conducted using haematoxylin and eosin staining, scanning and transmission electron microscopy.

RESULTS

We identified a homozygous missense variant (c.1237A > T: p.Ile413Phe) in the sperm-specific SLO3 in one Chinese patient with male infertility. This SLO3 variant was rare in human control populations and predicted to be deleterious by multiple bioinformatic tools. Sperm from the individual harbouring the homozygous SLO3 variant exhibited severe morphological abnormalities, such as acrosome hypoplasia, disruption of the mitochondrial sheath, coiled tails, and motility defects. The levels of SLO3 mRNA and protein in spermatozoa from the affected individual were reduced. Furthermore, the acrosome reaction, mitochondrial membrane potential, and membrane potential during capacitation were also afflicted. The levels of acrosome marker glycoproteins and PLCζ1 as well as the mitochondrial sheath protein HSP60 and SLO3 auxiliary subunit LRRC52, were significantly reduced in the spermatozoa from the affected individual. The affected man was sterile due to acrosome and mitochondrial dysfunction; however, intra-cytoplasmic sperm injection successfully rescued this infertile condition.

CONCLUSIONS

SLO3 deficiency seriously impact acrosome formation, mitochondrial sheath assembly, and the function of K⁺ channels. Our findings provided clinical implications for the genetic and reproductive counselling of affected families.

Homology modeling and molecular docking simulation of martentoxin as a specific inhibitor of the BK channel

Background

Large conductance calcium-activated potassium channel (BK channel) is gated by both voltage and calcium ions and is widely distributed in excitable and nonexcitable cells. BK channel plays an important role in epilepsy and other diseases, but BK channel subtype-specific drugs are still extremely rare. Martentoxin was previously isolated from the venom of members of Scorpionidae and shown to be composed of 37 amino acids. Research has shown that the pharmacological selectivity of martentoxin to the BK channel is higher than that to other potassium channels. Therefore, it is of great significance to study the mechanism of interaction between martentoxin and BK channels.

Methods

The three-dimensional structure of BK channel pore region was constructed by homologous modeling method, and the key amino acid sites of BK channel interaction with martentoxin were analyzed by protein-protein docking, molecular dynamic simulation and virtual alanine mutation.

Results

Based on homologous modeling of BK channel pore structure and protein-protein docking analysis, Phe1, Lys28 and Arg35 of martentoxin were found to be key amino acids in toxin BK channel interaction.

Conclusions

This study reveals the structural basis of martentoxin interaction with BK channel. These results will contribute to the design of BK channel specific blockers based on the structure of martentoxin.

Analysis of selection in protein-coding sequences accounting for common biases

The evolution of protein-coding genes is usually driven by selective processes, which favor some evolutionary trajectories over others, optimizing the subsequent protein stability and activity. The analysis of selection in this type of genetic data is broadly performed with the metric nonsynonymous/synonymous substitution rate ratio (dN/dS). However, most of the well-established methodologies to estimate this metric make crucial assumptions, such as lack of recombination or invariable codon frequencies along genes, which can bias the estimation. Here, we review the most relevant biases in the dN/dS estimation and provide a detailed guide to estimate this metric using state-of-the-art procedures that account for such biases, along with illustrative practical examples and recommendations. We also discuss the traditional interpretation of the estimated dN/dS emphasizing the importance of considering complementary biological information such as the role of the observed substitutions on the stability and function of proteins. This review is oriented to help evolutionary biologists that aim to accurately estimate selection in protein-coding sequences.

The influence of heterogeneous codon frequencies along sequences on the estimation of molecular adaptation

MOTIVATION

The nonsynonymous/synonymous substitution rate ratio (dN/dS) is a commonly used parameter to quantify molecular adaptation in protein-coding data. It is known that the estimation of dN/dS can be biased if some evolutionary processes are ignored. In this concern, common ML methods to estimate dN/dS assume invariable codon frequencies among sites, despite this characteristic is rare in nature, and it could bias the estimation of this parameter.

RESULTS

Here we studied the influence of variable codon frequencies among genetic regions on the estimation of dN/dS. We explored scenarios varying the number of genetic regions that differ in codon frequencies, the amount of variability of codon frequencies among regions and the nucleotide frequencies at each codon position among regions. We found that ignoring heterogeneous codon frequencies among regions overall leads to underestimation of dN/dS and the bias increases with the level of heterogeneity of codon frequencies. Interestingly, we also found that varying nucleotide frequencies among regions at the first or second codon position leads to underestimation of dN/dS while variation at the third codon position leads to overestimation of dN/dS. Next, we present a methodology to reduce this bias based on the analysis of partitions presenting similar codon frequencies and we applied it to analyze four real datasets. We conclude that accounting for heterogeneous codon frequencies along sequences is required to obtain realistic estimates of molecular adaptation through this relevant evolutionary parameter.

AVAILABILITY AND IMPLEMENTATION

The applied frameworks for the computer simulations of protein-coding data and estimation of molecular adaptation are SGWE and PAML, respectively. Both are publicly available and referenced in the study.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.