Anomalous behavior of membrane fluidity caused by copper-copper bond coupled phospholipids

Membrane fluidity, essential for cell functions, is obviously affected by copper, but the molecular mechanism is poorly understood. Here, we unexpectedly observed that a decrease in phospholipid (PL) bilayer fluidity caused by Cu2+ was more significant than those by Zn2+ and Ca2+, while a comparable reduction occurred in the last two ions. This finding disagrees with the placement in the periodic table of Cu just next to Zn and far from Ca. The physical nature was revealed to be an anomalous attraction between Cu+ cations, as well as the induced motif of two phospholipids coupled by Cu-Cu bond (PL-diCu-PL). Namely, upon Cu2+ ion binding to a negatively charged phosphate group of lipid, Cu2+ was reduced to Cu+. The attraction of the cations then caused one Cu+ ion simultaneously binding to two lipids and another Cu+, resulting in the formation of PL-diCu-PL structure. In contrast, this attraction cannot occur in the cases of Zn and Ca ions. Remarkably, besides lipids, the phosphate group also widely exists in other biological molecules, including DNA, RNA, ADP and ATP. Our findings thus provide a new view for understanding the biological functions of copper and the mechanism underlying copper-related diseases, as well as lipid assembly.

Proper fluidity of the biological membrane is critically essential for numerous cell functions, such as adapting to the thermal stress of the environment of the microorganism 1 , the binding of peripheral proteins associated at the lipid surface 2 , reaction rates of enzymes 1 , and even cell signaling and phagocytosis 3 . Both in-vivo and in-vitro evidences have indicated that copper, as a biologically trace element, plays an important role in the membrane, especially in regard to its fluidity [4][5][6][7][8] . However, the underlying mechanism is still far from being understood, partially because researches have been majorly devoted to the interactions of alkali and alkaline earth metal ions with phospholipids as well as the influences on the lipid bilayer [9][10][11][12][13][14][15][16][17][18] . Traditionally, the effect of metal ions on the membrane was majorly attributed to the electrostatic attraction with lipid headgroups 19 . This can explain the impact of divalent metal ions on membrane fluidity more than that of monovalent ones but cannot be applied to the differences of the influences of divalent ions, such as Ca 2+ , Mg 2+ , Zn 2+ and Cu 2+ . Recently, Cremer and his coworkers studied effects of the Cu 2+ ion on a bilayer comprised of both phosphatidylcholine (PC) and phosphatidylserine (PS), and proposed that the ion was specifically bound to PS 20 . Meanwhile, this binding was only stable under basic conditions, but not at acidic pH values. Further investigations suggested that a complex of Cu(PS) 2 formed upon Cu 2+ binding to PS molecules, which did not alter the net negative charge on the membrane 21 . This differed from the manner and impact of Ca 2+ or Mg 2+ binding. Very recently, it was determined that the cis isomer of the Cu(PS) 2 complex was preferred to the trans one 22 . Additionally, a synergetic effect of Cu 2+ with Ca 2+ were proposed, which potentially triggered the transition of PS membrane from fluid phase to soft solid phase 23 . Besides those, influences of the Cu 2+ ion on a bilayer consisting of PC and phosphatidylethanolamine (PE) have also been explored 24 . It was suggested that Cu 2+ could stably bind to the amine moieties of PE lipids, while other transition metal ions to PE bound in a similar manner. Noticeably, all these Cu 2+ -lipid interactions specifically relate to the existence of amine moiety in the headgroup of the lipid.
Here, we propose an amine-independent copper-phospholipid motif in the membrane and apply it to illuminate our measurements of the anomalous effect of copper on the fluidity of a bilayer composed of PC and phosphatidylglycerol (PG). We observed that the decrease of the PC/PG bilayer fluidity caused by Cu 2+ ions was much more significant than those induced by Zn 2+ and Ca 2+ ions, while a comparable reduction occurred in the last two cases. A model of two phospholipids coupled by a Cu-Cu bond (diCu) was built to explain the unexpected behavior of the bilayer fluidity induced by the Cu 2+ ion. Namely upon the interaction of two Cu 2+ ions with two phospholipids, one ion preferred simultaneously binding with the phosphate groups of the lipids and another ion after the Cu 2+ ions were reduced to Cu + . The underlying physics was then revealed to be an anomalous 3d 10 -3d 10 attraction between Cu + cations, which was resulted from a special 3d 10 closed shell of the outermost electron structure in Cu + . In contrast, this attraction cannot occur in the cases of Zn and Ca ions due to their electron structures. Moreover, Ångström-resolution atomic force microscope (AFM) imaging also supported the formation of diCu coupled to two lipids.

