Language Selection

Get healthy now with MedBeds!
Click here to book your session

Protect your whole family with Orgo-Life® Quantum MedBed Energy Technology® devices.

Advertising by Adpathway

         

 Advertising by Adpathway

A histidine-mediated supramolecular mechanism links temperature and pH activation of the sperm channel CatSper

1 week ago 8

PROTECT YOUR DNA WITH QUANTUM TECHNOLOGY

Orgo-Life the new way to the future

  Advertising by Adpathway

  • Loading metrics

Open Access

Peer-reviewed

Discovery Report

Discovery Reports capture exploratory initial findings centered around a valuable and interesting research question.

See all article types »

?

This is an uncorrected proof.

Abstract

The cation channel of sperm, CatSper, is a sperm-specific calcium channel essential for male fertility across metazoans. CatSper is activated by intracellular alkalinization, membrane depolarization, and elevated temperature, which together ensure maximal current generation. However, the structural mechanisms by which these physiological stimuli jointly regulate channel activation remain poorly understood, in part because CatSper is an unusually complex ~15-subunit assembly organized in zigzag supramolecular arrays along the sperm flagellum. Here, we combine comparative genomics across 47 species with AlphaFold3-based modeling and evolutionary sequence-structure analyses to investigate a mechanism underlying temperature and pH sensitivity of CatSper. Phylogenetic analysis shows that N-terminal length and histidine enrichment of the pore-forming subunit CatSper1 strongly correlate with species-specific fertilization temperatures, suggesting evolutionary adaptation of temperature and pH sensitivity. Structural modeling further revealed conserved surface-exposed histidine clusters that form predicted contact interfaces between adjacent CatSper assemblies, with this contact interface positioned near the dominant voltage-sensing module. Functional validation using electrophysiology and calcium imaging in mouse sperm shows that capacitation-associated partial removal of the CatSper1 N-terminus selectively impairs temperature-dependent activation. Together, these findings support a model in which temperature-dependent histidine deprotonation functionally couples neighboring CatSper complexes within the supramolecular assembly, thereby promoting synchronized channel activation.

Citation: Zhao B, Bhagwat S, Ferreira J, Ruff KM, Swain DK, Pappu RV, et al. (2026) A histidine-mediated supramolecular mechanism links temperature and pH activation of the sperm channel CatSper. PLoS Biol 24(10): e3003776. https://doi.org/10.1371/journal.pbio.3003776

Academic Editor: Raimund Dutzler, University of Zurich, SWITZERLAND

Received: March 28, 2026; Accepted: September 23, 2026; Published: October 9, 2026

Copyright: © 2026 Zhao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The code and data used in this study have been archived on Zenodo and can be accessed with the DOI https://doi.org/10.5281/zenodo.22716196. All other relevant data are within the paper and its Supporting information files.

Funding: This work was supported by BJC Investigator fund and NICHD 1R01HD119961 to P.V.L., NICHD 1R01HD119961 to C.S., startup funding from the McDonnell Center for Cellular and Molecular Neurobiology, the Department of Cell Biology & Physiology, the Center for the Investigation of Membrane Excitability Diseases at Washington University in St. Louis to Z.F, the US National Science Foundation (MCB-2227268) to R.V.P., and the WashU Center for Biomolecular Condensates (K.M.R. and R.V.P). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: P.L is a member of PLOS Biology’s Editorial Board.

Abbreviations: aa, amino acid; CatSper, cation channel of sperm; CIMED, Center for the Investigation of Membrane Excitability Diseases; cryo-ET, cryo-electron tomography; FCE, chain expanding; FCR, fraction of charged residues; HRC, Histidine-rich calcium-binding protein; HRG, histidine-rich glycoprotein; IDR, intrinsically disordered region; LLPS, liquid–liquid phase separation; NC, non-capacitating; NCPR, net charge per residue; NTD, N-terminal domain; RSA, Relative Solvent Accessibility; SASA, Solvent Accessible Surface Area; VSD, voltage-sensor domain.

Introduction

The cation channel of sperm (CatSper) is a sperm-specific calcium ion channel that plays an important role in sperm function across the evolutionary tree and is essential for male fertility in mice and humans [1–15]. Originally identified in mice [1], CatSper channels are indispensable for hyperactivation—an asymmetric flagellar bending pattern required for successful fertilization—which is triggered by calcium influx through this ion channel [1,8,16]. Because calcium entry must be tightly controlled, CatSper activity is regulated by multiple modalities that ensure it functions only under precise physiological conditions. CatSper is gated by warm temperatures [17], membrane depolarization [6,12,18,19], and intracellular alkalinization [6,12,19], while cytoplasmic acidification potently inhibits its function [6].

This multimodal regulation is supported by the structural complexity of CatSper, which forms a large multiprotein complex composed of at least 15 distinct subunits, each encoded by a separate gene [1,4,7,13,17,20–32], and is called CatSpermasome [31]. Each CatSpermasome is further assembled into highly ordered zigzag arrays that form four longitudinal rows running along the entire length of the sperm principal piece [21].

Comparative genomic studies have revealed that while CatSper genes are broadly conserved across vertebrates, they exhibit striking lineage-specific patterns of gain and loss, highlighting the channel’s evolutionary diversity [10,15,20,33]. In mammals such as mice and humans, CatSper is indispensable for male fertility; however, in certain avian and fish lineages, as well as mollusks, the channel has been lost, suggesting the evolution of alternative mechanisms for sperm activation [10,15,34].

The pore of the CatSper channel is formed by the relatively homologous subunits CatSper1, CatSper2, CatSper3, and CatSper4 [7]. All four are six-transmembrane domain polypeptides encoded by separate genes, and they share a conserved selectivity filter. While most other CatSper subunits are conserved, one in particular- CatSper1 – shows exceptional divergence across species, particularly in the length of its N-terminal domain highly enriched in histidine residues.

The activity of CatSper complex is also time-dependent. As sperm undergo their final maturation within the female reproductive tract, a process known as capacitation [35,36], the composition of the plasma membrane changes and sperm proteins undergo multiple modifications, including tyrosine phosphorylation [37,38] and targeted proteolytic processing [39]. These changes progressively prime sperm for optimal function. However, excessive or prolonged capacitation can ultimately compromise sperm performance. For example, murine sperm that undergo capacitation for more than 90 min show reduced functional performance [39]. Interestingly, the CatSper complex undergoes capacitation-dependent proteolytic cleavage that specifically removed the portion of N-terminal domain of CatSper1 [14,39]. The timing of this cleavage closely coincides with the reported loss of CatSper temperature sensitivity [17], suggesting that the CatSper1 N-terminal domain may play a central role in temperature-dependent channel regulation.

To explain the multimodal regulation of CatSper, such as its ability to integrate pH and thermal signals to fine-tune calcium entry during sperm hyperactivation, we combined comparative genomics with AlphaFold3 [40] prediction. This analysis produced a mechanistic model for the operation of this unique ion channel complex. According to this model, surface-exposed histidine residues in the N-terminal domain of CatSper1 engage in pH-dependent interactions, stabilizing a previously presumed disordered region and promote coupling among adjacent CatSper channels. Increasing temperature would favor histidine deprotonation, thereby strengthening these interactions and enhancing coupling across CatSper channel array. This cooperative mechanism could promote synchronized channel opening and provide a framework for understanding how CatSper integrates pH and temperature signals to regulate Ca2+ influx during sperm activation.

The importance of CatSper1 N-terminal domain for temperature sensitivity was further supported by electrophysiological and calcium imaging recording from murine sperm that underwent partial degradation of the N-terminal domain [14,17,39].

Results

Comparative genomics analysis revealed a correlation between N-terminal composition of CatSper1 and optimal fertilization temperature

The core of the CatSpermasome complex is formed by four pore-forming subunits: CatSper1, CatSper2, CatSper3, and CatSper4 [7], while four auxiliary subunits- CatSper β [27], δ [13], ε [23,41], and γ [25]—extend large extracellular domains that form a canopy over the channel pore [21,31]. All eight CatSper subunits are transmembrane proteins and are largely conserved across species, with the notable exception of CatSper1, which exhibits substantial variation in length (Figs 1A–1B and S1A) and shows a positive correlation between its overall length and species-specific fertilization temperatures (Figs 1C and S1B). The feature that sets CatSper1 apart is its unusually long N-terminal domain (NTD), whose length also scales up with fertilization temperature and is highly enriched in histidine (His, H) residues (Fig 1D). Interestingly, marine invertebrates such as corals and sea urchins that spawn their sperm directly into sea water, feature truncated NTDs barely reaching 120 amino acids (Fig 1A–1D and Table 1). This trend was shared with reptiles that presented a slightly elongated NTD around 127 amino acids. As the core body temperature increased, particularly for internal fertilizers, the length of the domain rose exponentially, reaching ~500 residues in species such as Lontra Canadensis (River Otter) and Oryctolagus Cuniculus (European Rabbit).

thumbnail

