key: cord-0991710-ermacmf8
authors: Verdiá-Báguena, Carmina; Aguilella, Vicente M.; Queralt-Martín, María; Alcaraz, Antonio
title: Transport mechanisms of SARS-CoV-E viroporin in calcium solutions: Lipid-dependent Anomalous Mole Fraction Effect and regulation of pore conductance
date: 2021-02-20
journal: Biochim Biophys Acta Biomembr
DOI: 10.1016/j.bbamem.2021.183590
sha: 42790cff264269c381ced646050d647d980dde5b
doc_id: 991710
cord_uid: ermacmf8

The envelope protein E of the SARS-CoV coronavirus is an archetype of viroporin. It is a small hydrophobic protein displaying ion channel activity that has proven highly relevant in virus-host interaction and virulence. Ion transport through E channel was shown to alter Ca(2+) homeostasis in the cell and trigger inflammation processes. Here, we study transport properties of the E viroporin in mixed solutions of potassium and calcium chloride that contain a fixed total concentration (mole fraction experiments). The channel is reconstituted in planar membranes of different lipid compositions, including a lipid mixture that mimics the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) membrane where the virus localizes within the cell. We find that the E ion conductance is non-monotonic with the total ionic concentration displaying an Anomalous Mole Fraction Effect (AMFE) only when charged lipids are present in the membrane. We also observe that E channel insertion in ERGIC-mimic membranes – including lipid with intrinsic negative curvature – enhances ion permeation at physiological concentrations of pure CaCl(2) or KCl solutions, with a preferential transport of Ca(2+) in mixed KCl-CaCl(2) solutions. Altogether, our findings demonstrate that the presence of calcium modulates the transport properties of the E channel by interacting preferentially with charged lipids through different mechanisms including direct Coulombic interactions and possibly inducing changes in membrane morphology.

The Covid-19 pandemic has sparked intense research into the structure of the SARS-CoV-2 coronavirus aimed to understand its mechanism of infection and hamper its virulence [1] [2] [3] . The SARS-CoV-2 RNA genome encodes several membrane proteins, one of which, the envelope protein E, is a small hydrophobic protein with ion channel activity localized at the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) . Importantly, the E protein of SARS-CoV-2 shares a high degree of similarity with the severe acute respiratory syndrome coronavirus (SARS-CoV) E protein studied here and their transmembrane domains -responsible for the ion channel activity -have identical amino acid sequences [1, 2, 4] . In addition, recent in silico studies predict analogous ion channel activity of both proteins [1] . Viroporins have been associated to viral replication and pathogenesis [5] . However, only recently the mechanism by which the SARS-CoV-E protein ion activity is linked to pathogenesis was elucidated [6] . The imbalance in ion concentration within cell organelles caused by channel activity trigger an immune response and activates the inflammasome and the subsequent pro-inflammatory response and lung damage. In addition, the E channel activity improves virus production, confers virulence in vivo, and inhibitory mutations in the E amino acid sequence induce compensatory mutations in the virus to restore channel activity [6] . All in all, this evidence makes the E protein a target to find therapies for SARS associated diseases and live attenuated vaccines.

The E viroporin is a 75-amino acid-long protein with a transmembrane domain (TM) of 30 residues (8-38) ( Figure 1A ). The 3D structure of E viroporin from SARS-CoV was solved using solid-state NMR spectroscopy in lipid micelles [7] ( Figure 1B ) and the 3D structure of its homologue in SARS-CoV-2 was also very recently obtained using the same technique in ERGICmimetic lipid bilayers [2] . It was also reported that that SARS-CoV E protein forms a Ca 2+ permeable channel in ERGIC-mimetic lipid bilayers. E viroporin channel activity was found to alter Ca 2+ homeostasis within cells and overstimulate inflammasomes leading to immunopathology [8] . The high calcium gradient concentration between the cytoplasm and ERGIC (100 nM in cytoplasm in contrast with 0.4 mM in ERGIC) points to a major role of calcium in different pathways during the viral cycle. It is known that regulation of intracellular calcium concentrations is essential for several cell functions [6] . The above preliminary study [8] aimed to prove the channel permeability and selectivity to Ca 2+ ions. However, solutions of pure CaCl 2 were used, which is an ionic environment somewhat different from the multiionic one existing within the cell organelles involved in Ca 2+ homeostasis. Here, we specifically address E viroporin permeability to calcium in the presence of another monovalent salt.