Results and Discussion
To avoid the effect of amine moiety in the specific lipid (PS and PE), we applied PC and PG to prepare a lipid bilayer. The fluorescence recovery after photobleaching (FRAP) method involves the production of a concentration gradient of fluorescent molecules by irreversibly bleaching a portion of fluorophores in the observed region. The disappearance of this gradient over time is an indicator of the mobility of the fluorophores in the membrane as the fluorophore diffuses from the adjacent unbleached regions of the membrane into the bleached zone. We chose to use PC/PG/NBD (73:25:2) supported lipid bilayers for the FRAP measurements due to their large size, which allows their visualization in the microscopic field. The supported lipid bilayers appear large and uniformly fluorescent when observed through the confocal microscope. The inset in Fig. 1a shows fluorescence images of a representative FRAP experiment performed on such lipid bilayers: the dark circular region represents the bleached spot immediately after bleaching (0 s, bleach) and after the recovery of fluorescence at 80 s and 160 s (postbleach). The scanning parameters for all FRAP experiments were optimized to ensure no significant fluorescence photobleaching due to repeated imaging. Nonlinear curve fitting analysis of NBD fluorescence recovery kinetics after bleaching from the experimental data using the equations described in previous studies 25,26 in the absence and presence of a 20 mM CuCl 2 /ZnCl 2 /CaCl 2 /ZnCl 2 treatment on PC/PG bilayer mixtures are shown in Fig. 1a. The corresponding lateral diffusion rates of the NBD probe evaluated from Fig. 1a are displayed in Fig. 1b. Surprisingly, in contrast to other metal ions, copper ions led to an anomalous behavior of membrane fluidity. In the DOPC/DOPG membrane system, the rate was 0.64 ± 0.03 for the sample incubated with CuCl 2 , 1.12 ± 0.03 for ZnCl 2 , 1.24 ± 0.03 for CaCl 2 , 1.083 ± 0.001 for MgCl 2 , and 1.80 ± 0.11 for the control (incubated with NaCl 2 ). To rule out the possible quenching effects of metal ions on the fluorophore molecules in the membranes, lipid bilayers containing tail-labeled NBD-PE are used to repeat the FRAP experiments, the similar results are obtained (Fig. S7).
To further determine whether copper ions affect membrane fluidity in a anomalous behavior of membrane fluidity, another technique fluorescence correlation spectroscopy (FCS) is applied. FCS operates on the principle that both the diffusion of fluorescent particles in and out of an area and the fluorescence intensity fluctuations over time can be recorded. After calculated by auto-correlation of the fluorescence intensity trace over time and fitted with established correlation function models, the diffusion coefficients of membrane are obtained (Fig. 1c,d). The results show that CuCl 2 treatment does significantly alter the membrane fluidity, contrast to ZnCl 2 and CaCl 2 treatments with a little decrease of the diffusion rate on the bilayer membrane.
These FRAP and FCS data indicate that the interaction of Cu 2+ /Mg 2+ /Ca 2+ /Zn 2+ ions with a PC/PG bilayer suppress the motion of the membrane. This can be attributed to the fact that divalent metal ions bind with the negatively charged lipid PG more stable than Na + (the control) 19 , hence decreasing the electrostatic repulsion of the lipids, resulting in a close packing of lipid molecules. Subsequently, the constraints imposed on the displacement of lipids due to their enhanced order in the presence of metal ions in membranes lead to a reduced rate of lateral diffusion.
Remarkably, an unexpected order of metal-ion impacts on the diffusion rate of lipids in the membranes was observed: Cu 2+ > Zn 2+ ~ Ca 2+ ~ Mg 2+ > Na + (control). This indicates that Cu 2+ plays a distinguished role in suppressing the mobility of lipids, more than the other ions, while a comparable effect occurs with the incubation of Zn 2+ , Ca 2+ and Mg 2+ . This result is obviously inconsistent with the placement in the periodic table of Cu just next to Zn and far from Mg and Ca.
To illustrate the mechanism under the anomalous influence of the Cu 2+ ion on the phospholipid membrane, we have studied the interactions of two Cu ions with two PG molecules. H 3 C-[PO 4 ] − -CH 3 , with a phosphate group as a large portion of its composition, was employed as a simplified model of the phospholipid (PL).
[Cu(H 2 O) 5 ] 2+ was applied to simulate the Cu ion because it was found to coordinate five water molecules in solution 27 .
First, two hydrated Cu 2+ ions interact with two phospholipids, respectively. We called the resultant state "State I" (Fig. 2c), denoted by [PL-Cu(aq)] + , in which the hydrated Cu cation bound the negatively charged oxygen atom in the phosphate group of the lipid. The label aq stands for the water molecules in the hydrated group. Using Eq. 1, the binding energy of a Cu 2+ ion in this state was calculated by an ab initio method based on density functional theory (DFT) with the solvation effect of the outer water environment,  30 . Further analyses showed that the property of Cu-O binding was similar to that in State I (seeing details in Fig. S1 and Table S1 of Supplementary Information). Third, two hydrated Cu 2+ ions simultaneously bind with two phospholipids, resulting in a diCu coupled lipid pair (Fig. 2e). The resulting conformation, referred to State II, had two positive charges (denoted by [PL-diCu(aq)-PL] 2+ ). The optimized Cu-O and Cu-Cu bond lengths in the resulting PL pair were 1.90 Å and 2.58 Å, respectively. The binding energy of a Cu 2+ ion in State II was calculated as follows, 14 kcal/mol in State I. Moreover, the Wiberg bond order was 0.420 for the Cu-Cu binding in State II (Fig. 2b), denoting that a chemical bond with a definite covalent characteristic is formed (more information shown in the following part with Fig. 4). All of these results suggest that the binding of Cu 2+ in State II is more stable than in States I and I′ with the presence of water. It should be noted that State II [PL-diCu-PL] 2+ can be taken as the coupling of two groups [PL-Cu] + (State I) (see Fig. 2c,e). The positive charges of the two groups would cause a Coulomb repulsion and subsequently hinder  the formation of State II. Hence, we performed an electron-structure analysis to reveal the physics underlying the formation of State II. The Cu ion in the structure [PL-diCu-PL] 2+ had an NBO charge of only +0.97 e, indicating that the Cu 2+ ion is reduced to Cu + with an electron configuration of [Ar]3d 10 4s 0 upon binding to the negatively charged phospholipid. Additionally, electrons were observed at the midpoint of two Cu ions with a density of 0.043 e/Å 3 (Fig. 4a left, and Table S2 in Supplementary Information). The pattern of electron localized function (ELF) was approximately square in the region close to the Cu ion (the yellow squares under the yellow balls in the lower-left region of Fig. 4b), suggesting that the outermost electrons majorly occupy the 3d orbital. The ELF value was 0.254 at the midpoint of two Cu ions (Fig. 4b upper-left), indicating that there exists an electron pair localized in this area with a probability of 25.4%. Further NBO analysis showed a chemical bond clearly occurring between the two Cu ions (Fig. 4c left). The bond was composed of 50% of the valance orbitals of each Cu ion, in which the ratio of Cu 3d orbital reached 84.3%. Thus, the 3d electrons of the Cu ions substantially contributed to the Cu-Cu bond. All of these results suggest that a strong 3d 10 -3d 10 attraction 31,32 occurs between the ions after the Cu 2+ is reduced, upon binding to the negatively charged phosphate group of lipid. Therefore, it is this attractive force that suppresses the Coulomb repulsion between two [PL-Cu] + , resulting in the formation of [PL-diCu-PL] 2+ .
For comparison, we have also studied the bindings of other metal ions M (M = Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ and Ag 2+ ) with H 3 C-[PO 4 ] − -CH 3 fragments. It was found that only Ag 2+ presented the same behavior as Cu 2+ upon binding to the phospholipids, while the others did not show this manner. We first calculated the binding energies of hydrated M ions with H 3 C-[PO 4 ] − -CH 3 in States I, I′ and II, respectively. The results are presented in Fig. 3a. For State I, the binding energies were −23.14 kcal/mol for Na + (aq), −22.37 kcal/mole for K + (aq), −33.86 kcal/mol for Mg 2+ (aq), −26.14 kcal/mol for Ca 2+ (aq), −33.29 kcal/mol for Zn 2+ (aq), and −38.94 kcal/mol for Ag 2+ (aq). For State I′, the binding energies were −20.84 kcal/mol for Na + (aq), −17.94 kcal/mole for K + (aq), −22.17 kcal/mol for Mg 2+ (aq), −22.97 kcal/mol for Ca 2+ (aq), −28.98 kcal/mol for Zn 2+ (aq), and −31.46 kcal/mol for Ag 2+ (aq). These results indicate that State I is more stable than State I′, and the binding of divalent ions is more stable than that of monovalent ions. The difference between States I and I′ can be attributed to the solvent effect 30 . For State II, negative values were observed only in the binding energies of the hydrated Ag 2+ ion (−45.10 kcal/ mol) and Zn 2+ ion (−20.14 kcal/mol), and not in the other ions. Moreover, a larger binding strength in State II than in State I occurred with the Ag 2+ ion but not with the Zn 2+ ion. Additionally, the Wiberg bond order of M-M binding was 0.298 for the Ag ions in State II and was less than 0.1 for the others, suggesting that a chemical bond only occurs in a PL-diAg-PL structure with a covalent characteristic, but not in the cases of other metal ions. Therefore, with the competition of State I, State II is stable only for Ag 2+ ions in the presence of water, but not for the other ions. The NBO charge of the Ag 2+ ion was +0.93 in State II, denoting that Ag 2+ is reduced to Ag + upon binding to the phosphate group of lipid, and then, a closed shell [Kr]4d 10 5s 0 of electron configuration occurs in the Ag ion. Therefore, the formation of the Ag + -Ag + chemical bond can be attributed to the 4d 10 -4d 10 attraction, similar to the case of [PL-diCu-PL] 2+ . We thus conclude that the similar behavior of Ag ion and Cu ion in State II results from the structure of valence electrons, which is the same for Group 11 metals. This special property consequently leads to the distinguished difference of these ions with Zn 2+ , Ca 2+ , Mg 2+ , K + and Na + ions upon binding to phospholipids.
Based on the previous DFT calculations, we have improved the CHARMM36 force field specifically for Cu-Cu and Cu-O bindings in PL-diCu-PL structure 33 (detailed information presented in Table S3 and the corresponding discussion of Supplementary Information, Section S4) and performed classical molecular dynamics simulations (MD) to explore the influences of the Cu-containing structure on lipid bilayer. A molecular model, consisting of two PG molecules coupled by two Cu + ions (Fig. 5a), was employed in the MD. The positive charge of the PL-diCu-PL structure was neutralized by Cl − ions in the presence of water molecules. Through a 1.0-μs simulation of each sample, the lipids assembled to form a pattern of stripes (Fig. 5b,c). A segment of assembled PL-diCu-PL structures are shown in Fig. 5d. Two Cl − ions were observed over and under the Cu-O plane, respectively, indicating that Cl − ions are involved in the assembly of the membrane due to an electrostatic attraction. To study the order of the degree of assembled lipids in a membrane, we applied three parameters: the angles of O-O direction in phospholipid along the x (θ) and z (ϕ) axes (upper-left inset of Fig. 5e) for the orders parallel to the membrane surface and the position (d z ) of the phosphorus atom in the membrane along the z-axis for the order perpendicular to the surface. The coordinate origin of the z-axis was set at the midpoint of two layers. First, in the distribution of lipids according to the parameters θ and ϕ with presence of Cu ions (the orange curves in Fig. 5e), two peaks were observed at θ = ±60.5°, and ϕ = ±42.5°. In the distribution according to d z , two peaks were located at approximately ±2.0 nm. The plus/minus symbols above resulted from two opposite O-O directions in the PL-diCu-PL structure (Fig. 5a) for the cases of the angles θ and ϕ and from the two layers of the membrane for d z . In contrast, without the presence of Cu ions, there was no peak for θ, and a broad peak for ϕ (at −90° or +90°, due to a periodic condition), clearly suggesting that the lipids are out of order in the direction parallel to the membrane surface. The full-width half-maximum (FWHM) of the peak was then employed to quantitatively study the order. The FWHM of the θ-related peak was 180° (no peak) without Cu ions and 42° with Cu ions. The reduction from 180° to 42° obviously indicates that Cu ions induce a change in lipids from a disordered state to an ordered state. Similar observations occurred for the other parameters. The FWHM of ϕ-related peak decreased from 63° (without Cu ions) to 25° (with Cu ions), and the FWHM of d z -related peak was reduced from 0.7 nm (without Cu ions) to 0.4 nm (with Cu ions). Moreover, the narrow peaks (FWHM = 25°) at ϕ = ±42.5° with Cu ions further suggested that for the assembled lipids, most of the Cu-O planes were parallel to each other, with an angle of approximately 42.5° (Fig. 5d). We thus conclude that Cu ions can clearly enhance the order of lipids in a membrane, especially in the directions parallel to the membrane surface. Roughness of the membrane surface has also been calculated by the mean square deviation (MSD) of phosphorus-atom displacements along the z-axis of the membrane. The values in the presence and absence of Cu 2+ ions were 0.20 ± 0.02 nm and 0.26 ± 0.04 nm, respectively (Fig. 5f). The significant difference between them obviously suggests that Cu ions suppress the roughness of membrane. Moreover, AFM measurements further supported the theoretical results above. As shown in Fig. 5g,h, a pattern of clear stripes was observed in the CuCl 2 -incubated PC/PG bilayer and not in the control (incubated with NaCl 2 ). The roughness of the membrane was 0.12 ± 0.01 nm for the presence of CuCl 2 and 0.11 ± 0.01 nm in its absence, and the P value between them was 0.015. The AFM data above indicate that Cu ions can significantly suppress the roughness of a lipid bilayer. It is noted that the difference of roughness values between AFM and MD results can be attributed to that a substrate for lipid bilayer was applied in AFM measurements but not in MD simulations.
From both the theoretical and experimental analyses, we conclude that a pattern of stripes can be assembled in the membrane with help of the special Cu-lipid interaction and the Cl ions involved. This significantly enhanced order then imposes a strong constraint on displacement of lipid molecules in the bilayer, clearly reducing the fluidity of membrane, which results in the anomalous effect of CuCl 2 previously observed by FRAP measurements.