Table 1. Species included in this study and their CatSper-related parameters. The lengths of the CatSper1 subunit and its N-terminal domain (NT) are reported as numbers of amino acid residues (aa). Additional columns include the theoretical isoelectric point (pI) of the CatSper1 N-terminus, fertilization temperature in °C, percentage (%) and absolute number of histidine residues (#) within the N-terminus, and the corresponding references for fertilization temperatures for each species. Colored area depicts cold-fertilizing species (6–23 °C; blue), intermediate temperatures (25–34.5 °C; green), and warm-fertilizing species (35–40 °C; red).

https://doi.org/10.1371/journal.pbio.3003776.t001

thumbnail

Fig 1. Comparative genomic analyses reveal temperature-dependent divergence of the CatSper1 N-terminus across CatSper-bearing species.

(Panel A) Phylogenetic tree of representative species that retain CatSper channels. Created in BioRender. Lishko, P. (2026) https://BioRender.com/a1bmcre. (Panel B) Comparison of the total amino acid length of the four CatSper pore-forming subunits across the species shown in (Panel A), highlighting pronounced divergence of CatSper1 relative to CatSper2–4. (Panel C) Total CatSper1 length plotted against fertilization temperature for the species shown in (A), revealing an exponential relationship. (Panel D) Length of the CatSper1 NTD plotted as a function of fertilization temperature for the same species, demonstrating a strong positive correlation. (Panel E) Linear regression analysis using 32 (n) internally fertilizing species with fertilization temperatures from ~30 °C and higher. Pearson correlation analysis shows significant positive correlation between the variables (Pearson’s r = 0.675, n = 32, P < 0.0001, R2 = 0.456. Dotted line shows x-intercept at 33.6 ˚C. (Panel F) The percentage of histidine residues within the CatSper1 N-terminal domain shows a strong positive correlation with fertilization temperature across CatSper-bearing species shown in Fig 1A. Outliers (blue) correspond to mole rat species, which inhabit subterranean, hypercapnic environments that result in chronically acidic systemic conditions due to elevated carbon dioxide levels. The unusually high apparent percentage of N-terminal histidines in corals (Paramuricea clavata) reflects the exceptionally short length of their CatSper1 N-terminal domain (9 residues). Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.g001

Because fertilization temperature spans a wide evolutionary range and NTD length varies substantially across taxa, we used Bayesian regression to quantify uncertainty and assess the robustness of the temperature–NTD length association (S2 Fig). This analysis supported a strong positive posterior association between fertilization temperature and CatSper1 NTD length, with the posterior distribution of the slope parameter overwhelmingly shifted toward positive values (P(β > 0 | data) > 0.99). Cold-fertilizing species, such as marine invertebrates and fishes (6–23 °C) exhibited short N termini (≈ 9–120 aa), whereas warm-fertilizing species (35–40 °C) displayed substantially expanded NTDs (≈ 200–550 aa), with intermediate temperatures showing intermediate NTD lengths. This relationship was consistent across phylogenetic distances (Figs 1D, S2A and Table 1).

To further evaluate the relationship between fertilization temperature and CatSper NTD length, we performed linear regression analysis on 34 internally fertilizing species with fertilization temperatures of approximately 30 °C or higher (Table 1). Pearson correlation analysis revealed a moderate positive correlation between the two variables (Pearson’s r = 0.510, n = 34 species, P = 0.002; R2 = 0.260; S2B Fig and Table 1). Because two mole rat species (Heterocephalus glaber and Fukomys damarensis) exhibited exceptionally low fertilization temperatures relative to the remaining mammals and represented clear biological outliers, we repeated the analysis after excluding these species. This analysis revealed a stronger positive correlation (r = 0.675, n = 32, P < 0.0001; R² = 0.456; Fig 1E). Interestingly, the estimated x-intercept of the regression was 33.6 ˚C, closely matching the reported threshold temperature for CatSper activation [17].

CatSper1 NTD’s extraordinary variability suggests it may serve as an evolutionary tuning module, potentially enabling species-specific modulation of CatSper activation in response to reproductive thermal environments. The presence and functional importance of the CatSper1 NTD, that features rich histidine content, has puzzled researchers since its original discovery [1,9]. While it was initially proposed to play a key role in CatSper pH gating, no experimental evidence has yet confirmed this function. The NTD of CatSper1 represents an extreme example of histidine enrichment, exhibiting the highest histidine density reported among pH-sensitive ion channels (Table 2). In mouse CatSper1, histidine accounts for 9.2% of the full-length protein, whereas the N-terminal domain alone contains ~17.2% histidine. This level of enrichment places the CatSper1 N-terminus within the range of classical histidine-rich proteins, such as histidine-rich calcium-binding protein (HRC) and histidine-rich glycoprotein (HRG), which contain ~11%–12.7% histidine (Table 2). This histidine enrichment of N-terminal region was absent from other pore forming subunits of CatSpermasome, such as CatSper2–4 across the entire evolutionary tree (Figs 1B and S3A).

thumbnail

Table 2. Histidine (His, H) abundance in known histidine-rich proteins. Histidine-rich proteins are defined by an unusually high abundance of histidine residues. In the average proteome, histidine constitutes approximately 2%–3% of amino acids; by contrast, histidine-rich proteins typically contain ≥5%–10% histidine, with some domains exhibiting substantially higher levels. Among all pH-sensitive ion channels (blue), only CatSper1 and CatSper N-terminus (CatSper1-NT; red) shows high His content, on par with other known His-rich proteins.

https://doi.org/10.1371/journal.pbio.3003776.t002

In many histidine-rich proteins, clusters of His residues confer sensitivity to pH, metal ions, and other environmental cues, prompting the hypothesis that CatSper1’s N-terminal His clusters contribute to temperature-dependent gating. To test this, we first quantified the total histidine content and found a strong positive correlation with species-specific fertilization temperatures (Fig 1F), whereas the histidine content of CatSper2–4 showed no correlation (S3A Fig). Interestingly, an inverse relationship was observed between fertilization temperature and the predicted isoelectric point of the CatSper1 N-terminus (S3B Fig), consistent with increasing histidine enrichment conferring enhanced pH sensitivity at higher temperatures.

Secondary structure analysis of the CatSper1 NTD indicates that it possesses hallmark features of an intrinsically disordered region (IDR). These are protein regions that are defined by sequence-specific conformational heterogeneity. Consistent with this prediction, recent high-resolution cryo-EM and in situ 3D cryo-ET structures of the murine CatSper complex [21,31], obtained at 25 °C, resolved the architecture of most CatSpermasome subunits but failed to visualize the CatSper1 NTD. The absence of defined density for this region further supports the notion that the NTD is intrinsically disordered.

IDRs are defined by “molecular grammars”. This refers to non-random amino acid compositions (the alphabet) and non-random patterning of distinct residue types with respect to one another (the syntax) [42,43]. Given the disordered nature of CatSper1 NTD, we performed an unbiased analysis of molecular grammars using the NARDINI+ algorithm. This algorithm takes an IDR sequence as input and quantifies 90 different z-scores for different compositional and patterning features. Specifically, NARDINI+ identifies non-random grammar features from z-scores calculated by comparing compositional features to the human IDRome and patterning features to 105 composition matched scrambles. When applied to the CatSper1 NTD sequences across evolution, the NARDINI+ analysis helped in identifying features of the IDR grammars that differed among warm-, intermediate, and cold-fertilizing species. We find that 24 grammar features have distinct z-score distributions between fertilization temperatures as assessed via the Kolmogorov–Smirnov test (Fig 2A–2C).

thumbnail

Fig 2. CatSper1 histidine content positively correlates with fertilization temperature.

(Panel A) Heatmap of NARDINI+ z-scores for CatSper1 N-terminal domains per species for grammar features highlighted in (Panel B). Green denotes features enriched and purple denotes features depleted relative to the human IDRome. Species were split into three fertilization temperature classes: cold (6–23 °C), intermediate (25–34.5 °C), and warm (35–40 °C). Colored boxes around the species names denote cold- (blue), intermediate- (green), and warm- (red) fertilizing species. (Panel B) Boxplots of z-scores extracted from NARDINI+ for grammar features within CatSper1 N-terminal domain that show significant differences in z-score distributions between warm-, intermediate-, and cold-fertilizing species (*: p-value < 0.05 and ≥ 0.001, **: p-value < 0.001). Z-scores greater than zero are enriched compared to the human IDRome, whereas z-scores less than zero are depleted compared to the human IDRome. The naked mole rat (Heterocephalus glaber) is shown as an outlier intermediate-fertilizing species by blue circles. (Panel C) Comparison of NARDINI+ z-score distributions of CatSper1 NTDs for grammar features between fertilization temperatures. The heatmap denotes the signed log10(p-value) as assessed via the Kolmogorov-Smirnov test. Only grammar features with p-values < 0.05 are shown. Here, positive values (blue-green) denote that the temperature class listed first had a mean z-score greater than or equal to the temperature class listed second. Negative values (brown) denote that the temperature class listed first had a mean z-score less than the temperature class listed second. (Panels D-G) AlphaFold3 modeling of the murine CatSper1 N-terminal region. (Panel D) The AlphaFold3 predicted structure of a lone murine CatSper1, shown in red. NTD is not predicted to fold into any secondary structures. (Panel E) AlphaFold3 predicted structure of the core subunits CatSper1-4, with CatSper1 shown in red. The NTD is still unstructured. (Panel F) The AlphaFold3 predicted structure of the transmembrane and cytoplasmic parts of the CatSper complex when including intracellular subunits and the transmembrane helices of the extracellular subunits. The NTD is predicted to sit in a region occupied by the cell membrane. (Panel G) The AlphaFold3 predicted structure of the transmembrane and cytoplasmic parts of the CatSper complex as in (Panel E), but with 50 oleic acid molecules added to the prediction. The free oleic acid clusters around the transmembrane domain, causing the NTD to adopt an intracellular position. Each structure is shown with a 100 Å scalebar. Dotted lines show plasma membrane bilayer. (Panel H) Integrative structural model of multiple adjacent murine CatSpermasomes fitted into the cryo-ET density of murine CatSpermasome from [21]. (Panel I) Same as (H) in the absence of Cryo-ET density. Single-particle cryo-EM structures and AlphaFold3-predicted models were combined to visualize the arrangement of pore-forming and auxiliary subunits within the array. Subunits are color-coded as indicated: CatSper1 (red), CatSper2 (green), CatSper3 (gray), CatSper4 (purple), CatSperβ (light pink), CatSperγ (yellow), CatSperδ (dark gray), and CatSperε (light green). Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.g002

The z-scores for the fraction of histidine residues were significantly different and more enriched for warm-fertilizing species compared to intermediate- and cold-fertilizing species (Fig 2A–2C). We find that intermediate-fertilizing species are generally enriched in hydrophobic residues, whereas warm-fertilizing species are generally depleted in residues deemed to be hydrophobic (Fig 2A–2C). This result suggests that species that fertilize at intermediate temperatures may exchange histidine residues for hydrophobic residues as the drivers of temperature-dependent gating. Specifically, the strength of hydrophobic interactions can increase with increasing temperature [44]. Lastly, cold-fertilizing species show unique grammar features including an increase in the fraction of charged residues (FCR) and a decrease in disorder promoting residues compared to warm- and intermediate-fertilizing species.

In addition, the CatSper1 NTD displays a high predicted droplet-promoting propensity (pDP score approaching 1.0), suggesting a strong likelihood of liquid–liquid phase separation behavior. High pDP values (~1.0) are characteristic of proteins capable of droplet formation and self-assembly into higher-order condensates, processes that are often sensitive to environmental factors such as temperature and ionic strength. Thus, phase separation mediated by the CatSper1 IDR may contribute to the distinctive regulatory properties of CatSper1 (S4 Fig).

Comparative genomics combined with AlphaFold3 modeling further supported a relationship between CatSper1 surface-exposed histidine residues and fertilization temperature

Next, to determine whether the CatSper1 IDR can adopt a stable three-dimensional structure under physiological conditions, we integrated AlphaFold3 predictions with available structural data. Specifically, we modeled the murine CatSper1 NTD in isolation, in the context of neighboring CatSpermasome components, and within a lipid membrane environment to assess whether molecular context promotes structural stabilization (Fig 2D–2I).

As shown in Fig 2D, initial modeling of CatSper1 NTD in isolation did not predict the formation of defined secondary structural elements, consistent with its intrinsically disordered nature. However, when the transmembrane and cytoplasmic components of the CatSper complex were included, comprising the intracellular subunits and the transmembrane helices of extracellular subunits, the NTD was positioned within a region corresponding to the plasma membrane interface (Fig 2E–2F). Upon further incorporation of free oleic acid molecules clustering around the transmembrane domain to better approximate a lipid environment, the NTD adopted an intracellular orientation (Fig 2G). This positioning is consistent with its proposed localization between adjacent CatSpermasomes [32,45,46] and enabled quantitative assessment of the density of surface-exposed histidine residues within the NTD. The calculated histidine surface density was then plotted against species-specific fertilization temperatures, revealing a strong positive correlation across mammals (S3C Fig).

Two exceptions- Heterocephalus glaber and Fukomys damarensis- possess extremely high N-terminal His content (~26%) yet intermediate core body temperatures (30–36 °C; Figs 1E–1F, 2A–2C, S2B, S3B–S3C, S5 and Table 1). These subterranean mole rats inhabit hypoxic, CO2-rich burrows, likely imposing unique selective pressures on sperm physiology. Such an environment may have contributed to lineage-specific adaptations in CatSper channel function and regulation, setting these species apart from other mammals in the evolutionary trajectory of sperm motility control.

To identify a potential interaction module between CatSper1 N-termini, we fitted the high-resolution cryo-EM structure of the CatSpermasome [31] into the density maps of CatSpermasome arrays previously resolved by cryo-electron tomography (cryo-ET [21]. The predicted intracellular and transmembrane architecture of the CatSper complex was then integrated into the resulting 3D reconstruction (Figs 2H–2I and 3).

thumbnail

Fig 3. Structural modeling identifies a putative CatSper1 N-terminal interaction interface between neighboring CatSpermasomes.

(Panels A–C) Integrative model of two adjacent CatSpermasomes viewed from the side (Panel A), bottom (Panel B), and top (Panel C). Models were generated by fitting the high-resolution single-particle cryo-EM structure of the CatSper complex from [31] into the in situ cryo-electron tomography (cryo-ET) density of CatSper arrays from [21], followed by incorporation of AlphaFold3-predicted transmembrane and intracellular regions. Pore-forming subunits are color-coded: CatSper1 (red), CatSper2 (green), CatSper3 (gray), and CatSper4 (purple). Dashed boxes indicate the region of closest approach between adjacent CatSpermasomes. Only first 290 amino acids of CatSper auxiliary subunits β, δ, ε, and γ are shown. (Panel D) Magnified view of the boxed region in (Panel B), highlighting paired α-helical elements from the CatSper1 NTDs positioned at the inter-complex interface. Histidine side chains are shown in blue, arginine side chains are shown in red. Dashed yellow lines indicate distances between conserved His-His residues and His-Arg residues across neighboring complexes: a, His15–His321 (~12.1 Å); b, Arg7–His321 (~4.3 Å); c, Arg6–His317 (~3.7 Å); d, His317–His317 (~14.2 Å); e, His321–Arg7 (~4.3 Å); and f, His321–His15 (~12.1 Å). (Panel E) Magnified view of the boxed region in (C), illustrating the spatial distribution of surface-exposed histidine residues along the interacting helices. Scale bars, 10 nm (Panels A–C) and 1 nm (Panels D, E).

https://doi.org/10.1371/journal.pbio.3003776.g003

This composite model revealed a previously unrecognized putative interaction interface formed by the proximal regions of CatSper1 N-terminal domains, spanning residues 6–321 in murine CatSper1. These segments are positioned at the junction between adjacent CatSpermasomes within the longitudinal arrays, suggesting a role in inter-complex coupling. AlphaFold3 predicts that this interaction module is composed of paired α-helical elements contributed by opposing CatSper1 N-termini. Each helix contains three conserved histidine residues (H15, H317, and H321; Figs 3 and S5). The spatial arrangement of these residues defines a potential interaction network, including contacts between H15 from one CatSpermasome and H321 from an adjacent complex, as well as reciprocal interactions between opposing H317 residues (Figs 3A–3C and S5).

Emerging evidence indicates that histidine residues within IDRs can mediate context-dependent stabilization [47–49]. In addition to π–π or His–His interactions, neutral histidines at pH values above ~6.8 can participate in strong cation–π interactions with positively charged residues such as arginine and lysine, promoting local structural ordering. Notably, two arginine residues (Arg6 and Arg7) are located within the predicted murine CatSper1 NTD interface. The modeled distances between these arginines and neighboring histidines (3.7–4.3 Å) fall within the range compatible with stabilizing cation–π interactions (3–6 Å). Such proximity suggests a potential mechanism for pH-dependent inter-CatSpermasome coupling mediated by Arg–His interactions (Figs 3A–3C and S5–S6).Our modeling is based on the available cryo-ET [21] map (Fig 2G) obtained at room temperature (25 °C), which is below the activation threshold of murine CatSper [17]. This likely explains the relatively long distances observed between putative interacting histidines (H317–H317 and H15–H321), which range from ~11.8 to 14.2 Å—too large to support direct π–π or hydrogen-bond interactions. We therefore propose that elevated temperature promotes histidine deprotonation, reducing electrostatic repulsion and enabling closer association of CatSper1 NTDs. Consistent with this idea, phylogenetic analysis reveals that H317 (in murine CatSper1) is the most conserved histidine residue within this interface, with notable divergence in cold-blooded internal fertilizers such as reptiles, supporting a role in temperature-dependent regulation (Figs 3 and S5). Together, these findings identify a histidine-arginine-based interface that may mediate pH- and temperature-dependent coupling between CatSpermasomes.

Histidine-rich proteins are known to serve diverse regulatory functions, including pH buffering, metal ion coordination, and modulation of protein stability. Importantly, the pKa of histidine (~6.3) [50] allows it to readily switch between charged and neutral states within a narrow physiological pH range (6.5–7.5), which coincides with the conditions that regulate CatSper activity.

In temperature-sensitive proteins, increased temperature can shift local pKa values, weakening hydrogen bonding networks and promoting histidine deprotonation. Conversely, under acidic conditions, histidine protonation leads to charge accumulation that can destabilize protein structure by introducing electrostatic repulsion and disrupting hydrophobic packing. In fact, as we have shown recently, cytoplasmic acidification to 6.0 profoundly attenuated CatSper’s temperature gating [17], which further supports CatSper N-terminal involvement in this process.

Experimental evidence indicates that removal of CatSper1 N-terminus results in the loss of temperature-dependent gating

In our recent work, we demonstrated that murine spermatozoa subjected to prolonged capacitation (>60 min) progressively lose temperature-dependent gating of CatSper [17]. This loss is accompanied by a time-dependent decline in CatSper-mediated currents for both divalent (Ba2+; Fig 4A–4C) and monovalent (Cs+) [17]. Using direct sperm patch-clamp recordings, we observed a progressive reduction in temperature-dependent divalent conductance through CatSper and significant decrease of its temperature sensitivity after prolonged sperm capacitation (Fig 4A–4C).

thumbnail

Fig 4. Sperm capacitation leads to removal of temperature sensitivity.

(Panel A) Decrease in current density of Ba2+ currents via CatSper with sperm capacitation (red; cap) after response to elevated temperatures. CatSper1−/− sperm (green, noncap) were used as negative control. Noncapacitated sperm (noncap) data are shown in blue. (Panel B) Representative divalent Ba2+ currents (IBa2+) recorded from wild type noncapacitated mouse spermatozoa (left) and from wild type capacitated sperm cell (right) in response to indicated voltage steps under two different temperatures: 24 °C and 38 °C. (Panel C) Exposure to 38 °C increases IBa2+ stronger in noncapacitated cells. Data are mean values ± S.E.M., n corresponds to the number of cells used; n = 3-6; total capacitation time more than 60–90 min. (Panel D) Representative cryo-electron tomography (cryo-ET) slice of the human CatSper channel array in the sperm flagellum from [21], viewed parallel to the principal piece axis. The dashed yellow box highlights a repeating CatSpermasome unit within the longitudinal zigzag array. (Panels E, F) Model for pH- and temperature-dependent coupling of CatSpermasomes mediated by the CatSper1 N-terminus. (Panel E) Schematic representation of CatSpermasome organization within the zigzag array under acidic intracellular pH and/or cooler temperature conditions. Under these conditions, protonation of histidine residues within the CatSper1 NTDs (1NT) is proposed to weaken 1NT interactions and reduce cooperative channel coupling: array uncoupled. Right insert: representative small CatSper current is shown. (Panel F) Corresponding schematic demonstrates array coupling under alkaline intracellular pH and warm temperature conditions. Deprotonation of histidine residues will ensure histidine-arginine coupling, stabilizing interactions between neighboring CatSpermasomes, and promoting cooperative coupling along the array. Representative enlarged CatSper current is shown on the right. Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.g004

To independently validate these findings, we performed calcium imaging experiments to assess CatSper activity in response to its threshold activation in non-capacitated versus capacitated sperm. In agreement with electrophysiological recordings, sperm subjected to prolonged capacitation (80–110 min) failed to exhibit temperature-dependent Ca2+ influx. CatSper1-deficient sperm were used as controls to exclude non-CatSper-mediated effects. As shown in S7A–S7B Fig, capacitation completely abolished temperature-dependent activation of murine CatSper. It must be noted that the optimal duration of capacitation differs among mammalian species and that mouse sperm undergo capacitation over a substantially shorter time scale than sperm from larger mammals, such as bovine or human sperm [8,51–53].

Prior studies from Chung and colleagues and Ded and colleagues [14,39] convincingly demonstrated that these capacitating conditions induce selective proteolytic removal of at least ~150 amino acids, (~20 kDa) fromCatSper1 N-terminus [14,39], temporally coinciding with the loss of temperature gating observed in our experiments (Figs 4A–4C and S7A–S7B). The cleaved region contains the majority of the histidine-rich N-terminal segment, as well as N-terminal arginines, thereby removing a substantial fraction of the residues that is proposed to contribute to temperature-sensitive interactions. Together, these findings support the conclusion that temperature gating is critical for maintaining CatSper channels organized within CatSpermasomes in a highly ordered and synchronized array, enabling their concerted opening. Disruption of this domain uncouples channels from the ordered CatSpermasome lattice, leading to asynchronous or inefficient channel opening. Loss of concerted gating reduces effective channel conductance and limits Ca2+ influx, thereby impairing the robust calcium signaling required to drive hyperactivated sperm motility.

Based on these principles, we propose that the high histidine content of the CatSper1 N-terminal domain paired with adjacent positively charged residues (Arg, Lys) confers sensitivity to temperature-dependent histidine deprotonation, stabilizing inter-CatSpermasome interactions in a manner analogous to alkaline pH. In contrast, intracellular acidification would promote histidine protonation, destabilizing the N-terminal interface and inhibiting channel activity (Fig 4D–4F). This model is fully consistent with experimental observations showing that CatSper is activated by warm temperatures (>34 °C) [17] and alkaline pH [6,53], while acidic pH potently suppresses channel function.

According to this model, pH- and temperature-dependent deprotonation of surface-exposed histidine residues within the CatSper1 N-termini, including His317 and His321, stabilizes cation–π interactions with Arg6 and Arg7, and promotes coupling between adjacent CatSpermasomes, thereby enabling synchronized channel activation at elevated temperatures (Fig 4D–4F).

Discussion

In addition to regulation by pH and temperature, CatSper activity is strongly modulated by membrane voltage, such that progressively greater depolarization is required to open the channel under acidic conditions [6,17,19,53,54]. Conversely, at warmer temperatures, CatSper opens at less depolarized potentials, whereas acidic pH completely abolishes temperature-dependent gating. This interdependence between temperature, pH, and voltage suggests a shared structural basis for polymodal regulation.

We propose that this integration is mediated by the spatial proximity of the CatSper1 N-terminal interaction module to the S4 segment of the CatSper1 voltage-sensor domain (VSD). The S4 segment of CatSper1 is the strongest voltage-sensing element among the pore-forming subunits, harboring an unusually high density of positively charged residues—four arginines (Arg; R) and two lysines (Lys; K)—that together form a classical voltage-gated motif (RxxRxxRxxR; S7C Fig and [9,31]. This charge density is predicted to confer heightened voltage sensitivity relative to CatSper2–4, whose S4 segments contain progressively fewer charged residues, with CatSper4 harboring only two Arg and one His. These structural differences point to a hierarchical organization of gating contributions, with CatSper1 acting as the dominant voltage sensor and the remaining subunits serving modulatory roles.

The close apposition of the CatSper1 VSD and NTD suggests that these elements function as a coupled gating unit, integrating voltage sensing with temperature- and pH-dependent regulation. Upon membrane depolarization, upward movement of the S4 helix is expected to mechanically displace the adjacent NTD. Based on the predicted organization of CatSpermasome arrays, neighboring CatSper channels may interact through their NTD, potentially via coiled-coil–like associations. In this configuration, depolarization-induced S4 movement would impose mechanical tension on the N-terminal interface, promoting coordinated conformational transitions and priming adjacent channels for synchronized opening.

Inter-CatSpermasome coupling interfaces have previously been identified within the canopy region and at the interface between the CatSper4 and CatSperβ subunits [45,46]. However, the interface formed by the CatSper1 NTD has remained unresolved and was largely considered an intrinsically disordered region with unknown functional properties. To our knowledge, this is the first report assigning a specific structural and functional role to the CatSper1 NTD in coordinating channel opening within longitudinal CatSper arrays.

Recent quantum chemical analyses by Calinsky and colleagues [47–49] have demonstrated that histidine residues within intrinsically disordered domains can mediate pH-dependent structural assembly through dual π–π and cation–π interactions. At physiologically relevant alkaline pH, deprotonated (neutral) histidine can form stabilizing cation–π interactions with positively charged residues such as arginine (Arg) and lysine (Lys). Notably, two arginine residues are present within the predicted interface between CatSper1 N-termini, supporting the possibility of such stabilizing interactions under alkaline conditions.

In contrast, at acidic pH, protonated histidine preferentially engages in cation–π interactions with aromatic residues (e.g., Tyr, Phe, Trp), interactions that are typically stronger than those formed by neutral histidine with cationic residues. Importantly, aromatic residues are absent from the predicted CatSper1 NTD interface, suggesting that acidification would not promote alternative stabilizing contacts.

Although the arginine residues predicted to mediate cation–π interactions within the murine CatSper1 NTD appear to be poorly conserved across warm-fertilizing species, their functional role may be substituted by lysine residues (S6 Fig). Lysines are more broadly conserved among these species and are positioned in close proximity to the predicted interface region, making them plausible alternative partners for stabilizing cation–π interactions with histidines. Nevertheless, high-resolution atomic structures of CatSper complexes from additional species will be required to directly test and validate this hypothesis.

An additional important variable that may influence histidine abundance in CatSper1 is pH of female reproductive tract pH that varies during estrous and menstrual cycles. Indeed, during the fertile period, cervical mucus becomes markedly more alkaline, reaching pH values as high as ~9.0 in some species, thereby facilitating sperm transport [53,55–57]. Such physiological pH fluctuations could influence the protonation state of histidine residues and thereby modulate the function of pH-sensitive proteins.

Histidine is a particularly effective physiological buffer because its imidazole side chain has a pKa of approximately 6.0–6.5, close to physiological pH. As environmental pH decreases, histidine residues become reversibly protonated, transiently binding excess H+ and thereby increasing buffering capacity. Consequently, proteins enriched in histidine residues generally exhibit greater tolerance to transient acidification, as reversible protonation helps preserve protein structure and function over a wider pH range.

In such cases, species that inhabit chronically hypoxic and hypercapnic underground burrows, such as naked mole rat (Heterocephalus glaber), could present a condition that is expected to promote tissue acidification, particularly of the female reproductive tract. Combined with reduced thermogenesis and lower core body temperature reported for these species, these environmental pressures may have favored proteins with enhanced buffering capacity. We therefore suggest that the increased abundance of histidine residues in naked mole rat CatSper1 could provide greater local buffering, helping to stabilize the channel and preserve its function under conditions of fluctuating pH. At present, however, this hypothesis cannot be directly evaluated because reproductive tract pH has not been measured in naked mole rats.

Together, these findings identify a histidine-dependent coupling interface positioned adjacent to the dominant voltage-sensing element of the CatSpermasome, providing a structural framework for the coordination between temperature, pH, and voltage in CatSper gating. In this model, acidic pH is predicted to disrupt N-terminal coupling by protonating histidine residues within the interface, increasing electrostatic repulsion and weakening inter-channel interactions. This would reduce coupling efficiency, leading to asynchronous channel opening and diminished cumulative Ca2+ influx. Conversely, elevated temperature and alkaline pH would favor histidine deprotonation, stabilize N-terminal interactions, and enhance cooperative and, moreover, synchronized gating within CatSper arrays (Fig 4E – 4F). More broadly, this work establishes CatSper as a new model for how evolutionary pressures sculpt ion-channel architecture to integrate environmental and physiological signals essential for optimal sperm performance.

As we recently reported [58], the estimated unitary conductance of mammalian CatSper appears to be remarkably low. This estimate is based on the ~6,000 CatSper complexes present along the murine sperm flagellum, as determined by cryo-EM and cryo-electron tomography analyses [21,31,46,58]. Assuming this channel density, the calculated unitary conductance of an individual CatSpermasome ranges from ~121 fS to 1.2 pS, indicating that CatSper is an exceptionally low-conductance, yet highly calcium-selective, ion channel. Given such a small unitary conductance, the large macroscopic CatSper currents routinely recorded from mammalian sperm are unlikely to result from the stochastic opening of individual channels. Instead, these currents likely require the coordinated activation of thousands of CatSper complexes assembled into supramolecular arrays along the flagellum. It is therefore feasible to assume that cooperative coupling between neighboring CatSpermasomes promotes their synchronized opening, providing a mechanistic explanation for how a sperm cell can generate robust macroscopic calcium currents despite the extremely low conductance of individual CatSper channels.

The model presented here is grounded in comparative genomics, structural modeling, and integration of high-resolution cryo-EM and cryo-ET data with functional analysis of murine sperm cells. However, definitive validation of the proposed histidine-dependent coupling mechanism will require atomic resolution of CatSper conformational dynamics across defined pH, temperature, and voltage conditions. In particular, high-resolution structural snapshots captured under physiologically relevant temperatures and targeted mutagenesis of conserved histidine residues will be essential to resolve the precise molecular rearrangements underlying synchronized gating. Such experiments are beyond the scope of the present study but are a natural and necessary extension of this work, and we anticipate that future experimental structures will further refine and strengthen the mechanistic model proposed here.

Methods

Comparative genomics and sequence alignment

Protein sequences for CatSper subunits and related proteins across a wide range of species were retrieved using Python Scripts accessing the UniProtKB, UniParc, and NCBI APIs. Data extraction was done in June 2025. Sequences labeled as partial or containing the names of other subunits were removed. In cases where a single species had multiple isoforms, the longest sequence was selected. Suspected misannotated sequences were visually identified in an alignment and replaced using NCBI Protein BLAST searches with the corresponding mouse subunit. Species with at least 4 pore-forming subunits (Catsper1–4) conserving pore-forming subunits (Catsper1–4) were chosen for sequence alignment. An evolutionary tree was built using NCBI Common Tree [PMID: 32761142].

Definition of the CatSper1 N-terminal region

To define the boundaries of the Catsper1 N-terminal region, all full-length Catsper1 protein sequences from the dataset were aligned using MUSCLE (v5.2) [102] with default parameters. The resulting alignment was used to identify the start of the conserved pore-forming region, defined as the first position in the alignment containing a block of at least 10 consecutive columns with greater than 80% occupancy. For each individual sequence, the non-gap characters in the alignment preceding this position were used as the N-terminal region. The N-terminal sequence spanned from the beginning of the sequence to the beginning of the short 20 amino acid (murine) linker region connecting the N-terminal domain to the S1 transmembrane domain.

Quantification of N-Terminal histidine content

To calculate N-terminal histidine abundance, the number of histidine residues in an N-terminal sequence was normalized to the total number of amino acids in the N-terminus.

Computation of physicochemical properties of the CatSper1 N-terminus

Theoretical isoelectric points for each N-terminal sequence were calculated using the Expasy ProtParam web server [103]. Liquid-Liquid Phase Separation (LLPS) probabilities were calculated using the FuzDrop [104] web server, using the mouse N-terminal sequence as a representative model. A pDP score greater than 0.6 was used as the threshold for identifying likely LLPS regions.

Quantification of grammar differences between warm-, intermediate, and cold-fertilizing species

To perform an unbiased assessment of grammar differences between species that fertilize at different temperatures, we ran NARDINI+ utilizing NARDINI 2.0 (May 2026) [42,43]. NARDINI+ quantifies z-scores of 54 compositional and 36 patterning grammar features in order to determine non-random grammars within an IDR sequence. This analysis was only performed on CatSper1 NTDs of length ≥ 30 to avoid erroneous interpretations from short sequences when calculating features based on fractions of a given residue or residue type.

The 54 compositional features are: the fraction of each of the 20 amino acids (20 features); patch parameters defined as the fraction of the IDR comprised of a given residue or RG stretches (20 features); the fraction of polar (Asn, Cys, Gln, Gly, His, Ser, Thr), aliphatic (Ala, Ile, Leu, Met, Val), aromatic (Phe, Trp, Tyr), K + R (Lys, Arg), D + E (Asp, Glu), charged (FCR), chain expanding (FCE), and disorder promoting residues (8 features); the ratios of Arg to Lys (R/K Ratio) and Glu to Asp (E/D ratio) (2 features); and the net charge per residue (NCPR), isoelectric point, hydrophobicity, and polyproline II propensity (4 features) [105]. Z-scores for these features are calculated by comparing the feature value of the given CatSper1 NTD to the feature distribution from the human IDRome. The human IDRome contains 24,508 IDRs of length ≥ 30 and was extracted using MobiDB on the human proteome (Swiss-Prot database) [43,106,107].In general, for compositional features, a z-score value ≥ 1 implies that feature is non-randomly enriched and a z-score value ≤ −1 implies that feature is non-randomly depleted compared to the human IDRome.

The 36 patterning features consider the patterning between eight types of residues: pol = (S, T, N, Q, C, H), hyd = (I, L, M, V), pos = (K, R), neg = (E, D), aro = (F, W, Y), ala = (A), pro = (P), and gly = (G). To calculate non-random grammars of patterning features the sequence of each CatSper1 NTD studied was scrambled 105 times. The composition matched scrambles form the prior distribution of patterning of two types of residues. For patterning features, z-scores less than zero imply the input sequence is more well-mixed than expected and z-scores greater than zero imply the input sequence is blockier than expected.

To determine if CatSper1 NTDs from warm (35–40 °C), intermediate (25–34.5 °C), and cold (6–23 °C) fertilizing species house different grammars, we compared the z-score distributions for all 90 features between warm-intermediate, warm-cold and intermediate-cold fertilizing species using the Kolmogorov-Smirnov test. We identified 24 grammar features have distinct z-score distributions between fertilization temperatures (p-value < 0.05). All of the significantly different features were compositional features.

Catspermasome structural modeling and bioinformatic analysis

Structural models of the transmembrane and intracellular parts of the CatSper complex were generated using AlphaFold3 [40] In cases where model inference failed, auxiliary subunits distant from the N-terminal region were excluded to improve model inference. The last 280 residues of each of the extracellular subunit sequences, containing the transmembrane helices, were included in the AlphaFold3 prediction as separate chains. Fifty oleic acid ligands were added to the Alphafold3 prediction to mimic the membrane environment. The N-terminal domain was only predicted to fold into secondary structures when the transmembrane helices of the extracellular subunits were included in the prediction (S5 Fig). Five models were generated for each complex. For each species, the model with the highest-ranking score as calculated by Alphafold3 was used for all further analysis.

Single-particle cryo-EM maps (EMD-31076) [31] were fitted as rigid bodies into the cryo-ET structure for the in situ arrangement of the CatSper complex (EMD-24210) [21] using ChimeraX [108]. Alphafold3 predicted structures for the transmembrane and intracellular regions of the CatSper complex were fitted to these cryo-EM maps to visualize the relative orientations of complexes. A residue was defined as being on the surface if its Relative Solvent Accessibility (RSA) was greater than 25%. RSA was calculated for each residue using the formula RSA = (SASA / MaxSASA) * 100. Only residues in the N-terminal region were accounted for in the calculation. The Solvent Accessible Surface Area (SASA) was calculated for each PDB model using the Shrake-Rupley algorithm as implemented in the Bio.PDB.SASA package in Biopython (v1.85) [109]. Maximum SASA values (MaxSASA) for each amino acid were taken from the reference values published by Tien and colleagues (2013) [110].

Bayesian regression analysis

To assess the relationship between fertilization temperature and CatSper1 N-terminal (NT) length across species, we performed a Bayesian regression analysis using data summarized in Table 1 and Fig 1d. NT length was modeled as a shifted exponential function (NTi) of fertilization temperature:

where represents the standardized temperature , α represents the baseline NT length, β represents the scaling parameter of exponential growth, γ represents the growth rate parameter, and σ represents the residual variance. Priors were chosen to be weakly informative and biologically plausible, reflecting observed NT length ranges:

α ~ N (100,50),

β ~ HalfNormal (50),

γ ~ LogNormal (0,1),

σ ~ HalfNormal (50).

Posterior distributions were estimated by Markov chain Monte Carlo sampling using the No-U-Turn Sampler implemented in PyMC [111]. Four chains were run for 3,000 iterations, each following 2,000 tuning steps with a target acceptance rate of 0.95. Inference was based focused on the posterior probability distributions of all parameters.

Linear regression analysis

Linear regression was performed to evaluate the relationship between fertilization temperature and CatSper NTD length across internal fertilizers. Pearson's correlation coefficient (r), the coefficient of determination (R2), and two-tailed P values were calculated to assess the strength and significance of the association. Regression analyses were performed using Origin 9.0 Software. Outlier-excluded analyses were performed after removal of the two mole rat species (Heterocephalus glaber and Fukomys damarensis), which exhibit exceptionally low fertilization temperatures compared with the other internally fertilizing mammals analyzed.

Data analysis

All the data are presented as mean values ± S.E.M., unless otherwise specified, and the number of individual cells is shown as (n). For current densities calculations, the CatSper Ba2+ amplitudes were normalized to capacitance (Cm, pF); during data analysis capacitance artifacts are graphically removed. Post-hoc Tukey’s test was used for comparison of mean values and statistical significance.

Animals and ethics statement

All mouse procedures were performed according to the National Institutes of Health Guiding Principles for the care and use of laboratory animals. These procedures were reviewed and approved by the Institutional Animal Care and Use Committee of Washington University in St. Louis (St. Louis, MO, USA) (protocol # 20-0126 and protocol 22-0251). C57BL/6, C57BL/6NJ, and CatSper1 knock-out male mice (CatSper KO) (60–90 days old) were purchased from Charles River and Jackson Laboratory (JAX) and kept at a constant temperature of 22 ± 2 °C under a 12/12 h dark/light cycle with free access to food and water at the barrier animal facility of WashU School of Medicine. The animals were fed a standard chow diet. There was no difference in the sperm quality between mouse lines. CatSper1−/− males were born at a mendelian ratio by breeding CatSper1−/− females with CatSper1+/− males. Genotypes were continuously monitored by TransnetYX (TN) and verified by in-house genotyping as reported in [17].

Mouse sperm collection and capacitation

Sexually mature male mice (60–90 days old) were euthanized using CO2 followed by cervical dislocation, and sperm were isolated from the cauda epididymis. Sperm were allowed to swim-up in non-capacitating (NC) TYH media buffered with HEPES (in mM: 119.3 NaCl, 4.7 KCl, 1.71 CaCl2.2H2O, 1.2 KH2PO4, 1.2 MgSO4.7H2O, 25.1 NaHCO3, 5.56 glucose, 0.51 sodium pyruvate, 10 HEPES) at pH 7.4 for 15–20 min at 37 °C, unless specified. The motile fraction of the sample was removed from the tube. To induce capacitation, sperm were incubated in capacitating (CAP) TYH medium supplemented with 5 mg/mL BSA and 15 mM NaHCO3 at 37 °C for 60–90 min. NC sperm were incubated in TYH without BSA and NaHCO3.

Sperm intracellular calcium concentration measurement

After sperm swim-up in NC or CAP conditions for 20 min, motile sperm were collected and incubated with 10 μM Fluo-4 AM and 0.05% Pluronic F-127 at 37 °C for 60−90 min, under NC or CAP conditions. Sperm cells were centrifuged at 300−400 g for 10 min and further resuspended in the corresponding media. Overall, sperm cells were capacitated for a total time of 80−110 min. Cells (5 X 106/well) were seeded in a 384-well black, clear-bottom plate (Greiner Bio-One) and sperm [Ca2+]i was recorded at 25, 30, and 38 °C using a CLARIOstar plate reader (BMG LABTECH). The instrument was set to bottom-reading mode with an excitation wavelength of 483−14 nm and emission of 530−30 nm. Acquisition parameters used were: 10 flashes per well with a cycle time of 30 s per read. Gain was adjusted to 50% on the highest signal well. The protocol was set to double orbital shaking at 300 rpm for 10 s before each measurement. Data were collected with MARS Data Analysis software. For each condition, an initial 10-minute baseline reading at 25 °C was followed by sequential 10-min recordings at 30 and 38 °C, followed by a final perfusion with 10 μM Ionomycin (Iono) at the end of every experiment. Clampfit 10 (Molecular Devices), and GraphPad were used for data analysis. To compare fluorescence across different samples and experimental conditions, Fluo-4 fluorescence (F) in response to different temperatures for each sperm was normalized to its respective FIono (F/Fiono).

Supporting information

S1 Fig. Comparative genomic analyses of CatSpermasome transmembrane subunits.

(Panel A) Comparison of the total amino acid length of eight CatSpermasome core-forming subunits across species shown on Fig 1A. (Panel B) CatSper1–4 entire length plotted against fertilization temperature for the species shown in Fig 1A. Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.s001

(TIFF)

S2 Fig. Bayesian regression analysis supports strong positive posterior association between fertilization temperature and CatSper1 N-terminal length.

(Panel A) Priors were chosen to be weakly informative and biologically plausible, reflecting observed NT length ranges: α ~ N (100,50); β ~ HalfNormal(50); γ ~ LogNormal(0,1), σ ~ HalfNormal(50). Parameters were inferred using Monte Carlo sampling. Four chains were run for 3,000 iterations, each following 2,000 tuning steps with a target acceptance rate of 0.95. Trace and posterior plots show convergence of the chains and the distribution of model parameters. (Panel B) Linear regression analysis using 34 internally fertilizing species with fertilization temperatures from 30 °C and higher. Pearson correlation analysis demonstrated a moderate positive correlation between the two variables (Pearson’s r value = 0.510, n (number of species) = 34, p = 0.002; R² = 0.260). The analysis includes two outliers (two species of mole rats Heterocephalus glaber and Fukomys damarensis). Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.s002

(TIFF)

S3 Fig. Correlation between CatSper1 N-terminal histidine content, pI and fertilization temperatures.

(Panel A) The percentage of histidine residues within entire CatSper1–4 subunits plotted against fertilization temperature. Only CatSper1 shows positive correlation. (Panel B) Theoretical isoelectric point (pI) calculated for N-terminus domain of CatSper1 plotted against fertilization temperature for mammalian species shown in Fig 1a. (Panel C) Surface His density correlates with fertilization temperatures. Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.s003

(TIFF)

S4 Fig. The CatSper1 N-terminal domain uniquely exhibits high predicted liquid-liquid phase separation propensity.

Predicted liquid–liquid phase separation (LLPS) propensity of CatSper pore-forming subunits. pDP scores calculated using FuzDrop are plotted as a function of residue position for CatSper1, CatSper2, CatSper3, and CatSper4. Domain boundaries are indicated schematically above each plot: N-terminal domain (NTD, peach), transmembrane domain (TMD, gray), and C-terminal domain (CTD, peach). The red dashed line denotes the pDP threshold of 0.6, above which regions are predicted to have a high propensity for droplet formation.CatSper1 displays an extended N-terminal region with consistently high pDP scores exceeding the LLPS threshold, whereas CatSper2–4 lack comparable LLPS-prone regions in their N-terminal domains and instead show elevated pDP scores primarily within C-terminal regions. These results indicate that a predicted intrinsically droplet-forming region is a unique feature of the CatSper1 N-terminus among the pore-forming CatSper subunits. Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.s004

(TIFF)

S7 Fig. Sperm capacitation leads to removal of temperature sensitivity as shown by calcium imaging.

(Panel A) Average traces of normalized Fluo4-AM fluorescence in wild-type noncapacitated mouse spermatozoa (blue) and from capacitated sperm (red; total capacitation time 80–110 min) in response to increasing temperatures: 25 °C, 30 °C, and 38 °C. (Panel B) Average traces of normalized Fluo4-AM fluorescence in noncapacitated wild-type (blue, same data as above in (Panel A)) and CatSper1−/− (green) mouse spermatozoa in response to increasing temperatures: 25 °C, 30 °C, and 38 °C. Each trace was normalized to its respective ionomycin (Ionomycin, 10 μM) response. Data are presented as mean and standard deviation (n = 3–4 biological replicates). *P < 0.050, **P < 0.010, ****P < 0.00010, ns = non-significant, by unpaired t test. (Panel C) Hierarchical organization of CatSper voltage-gated module located in S4 regions. Sequence alignment between S4 regions of murine CatSper1–4. An unusually high density of positively charged residues, four R and two K, is located in the S4 region of CatSper1, which is predicted to enable enhanced voltage sensitivity. CatSper2–4, however, contain progressively fewer charged residues, with CatSper4 harboring only two R. Conserved positively charged residues are indicated by green triangles. Blue triangles indicate additional residues unique to CatSper1. Underlying data can be found in S1 Data file.

https://doi.org/10.1371/journal.pbio.3003776.s007

(TIFF)

References

  1. 1. Ren D, Navarro B, Perez G, Jackson AC, Hsu S, Shi Q, et al. A sperm ion channel required for sperm motility and male fertility. Nature. 2001;413(6856):603–9. pmid:11595941
  2. 2. Quill TA, Ren D, Clapham DE, Garbers DL. A voltage-gated ion channel expressed specifically in spermatozoa. Proc Natl Acad Sci U S A. 2001;98(22):12527–31. pmid:11675491
  3. 3. Carlson AE, Westenbroek RE, Quill T, Ren D, Clapham DE, Hille B, et al. CatSper1 required for evoked Ca2+ entry and control of flagellar function in sperm. Proc Natl Acad Sci U S A. 2003;100(25):14864–8. pmid:14657352
  4. 4. Lobley A, Pierron V, Reynolds L, Allen L, Michalovich D. Identification of human and mouse CatSper3 and CatSper4 genes: characterisation of a common interaction domain and evidence for expression in testis. Reprod Biol Endocrinol. 2003;1:53. pmid:12932298
  5. 5. Avidan N, Tamary H, Dgany O, Cattan D, Pariente A, Thulliez M, et al. CATSPER2, a human autosomal nonsyndromic male infertility gene. Eur J Hum Genet. 2003;11(7):497–502. pmid:12825070
  6. 6. Kirichok Y, Navarro B, Clapham DE. Whole-cell patch-clamp measurements of spermatozoa reveal an alkaline-activated Ca2+ channel. Nature. 2006;439(7077):737–40. pmid:16467839
  7. 7. Qi H, Moran MM, Navarro B, Chong JA, Krapivinsky G, Krapivinsky L, et al. All four CatSper ion channel proteins are required for male fertility and sperm cell hyperactivated motility. Proc Natl Acad Sci U S A. 2007;104(4):1219–23. pmid:17227845
  8. 8. Suarez SS. Control of hyperactivation in sperm. Hum Reprod Update. 2008;14(6):647–57. pmid:18653675
  9. 9. Navarro B, Kirichok Y, Chung J-J, Clapham DE. Ion channels that control fertility in mammalian spermatozoa. Int J Dev Biol. 2008;52(5–6):607–13. pmid:18649274
  10. 10. Cai X, Clapham DE. Evolutionary genomics reveals lineage-specific gene loss and rapid evolution of a sperm-specific ion channel complex: CatSpers and CatSperbeta. PLoS One. 2008;3(10):e3569. pmid:18974790
  11. 11. Hildebrand MS, Avenarius MR, Smith RJH. CATSPER-related male infertility; 2009.
    • 12. Lishko PV, Kirichok Y. The role of Hv1 and CatSper channels in sperm activation. J Physiol. 2010;588(Pt 23):4667–72. pmid:20679352
    • 13. Chung J-J, Navarro B, Krapivinsky G, Krapivinsky L, Clapham DE. A novel gene required for male fertility and functional CATSPER channel formation in spermatozoa. Nat Commun. 2011;2:153. pmid:21224844
    • 14. Chung J-J, Shim S-H, Everley RA, Gygi SP, Zhuang X, Clapham DE. Structurally distinct Ca(2+) signaling domains of sperm flagella orchestrate tyrosine phosphorylation and motility. Cell. 2014;157(4):808–22. pmid:24813608
    • 15. Cai X, Wang X, Clapham DE. Early evolution of the eukaryotic Ca2+ signaling machinery: conservation of the CatSper channel complex. Mol Biol Evol. 2014;31(10):2735–40. pmid:25063443
    • 16. Quill TA, Sugden SA, Rossi KL, Doolittle LK, Hammer RE, Garbers DL. Hyperactivated sperm motility driven by CatSper2 is required for fertilization. Proc Natl Acad Sci U S A. 2003;100(25):14869–74. pmid:14657366
    • 17. Swain DK, Vergara C, Castro-Arnau J, Lishko PV. The essential calcium channel of sperm CatSper is temperature-gated. Nat Commun. 2025;16(1):3657. pmid:40246904
    • 18. Lishko PV, Botchkina IL, Fedorenko A, Kirichok Y. Acid extrusion from human spermatozoa is mediated by flagellar voltage-gated proton channel. Cell. 2010;140(3):327–37. pmid:20144758
    • 19. Lishko PV, Botchkina IL, Kirichok Y. Progesterone activates the principal Ca2+ channel of human sperm. Nature. 2011;471(7338):387–91. pmid:21412339
    • 20. Hwang JY, Chung JJ. CatSper Calcium Channels: 20 Years On. Physiol (Bethesda). 2023;38(0).
    • 21. Zhao Y, Wang H, Wiesehoefer C, Shah NB, Reetz E, Hwang JY, et al. 3D structure and in situ arrangements of CatSper channel in the sperm flagellum. Nat Commun. 2022;13(1):3439. pmid:35715406
    • 22. Hwang JY, et al. Dual sensing of physiologic pH and calcium by EFCAB9 regulates sperm motility. Cell. 2019;177(e1419):1480–94.
    • 23. Chung J-J, Miki K, Kim D, Shim S-H, Shi HF, Hwang JY, et al. CatSperζ regulates the structural continuity of sperm Ca2+ signaling domains and is required for normal fertility. Elife. 2017;6:e23082. pmid:28226241
    • 24. Brenker C, Goodwin N, Weyand I, Kashikar ND, Naruse M, Krähling M, et al. The CatSper channel: a polymodal chemosensor in human sperm. EMBO J. 2012;31(7):1654–65. pmid:22354039
    • 25. Wang H, Liu J, Cho K-H, Ren D. A novel, single, transmembrane protein CATSPERG is associated with CATSPER1 channel protein. Biol Reprod. 2009;81(3):539–44. pmid:19516020
    • 26. Ho K, Wolff CA, Suarez SS. CatSper-null mutant spermatozoa are unable to ascend beyond the oviductal reservoir. Reprod Fertil Dev. 2009;21(2):345–50. pmid:19210926
    • 27. Liu J, Xia J, Cho K-H, Clapham DE, Ren D. CatSperbeta, a novel transmembrane protein in the CatSper channel complex. J Biol Chem. 2007;282(26):18945–52. pmid:17478420
    • 28. Xia J, Reigada D, Mitchell CH, Ren D. CATSPER channel-mediated Ca2+ entry into mouse sperm triggers a tail-to-head propagation. Biol Reprod. 2007;77(3):551–9. pmid:17554080
    • 29. Jin J, Jin N, Zheng H, Ro S, Tafolla D, Sanders KM, et al. Catsper3 and Catsper4 are essential for sperm hyperactivated motility and male fertility in the mouse. Biol Reprod. 2007;77(1):37–44. pmid:17344468
    • 30. Carlson AE, Quill TA, Westenbroek RE, Schuh SM, Hille B, Babcock DF. Identical phenotypes of CatSper1 and CatSper2 null sperm. J Biol Chem. 2005;280(37):32238–44. pmid:16036917
    • 31. Lin S, Ke M, Zhang Y, Yan Z, Wu J. Structure of a mammalian sperm cation channel complex. Nature. 2021;595(7869):746–50. pmid:34225353
    • 32. Zhao Q, Lin S, Kang H, Ru Y, Xu Q, Yu Z, et al. ARMH2 is a cytosolic component of CatSper crucial for sperm function. Nat Commun. 2025;16(1):10243. pmid:41271765
    • 33. Huang X, Miyata H, Wang H, Mori G, Iida-Norita R, Ikawa M, et al. A CUG-initiated CATSPERθ functions in the CatSper channel assembly and serves as a checkpoint for flagellar trafficking. Proc Natl Acad Sci U S A. 2023;120(39):e2304409120. pmid:37725640
    • 34. Cai X, Clapham DE. Ancestral Ca2 signaling machinery in early animal and fungal evolution. Mol Biol Evol. 2011.
    • 35. Austin CR. Observations on the penetration of the sperm in the mammalian egg. Aust J Sci Res B. 1951;4(4):581–96. pmid:14895481
    • 36. Chang MC. Fertilizing capacity of spermatozoa deposited into the fallopian tubes. Nature. 1951;168(4277):697–8. pmid:14882325
    • 37. Visconti PE, Moore GD, Bailey JL, Leclerc P, Connors SA, Pan D, et al. Capacitation of mouse spermatozoa. II. Protein tyrosine phosphorylation and capacitation are regulated by a cAMP-dependent pathway. Development. 1995;121(4):1139–50. pmid:7538069
    • 38. Visconti PE, Krapf D, de la Vega-Beltrán JL, Acevedo JJ, Darszon A. Ion channels, phosphorylation and mammalian sperm capacitation. Asian J Androl. 2011;13(3):395–405. pmid:21540868
    • 39. Ded L, Hwang JY, Miki K, Shi HF, Chung J-J. 3D in situ imaging of the female reproductive tract reveals molecular signatures of fertilizing spermatozoa in mice. Elife. 2020;9:e62043. pmid:33078708
    • 40. Abramson J, Adler J, Dunger J, Evans R, Green T, Pritzel A, et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature. 2024;630(8016):493–500. pmid:38718835
    • 41. Brown SG, Miller MR, Lishko PV, Lester DH, Publicover SJ, Barratt CLR, et al. Homozygous in-frame deletion in CATSPERE in a man producing spermatozoa with loss of CatSper function and compromised fertilizing capacity. Hum Reprod. 2018;33(10):1812–6. pmid:30239785
    • 42. Cohan MC, Shinn MK, Lalmansingh JM, Pappu RV. Uncovering non-random binary patterns within sequences of intrinsically disordered proteins. J Mol Biol. 2022;434(2):167373. pmid:34863777
    • 43. Ruff KM, King MR, Ying AW, Liu V, Pant A, Lieberman WE, et al. Molecular grammars of predicted intrinsically disordered regions that span the human proteome. Cell. 2026;189(1):323–42.e17. pmid:41232529
    • 44. Durell SR, Ben-Naim A. Temperature dependence of hydrophobic and hydrophilic forces and interactions. J Phys Chem B. 2021;125(48):13137–46. pmid:34850632
    • 45. Hwang JY, Wang H, Oh J-N, Kim J, Finnegan SF, Skakkebaek NE, et al. CATSPERε extracellular domains are essential for sperm calcium channel assembly and activity modulation. Sci Adv. 2026;12(5):eadw3414. pmid:41604475
    • 46. Xu Q, Lin S, Zhao Q, Ru Y, Kang H, Zeng X, et al. Molecular basis of the higher-order assembly of CatSper. Proc Natl Acad Sci U S A. 2026;123(2):e2510754123. pmid:41490491
    • 47. Calinsky R, Levy Y. Aromatic residues in proteins: re-evaluating the geometry and energetics of π-π, Cation-π, and CH-π Interactions. J Phys Chem B. 2024;128(36):8687–700. pmid:39223472
    • 48. Calinsky R, Levy Y. Histidine in proteins: pH-dependent interplay between π-π, Cation-π, and CH-π Interactions. J Chem Theory Comput. 2024;20(15):6930–45. pmid:39037905
    • 49. Calinsky R, Levy Y. A pH-dependent coarse-grained model for disordered proteins: histidine interactions modulate conformational ensembles. J Phys Chem Lett. 2024;15(37):9419–30. pmid:39248414
    • 50. Bashford D, Karplus M. pKa’s of ionizable groups in proteins: atomic detail from a continuum electrostatic model. Biochemistry. 1990;29(44):10219–25. pmid:2271649
    • 51. Alvau A, Battistone MA, Gervasi MG, Navarrete FA, Xu X, Sánchez-Cárdenas C, et al. The tyrosine kinase FER is responsible for the capacitation-associated increase in tyrosine phosphorylation in murine sperm. Development. 2016;143(13):2325–33. pmid:27226326
    • 52. Sumigama S, Mansell S, Miller M, Lishko PV, Cherr GN, Meyers SA, et al. Progesterone accelerates the completion of sperm capacitation and activates CatSper channel in spermatozoa from the rhesus macaque. Biol Reprod. 2015;93(6):130. pmid:26490839
    • 53. Lishko PV, Kirichok Y, Ren D, Navarro B, Chung J-J, Clapham DE. The control of male fertility by spermatozoan ion channels. Annu Rev Physiol. 2012;74:453–75. pmid:22017176
    • 54. Lishko P, Clapham DE, Navarro B, Kirichok Y. Sperm patch-clamp. Methods Enzymol. 2013;525:59–83. pmid:23522465
    • 55. Eggert-Kruse W, Köhler A, Rohr G, Runnebaum B. The pH as an important determinant of sperm-mucus interaction. Fertil Steril. 1993;59(3):617–28. pmid:8458467
    • 56. Maas DH, Stein B, Metzger H. PO2 and pH measurements within the rabbit oviduct following tubal microsurgery: reanastomosis of previously dissected tubes. Adv Exp Med Biol. 1984;169:561–70. pmid:6731111
    • 57. Maas DH, Storey BT, Mastroianni L Jr. Hydrogen ion and carbon dioxide content of the oviductal fluid of the rhesus monkey (Macaca mulatta). Fertil Steril. 1977;28(9):981–5. pmid:19307
    • 58. Niu Q, Lishko PV. Ion channels as gatekeepers of fertility: from uterine Kir7.1 to sperm CatSper. Physiology (Bethesda). 2026:10.1152/physiol.00015.2026. pmid:42262737
    • 59. Dodge DP, Maccrimmon HR. Environmental influences on extended spawning of rainbow trout (Salmo gairdneri). Trans Am Fish Soc. 1971;100(2):312–8.
    • 60. Peterson RH, Spinney HCE, Sreedharan A. Development of Atlantic salmon (Salmo salar) eggs and alevins under varied temperature regimes. J Fish Res Bd Can. 1977;34(1):31–43.
    • 61. Southward EC, Andersen AC, Hourdez S. Lamellibrachia anaximandrin. sp., a new vestimentiferan tubeworm (Annelida) from the Mediterranean, with notes on frenulate tubeworms from the same habitat. Zoosystema. 2011;33(3):245–79.
    • 62. Clutton EA, Alurralde G, Repolho T. Early developmental stages of native populations of Ciona intestinalis under increased temperature are affected by local habitat history. J Exp Biol. 2021;224(Pt 5):jeb233403. pmid:33472872
    • 63. Hammond LM, Hofmann GE. Thermal tolerance of Strongylocentrotus purpuratus early life history stages: mortality, stress-induced gene expression and biogeographic patterns. Mar Biol. 2010;157(12):2677–87. pmid:24391252
    • 64. Yuan X, Yang H, Wang L, Zhou Y, Gabr HR. Bioenergetic responses of sub-adult sea cucumber Apostichopus japonicus (Selenka) (Echinodermata: Holothuroidea) to temperature with special discussion regarding its southernmost distribution limit in China. J Thermal Biol. 2009;34(6):315–9.
    • 65. Stefanik DJ, Friedman LE, Finnerty JR. Collecting, rearing, spawning and inducing regeneration of the starlet sea anemone, Nematostella vectensis. Nat Protoc. 2013;8:916–23.
    • 66. Gómez-Gras D, Bensoussan N, Ledoux JB, López-Sendino P, Cerrano C, Ferretti E, et al. Exploring the response of a key Mediterranean gorgonian to heat stress across biological and spatial scales. Sci Rep. 2022;12(1):21064. pmid:36473926
    • 67. Van Dooremalen C, Koekkoek J, Ellers J. Temperature-induced plasticity in membrane and storage lipid composition: thermal reaction norms across five different temperatures. J Insect Physiol. 2011;57(2):285–91. pmid:21115015
    • 68. Putnam HM, Edmunds PJ, Fan T-Y. Effect of temperature on the settlement choice and photophysiology of larvae from the reef coral Stylophora pistillata. Biol Bull. 2008;215(2):135–42. pmid:18840774
    • 69. Emig CC. Ecology of inarticulated brachiopods. In: Treatise on Invertebrate Paleontology; 1997. pp. 473–526.
      • 70. Manullang C, Hanahara N, Tarigan AI, Abe Y, Furukawa M, Morita M. Slight thermal stress exerts genetic diversity selection at coral (Acropora digitifera) larval stages. BMC Genomics. 2025;26(1):36. pmid:39810102
      • 71. dela Cruz DW, Harrison PL. Optimising conditions for in vitro fertilization success of Acropora tenuis, A. millepora and Favites colemani corals in northwestern Philippines. J Exp Mar Biol Ecol. 2020;524:151286.
      • 72. Service, U.S.F.W. Chinese softshell turtle (Pelodiscus sinensis) ecological risk screening summary (ed. Web Version); 2018.
        • 73. Ferguson MW, Joanen T. Temperature of egg incubation determines sex in Alligator mississippiensis. Nature. 1982;296(5860):850–3. pmid:7070524
        • 74. Riccio AP, Goldman BD. Circadian rhythms of body temperature and metabolic rate in naked mole-rats. Physiol Behav. 2000;71(1–2):15–22. pmid:11134680
        • 75. Smyth DM. Temperature regulation in the platypus, Ornithorhynchus anatinus (Shaw). Comp Biochem Physiol A Comp Physiol. 1973;45(3):705–15. pmid:4147924
        • 76. Spears S, Pettit C, Berkowitz S, Collier S, Colwell C, Livingston EH, et al. Lizards in the wind: the impact of wind on the thermoregulation of the common wall lizard. J Therm Biol. 2024;121:103855. pmid:38648702
        • 77. Descovich KA, Johnston S, Lisle A, Nicolson V, Janssen T, Brooks P, et al. Long-term measurement of body temperature in the southern hairy-nosed wombat (Lasiorhinus latifrons). Aust Mammal. 2016;39(1):48–55.
        • 78. Busse S, Lutter D, Heldmaier G, Jastroch M, Meyer CW. Torpor at high ambient temperature in a neotropical didelphid, the grey short-tailed opossum (Monodelphis domestica). Naturwissenschaften. 2014;101(11):1003–6. pmid:25142634
        • 79. Harlow HJ, Purwandana D, Jessop TS, Phillips JA. Body temperature and thermoregulation of Komodo dragons in the field. J Thermal Biol. 2010;35(7):338–47.
        • 80. Boyles JG, Smit B, Sole CL, McKechnie AE. Body temperature patterns in two syntopic elephant shrew species during winter. Comp Biochem Physiol A Mol Integr Physiol. 2012;161:89–94.
        • 81. Jones ME, Grigg GC, Beard LA. Body temperatures and activity patterns of Tasmanian devils (Sarcophilus harrisii) and eastern quolls (Dasyurus viverrinus) through a subalpine winter. Physiol Zool. 1997;70(1):53–60. pmid:9231376
        • 82. Streicher S, Boyles JG, Oosthuizen MK, Bennett NC. Body temperature patterns and rhythmicity in free-ranging subterranean Damaraland mole-rats, Fukomys damarensis. PLoS One. 2011;6(10):e26346. pmid:22028861
        • 83. Adam D, Johnston SD, Beard L, Nicolson V, Gaughan JB, Lisle AT, et al. Body temperature of free-ranging koalas (Phascolarctos cinereus) in south-east Queensland. Int J Biometeorol. 2020;64(8):1305–18. pmid:32388686
        • 84. Liu R, Z X. Artificially rearing the Baiji, Lipotes vexillifer, for twenty-one years. China Acad J Electron Publish House. 2018.
        • 85. Mole MA, Rodrigues DÁraujo S, van Aarde RJ, Mitchell D, Fuller A. Savanna elephants maintain homeothermy under African heat. J Comp Physiol B. 2018;188(5):889–97. pmid:30008137
        • 86. Snyder GK, Baudinette RV, Gannon BJ. Oxygen transport and acid-base balance during exercise in the tammar wallaby. Respir Physiol. 1999;117(1):41–51. pmid:10505478
        • 87. Melero M, Rodríguez-Prieto V, Rubio-García A, García-Párraga D, Sánchez-Vizcaíno JM. Thermal reference points as an index for monitoring body temperature in marine mammals. BMC Res Notes. 2015;8:411. pmid:26338544
        • 88. Baak NA, Cantineau AE, Farquhar C, Brison DR. Temperature of embryo culture for assisted reproduction. Cochrane Database Syst Rev. 2019;9(9):CD012192. pmid:31529804
        • 89. Walters EA, Brown JL, Krisher R, Voelkel S, Swain JE. Impact of a controlled culture temperature gradient on mouse embryo development and morphokinetics. Reprod Biomed Online. 2020;40(4):494–9. pmid:32171708
        • 90. Taffe MA. A comparison of intraperitoneal and subcutaneous temperature in freely moving rhesus macaques. Physiol Behav. 2011;103(5):440–4. pmid:21443893
        • 91. McGinnis SM, Southworth TP. Body temperature fluctuations in the northern elephant seal. J Mammal. 1967;48(3):484.
        • 92. Silva RG. Assessment of body surface temperature in cetaceans: an iterative approach. Braz J Biol. 2004;64(3B):719–24. pmid:15620012
        • 93. Holyoak GW, Stones RC. Temperature regulation of the little brown bat, Myotis lucifugus after acclimation at various ambient temperatures. Comp Biochem Physiol A Comp Physiol. 1971;39(3):413–20. pmid:4150224
        • 94. D R. Temperature regulation and thermal environment. Dukes’ physiology of domestic animals. Cornell University Press; 2004.
          • 95. Shah J, et al. The brown bear and hibernating mammals as a translational model for human resilience: insights for space medicine, critical care, and austere environments. Biology. 2025;14:1434.
          • 96. North American River Otter. (ed. Institute, S.s.N.Z.C.B.) Smithsonian's National Zoo&Conservation Biology Institute; 2025.
            • 97. LafeberVet. The resource for exotic animal veterinary professionals; 2010.
              • 98. Liles HL, Schneider LG, Pohler KG, Oliveira Filho RV, Neal Schrick F, Payton RR, et al. Positive relationship of rectal temperature at fixed timed artificial insemination on pregnancy outcomes in beef cattle. J Anim Sci. 2022;100(7):skac100. pmid:35772756
              • 99. Jass CN, Mead JI. Capricornis crispus. Mammalian Species. 2004;750:1–10.
              • 100. Namiki M, Fukayama T, Suzuki T, Masaiwa A. Relationship between ear temperature, behaviour and stress hormones in guinea pigs (Cavia porcellus) during different interactive activities in zoos. Animals (Basel). 2024;14(7):1111. pmid:38612351
              • 101. MacArthur RA, Wang LC. Physiology of thermoregulation in the pika, Ochotona princeps. Can J Zool. 1973;51(1):11–6. pmid:4696076
              • 102. Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004;5:113. pmid:15318951
              • 103. Wilkins MR, Gasteiger E, Bairoch A, Sanchez JC, Williams KL, Appel RD, et al. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol. 1999;112:531–52. pmid:10027275
              • 104. Hardenberg M, Horvath A, Ambrus V, Fuxreiter M, Vendruscolo M. Widespread occurrence of the droplet state of proteins in the human proteome. Proc Natl Acad Sci U S A. 2020;117(52):33254–62. pmid:33318217
              • 105. Holehouse AS, Das RK, Ahad JN, Richardson MOG, Pappu RV. CIDER: resources to analyze sequence-ensemble relationships of intrinsically disordered proteins. Biophys J. 2017;112(1):16–21. pmid:28076807
              • 106. UniProt C. UniProt: the universal protein knowledgebase in 2025. Nucleic Acids Res. 2025;53:D609–17.
              • 107. Potenza E, Di Domenico T, Walsh I, Tosatto SCE. MobiDB 2.0: an improved database of intrinsically disordered and mobile proteins. Nucleic Acids Res. 2015;43(Database issue):D315–20. pmid:25361972
              • 108. Goddard TD, Huang CC, Meng EC, Pettersen EF, Couch GS, Morris JH, et al. UCSF ChimeraX: meeting modern challenges in visualization and analysis. Protein Sci. 2018;27(1):14–25. pmid:28710774
              • 109. Cock PJA, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009;25(11):1422–3. pmid:19304878
              • 110. Tien MZ, Meyer AG, Sydykova DK, Spielman SJ, Wilke CO. Maximum allowed solvent accessibilites of residues in proteins. PLoS One. 2013;8(11):e80635. pmid:24278298
              • 111. Abril-Pla O, Andreani V, Carroll C, Dong L, Fonnesbeck CJ, Kochurov M, et al. PyMC: a modern, and comprehensive probabilistic programming framework in Python. PeerJ Comput Sci. 2023;9:e1516. pmid:37705656
              Read Entire Article

                       

                      

              Start the new Vibrations with a Medbed Franchise today!  

              Protect your whole family with Quantum Orgo-Life® devices

                Advertising by Adpathway