Mechanisms of permselectivity of connexin hemichannels to small molecules

preprint OA: gold CC-BY-4.0
📄 Open PDF Full text JSON View at publisher
Full text 73,792 characters · extracted from oa-pdf · 9 sections · click to expand

Keywords

Connexins, hemichannels, ATP, glutamate, lactate, N-terminus, permeability .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint

Abstract

Connexins can either that act as hemichannels, to facilitate ion and small molecule movement from the cytosol to the extracellular space or as gap junction channels to provide a pathway for solute exchange between cells. Connexins are ubiquitously expressed throughout the body and are implicated in a wide range of physiological and pathological processes. The permselectivity of connexin hemichannels for small neurochemicals remains poorly understood . By coexpressing genetically encoded fluorescent sensors for ATP, glutamate and lactate with a range of connexins, we examined the ability of different hemichannels to permit release of these compounds under physiological conditions and in response to physiological stimuli ( small changes in PCO 2 and transmembrane depolarisation). We found that some connexins were relatively non-selective (Cx26, Cx32, Cx43, Cx31.1) allowing passage of ATP, glutamate and lactate. By contrast other connexins (Cx36, Cx46 and Cx50) were highly selective. Cx36 and Cx46 allowed release of ATP, but not glutamate or lactate. This shows that size of t he permeating molecule is not the sole determinant of permselectivity. By contrast, Cx50 permitted the release of lactate and glutamate but not ATP. We also found that the nature of the opening stimulus could alter the permselectivity of the hemichannel -for some of the relatively non-selective connexins, hemichannel opening via depolarisation was ineffective at allowing release of lactate. By comparing the differential selectivity of the closely related Cx46 and Cx50, we showed that the charge on the N -terminus and N -terminus-TM2 interactions are key contributors to permselectivity. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint

Introduction

There are 21 connexin genes in the human genome (Lucaciu et al., 2023) . Connexins form hexamers that, if unopposed , can act as a plasma membrane hemichannel that opens to the extracellular space. However, hemichannels of closely apposed cells can also dock together to form gap junction channels to provide an aqueous passageway between cells. The structure of the hemichannel is highly conserved in all 21 isoforms; each connexin has an N -terminal helix, 4 transmembrane helices, a cytoplasmic loop, 2 extracellular loops and a cytoplasmic C- terminus (Maeda et al., 2009; Myers et al., 2018; Flores et al., 2020; Brotherton et al., 2022; Lee et al., 2023; Qi et al., 2023; Brotherton et al., 2024) . 6 subunits then co -assemble to form a hemichannel with a central pore, spanning ~1.2 nm that is permeable to small molecules up to a molecular weight of about 1000. Major differences in structure between isoforms lie within the cytoplasmic loop and C -terminus, which vary in sequence and length (Mese et al., 2007) . The 6 N -terminal helices line the hemichannel pore, to form the narrowest part of the permeation pathway, suggesting that the N-terminus may be an important for determining the permeability of the channel (Oshima et al., 2007, 2008; Maeda et al., 2009; Nielsen et al., 2019; Yue et al., 2021). Connexin hemichannels have been documented to release small molecules such as ATP under physiological conditions (Weissman et al., 2004; Pearson et al., 2005; Huckstepp et al., 2010b; Chever et al., 2014; van de Wiel et al., 2020) . Yet the mechanisms that control hemichannel permeability to different molecules , and whether there is specificity to which molecules may permeate is still unclear. Traditionally, this has been investigated using various fluorescent dyes such as ethidium bromide, and the size and charge of dyes provided some evidence for selectivity to release of larger molecules (Li et al., 1996; Saez et al., 2010; Johnson et al., 2016). Investigation of hemichannel permeability via dye fluxes, while valuable , may differ from how physiological metabolites such as ATP, glutamate or lactate permeate these channels. Traditionally, connexin permeability studies have used the removal of extracellular divalent cations to unblock the hemichannels (Hansen et al., 2014; Nielsen et al., 2019). The mechanism was defined in Cx26 and Cx32 and involves a ring of 12 aspartate residues within the extracellular loop that provide a carboxylate cluster able to bind Ca2+ ions with millimolar affinity (Gomez-Hernandez et al., 2003; Bennett et al., 2016; Lopez et al., 2016) . While hemichannels are essentially blocked at Ca2+ concentrations over 1 mM, there are very few if any physiological conditions in which extracellular Ca2+ is lower than 1 mM. Thus, unblocking of hemichannels via Ca2+ removal may open a permeation pathway that is not representative of physiological gating of connexin hemichannels. Nevertheless, this method has been used to suggest differential permeability of Cx30 and Cx43 to a variety of small molecules (Hansen et al., 2014; Nielsen et al., 2019). Connexin hemichannels can be opened under physiological conditions by other gating stimuli. As the N-termini of connexins have charged residues and are within the membrane electric field, almost all c onnexins can be opened by sufficient depolarisation, without the need to lower extracellular Ca 2+ (Pinto et al., 2016) . A subset of connexins is directly sensitive to the concentration of gaseous CO2 and can be opened by relatively small changes in the partial pressure of CO2 (PCO2) around the physiological norm. CO2-dependent gating was elucidated in Cx26 (Meigh et al., 2013): K125 is carbamylated by CO2, which facilitates the formation of a salt bridge with R104 of the adjacent subunit to bias the channel into an open conformation. This .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint mechanism involves movements of the N -terminus, along with related movements of the transmembrane helices (particularly TM2) (Brotherton et al., 2022; Brotherton et al., 2024). CO2 gating has subsequently been discovered in a subset of connexins: Cx30, Cx32 (Huckstepp et al., 2010a), Cx43 (Nijjar et al., 2025) and Cx50, since they all have the carbamylation motif that is required for CO2-mediated opening. Their CO2 sensitivity lies close to the physiological range of PCO2; Cx26, Cx43 and Cx50 are maximally open at ~55 mmHg PCO2, and Cx32 at ~70 mmHg. To permit the study of the permeation of connexin hemichannels, for physiological molecules under their normal electrochemical gradients and with physiological gating stimuli (voltage and PCO2), we have developed an assay to allow real-time imaging of analyte release at single cell resolution. We utilised the genetically encoded sensors GRABATP (Wu et al., 2022) , iGluSnFR (Marvin et al., 2013) and eLACCO1.1 (Nasu et al., 2021) to measure the release of ATP, glutamate and lactate respectively via coexpressed connexin hemichannels . We find that connexins can be divided essentially into relatively non-selective and highly selective categories. Connexins such as Cx26, Cx32 and Cx43 fall into the relatively non-selective category, whereas Cx36, Cx46 and Cx50 are highly selective. By comparing the highly homologous Cx46 and Cx50, we have shown that key residues in the N -terminus and in the interacting portion of TM2 determine the permeability profile of the hemichannel.

Results