Transmembrane domain (TM) sequence of E protein (accession number AYV99820). Blue and red letters represent acidic and basic residues, respectively, located at opposite ends of the TM. Hydrophobic amino acids within the TM are highlighted in yellow. (B) Top and side cartoon representations of E protein NMR pentameric structure of residues 8-65 (PDB code 5X29). Residues 40-65 are excluded in the side view for a better representation of the TM. Each monomer is shown in a different color. (C) Cartoon sketching a putative arrangement of five E protein monomers and lipid molecules to form a proteolipidic toroidal pore. (D) Representative E protein-induced channel current recording (left) and histogram of the current jump amplitudes (right) obtained in 1 M KCl at pH 6 with neutral PC membranes. Current recordings were digitally filtered at 1 kHz using a low-pass Bessel (8-pole) filter. Solid line in the histogram represents a single-Gaussian fitting with mean value and standard deviation as indicated.

In previous investigations we showed that the membrane composition has an important impact on ionic transport through E channels, what lead us to hypothesize that E TMs assemble with lipid molecules to form some sort of proteolipidic pore structures [9] [10] [11] ( Figure  1C ). The two abovementioned 3D structural models for E viroporin [2, 7] agree in a pentameric organization of E TMs stabilized by lipids (see Figure 1B ), with small differences between the two models in the -helices tilting and in the resulting pore size. Rearrangements between different conformations would likely involve lipid-protein interactions and require notable internal flexibility. Having this in mind, the lipid molecules could participate in the pore structure in different forms. In the so-called "barrel stave" pores the lipids are located only in J o u r n a l P r e -p r o o f Journal Pre-proof the pore mouths while the hydrophilic regions of the protein are the ones forming the channel walls across the membrane in parallel to the membrane normal (Alamethicin seems to be the only peptide of this category) [12, 13] . Alternatively, if the length of proteins does not match the hydrophobic length of the bilayer, arrangements in form of "toroidal" pores should appear involving either protein tilting or membrane deformation [12, [14] [15] [16] [17] . This could provoke that lipid head groups may be interspersed between the proteins as shown in Figure 1C .

In this scenario, calcium could modulate the transport properties of the E channel by a number of different mechanisms, such as direct electrostatic interactions with protein residues and lipid molecules, generating hydrophobic effects or inducing membrane deformations, among many others [15, 18] . The affinity of a membrane channel for a particular ionic species can be assessed in different ways, being reversal potential (the potential needed to obtain zero current under a concentration gradient) and subsequent calculations of permeability ratios the most usual by far [19] [20] [21] [22] [23] . Here, we focus on another type of experiment that resembles more the multiionic environment of the E channel activity. It consists in measuring channel conductance in the presence of a mixture of two different species whose individual concentration is changed while the total concentration is kept constant [19, 24] . When conductance changes monotonically between the two endpoints corresponding to each pure electrolyte, it is inferred that the channel has no affinity for any of the considered ions that behave independently [19, 24] . However, in several channels conductance vs mole fraction curves do not change monotonically but show a marked minimum, being the calcium channels particularly prolific in this sense [19, [24] [25] [26] [27] [28] . Such non-monotonic behavior is known as Anomalous Mole Fraction Effect (AMFE) and it is considered a measure of preferential selectivity [24] because the presence of a second ion inhibits the conductance of the first one.

1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC), 1,2-dipalmitoyl-sn-glycero-3-phospho-L-serine (DPhPS), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2-dioleoyl-sn-glycero-3phosphoethanolamine (DOPE), and 1,2-dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) were purchased from Avanti Polar Lipids (Alabaster, AL). Full-length SARS-CoV E protein was kindly provided by Dr. Jaume Torres and synthesized and purified as previously described [9] .