Conclusions
In summary, we proposed a motif of diCu coupled phospholipids, which can illuminate the anomalous decrease of membrane fluidity caused by Cu 2+ , compared to those by Zn 2+ , Ca 2+ and Mg 2+ . The mechanism under the motif formation was further revealed. Upon the Cu 2+ ion interacting with the lipid, Cu 2+ was reduced to Cu + . After that, one Cu + ion preferred simultaneous binding to two phospholipids and another Cu + , due to the anomalous 3d 10 -3d 10 attraction between the metal ions. In contrast, this attraction cannot occur in the cases of Zn, Ca and Mg ions due to their electronic structures. It is worth noting that besides lipids, the phosphate group also widely exists in other biological molecules, such as DNA, RNA, ADP, ATP, and enzymes. Therefore, as a kernel of the motif, the structure of diCu coupled phosphate groups and the anomalous Cu-Cu attraction (see Fig. 2d) will provide a new direction for understanding the biological function of copper as a trace element essential to our life, as well as to the mechanisms of copper-related clinical diseases.

Confocal microscope experimental setup and data acquisition.
A commercial confocal microscope (Leica TCS SP5) was used for the fluorescence recovery after photobleaching (FPAR) measurements. A 488-nm Argon laser was used as the excitation source. The sample was illuminated and the fluorescence emission was collected by a 63× oil immersion objective. The dimensions of the acquired regions were 180 μm × 180 μm. The bleaching pulse was applied by rapidly scanning a focused laser beam over an area with a dimension of 8 μm × 8 μm for 30 s with an interval of 1 s at full laser power. Immediately after bleaching, the region of 180 μm × 180 μm was recorded for 300 s with an interval of 5 s at low excitation energy.