We first transfected HeLa cells with the genetically encoded sensors on their own to ensure that the sensors had no intrinsic responses to CO2 or high KCl solutions, or that parental HeLa cells exhibited ATP, glutamate or lactate release in the absence of connexin expression. The median change in normalised fluorescence (∆F/F0) for GRABATP was -0.007 (95% CI: 0.0043, -0.013) and -0.004 (0.0017, -0.0064) for 55 mmHg and 50 mM KCl respectively. GRABATP was functional as it gave a median response of 0.2 (0.25, 0.16) to 3 µM ATP (Fig.1). The median change in ∆F/F 0 for iGluSnFR with 55 mmHg was 0 (0.0008, -0.0026), for 50 mM KCl it was -0.014 (-0.011, -0.030), and for 3 µM glutamate was 0.1 (0.13, 0.05). Finally, the median ∆F/F0 for eLACCO1.1 (modified by insertion into iGluSnFR backbone) was 0.005 (0.0073, -0.0046), -0.002 (0.0048, -0.0069) and 0.04 (0.046, 0.032) for 55 mmHg, 50 mM KCl and 3 µM lactate, respectively. The n egative recorded values are an artefact of photobleaching. As there were no responses of any of the genetically sensors to a change in PCO2 or membrane depolarisation, we conclude that parental HeLa cells do not express any channels capable of releasing ATP, glutamate or lactate to these stimuli. Cx26, Cx32,Cx43 and Cx31.3 are permeable to ATP, glutamate, lactate We next co-transfected HeLa cells with Cx26 , Cx32, Cx43 or Cx31.3 and one of the genetically encoded fluorescent sensors. Throughout all of the assays, we selected cells that co-expressed the connexin and the sensor for measurement and analysis (Fig. 2 figure supplement 1) . Cx26 hemichannels are highly permeable to ATP, glutamate and lactate (Fig. 2). The median ATP release to hypercapnic stimuli was 1.5 µM (1.91, 1.32) and with a depolarising stimulus was 2.5 µM (3.67, 2.30). This was significantly less than for glutamate and lactate: the median glutamate release was 5.2 µM (5.71, 4.36) and 1.1 µM (6.10, 0.74) when opened with hypercapnia and a depolarising stimulus, respectively. Lactate release evoked by hypercapnia was comparable to glutamate: median release of 4.7 µM (6.26, 2.52). However when stimulated by depolarisation .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint no lactate release was evident : median -0.1 µM (0.23, -1.04). For voltage -gated glutamate release it appeared that in some cells the depolarising stimulus was ineffective at evoking release. This suggests that KCl evoked depolarisation (estimated to be about 70 mV from the Nernst equation) is a less reliable gating mechanism than CO2 for Cx26. We determined that Cx32 was also permeable to all three tested analytes . Significantly more glutamate and lactate was released compared to ATP regardless of the stimulus (Fig. 3). However the nature of the stimulus did alter the relative amounts of release of the three metabolites. The median release of ATP was 1. 3 (1.51, 1.16) and 3.0 (3.12, 2.86) µM whe n stimulated by hyper capnia and voltage, respectively. Glutamate releas ed via Cx32 by hypercapnia and depolarisation was respectively 3.9 (4.64, 3.41) µM and 5.2 (7.44, 4.62) µM. Lactate release via Cx32 by hypercapnia and depolarisation was respectively 7.2 (8.32, 5.31) µM and 5.5 (6.84, 4.52) µM. With hypercapnia, the relative release of both glutamate and lactate compared to ATP was considerably more than might be expected from the calculated electrochemical driving force (Tables 1 and 2) . This suggests that Cx32 may have enhanced permeability for these molecules when opened by hypercapnia. With a hypercapnic stimulus Cx43 was permeable to ATP, glutamate an d lactate (Fig. 4). The median ATP, glutamate and lactate release were repsectively 2.2 (2.34, 1.49) µM, 5.3 (5.71, 4.36) µM and 4. 5 (5.84, 2.46) µM (Fig. 4). This amount of release is roughly proportional to the electrochemical driving force on these molecules ( Tables 1 and 2 ) suggesting no selectivity. However during membrane depolarisation evoked by 50 mM KCl, the permeation profile was different. Whereas median ATP release was 2.5 (2.88, 2.22) µM (similar to that evoked by CO2), the release of glutamate was reduced, median 1.0 (6.10, 0.74) µM, and lactate was not released at all during this stumulus, median -1.0 (-0.069, -1.49) µM. The gating mechanism of Cx43 thus alters the permeability profile of the hemichannel to small molecules. While permeability of Cx43 hemichannels to ATP has been elegantly demonstrated (Kang et al., 2008), previous reports suggest that Cx43 hemichannels are apparently not permeable to glutamate or lactate (Hansen et al., 2014; Nielsen et al., 2019). However these previous studies used removal of extracellular Ca 2+ to unblock the channel, and the channel might thus have a different permeability profile. Whereas zero [Ca2+]ext is often achieved by the use of chelators such as EGTA, the genetically encoded fluorescent sensors have some degree of Ca2+- dependency. To ensure compatibility with the sensors, we omitted EGTA and Ca 2+ to lower but not eliminate extracellular Ca2+. To ensure this was still sufficient to open the hemichannels, we employed a dye loading assay using the hemichannel-permeable dye , FITC. Under control conditions (PCO2 20 mmHg, 2 mM [Ca2+]ext) the median pixel intensity was 15.7 (24.72, 14.04). In low [Ca2+]ext (PCO2 20 mmHg), the median pixel intensity was 54.3 (86,74, 51.82) (Fig. 4, figure supplement 1). For a comparison to a reliable opening stimulus, we use a PCO2 of 55 mmHg to open the hemichannels and permit loading with FITC . This yielded a median pixel intensity of 48.0 (54.83, 40.66) and demonstrates that the low [Ca2+]ext solution was an effective stimulus to open Cx43 hemichannels. Having established the efficacy of low [Ca2+]ext at opening Cx43 hemichannels, we assessed the permeation of ATP, glutamate and lactate during this stimulus and compared it to permeation in response to hypercapnia in the same cells. Consistent with our previous findings, hypercapnia evoked the release of ATP, glutamate and lactate. However, low [Ca 2+]ext .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint significantly reduced the analyte release through Cx43 hemichannels, with a median release of 0.1 (0.86, 0.08) µM ATP, 0.2 (0.76, 0.004) µM glutamate and 0. 5 (2.11, -0.25) µM lactate (Fig. 4, figure supplement 2 ). The nature of the gating mechanism appears to change channel permeability and may account for why our results differ from previous release studies. Cx31.3 (also called Cx29) is mainly expressed in myelinating cells (Cisterna et al., 2019). There has been some suggestion that Cx31.3 predominantly forms hemichannels and mediates ATP efflux from cells that expressing this isoform (Sargiannidou et al., 2008; Liang et al., 2011) . We were able to demonstrate a median release of 2.6 µM (95% CI 2.91, 2.26) ATP, 2.3 µM (2.58, 1.91) glutamate and 2.8 µM lactate (95% CI 2.92, 2.24) (Fig. 5). While relatively non-selective, Cx31.3 nevertheless shows some preference for ATP as the relative permeation of glutamate and lactate is less than that predicted by the electrochemical driving force on these molecules (Tables 1 and 2). Cx36, Cx46 and Cx50 hemichannels have highly specific permeability profiles Cx36 acts predominantly as a major neuronal connexin (Condorelli et al., 1998) and forms the gap junctions or electrical synapses that facilitate fast synaptic transmission and synchronous neuronal firing (Srinivas et al., 1999; Deans et al., 2001; Buhl et al., 2003) . We assayed the release of ATP, glutamate and lactate from HeLa cells transfected with Cx36. Because Cx36 is insensitive to CO 2 (Huckstepp et al., 2010a) , we used the depolarising stimulus to gate the channel. Surprisingly, we found that of the three analytes, Cx36 was only permeable to ATP (Fig. 6). The median release of ATP evoked by 50 mM KCl was 2.6 (3.27, 2.41) µM, compared to -0.3 (0.22, -0.68) and -0.2 ( -0.05, -0.21) µM glutamate and lactate respectively. As t he electrochemical driving force for release of glutamate and lactate is about double that of ATP (Tables 1 and 2), and these molecules are smaller than ATP, the differential permeability of Cx36 suggests the existence of a selectivity filter within the pore. Human Cx46 also lacks the carbamylation motif and is not sensitive to CO2. We therefore used the high K + stimulus to open Cx46 hemichannels. Like Cx36, Cx46 was only permeable to ATP giving a median release of 2.6 µM (2.74, 2.34) (Fig. 7). No glutamate or lactate was released via Cx46 (median release -0.07 µM (0.11, -0.28) and 0.2 µM (0.41, -0.28) for glutamate and lactate respectively, Fig. 7). By contrast, Cx50 was readily permeable to glutamate (median release 1.5 µM (1.57, 1.32) ) when stimulated by hypercapnia and lactate when stimulated by either hypercapnia or depolarisation (median release 5.8 µM (7.36, 2.55) and 3. 4 µM (3.78, -0.55) respectively, Fig 7). However, no ATP could permeate Cx50 (Fig. 7). Mutational analysis of the differential permeability of Cx46 and Cx50 Cx46 and Cx5 0 are structurally quite similar yet have markedly different permeability profiles (Myers et al., 2018; Flores et al., 2020; Yue et al., 2021). These two connexins therefore offer an opportunity to explore the mechanistic basis of differential permeability. The N-termini fold into the gap junction pore to form the narrowest point, with hydrophobic residues anchoring the helix to transmembrane regions 1 and 2 (TM1/2) to stabilise the open state (Myers et al., 2018). Aligning the N -termini sequences of Cx46 and Cx50 shows that they have a difference in the overall net charge (respectively 0 and -2 for Cx46 and Cx50, Fig. 8). The first divergence is at position 9, where Cx46 has a positively charged arginine and Cx50 has an asparagine. At position .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint 13, Cx50 has a negatively charged glutamate residue where Cx46 has a neutrally charged asparagine. To explore the possible roles of the difference in net charge, we introduced the single mutations N9R, E13N and N9 K into Cx50 to make it more like Cx46. These mutations gave a gain of ATP permeability to Cx50, though this was still significantly less than that of Cx46 WT (Fig. 8). The double mutation, Cx50 N9R E13N , gave an increase in ATP release that matched that of Cx46 WT release (for a depolarising stimulus) , with a median release of 2.4 (3.37, 2.07) µM from the mutant connexin compared to 2.5 (2.74, 2.34) for Cx46WT. We have therefore shown that the net charge of the N -terminus, and the specific residues N9 and E13, act to regulate the ATP permeability of Cx50. We next examined whether, if we made the N-terminus of Cx46 more like that of Cx50, we could change the permeability profile of Cx46 to match that of Cx50. We therefore made the double mutation R9N ,N13E in Cx46. However, this d id not diminish ATP permeation via Cx46 hemichannels (Fig. 9). The median release from Cx46R9N,N13E was 2.8 (3.10,2.47), compared with 2.5 (2.74, 2.34) µM for Cx46WT. This indicates that the regulation of permeability of Cx46 to ATP is more complex than just the net charge of the N-terminus. We considered the possibility that interactions between the N -terminus and TM2 may differ between Cx46 and Cx50 and this could permit ATP permeation even if the net charge of t he N- termini were negative. Both Cx46 and Cx50 have hydrophobic residues at position 14: Cx50 has a valine, but Cx46 has an alanine. We therefore used experimentally determined structures of these connexins (Jaradat et al., 2022) to identify possible interacting residues in TM2 . This highlighted residue 89 as potentially important: in Cx50 this is serine, but in Cx46 it is threonine. We sought to make these interactions in Cx46 more similar to those occurring in Cx50 by introducing the double mutation A14V,T89S . We note that simple introduction of the Cx50 N - terminal helix into Cx46 has been reported as resulting in non-functional gap junction channels due to a steric clash between V14 and T89 in the chimaeric channel (Yue et al., 2021). By making two mutations in Cx46, A14V and T89S, we have avoided this clash. The double mutation A14V, T89S did not by itself alter ATP permeation (Fig. 9), indicating that the Cx46A14V,T89S hemichannel gated normally to voltage. However, when A14V and T89S were then combined with R9N and N13E, the quadruply mutated Cx46R9N,N13E,A14V,T89S hemichannels (Cx46QM) exhibited significantly reduced ATP permeation compared to the wild-type, with a median release of 0.54 (1.53, 0.34) µM (Fig. 9). One possible interpretation of this result is that the quadruple mutation simply reduces overall hemichannel permeability rather than having a specific effect on that of ATP . To evaluate this , we examined depolarisation dependent loading of FITC into HeLa cells expressing either Cx46WT or Cx46QM (Fig. 9, figure supplement 1). We found that FITC still permeated into cells expressing Cx46QMand was not significantly different from the permeation observed in those expressing Cx46WT. This suggests that the quadruple mutant has a somewhat selective effect on ATP permeation. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint

Discussion

This study explored whether there is selectivity to the release of small molecules from connexin hemichannels expressed in HeLa cells . To assess their relative permeability via the different connexins, we need to understand the electrochemical driving force on the three metabolites . As these metabolites are charged , we can use the Nernst equation , along with typical concentrations of these molecules in HeLa cells (Piva and McEvoy-Bowe, 1998; Imamura et al., 2009; San Martin et al., 2013) to calculate the equilibrium potential (Table 1). If we assume that the concentrations of these metabolites are scattered around a mean value that is consistent across all HeLa cells , sufficient recordings of release should statistically reflect these transmembrane concentrations. Based on this analysis, if release through connexins were to be non-selective, the relative proportions of release should corre spond to the relative electrochemical driving forces. Thus, for a completely non-selective connexin we would expect to see ATP, glutamate and lactate released approximately in the proportion 1:1.8:1.4 (assuming a resting potential of -60 mV) to hypercapnia and in the proportion 1:2.1:1.6, for the depolarising stimulus. Any notable deviation from these proportions would indicate that there is some selectivity to release of small molecules through connexin hemichannels (Table 2). Our study is the first to produce a comprehensive permeability profile of a wide range of connexin isoforms for release of physiological metabolites with physiological opening stimuli . Connexin hemichannels fall into two b road categories : relatively non-selective (Cx26, Cx32, Cx43, Cx31.3); and highly selective (Cx36, Cx46 and Cx50). For Cx26 and Cx32 , the release of glutamate and lactate relative to ATP is more than predicted by the electrochemical driving force, suggesting that the smaller molecules may permeate more readily than ATP (Table 2) . Cx32 in particular shows enhanced permeability to lactate (nearly 4 times that predicted by driving force). Interestingly, when Cx32 was opened by depolarisation , the relative release of analytes followed that predicted by the electrochemical driving force more closely. For Cx43, when opened by hypercapnia, the relative release of ATP, glutamate and lactate follows very closely the pattern predicted by the electrochemical driving force (Table 2). Cx31.3 shows some preference for ATP over glutamate and lactate but is nevertheless permeable to these smaller molecules (Table 2). Rather surprisingly, Cx50 hemichannels are impermeable to ATP . To our knowledge this is the only connexin hemichannel t hat ATP cannot permeate. Traditionally, connexin permeability studies have studied selectivity by increasing the size of fluorescent dyes (Weber et al., 2004; Harris, 2007), with the idea being that any molecule that was below the limiting pore diameter (around 12 Å) should permeate. This would suggest minimal selectivity between permeants, leaving the major driving force as intracellular concentration. Our resul ts modify this idea . Although the smaller molecules glutamate and lactate do permeate the relatively non-selective channels more easily, lactate being the smallest molecule should permeate all connexins : but it cannot permeate Cx36 or Cx46, whereas ATP, a much larger molecule can. We also find that the permeability profile of the hemichannel can alter with the nature of the gating stimulus. For those connexins that are directly CO 2 sensitive, opening the hemichannel by hypercapnia seems to give greater release than depolarisation (summarized in Table 2). One possibility is that 50 mM KCl , while predicted to depolarise the cell by about 70 mV , may not sufficiently depolarise the membrane to obtain full channel opening. We also find that lowered .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint [Ca2+]ext seems to be the least effective stimulus for release. While this manipulation may unblock the channel (by removing the ring of bound Ca 2+ ions), the N-termini may still partially block the channel and alter the permeation pathway in a way that may not happen with a more physiological stimulus. The selective permeability of connexins may match their physiological roles That there are 21 connexin genes in the human genome, indicates distinct fundamental roles in physiological processes. The need for so many isoforms suggest functional specialisation - combinations of properties that match connexin to function. Sensitivity to gating stimuli, and permeability to small molecules are two important properties that may determine which functional roles connexins are suited to. The relatively non -selective connexins, Cx26, Cx32 and Cx43, have different CO 2 sensitivity profiles. Cx26 is suited to detection of systemic PCO 2 and has a role in the control of breathing (Huckstepp et al., 2010b; van de Wiel et al., 2020; Dale, 2021). Cx32 requires much higher levels of PCO2 to open and may be more suited to detecting local CO 2 production (Huckstepp et al., 2010a; Dospinescu et al., 2019; Butler and Dale, 2023) . It is interesting the Cx32 is highly permeable to lactate when opened by hypercapnia. This lends support to an attractive hypothesis that Cx32 may detect hotspots of metabolic activity and, by opening and permitting release of lactate, could provide metabolic support for highly active cells (Barros, 2013). Cx43 with its extensive C-terminus interacts with many other proteins (Iacobas et al., 2003; Iacobas et al., 2007) , yet is also CO 2 sensitive (Dospinsecu et al., 2025) and is partially open under physiological conditions (Chever et al., 2014; Turovsky et al., 2020). We also discovered a group of connexins with highly selective permeability profiles: Cx36, Cx46 and Cx50. Cx36 is expressed predominantly in neurons and is impermeable to glutamate and lactate. The lack of glutamate permeability may be functionally significant as most excitatory neurons use this as their main neurotransmitter. This lack of permeability may ensure that glutamate release is tightly regulated under ph ysiological conditions and occurs mainly via vesicular exocytosis at synaptic sites. As lactate is an effective metabolite for oxidative phosphorylation in neurons , the lack of permeability of Cx36 hemichannels to lactate would prevent unregulated efflux of lactate from neurons. Cx46 and Cx50 are expressed almost exclusively in the lens of the eye (Mathias et al., 2010; Berthoud and Ngezahayo, 2017). Cx46 and Cx50 form gap junctions between the lens epithelial and fibre cells, and between lens fibre cells. Cx50 is also present as hemichannels in lens fibre cells. As lens fibre cells mature , they lose their intracellular organelles including mitochondria (Bassnett, 2002) and thus the ability to make ATP via oxidative phosphorylation. One can speculate that having ATP permeable C x46 gap junction channels may be valuable in allowing diffusion of this key metabolite from the metabolically active cells in the lens into the relatively inactive lens fibre cells. Equally, the lack of ATP permeability of Cx50 hemichannels may permit preservation of this scarce resource in the lens fibre cells. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint The N-terminus is a fundamental determinant of selectivity Considerable evid ence suggest that t he N -terminus is part of the gating mechanism of connexins. Depending on the isoform , the N -terminus projects into the pore to form a n occluding plug or forms a cap at the cytoplasmic vestibule also to close it. In the open state the N-termini line the pore and thus experience the membrane electric field. This may explain why increasing the charge of the Cx50 N -terminus appears to enhance the voltage -gating of the mutants (Figs 7 & 8 ) (Kalmatsky et al., 2009) . As the N-terminus lines the entrance to the pore (Maeda et al., 2009; Myers et al., 2018; Flores et al., 2020; Yue et al., 2021; Brotherton et al., 2022; Jaradat et al., 2022; Brotherton et al., 2024) it may regulate selectivity to small molecules. Our evidence suggests that this is the case: by manipulating the charge on the N -terminus of Cx50 to make it more like Cx46, we were able to create a mutant channel that gained ATP sensitivity over the wild type Cx50 hemichannel. We note that the mutation Cx50N9R increase the permeability to ATP through the hemichannel, despite this mutation being reported to decrease the conductance of the Cx50 gap junction channel (Yue et al., 2021) . However, the reverse set of mutations -changing charge on the N -terminus of Cx46 to make it more Cx50 -like did not diminish ATP permeability of the mutant Cx46 hemichannel. This shows that other mechanisms can be important. Our data suggests that residues that alter the interaction between the N-terminus and the lining of the pore -notably TM2 – play a role in selectivity (Brotherton et al., 2022; Brotherton et al., 2024). We identified Ala14 and Thr89 of Cx46 as residues that might be important in this interaction and mutated those to their equivalent in Cx50 ( A14V,T89S). Thr89 in Cx46 is equivalent to Ala88 of Cx26, which interacts with Val13 on the N terminus (Brotherton et al., 2022). Substitution of larger residues (A88V) change Cx26 channel function and are pathological (Koval, 2013; Meigh et al., 2014) . However, the double mutation in which a smaller residue is substituted (T89S) avoids the steric clash between V14 and T89 previously reported in Cx46-50 chimaeric channels which were non -functional (Yue et al., 2021) . The double mutation s A14V and T89S or R9N and N13E by themselves did not alter ATP permeability. However, when all four mutations were combined to create Cx46QM there was a considerable reduction of ATP release compared to Cx46WT. The effect of the quadruple mutation FITC permeation was not statistically significant. This suggests that rather than simply resulting in a poorly permeable channel, the quadruple mutation has a somewhat selective effect on the ATP permeation pathway. Our study not only highlights the importance of the N -terminus in determining selectivity of connexin hemichannels to small molecules, but also that even in two very closely related connexins there are other structural elements that control permselectivity.