Planar bilayers were formed by apposition of two monolayers prepared from a solution 5 mg/ml in pentane of pure DPhPC, pure DPhPS, or a mixture of DOPC, DOPS and DOPE at different ratios as indicated in the text. Lipids were added on 100 μm diameter orifices in the 15 μm-thick Teflon partition that separated two identical chambers [29, 30] . The orifices were pretreated with a 1% solution of hexadecane in pentane. Aqueous solutions of KCl were buffered with 5 mM HEPES and the pH was fixed at pH = 6 except otherwise noticed. All measurements were performed at room temperature (23 ± 1 °C). Ion channel insertion was achieved by adding 0.5-1 μl of a 300 μg/ml solution of synthetic protein in acetonitrile:isopropanol (40:60) to one side of the chamber (cis side). An electric potential was applied using Ag/AgCl electrodes in 2 M KCl, 1.5% agarose bridges assembled within standard 250 μl pipette tips. The potential was defined as positive when it was higher on the side of protein addition (cis side), whereas the trans side was set to ground. An Axopatch 200B amplifier (Molecular Devices, Sunnyvale, CA) in the voltage-clamp mode was used to measure the current and the applied potential. Data were filtered by an integrated low pass 8-pole Bessel filter at 10 kHz, saved with a sampling frequency of 50 kHz and analyzed using pClamp 10.7 software (Molecular Devices, Sunnyvale, CA). The chamber and the head stage were isolated from external noise sources with a double metal screen (Amuneal Manufacturing Corp., Philadelphia, PA).

The channel conductance was obtained from current measurements under an applied potential of +100 mV in symmetrical salt solutions. Current was calculated from a single-Gaussian fitting of the histogram of current jump amplitudes generated after addition of SARS-CoV E protein to the chamber. The standard deviation of the data corresponds to the sigma obtained from the Gaussian fit. Histograms contained at least 300 events from a minimum of 3 independent experiments. The reversal potential (the voltage corresponding to zero current) was measured on one or several channels under a salt concentration gradient of 500/50 mM KCl with or without a symmetrical addition of 15 mM CaCl 2 . The measured value was corrected by the liquid junction potential from Henderson's equation [31] to obtain the final reversal potential. Permeability ratios between K + and Cl − , P + /P -, where calculated according to the Goldman-Hodgkin-Katz equation [19] . In experiments with CaCl 2 , the presence of 15 mM CaCl 2 at both sides of the membrane was taken into account to calculate the permeability ratios.

Here, we investigate the interaction between calcium and E-induced pores by analyzing elementary current jumps in different mixtures of KCl and CaCl 2 (as shown in Figure 1D , left panel). Conductance values (G) correspond to the mean value obtained from a single-Gaussian fit of the current jump amplitudes (see Figure 1D ) divided over the applied potential (G = I/V). Figure 2A shows measurements performed in negatively charged PS membranes in 50 mM solutions. The curve in Figure 2A shows a minimum at 0.15 Ca 2+ mole fraction, which corresponds to 7 mM CaCl 2 and to a 43% reduction in the K + current relative to a pure KCl solution. As far as we know, this effect has not been described up to now in viroporins. Channel conductance is measured at +100 mV applied potential. Here and elsewhere, each data point is calculated from a single-Gaussian fitting of the histogram of current jump amplitudes as shown in Figure 1D .