Ab initio calculations.
Our ab initio calculations based on the density functional theory (DFT) as well as the electron structure analyses were implemented in the Gaussian09 package 34 . The geometry optimizations and vibrational frequencies of all compounds were carried out at DFT level, employing the M06L functional 35  We performed the MD simulations for the system relaxation in an NPT ensemble at high temperature 320 K for 600 ns. After that, we selected one conformation per 20 ns in the time interval from 500 ns to 600 ns, and obtained five samples as initial structures. Finally, The MD simulations were performed in an NPT ensemble at the temperature of 303 K with 1.0-μs for each sample.
All simulations were performed using GROMACS 5.1 38 with a time step of 2 fs. The CHARMM36 force field for lipids 39,40 and the CHARMM TIP3P water model 41 were used. The particle mesh Ewald (PME) method 42,43 was used to treat long-range electrostatic interactions, whereas the van der Waals interactions were treated with a 1.0 nm-1.2 nm force-based switching function 44 . The temperature was maintained at 303 K using Nose-Hoover thermostat 45,46 with a coupling constant of 1 ps, and the pressure was kept constant at 1 bar using semi-isotropic Parrinello-Rahman barostat 47 with a coupling constant of 5 ps and a compressibility of 4.5 × 10 −5 bar −1 . Periodic boundary conditions were applied in the three directions. After a series of minimization and equilibration steps suggested by CHARMM-GUI 44 , the data were collected every 2 ps during the next 1 μs production run.
Preparation of lipid bilayers for AFM imaging. The negatively charged lipid, 1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol (DMPG), and neutral lipid, 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) were purchased from Avanti Polar Lipids (Alabaster, AL) and used without further purication. Lipid bilayers on freshly cleaved mica surface were prepared following the vesicle fusion method. Briefly, DMPC/DMPG (3:1) mixtures were first dissolved in chloroform, followed by evaporation of the solvent under nitrogen. After that lipid mixtures were dissolved in 50 mM Tris, 150 mM NaCl to a concentration of 1.5 mg/ml and sonicated in a bath sonicator until clear to form small vesicles. A 20 μL droplet of the vesicle solution was then applied to a freshly cleaved fragment of mica, incubated for about 2 h at room temperature, and then the sample was incubated at 35 °C for 40 min to fluidize the lipid, a necessary step to form the bilayer.
AFM imaging. The sample was placed within the AFM (Nano III, Veeco) and imaged in the contact mode using DNP tips (Bruker), with a spring constant of 0.06 N/m. The scan rate was 10 Hz and the applied force was minimized to about 0.1 nN.