Materials and methods

Connexin mutagenesis cDNAs for the Cx26 and Cx32 genes were synthesised by Genscript and for the Cx46, Cx50, Cx36, Cx43 and Cx 31.3 genes by IDT. These were subsequently subcloned into pCAG -GS- mCherry vector prior to transfection. Point mutations were introduced using Gibson assembly. Overlapping fragments both containing the desired mutation were PCR amplified with primers (IDT). Successful mutagenesis was confirmed using Sanger sequencing (GATC Biotech). Double .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint mutants (Cx50 N9R E13N , Cx46R9N N13E ) were cloned using successive Gibson assemblies. All Cx constructs were inserted upstream of an mCherry tag, linked via a 12 AA linker (GVPRARDPPVAT). Cell culture Parental HeLa DH cells (ECACC 96112022 RRID:CVCL_2483) were cultured with low-glucose DMEM (Merck Life Sciences UK Ltd, CAT# D6046) supplemented with 10% foetal bovine serum (Labtech.com, CAT# FCS-SA) and 5% penicillin/streptomycin. Cells were seeded onto coverslips at a density of 4x104 cells per well. Cells were transiently transfected to co-express one connexin isoform and one genetically encoded fluorescent sensor: pDisplay-GRAB_ATP1.0-IRES-mCherry-CAAX was a gift from Yulong Li (Addgene plasmid # 167582 ; http://n2t.net/addgene:167582 ; RRID:Addgene_167582) (Wu et al., 2022). pAEMXT-eLACCO1.1 was a gift from Robert Campbell (Addgene plasmid # 167946 ; http://n2t.net/addgene:167946 ; RRID:Addgene_167946) (Nasu et al., 2021) . To improve expression of eLACCO1.1, this construct was subcloned into the iGluSnFR exxpression vector backbone. Sequences were verified with Sanger sequencing (GATC). pCMV(MinDis).iGluSnFR was a gift from Loren Looger (Addgene plasmid # 41732 ; http://n2t.net/addgene:41732 ; RRID:Addgene_41732) (Marvin et al., 2013). A mixture of 1 µg of DNA from pCAG -Cx-mCherry construct and 1 µg sensor with 3 µg PEI was added to cells for 4-8h. Cells were imaged 48 hours after transfection. aCSF Control (20 mmHg pCO2) – 140 mM NaCl, 10 mM NaHCO3, 1.25 mM NaH2PO4, 3mM KCl, 1 mM MgSO4. Control (35 mmHg pCO2) – 124 mM NaCl, 26 mM NaHCO 3, 1.25 mM NaH2PO4, 3mM KCl, 1 mM MgSO4. Hypercapnic (55 mmHg pCO2) - 100 mM NaCl, 50 mM NaHCO3, 1.25 mM NaH2PO4, 3 mM KCl, 1 mM MgSO4. Hypercapnic (70 mmHg pCO2) - 70 mM NaCl, 80 mM NaHCO 3, 1.25 mM NaH2PO4, 3 mM KCl, 1 mM MgSO4. High K+ (20 mmHg) – 47 mM NaCl, 10 mM NaHCO3, 1.25 mM NaH2PO4, 50mM KCl, 1 mM MgSO4. For 35 mmHg, the recipe is the same but with 77 mM NaCl. All solutions ha d 10 mM D -glucose and 2 mM CaCl 2 (or MgCl2 where [Ca2+]0 solution desired) added just before use and saturated with 98% O2/2% CO 2 (20 mmHg), 95% O 2/5% CO 2 (carbogen) (35 mmHg), or carbogen plus CO2 (55 mmHg and 70 mmHg). CO2 in all solutions was adjusted to give a pH of ~7.4. Live cell fluorescence imaging and analysis .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Cells were transiently transfected with one pCAG -Cx-mCherry construct and one genetically encoded fluorescent sensor 48 hours prior to imaging. Cells were perfused with control aCSF until a stable baseline was reached, before perfusion with either hypercapnic or high K+ aCSF. Once a stable baseline was reached after solution change, cells were again perfused with control aCSF and when a stable baseline reached, recordings were calibrated by direct application of 3 µM of the corresponding analyte. All cells were imaged by epifluorescence ( Scientifica Slice Scope, Cairn Research OptoLED illumination, 60x water Olympus immersion objective, NA 1.0, Hamamatsu ImagEM EM -SSC camera, Metafluor software). cpGFP in the sensors were excited by a 470 nm LED, with emission captured between 504 -543 nm. Connexin constructs have a C -terminal mCherry tag, which is excited by a 535 nm LED and emission captured between 570 -640 nm. Only cells expressing both cpGFP and mCherry were selected for recording, with cpGFP images acquired ev ery 4 seconds. For each condition, at least 3 independent transfections were performed with at least 2 coverslips per transfection. Analysis of all experiments was carried out in ImageJ. Images were opened as a stack and stabilised (Li, 2008). ROIs were drawn around cells co -expressing both sensor and connexin. Median pixel intensity was plotted as normalised fluorescence change (∆F/F 0) versus time to give traces of fluorescence change. The a mount of analyte release was quantified as concentration by normalising to the ∆F/F 0 caused by application of 3 µM of analyte, which was within the linear portion of the dose response curve for each sensor . Release froma single cell was considered to be a statistical replicate. Dye loading Cells were transiently transfected with pCAG-Cx43-mCherry, pCAG-Cx46WT-mCherry or pCAG- Cx43QM-mCherry constructs 48 hours prior to experiments. After washing in 20 mmHg aCSF , cells were perfused with one of the following solutions: 20 mmHg aCSF, 55 mmHg aCSF (Cx43 only), Ca2+-free aCSF (Cx43 only) or 50 mM KCl aCSF (Cx46WT and Cx46QM only) each containing 50 µM fluor escein isothiocyanate (FITC) for 10 mins. The cells were then washed, first by perfusion with 20 mmHg in the absence of FITC, before being placed in serial washes of 20 mmHg aCSF. Cells were fixed in 4% paraformaldehyde for 30 mins before being washed 3 times with PBS. Coverslips were mounted inverted on a microscope slide using Fluorshield™ with DAPI mounting medium (Sigma-Aldrich, Cat# F6057). Images were taken using the Zeiss-880 or Zeiss- 980 confocal LSM, specifically using the 488 and 561 nm lasers. Subsequent analysis was done using the FIJI software. For Figure 9, figure supplement 1 , the median pixel intensity for each transfection was considered to be a statistical replicate. Statistical analysis All quantitative data is presented as box and whisker plots, where the line represents the median, the box is the interquartile range, and the whiskers are the range, with all individual data points included. Pairwise Mann Whitney U-tests were carried out with corrections for multiple comparisons using the false discovery method, with a maximum false discovery set to 0.05 (Curran-Everett, 2000). In text, all data is presented as median (95% CI, upper, lower limit). All calculations were performed on GraphPad Prism. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Author Contributions: AL: Conceptualisation, data curation, investigation, writing – original draft, review and editing, formal analysis . JB: Investigation (dye loading and confocal microscopy, sensor optimisation). ND: Conceptualisation, supervision, writing – review and editing. Funding: JB was supported by the Biotechnology and Biological Sciences Research Council (BBSRC) and University of Warwick funded Midlands Integrative Biosciences Training Partnership (MIBTP) grant number BB/T00746X/1. AL was funded by the Medical Research Council through the University of Warwick Doctoral Training Partnership, grant number MR/N014294/1 Conflicts of interest: The authors declare that there are no conflicts of interest.

Acknowledgements

We thank Prof Alexander Cameron for reading a draft of this paper.

References

Barros LF (2013) Metabolic signaling by lactate in the brain. Trends Neurosci 36:396-404. Bassnett S (2002) Lens organelle degradation. Experimental eye research 74:1-6. Bennett BC, Purdy MD, Baker KA, Acharya C, McIntire WE, Stevens RC, Zhang Q, Harris AL, Abagyan R, Yeager M (2016) An electrostatic mechanism for Ca(2+)-mediated regulation of gap junction channels. Nat Commun 7:8770. Berthoud VM, Ngezahayo A (2017) Focus on lens connexins. BMC cell biology 18:6. Brotherton DH, Savva CG, Ragan TJ, Dale N, Cameron AD (2022) Conformational changes and CO2-induced channel gating in connexin26. Structure 30:697-706. Brotherton DH, Nijjar S, Savva CG, Dale N, Cameron AD (2024) Structures of wild -type and a constitutively closed mutant of connexin26 shed light on channel regulation by CO2. eLife 13:RP93686. Buhl DL, Harris KD, Hormuzdi SG, Monyer H, Buzsaki G (2003) Selective impairment of hippocampal gamma oscillations in connexin -36 knock -out mouse in vivo. J Neurosci 23:1013-1018. Butler J, Dale N (2023) X -linked Charcot Marie Tooth mutations alter CO(2) sensitivity of connexin32 hemichannels. Front Cell Neurosci 17:1330983. Chever O, Lee CY, Rouach N (2014) Astroglial connexin43 hemichannels tune basal excitatory synaptic transmission. J Neurosci 34:11228-11232. Cisterna BA, Arroyo P, Puebla C (2019) Role of Connexin -Based Gap Junction Channels in Communication of Myelin Sheath in Schwann Cells. Front Cell Neurosci 13:69. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Condorelli DF, Parenti R, Spinella F, Salinaro AT, Belluardo N, Cardile V, Cicirata F (1998) Cloning of a new gap junction gene (Cx36) highly expressed in mammalian brain neurons. European Journal of Neuroscience 10:1202-1208. Curran-Everett D (2000) Multiple comparisons: philosophies and illustrations. Am J Physiol Regul Integr Comp Physiol 279:R1-8. Dale N (2021) CO2 sensing by connexin26 and its role in the control of breathing. Interface Focus 11:20200029. Deans MR, Gibson JR, Sellitto C, Connors BW, Paul DL (2001) Synchronous activity of inhibitory networks in neocortex requires electrical synapses containing connexin36. Neuron 31:477-485. Dospinescu V-M, Nijjar S, Spanos F, Cook J, de Wolf E, Biscotti MA, Gerdol M, Dale N (2019) Structural determinants of CO2 -sensitivity in the β connexin family suggested by evolutionary analysis. Communications Biology 2:331. Dospinsecu VM, Mascarenhas A, Butler J, Nijjar S, de Oliveira Taborda K, Connors S, Huang L, Dale N (2025) CO2 -dependent opening of Connexin 43 hemichannels. bioRxiv:2025.2001.2013.632711. Flores JA, Haddad BG, Dolan KA, Myers JB, Yoshioka CC, Copperman J, Zuckerman DM, Reichow SL (2020) Connexin -46/50 in a dynamic lipid environment resolved by CryoEM at 1.9 A. Nat Commun 11:4331. Gomez-Hernandez JM, de Miguel M, Larrosa B, Gonzalez D, Barrio LC (2003) Molecular basis of calcium regulation in connexin -32 hemichannels. Proc Natl Acad Sci U S A 100:16030 - 16035. Hansen DB, Braunstein TH, Nielsen MS, MacAulay N (2014) Distinct permeation profiles of the connexin 30 and 43 hemichannels. FEBS Lett 588:1446-1457. Harris AL (2007) Connexin channel permeability to cytoplasmic molecules. Prog Biophys Mol Biol 94:120-143. Huckstepp RT, Eason R, Sachdev A, Dale N (2010a) CO2-dependent opening of connexin 26 and related beta connexins. J Physiol 588:3921-3931. Huckstepp RT, id Bihi R, Eason R, Spyer KM, Dicke N, Willecke K, Marina N, Gourine AV, Dale N (2010b) Connexin hemichannel-mediated CO2-dependent release of ATP in the medulla oblongata contributes to central respiratory chemosensitivity. J Physiol 588:3901-3920. Iacobas DA, Iacobas S, Spray DC (2007) Connexin -dependent transcellular transcriptomic networks in mouse brain. Prog Biophys Mol Biol 94:169-185. Iacobas DA, Urban-Maldonado M, Iacobas S, Scemes E, Spray DC (2003) Array analysis of gene expression in connexin-43 null astrocytes. Physiol Genomics 15:177-190. Imamura H, Nhat KP, Togawa H, Saito K, Iino R, Kato -Yamada Y, Nagai T, Noji H (2009) Visualization of ATP levels inside single living cells with fluorescence resonance energy transfer-based genetically encoded indicators. Proc Natl Acad Sci U S A 106:15651 - 15656. Jaradat R, Li X, Chen H, Stathopulos PB, Bai D (2022) The Hydrophobic Residues in Amino Terminal Domains of Cx46 and Cx50 Are Important for Their Gap Junction Channel Ion Permeation and Gating. Int J Mol Sci 23. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Johnson RG, Le HC, Evenson K, Loberg SW, Myslajek TM, Prabhu A, Manley AM, O'Shea C, Grunenwald H, Haddican M, Fitzgerald PM, Robinson T, Cisterna BA, Saez JC, Liu TF, Laird DW, Sheridan JD (2016) Connexin Hemichannels: Methods for Dye Uptake and Leakage. J Membr Biol 249:713-741. Kalmatsky BD, Bhagan S, Tang Q, Bargiello TA, Dowd TL (2009) Structural studies of the N - terminus of Connexin 32 using 1H NMR spectroscopy. Arch Biochem Biophys 490:9-16. Kang J, Kang N, Lovatt D, Torres A, Zhao Z, Lin J, Nedergaard M (2008) Connexin 43 hemichannels are permeable to ATP. J Neurosci 28:4702-4711. Koval M (2013) Drowning out communication. Focus on "The human Cx26-D50A and Cx26-A88V mutations causing keratitis -ichthyosis-deafness syndrome display increased hemichannel activity". Am J Physiol Cell Physiol 304:C1129-1130. Lee SN, Cho HJ, Jeong H, Ryu B, Lee HJ, Kim M, Yoo J, Woo JS, Lee HH (2023) Cryo-EM structures of human Cx36/GJD2 neuronal gap junction channel. Nat Commun 14:1347. Li H, Liu TF, Lazrak A, Peracchia C, Goldberg GS, Lampe PD, Johnson RG (1996) Properties and regulation of gap junctional hemichannels in the plasma membranes of cultured cells. J Cell Biol 134:1019-1030. Li K (2008) The image stabilizer plugin for ImageJ. In. Liang WG, Su CC, Nian JH, Chiang AS, Li SY, Yang JJ (2011) Human connexin30.2/31.3 (GJC3) does not form functional gap junction channels but causes enhanced ATP release in HeLa cells. Cell Biochem Biophys 61:189-197. Lopez W, Ramachandran J, Alsamarah A, Luo Y, Harris AL, Contreras JE (2016) Mechanism of gating by calcium in connexin hemichannels. Proc Natl Acad Sci U S A 113:E7986-E7995. Lucaciu SA, Leighton SE, Hauser A, Yee R, Laird DW (2023) Diversity in connexin biology. J Biol Chem 299:105263. Maeda S, Nakagawa S, Suga M, Yamashita E, Oshima A, Fujiyoshi Y, Tsukihara T (2009) Structure of the connexin 26 gap junction channel at 3.5 A resolution. Nature 458:597-602. Marvin JS, Borghuis BG, Tian L, Cichon J, Harnett MT, Akerboom J, Gordus A, Renninger SL, Chen TW, Bargmann CI, Orger MB, Schreiter ER, Demb JB, Gan WB, Hires SA, Looger LL (2013) An optimized fluorescent probe for visualizing glutamate neurotransmission. Nat

Methods

10:162-170. Mathias RT, White TW, Gong X (2010) Lens gap junctions in growth, differentiation, and homeostasis. Physiol Rev 90:179-206. Meigh L, Hussain N, Mulkey DK, Dale N (2014) Connexin26 hemichannels with a mutation that causes KID syndrome in humans lack sensitivity to CO2. eLife 3:e04249. Meigh L, Greenhalgh SA, Rodgers TL, Cann MJ, Roper DI, Dale N (2013) CO2 directly modulates connexin 26 by formation of carbamate bridges between subunits. eLife 2:e01213. Mese G, Richard G, White TW (2007) Gap junctions: basic structure and function. J Invest Dermatol 127:2516-2524. Myers JB, Haddad BG, O'Neill SE, Chorev DS, Yoshioka CC, Robinson CV, Zuckerman DM, Reichow SL (2018) Structure of native lens connexin 46/50 intercellular channels by cryo- EM. Nature 564:372-377. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Nasu Y, Murphy-Royal C, Wen Y, Haidey JN, Molina RS, Aggarwal A, Zhang S, Kamijo Y, Paquet ME, Podgorski K, Drobizhev M, Bains JS, Lemieux MJ, Gordon GR, Campbell RE (2021) A genetically encoded fluorescent biosensor for extracellular L -lactate. Nat Commun 12:7058. Nielsen BS, Zonta F, Farkas T, Litman T, Nielsen MS, MacAulay N (2019) Structural determinants underlying permeant discrimination of the Cx43 hemichannel. J Biol Chem 294:16789 - 16803. Nijjar S, Brotherton DH, Butler J, Dospinescu V-M, Gannon H, Linthwaite VL, Cann MJ, Cameron A, Dale N (2025) Multiple carbamylation events are required for differential modulation of Cx26 hemichannels and gap junctions by CO2. The Journal of physiology. Oshima A, Tani K, Hiroaki Y, Fujiyoshi Y, Sosinsky GE (2007) Three -dimensional structure of a human connexin26 gap junction channel reveals a plug in the vestibule. Proc Natl Acad Sci U S A 104:10034-10039. Oshima A, Tani K, Hiroaki Y, Fujiyoshi Y, Sosinsky GE (2008) Projection structure of a N-terminal deletion mutant of connexin 26 channel with decreased central pore density. Cell Commun Adhes 15:85-93. Pearson RA, Dale N, Llaudet E, Mobbs P (2005) ATP released via gap junction hemichannels from the pigment epithelium regulates neural retinal progenitor proliferation. Neuron 46:731- 744. Pinto BI, Garcia IE, Pupo A, Retamal MA, Martinez AD, Latorre R, Gonzalez C (2016) Charged Residues at the First Transmembrane Region Contribute to the Voltage Dependence of the Slow Gate of Connexins. J Biol Chem 291:15740-15752. Piva TJ, McEvoy -Bowe E (1998) Oxidation of glutamine in HeLa cells: role and control of truncated TCA cycles in tumour mitochondria. J Cell Biochem 68:213-225. Qi C, Acosta Gutierrez S, Lavriha P, Othman A, Lopez-Pigozzi D, Bayraktar E, Schuster D, Picotti P, Zamboni N, Bortolozzi M, Gervasio FL, Korkhov VM (2023) Structure of the connexin - 43 gap junction channel in a putative closed state. eLife 12. Saez JC, Schalper KA, Retamal MA, Orellana JA, Shoji KF, Bennett MV (2010) Cell membrane permeabilization via connexin hemichannels in living and dying cells. Exp Cell Res 316:2377-2389. San Martin A, Ceballo S, Ruminot I, Lerchundi R, Frommer WB, Barros LF (2013) A genetically encoded FRET lactate sensor and its use to detect the Warburg effect in single cancer cells. PloS one 8:e57712. Sargiannidou I, Ahn M, Enriquez AD, Peinado A, Reynolds R, Abrams C, Scherer SS, Kleopa KA (2008) Human oligodendrocytes express Cx31.3: function and interactions with Cx32 mutants. Neurobiol Dis 30:221-233. Srinivas M, Rozental R, Kojima T, Dermietzel R, Mehler M, Condorelli DF, Kessler JA, Spray DC (1999) Functional properties of channels formed by the neuronal gap junction protein connexin36. J Neurosci 19:9848-9855. Turovsky EA, Braga A, Yu Y, Esteras N, Korsak A, Theparambil SM, Hadjihambi A, Hosford PS, Teschemacher AG, Marina N, Lythgoe MF, Haydon PG, Gourine AV (2020) Mechanosensory Signaling in Astrocytes. J Neurosci 40:9364-9371. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint van de Wiel J, Meigh L, Bhandare A, Cook J, Nijjar S, Huckstepp R, Dale N (2020) Connexin26 mediates CO2-dependent regulation of breathing via glial cells of the medulla oblongata. Communications Biology 3:521. Weber PA, Chang HC, Spaeth KE, Nitsche JM, Nicholson BJ (2004) The permeability of gap junction channels to probes of different size is dependent on connexin composition and permeant-pore affinities. Biophys J 87:958-973. Weissman TA, Riquelme PA, Ivic L, Flint AC, Kriegstein AR (2004) Calcium waves propagate through radial glial cells and modulate proliferation in the developing neocortex. Neuron 43:647-661. Wu Z, He K, Chen Y, Li H, Pan S, Li B, Liu T, Xi F, Deng F, Wang H, Du J, Jing M, Li Y (2022) A sensitive GRAB sensor for detecting extracellular ATP in vitro and in vivo. Neuron 110:770-782 e775. Yue B, Haddad BG, Khan U, Chen H, Atalla M, Zhang Z, Zuckerman DM, Reichow SL, Bai D (2021) Connexin 46 and connexin 50 gap junction channel properties are shaped by structural and dynamic features of their N-terminal domains. J Physiol 599:3313-3335. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figures and Legends Figure 1. HeLa cells expressing the g enetically encoded fluorescent sensors GRAB ATP, iGluSnFr and eLACCO1.1 alone do not respond to connexin gating stimuli. A. Representative images of cells at 55 mmHg PCO2 (inset is 20 mmHg control), 50 mM KCl, and after application of 3 µM corresponding analyte. Scale bar represents 20 µm. B. Representative traces showing sensor responses to 50 mM KCl (green bar), 55 mmHg (blue bar), and 3 µM corresponding analyte (red bar). C. Summary data showing median ∆F/F0 for ATP (n = 16 cells), glutamate (n = 8) and lactate (n = 8). 200s 10% ΔF/F0 100s 10% ΔF/F0 200s ATPGlutamateLactate 55 mmHg 50 mM KCl20 mmHg 3 µM A B C 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 200s 200s 200s CO 2 50 mM KCl 3 µMCO 2 50 mM KCl 3 µMCO 2 50 mM KCl 3 µM -0.1 0.0 0.1 0.2 0.3 0.4 0.5 ΔF/F0 ATP LactateGlutamate .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 2. Cx26 mediates release of ATP, glutamate and lactate. A. Representative images showing cells from each sensor under each condition. Scale bar represents 20 µm. B. Representative traces of normalised fluorescence changes in response to 55 mmHg PCO2 (blue bar), 50 mM KCl (green bar) and a 3 µM calibration of the corresponding analyte (red bar). C. Summary data showing the median release for ATP (n = 15 cells), glutamate (n = 19) and lactate (n = 27) through Cx26 hemichannels stimulated either by 55 mmHg (filled circles), or 50 mM KCl (open circles). Mann Whitney U-test, p < 0.0001 (55 mmHg ATP vs glutamate, 55 mmHg ATP vs lactate). The order of stimulus (CO2 or KCl) was regularly reversed between recordings to avoid any potential depletion of ATP release. Data present is from at least 3 independent transfections. 55 mmHg 50 mM KCl20 mmHgATPGlutamateLactate 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 100s 200s 100s A B C -5 0 5 10Concentration (µM) ATP Glutamate Lactate <0.0001 <0.0001 10% ΔF/F0 200s -5 0 5 10Concentration (µM) ATP Glutamate Lactate <0.0001 <0.0001 .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 2, figure supplement 1. Co-expression of Cx-mCherry and sensor (GFP). Using Cx36 as an example, the images demonstrate that cells were selected only if they co-expressed both the connexin construct, which is tagged with mCherry (top row), and the sensor, which is tagged with GFP (bottom row). Scale bar represents 20 µm. Cx36-mCherry GRABATP iGluSnFR eLACCO1.1 .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 3. Cx32 mediates release of ATP, glutamate and lactate. A. Representative images showing cells from each sensor under each condition. Scale bar represents 20 µm. B. Representative traces of normalised fluorescence changes in response to 70 mmHg PCO2 (blue bar), 50 mM KCl (green bar) and a 3 µM calibration of the corresponding analyte (red bar). C. Summary data showing the median release of: ATP (n = 18), glutamate (n = 27 for CO2, n = 14 for KCl), and lactate (n = 15), through Cx32 showing changes in release depending on whether the channel was opened by CO 2 (open circles), or a depolarising stimulus (closed circles). Data presented is from at least 3 independent transfections. Mann Whitney U-test, p < 0.0001 (CO 2 ATP vs glutamate), p < 0.0001 (CO2 ATP vs lactate), p < 0.0001 (50 mM KCl ATP vs glutamate), p < 0.0001 (50 mM KCl ATP vs lactate). 70 mmHg 50 mM KCl35 mmHg ATPGlutamateLactate 10% ΔF/F0 200s 20% ΔF/F0 200s 10% ΔF/F0 100s 10% ΔF/F0 200s 20% ΔF/F0 200s A B C 0 4 8 12Concentration (µM) ATP Glutamate Lactate <0.0001 <0.0001 .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 4. Cx43 mediates release of ATP, glutamate and lactate. A. Representative images showing cells from each sensor under each condition. Scale bar represents 20 µm. B. Representative traces of normalised fluorescence changes in response to 55 mmHg PCO2 (blue bar), 50 mM KCl (green bar) and a 3 µM calibration of the corresponding analyte (red bar). C. Summary data showing the median release of analytes through Cx43: ATP (n = 33 (CO 2), n = 17 (high KCl)), glutamate (n = 19), lactate (n = 18). Data shows CO 2-dependent release (open circles), and release to a depolarising stimulus (closed circles). Data presented is from at least 3 independent transfections. Mann Whitney U-test, p < 0.0001 (CO 2 ATP vs glutamate), p < 0.0001 (CO2 ATP vs lactate). 200s 10% ΔF/F0 200s 20% ΔF/F0 200s 55 mmHg 50 mM KCl20 mmHgATPGlutamateLactate 10% ΔF/F0 200s 100s 20% ΔF/F0 20% ΔF/F0 200s A B C -5 0 5 10Concentration (µM) <0.0001 <0.0001 ATP Glutamate Lactate .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 4 , Figure supplement 1. Low [Ca 2+]ext is sufficient to open Cx43 hemichannels. A. Representative images of FITC loaded into mCherry-tagged Cx43 positive cells. Scale bar 20 µm. B. Summary data showing the median pixel intensity under control (20 mmHg) (n = 19), hypercapnia (55 mmHg) (n = 18) and low [Ca2+]ext (n = 17). Mann Whitney U-test, p < 0.0001 (low [Ca2]ext vs 20 mmHg). 20 mmHg 55 mmHg Low [Ca2+]ext 20 mmHg55 mmHg low [Ca 2 ]ext 0 30 60 90 Pixel Intensity <0.0001 20mmHg55 mmHglow [Ca 2+]ext .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 4, figure supplement 2 . Removal of extracellular Ca 2+ changes the permeability of Cx43 to physiological molecules . A. Representative images showing cells from each sensor under each condition. Scale bar represents 20 µm. B. Representative traces of normalised fluorescence changes in response to 55 mmHg pCO 2 (blue bar), low [Ca2+]ext (yellow bar) and a 3 µM calibration of the corresponding analyte (red bar). C. Summary data showing the median release of analytes through Cx43 and how it differs depending on the stimulus: ATP (n = 33 for CO2, n = 17 for high KCl), glutamate (n = 19), lactate (n = 18). Data presented is from a t least 3 independent transfections. For 55 mmHg vs low [Ca2+]ext for ATP, glutamate and lactate release, Mann Whitney U-test was performed, p < 0.0001. 55 mmHg Low [Ca2+]ext ATP GlutamateLactate 20 mmHg 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 200s A B C ATP CO2ATP 0 Ca Glutamate CO2Glutamate 0 Ca Lactate CO2Lactate 0 Ca -5 0 5 10 15Concentration (µM) ATP GlutamateLactate +55 mmHg Low [Ca2+]ext + + + + + <0.0001 <0.0001 <0.0001 ATP .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 5. Cx31.3 mediates release of ATP, glutamate and lactate. A. Representative images showing cells from each sensor under each condition. Scale bar represents 20 µm. B. Representative traces of normalised fluorescence changes in response to 50 mM KCl (green bar), and 3 µM analyte (red bar). C. Summary data showing the median ATP (n = 32) , glutamate (n = 19) and lactate (n= 50) release from Cx31.3 hemichannels in response to a depolarising stimulus. Data presented is from at least 3 independent transfections. Mann Whitney U-test, p = 0.0371 (ATP vs glutamate) and p=0.004 (lactate vs glutamate). 20% ΔF/F0 10% ΔF/F0 100s 200s 20 mmHg 50 mM KCl ATPGlutamate A B C Lactate 10% ΔF/F0 100s ATP Glutamate Lactate 0 1 2 3 4Concentration (µM) 0.03 0.004 .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 6. Cx36 hemichannels mediate release of ATP but are impermeable to glutamate and lactate. A. Representative images showing cells from each sensor under each condition: ATP (n = 24), glutamate (n = 22), lactate (n = 15). Scale bar represents 20 µm. B. Representative traces of normalised fluorescence changes in response to 50 mM KCl (green bar) and 3 µM analyte (red bar). C. Summary data depicting ATP release from Cx36 and no permeability to glutamate or lactate. Mann Whitney U-test, p < 0.0001 (ATP vs glutamate, ATP vs lactate). ATPGlutamateLactate 10% ΔF/F0 100s 10% ΔF/F0 100s 10% ΔF/F0 100s 100s 100s 10% ΔF/F0 10% ΔF/F0 A B C 3 mM KCl50 mM KCl ATP Glutamate Lactate -4 -2 0 2 4 6 Cx36 Concentration (µM) .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 7. Cx46 and Cx50 have opposing permeability profiles. A. Representative traces showing normalised fluorescence changes in response to 55 mmHg (blue bar), 50 mM KCl (green bar) and 3 µM (analyte). B. Summary data depicting analyte release under each stimulus: Cx46 ATP (n = 40), Cx46 glutamate (n = 25), Cx46 lactate (n = 16), Cx50 ATP (n = 28), Cx50 glutamate (n = 12), Cx50 lactate (n = 11). Cx50 is CO 2 sensitive, and thus release is stimulated by 55 mmHg pCO 2 (open circles), and 50 mM KCl (closed circles). Data presented is from 3 independent transfections. 200s 20% ΔF/F0 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 10% ΔF/F0 5% ΔF/F0 100s 100s 100s 200s 100s 100s Cx50 ATP Cx50 glutamate Cx50 lactate Cx46 ATP Cx46 glutamate Cx46 lactate A B -2 0 2 4 6 8 Concentration (µM) Cx46 Cx50 ATPGlutamateLactate ATP Glutamate Lactate .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 8. N -terminal mutations make Cx50 permeable to ATP. A. Sequence alignments for residues 1-26 of Cx46 and Cx50, to show the residues in Cx50 that have been mutated to those in Cx46 (purple boxes). B. Representative images showing cells from each sensor under each condition. Scale bar represents 20 µm. C. Representative traces of normalised fluorescence changes in response to 55 mmHg (blue bar), 50 mM KCl (green bar) and 3 µM ATP (red bar). D. Summary data showing ATP release in response to depolarisation from Cx50 mutants: N9R (n = 26), N9K (n = 13), E13N (n = 18), N9R E13N ( n = 11). Data shows CO 2-dependent release (open circles), and release to a depolarising stimulus (closed circles). The dotted lines show the median and 95% confidence limits for ATP release from Cx46 in response to depolarisation. Mann Whitney U-test, p = 0.0018, p = 0.0155 and p = 0.0004 (N9R, N9K and E13N respectively, where release is evoked by 50 mM KCl, compared to Cx46 WT). Mann Whitney U-test, p = 0.96 (Cx46WT vs Cx50N9R E13N). E. Summary data demonstrating CO2-evoked ATP release form the Cx50 mutant channels N9R (n = 26), N9K (n = 13), E13N (n = 18), N9R E13N (n = 11). MGDWSFLGRLLENAQEHSTVIGKVWL MGDWSFLGNILEEVNEHSTVIGRVWLCx50/1-26 Cx46/1-26A N9R 55 mmHg 50 mM KCl N9KE13NN9R E13N 20 mmHg 200s 10% ΔF/F0 200s 10% ΔF/F0 200s 10% ΔF/F0 100s 10% ΔF/F0 200s 10% ΔF/F0 200s 10% ΔF/F0 200s 10% ΔF/F0 100s 10% ΔF/F0 ED B C N9R N9K E13N N9R E13N 0 1 2 3[ATP] release (µM) 55 mmHg N9R N9K E13N N9R E13N 0 1 2 3 4[ATP] release (µM) 50 mM KCl N9R N9K E13N N9R E13N 0 1 2 3 4[ATP] release (µM) 50 mM KCl N9R N9K E13N N9R E13N 0 1 2 3[ATP] release (µM) 55 mmHg .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 9. The permeability of Cx46 to ATP appears to be reliant on N-terminal charge and interactions between the N-terminus and TM2. A. Sequence alignments from 1-26 and 80-89 for Cx50 and Cx46, to show the residues in Cx46 which have been mutated to correspond to those in Cx50 (red). B. Representative images showing cells expressing Cx46WT, Cx46R9N N13E and Cx46QM under each condition. Scale bar represents 20 µm. C. Representative traces showing normalised change in fluorescence to 50 mM KCl (green bar), and 3 µM ATP (red bar). D. Summary data depicting the median ATP release from Cx46WT (the dotted lines show the median and 95% confidence interval for Cx46WT), Cx46R9N N13E (n = 7) and Cx46QM (n = 9). Mann Whitney U- test, p < 0.0001 (Cx46WT vs Cx46QM). .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Figure 9, figure supplement 1. Cx46QM remains permeable to FITC. A, Representative images showing cells expressing Cx46WT or Cx46QM loaded with FITC under each condition. Scale bar represents 15 µm. B, Summary data depicting the median change in fluorescence pixel intensity between 20 mmHg and 50 mM KCl for cells expressing Cx46WT or Cx46QM from 5 independent transfections. There is no significant difference in FITC loading into the Cx46 QM compared to Cx46WT (Mann Whitney U Test). For Cx46WT, 43 cells in 20 mmHg and 35 cells in 50 mM KCl were analyzed. For Cx46QM, 79 cells in 20 mmHg and 79 cells for Cx46QM in 50 mM KCl were analyzed. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Table 1. Calculations of the Nernstian equilibrium potential for ATP, glutamate and lactate in HeLa cells. We assigned the valence of ATP under the assumption that it is chelated to Mg2+. Driving force is calculated assuming a resting potential of -60 mV. The intracellular concentrations were taken from studies that investigated HeLa cells: a, (Imamura et al., 2009); b, (Piva and McEvoy-Bowe, 1998); c, (San Martin et al., 2013). Analyte Valence [Analyte]I (mM) [Analyte]o (mM) Veq (mV) Driving force at rest ATP -2 10a 10-4 145 205 Glutamate -1 20b 10-4 307 367 Lactate -1 1c 10-4 232 292 .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint Connexin ATP Glu Lac Cx26 1 3.5 (0.44) 3.1 (0) Cx32 1 3.0 (1.73) 5.5 (1.83) Cx43 1 2.4 (0.4) 2.0 (0) Cx31.3 1 0.9 1.1 Cx36 1 0 0 Cx46 1 0 0 Cx50 0 1.5 (0) µM 5.8 (3.4) µM Table 2. Relative release of ATP, glutamate and lactate of connexin hemichannels normalised to amount of ATP release. The numbers given are for the hypercapnic stimulus, except for Cx31.3, Cx36 and Cx46 which are insensitive to CO2. For those channels sensitive to CO2, the numbers in brackets are normalised values for the depolarising stimulus. For Cx50, which is not permeable to ATP, the numbers are absolute concentration in µM. The ratio of the Nernstian driving force for each analyte (normalised to that of ATP) during the hypercapnic stimulus is 1: 1.8:1.4, during the depolarising stimulus it is (1:2.1:1.6) assuming that the membrane depolarises to 0 mV. .CC-BY 4.0 International licensemade available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is The copyright holder for this preprintthis version posted March 12, 2025. ; https://doi.org/10.1101/2025.03.12.642803doi: bioRxiv preprint

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source — PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

Ask this paper AI returns verbatim quotes from the full text · source: oa-pdf

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. This is a recent paper (2025) — citers typically take a year or two to land, and the OpenAlex reference graph may still be filling in.

Source provenance

europepmc
last seen: 2026-05-20T01:45:00.602351+00:00
unpaywall
last seen: 2026-05-21T05:10:58.409756+00:00
License: CC-BY-4.0