To probe whether membrane charge has an impact on the observed AMFE in Figure 2A , we carried out similar experiments in neutral PC membranes. We found in this case that channel conductance seems to decrease monotonically with Ca 2+ mole fraction ( Figure 2B ). Even though data error bars are too large to be categorical, the trend shows no conclusive evidence on AMFE unlike that observed in Figure 2A . The conventional explanation of AMFE points to the correlated movement of ions in a single file fashion [19] characteristic of very narrow channels like gramicidin A [32] . However, E channel unitary conductance is similar to another multiionic proteolipidic pore, Alamethicin, in its lowest conductance state [9, [33] [34] [35] , which implies that it allows the transport of fully hydrated ions, water molecules and small molecules [22] . Since single-file transport of dehydrated ions seems highly improbable in multi-ionic pores such as synthetic nanopores where AMFE has already been reported [25] , it has been suggested that AMFE may appear also because of localized, ion-specific binding within the pore [25] . According to this explanation, in the presence of different permeant electrolytes, one of the ionic species binds tighter within the channel, excluding the other from the adjacent region. This creates a depletion zone that results in a decrease of channel conductance, given that the conductance for each ion is mostly determined by the region where it is present in the lowest concentration [36] . Thus, the reported AMFE (discernible in Figure 2A and Figure Figure 2B and Figure 2D ) points to a preferential interaction of Ca 2+ over K + ions in localized sites within the E proteolipidic channel.

This observation of AMFE only when the E viroporin is reconstituted in negatively charged membranes hints at the fact that Ca 2+ ions interact preferably with lipid charges rather than with the acidic residues of the E protein TM. One distinctive trait of calcium-membrane interactions reported in GUVs is that increasing concentration of monovalent cations affects the binding affinity of Ca 2+ [15] . To explore this effect in the E viroporin, we performed AMFE experiments in concentrated solutions. Figures 2C and 2D display the change in conductance vs Ca 2+ mole fraction when the total salt concentration is 1 M. Similar to what is shown in Figure 2A -B for 50 mM solutions, AMFE is observed in PS membranes ( Figure 2C ), but not in neutral ones ( Figure 2D ) in 1 M solutions. However, the lowest conductance appears at higher Ca 2+ mole fraction than in the respective low-concentration experiment (0.3 in Figure 2C versus 0.15 in Figure 2A ) and the reduction in K + current is also smaller (26% in Figure 2C vs 43% in Figure 2A ). Thus, in agreement with experiments in GUVS, high salt concentration yields screening effects that reduce the binding affinity of Ca 2+ to negatively charged phospholipids [15] .

An alternative route to probe the role of membrane charge in the interaction of Ca 2+ ions with the pore is to assess the effect of solution acidification on AMFE experiments. pH acidification causes a charge neutralization of both the lipid and protein residues within the pore that is reflected in the channel ionic selectivity, which changes from cationic at neutral pH to anionic at acidic pH [8, 10] . This effect is quantified in Figure 3A by means of the permeability ratio P K+ /P Cl-in E pores as a function of pH for neutral (PC) and charged membranes (PS) in KCl solutions, and for charged membranes in KCl+CaCl 2 solutions (PS + CaCl 2 ). It is clearly seen that under very acidic conditions (pH 1.5) the channel selectivity in KCl is very similar in neutral and charged membranes (P K+ /P Cl-= 0.43 ± 0.06 and P K+ /P Cl-= 0.54 ± 0.07, respectively). Hence, at this acidic pH the negative charges of the PS lipid are titrated by protons, so the lipid barely influences channel selectivity. To test the existence of AMFE in these conditions, we performed channel conductance experiments for varying mixtures of CaCl 2 and KCl at pH 1.5 (total salt concentration of 50 mM) ( Figure 3B ). The monotonic change of conductance with Ca 2+ molar fraction demonstrates that there is no AMFE when membrane lipid charge is neutralized, confirming that lipid charges are largely responsible for the calcium binding within this proteolipidic channel. Furthermore, the interaction between Ca 2+ ions and negative lipid charges explains why addition of millimolar CaCl 2 in selectivity experiments with PS membranes ( Figure 3A , grey squares) reduces the channel preference for cations, which becomes similar to that obtained with PC membranes. Interestingly, measurements of P K+ /P Clin PS with and without Ca 2+ ions ( Figure 3A , grey squares and blue triangles) almost overlap at pH 1.5 (P K+ /P Cl-= 0.55 ± 0.07 and P K+ /P Cl-= 0.54 ± 0.07, respectively), implying that the interaction between lipid charges and Ca 2+ vanishes once the lipid is titrated. Still, an interaction of Ca 2+ ions with the E protein residues cannot be completely excluded, given that there is a small pKa shift in the titration experiments between neutral PC and neutralized PS + CaCl 2 membranes.

Journal Pre-proof Next, we consider the more physiologically relevant lipid mixture PC/PE/PS 3:1:1 (w/w), which mimics the phospholipid headgroup composition of ERGIC membranes, the environment where E protein is localized [8] . Figure 4A shows conductance vs Ca 2+ mole fraction in ERGIClike membranes (labelled as ERGIC) in 50 mM solutions of KCl+CaCl 2 . Interestingly, AMFE appears in ERGIC membranes with a minimum at 0.4 mole fraction (corresponding to 20 mM Ca 2+ ) and a 44% reduction of K + current. To facilitate direct comparison, curves are shown together with measurements in PC and PS under the same conditions (50 mM). At each mole fraction, conductance is lower in neutral PC membranes than in charged PS or ERGIC membranes, even though conductance in PC is monotonic and it is not in PS and ERGIC. Of note, for solutions with predominant amount of KCl or CaCl 2 (mole fractions close to 0 and 1, respectively) the conductance in ERGIC membranes (with 20% charged lipid) is higher than in PS. Thus, an interesting observation is that in the native-like lipid composition E pores present an optimal combination of curvature and charge that yields maximum ion conduction in both pure KCl and CaCl 2 . This implies that the intrinsic negative curvature of PE lipid (present in the ERGIC mixture with a 20% percentage) may be decisive to optimize the lipid packing, either inducing a membrane thinning or modifying the hydrophobic mismatch between the protein TM and the lipids [37] . Similar experiments performed at high electrolyte concentrations give a different scenario. Figure 4B shows G vs Ca 2+ mole fraction for PC, PS and ERGIC membranes at a total mixture concentration of 1 M. AMFE does exist in ERGIC membranes in these conditions, but it is barely visible. Importantly, the curve with lower conductance values is that of the fully charged membrane PS, being the neutral PC the one with higher measured conductance for all Ca 2+ mole fractions. ERGIC (20% of PS) conductance is now in-between PC and PS, as could be expected if only the percentage of charged lipid is considered.

To examine in more detail the effect of lipid charge in ERGIC membranes, we performed channel conductance experiments in lipid mixtures with a fixed ratio of PE (20% w/w, as in ERGIC membranes) and a varying ratio of PC and PS ( Figure 4B ). In line with the results in Figure 4A , opposing trends are observed at low and high concentrations in Figure 4B . At 100 mM CaCl 2 , conductance increases with the relative amount of PS, which could be explained by the enhanced accumulation of ions available for conduction due to the increase of charge. However, in 1 M CaCl 2 channel conductance decreases with the increasing proportion of charged lipids, suggesting that the binding of cations to the lipid polar heads is so intense that it has an impact on ion mobility. The stronger the interaction, the slower the movement of cations through the channel [19] . Next, we examine how monovalent and divalent cations interact with E pores by focusing on pure electrolytes, either KCl or CaCl 2 . Figure 5A summarizes conductance experiments performed in 50 mM and 1 M solutions of KCl or CaCl 2 , including now additional measurements in 1:1 PC:PS membranes. At low salt concentration, channel conductance in PC:PS (50% of charged lipid) falls in between pure PS and ERGIC (20% of charged lipid), demonstrating that the increased ion conduction observed in ERGIC is due to curvature effects (presence of PE) and not to a non-linear combination of PC and PS. Interestingly, at high concentration the channel conductance in ERGIC is halfway between PC and PC:PS 1:1, i.e. changing monotonically with the charged lipid content of the membrane and proving that curvature effects are no longer decisive. Figure 5A also shows that E channel conductance in KCl solutions is in general higher than in CaCl 2 , which seems counterintuitive given that bulk solutions of CaCl 2 are more conductive J o u r n a l P r e -p r o o f than KCl ones (both for 50 mM and 1 M solutions). To explore this effect, we calculated the ratio between channel conductance in CaCl 2 and KCl for the same conditions of salt concentration and membrane composition. Results are presented in Figure 5B , including the corresponding ratio between bulk conductivities as well. Figure 5B shows that the ratio between CaCl 2 and KCl conductance is lower than that of bulk conductivities for all membrane compositions tested. This implies that the interaction with the E proteolipidic pore is in general stronger for Ca 2+ than for K + ions, so that Ca 2+ effective mobility is further reduced compared to bulk [10, 22] . In addition, Figure 5B demonstrates that the less charged the membrane, the lower the ratio between CaCl 2 and KCl conductance (except for ERGIC at 50 mM salt concentration). This corresponds to a larger difference in binding strength between K + and Ca 2+ in PC than in PS membranes that causes a relatively higher decrease of Ca 2+ than K + ion mobility in PC, compared to PS. This seems reasonable because, while K + ions interact very weakly with PC membranes, Ca 2+ ions still bind to this neutral lipid [18, 38] , although to a lesser extent, given that no AMFE is observed in PC ( Figure 2 ). As commented before, the conductance ratio in ERGIC membranes at low salt concentration depicted in Figure 5B is the only one that deviates from the general trend, pointing out again to the paramount importance of membrane curvature at physiologically low concentrations, where the CaCl 2 /KCl conductance ratio approaches that of bulk conductivities. To rationalize our findings, we may consider some general features of the interaction between cations and phospholipid membranes. The binding of Ca 2+ to lipid membranes has been studied thoroughly [15, 18, 39, 40] , with a large scatter of measured calcium binding constants that may actually reflect a variety of different phenomena such as condensation effect of Ca 2+ ions, phase separation of lipids or shape remodeling of membranes [15, 18] . Ca 2+ is thought to interact with lipid membranes spontaneously in an entropy-driven process [18] and occurs in a comparable fashion for negatively charged lipids and for zwitterionic lipids. The higher binding constants for negative lipids are potentially caused by the electrostatic enrichment of Ca 2+ ions in the vicinity of the binding sites [18] and explains why we observe AMFE in charged PS but not in neutral PC (Figure 2 ). Regarding the interaction between monovalent cations and lipid membranes, literature data also show considerable dispersion, but binding constants are smaller than for Ca 2+ [38, 39, 41] . The interaction of monovalent cations with charged membranes is usually slightly stronger than with neutral ones, but the difference is not as marked as for divalent cations [41] .

The most direct consequence of the binding process could be that bound cations might alter the dimensions of the pore eyelet via mutual repulsion or by other electrostatic interactions with the protein TM residues. However, more subtle mechanisms could operate as well. As mentioned before, notable changes in membrane morphology have been reported due to presence of divalent Ca 2+ but also for monovalent cations like K + or Na + [15] . Interestingly, those changes seem to operate in opposite directions. On the one side, experiments in GUVs show that interactions of Ca 2+ with charged lipids induce negative curvature via membrane invaginations resulting in a decrease of the membrane area [15] . In contrast, similar experiments in presence of Na + suggest that these cations could yield positive curvature [15] in agreement with simulations showing that the membrane curves away from monovalent cations to increase the separation between them [42] . Although it is unclear how those findings in GUVs could be directly used to interpret experiments in planar membranes, it is indisputable that the presence of cations near the lipid headgroups provokes different adaptation mechanisms [37] . The first one implies that bound cations could induce a change in the hydrophobic thickness of proteins and lipids. The second is based on membrane stretching that cause changes in lateral pressure profiles [37] . Hydrophobic interactions may be also an additional player to understand some of the counterintuitive findings relative to KCl and CaCl 2 conductance. According to the latest atomic structure of the SARS-CoV-2 E viroporin [2] , the narrowest section of the E pore (2 Å radius) is flanked by four layers of hydrophobic residues that would force ion dehydration to pass through. Given that Ca 2+ hydration enthalpy is several times greater than K + one, this would be consistent with the higher conductance of E channel in KCl solutions than in CaCl 2 in low concentration solutions ( Figure 5B ).

The hitherto accepted view of the E viroporin as a slightly cation selective multiionic pore with five transmembrane helices stabilized by the lipid and poor specific selectivity towards any particular ionic species seems insufficient to explain its calcium permeation properties in J o u r n a l P r e -p r o o f solutions where there are other monovalent ions. We have demonstrated that the E channel conductance in solutions of CaCl 2 +KCl does not change monotonically with the mole fraction of each salt but exhibits a minimum, an effect known in the literature as AMFE and characteristic of several calcium channels.

Our experiments suggest that the ionic transport properties of the SARS-CoV-E viroporin in solutions containing both monovalent and divalent cations are strongly dependent on the membrane lipid composition, the ratio between CaCl 2 and KCl concentrations, as well as on the total salt concentration. The influence of the membrane lipid composition has two aspects, probably intertwined: on the one hand the interaction of cations with negatively charged lipid polar heads; on the other hand, the effect of lipids with intrinsic negative curvature as DOPE on the conformation of the five E transmembrane helices. This is consistent with the reported effect of membrane environment on the protein helical conformation and the established proteolipidic character of the SARS-CoV-E viroporin.

Both monovalent and divalent cations interact strongly with negatively charged lipids yielding two main effects: an accumulation of ions due to coulombic screening that enhances channel conductance at low concentration, and specific binding that becomes dominant at high concentration and induces relatively slower permeation. These two effects are clearly visible when comparing conductance in charged and neutral lipids: at low salt concentrations, conductance G(PC) < G(PS) because there is hardly any ion accumulation in neutral PC membranes. At high electrolyte concentration, G(PC) > G(PS) because the interaction between cations and neutral lipids is milder than with charged PS and has a more limited impact on ion permeation. As regards the comparison between Ca 2+ and K + binding, we may conjecture that divalent cations have a higher affinity for lipids than monovalent ones. We believe this is the cause of the reported AMFE, which occurs at Ca 2+ mole fractions lower than 0.5, consequently causing a larger difference between channel conductance and bulk conductivity as salt concentration increases. E channels reconstituted in ERGIC-mimic membranes display effects that cannot be easily attributed to the small lipid charge but are probably linked to a different protein-lipid arrangement induced by DOPE. At physiological CaCl 2 or KCl concentrations, these membranes enhance ion permeation in pure salts and optimize Ca 2+ over K + transport in mixed KCl+CaCl 2 solutions. It is remarkable that these effects arise just in lipid bilayers resembling the natural environment of the E protein channel activity.

Altogether, this work provides novel insights on SARS-CoV E protein function and demonstrates that the presence of calcium modulates the transport properties of the E channel by interacting preferentially with charged lipids within the pore through different interaction mechanisms including direct Coulombic interactions and possibly changes in membrane morphology.

Computational Study of the Ion and Water Permeation and Transport Mechanisms of the SARS-CoV-2 Pentameric E Protein Channel

Structure and drug binding of the SARS-CoV-2 envelope protein transmembrane domain in lipid bilayers

The deadly coronaviruses: The 2003 SARS pandemic and the 2020 novel coronavirus epidemic in China

Structural insight into the role of novel SARSCoV-2 E protein: A potential target for vaccine development and other therapeutic strategies

Relevance of Viroporin Ion Channel Activity on Viral Replication and Pathogenesis., Viruses

Severe Acute Respiratory Syndrome Coronavirus Envelope Protein Ion Channel Activity Promotes Virus Fitness and Pathogenesis

Structural model of the SARS coronavirus E channel in LMPG micelles

Severe acute respiratory syndrome coronavirus E protein transports calcium ions and activates the NLRP3 inflammasome

Coronavirus E protein forms ion channels with functionally and structurallyinvolved membrane lipids

Analysis of SARS-CoV E protein ion channel activity by tuning the protein and lipid charge

Lipid charge regulation of non-specific biological ion channels

Fluctuations of Ionic Current Through Lipid Bilayers at the Onset of Peptide Attacks and Pore Formation

Extending the Hydrophobic Mismatch Concept to Amphiphilic Membranolytic Peptides

Hydrophobic mismatch demonstrated for membranolytic peptides, and their use as molecular rulers to measure bilayer thickness in native cells

Cations induce shape remodeling of negatively charged phospholipid membranes

Single-molecule conformational dynamics of viroporin ion channels regulated by lipid-protein interactions

Ion channel activity of the CSFV p7 viroporin in surrogates of the ER lipid bilayer

Binding of calcium to phosphatidylcholinephosphatidylserine membranes

Ion Channels of Excitable Membranes

Theoretical description of the ion transport across nanopores with titratable fixed charges: analogies between ion channels and synthetic pores

Critical assessment of OmpF channel selectivity: merging information from different experimental protocols

Insights on the permeability of wide protein channels: measurement and interpretation of ion selectivity

Physical descriptions of experimental selectivity measurements in ion channels

The anomalous mole fraction effect in calcium channels: A J o u r n a l P r e -p r o o f measure of preferential selectivity

Synthetic nanopores as a test case for ion channel theories: the anomalous mole fraction effect without single filing

Non-selective conductance in calcium channels of frog muscle: calcium selectivity in a single-file pore

A non-selective cation conductance in frog muscle membrane blocked by micromolar external calcium ions

Energetics of divalent selectivity in a calcium channel: The ryanodine receptor case study

Probing alamethicin channels with water-soluble polymers. Effect on conductance of channel states

Formation of bimolecular membranes from lipid monolayers and a study of their electrical properties

Diffusion, exclusion, and specific binding in a large channel: a study of OmpF selectivity inversion

The gramicidin ion channel: A model membrane protein

Scaling Behavior of Ionic Transport in Membrane Nanochannels

Alamethicin channel conductance modified by lipid charge

Channel-forming activity of alamethicin: Effects of covalent tethering

Anomalous mole fraction effect, electrostatics, and binding in ionic channels

Regulation of membrane calcium transport proteins by the surrounding lipid environment

Effect of NaCl and KCl on Phosphatidylcholine and Phosphatidylethanolamine Lipid Membranes: Insight from Atomic-Scale Simulations for Understanding Salt-Induced Effects in the Plasma Membrane

Ca2+, Mg2+, Li+, Na+, and K+ distributions in the headgroup region of binary membranes of phosphatidylcholine and phosphatidylserine as seen by deuterium NMR

The complex nature of calcium cation interactions with phospholipid bilayers

Binding of monovalent alkali metal ions with negatively charged phospholipid membranes

Spontaneous curvature of bilayer membranes from molecular simulations: Asymmetric lipid densities and asymmetric adsorption

Authors acknowledge financial support by the Spanish Agencia Estatal de Investigación (Projects ID2019-108434GB-I00/AEI and IJC2018-035283-I/AEI), Generalitat Valenciana (Project AICO/2020/066), and Universitat Jaume I (Projects UJI-B2018-53 and UJI-A2020-21). 

☐ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.