Preparation of Lipid Bilayers for FCS experiment.
The supported lipid bilayers were prepared from the negatively charged lipid DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), neutral lipid DMPC (1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine) and tail labeled NBD-PE(1-myristoyl-2-(12-((7-nitro-2-1,3-benzoxadiazol-4-yl)amino)dodecanoyl)-sn-glycero-3-phosphoethanolamine). DMPC, DMPG and NBD-PE solutions in chloroform were mixed to achieve a DMPC/DMPG/NBD-PE molar ratio of 74.9:25:0.1; the solvent was evaporated under nitrogen and the dried lipid film was suspended in pure water to a concentration of 2.5 mM. The lipid suspension was then sonicated to clarity, yielding a suspension of small unilamellar vesicle liposomes. The small unilamellar vesicle suspension was then exposed to a clean glass surface (a microscope petri dishes first was etched by plasma-(Harrick, PDC-32G, 4 minutes in air/vacuum) then was cleaned by %1 hydrofluoric acid and thoroughly rinsed with deionized water and dried under nitrogen) and incubated for 1 hour at room temperature to form lipid bilayers. The excessive unfused liposomes were removed from rinsing with excess of the buffer.
SCIEnTIfIC RePoRTS | (2018) 8:14093 | DOI:10.1038/s41598-018-32322-4 FCS measurements. The FCS experiment was performed on an inverted confocal microscope same as FLIM experiment. The sample was excited at a wave length of 470 nm and an intensity of 30 μW by a supercontinuum fiber laser (Fianium ultrafast fiber laser SC-400-4-PP) equipped with an acousto-optical tunable filter AOTF (Fianium AOTF PX-00027).The excitation light was reflected into a 60x immersion oil objective (Olympus Plan Apo 60X; Cargille Type DF Immersion Oil) by a 488 nm dichroic mirror (Semrock FF488-SDi01). The high speed single-photon counting modules (Becker &Hickl HPM-100-40) equipped with Semrock 550/88 bandpass filter is used for fluorescence emission collection. The resulting photon count was recorded with a counter timer card (Magma Express Card/34) controlled by a TCSPC software suite (Becker & Hickl). All measurements were carried out at room temperature (25 °C).
FCS data fitting. Two models are used for fitting out the diffusion times under different ion conditions. Firstly, the correlation function for a membrane-localized fluorophores is, where N is the average number of molecules in the observed volume. τ and τ d are the correlation and diffusion time, respectively. In order to correct the intersystem crossing effects, such as the photon blinking and bleaching phenomena, the second model containing the triplet decay is also employed for data fitting: where f is the fraction of molecules in triplet state and τ T is the relaxation time for singlet-triplet relaxation.