Structural insights into a cooperative switch between one and two FimH bacterial adhesins binding pauci- and high-mannose type N-glycan receptors.

The FimH type-1 fimbrial adhesin allows pathogenic Escherichia coli to adhere to glycoproteins in the epithelial linings of human bladder and intestinal tract, by using multiple fimbriae simultaneously. Pauci- and high-mannose type N-glycans are natural FimH receptors on those glycoproteins. Oligomannose-3 and oligomannose-5 bind with the highest affinity to FimH by using the same Manα1,3Man branch. Oligomannose-6 is generated from oligomannose-5 in the next step of the biogenesis of high-mannose N-glycans, by the transfer of a mannose in α1,2-linkage onto this branch. Using serial crystallography and by measuring the kinetics of binding, we demonstrate that shielding the high-affinity epitope drives the binding of multiple FimH molecules. First, we profiled FimH glycan binding on a microarray containing paucimannosidic N-glycans and in a FimH LEctPROFILE assay. To make the transition to oligomannose-6, we measured the kinetics of FimH binding using paucimannosidic N-glycans, glycoproteins and all four α-dimannosides conjugated to bovine serum albumin. Equimolar mixed interfaces of the dimannosides present in oligomannose-6 and molecular dynamics simulations suggest a positive cooperativity in the bivalent binding of Manα1,3Manα1 and Manα1,6Manα1 dimannosides. The binding of core α1,6-fucosylated oligomannose-3 in cocrystals of FimH is monovalent but interestingly the GlcNAc1-Fuc moiety retains highly flexibility. In cocrystals with oligomannose-6, two FimH bacterial adhesins bind the Manα1,3Manα1 and Manα1,6Manα1 endings of the second trimannose core (A-4'-B). This cooperative switch towards bivalent binding appears sustainable beyond a molar excess of oligomannose-6. Our findings provide important novel structural insights for the design of multivalent FimH antagonists that bind with positive cooperativity.

Polyvalent Transition-State Analogues of Sialyl Substrates Strongly Inhibit Bacterial Sialidases*.

Bacterial sialidases (SA) are validated drug targets expressed by common human pathogens such as Streptococcus pneumoniae, Vibrio cholerae, or Clostridium perfringens. Noncovalent inhibitors of bacterial SA capable of reaching the submicromolar level are rarely reported. In this work, multi- and polyvalent compounds are developed, based on the transition-state analogue 2-deoxy-2,3-didehydro-N-acetylneuraminic (DANA). Poly-DANA inhibits the catalytic activity of SA from S. pneumoniae (NanA) and the symbiotic microorganism B. thetaiotaomicron (BtSA) at the picomolar and low nanomolar levels (expressed in moles of molecules and of DANA, respectively). Each DANA grafted to the polymer surpasses the inhibitory potential of the monovalent analogue by more than four orders of magnitude, which represents the highest multivalent effect reported so far for an enzyme inhibition. The synergistic interaction is shown to operate exclusively in the catalytic domain, and not in the flanked carbohydrate-binding module (CBM). These results offer interesting perspectives for the multivalent inhibition of other SA families lacking a CBM, such as viral, parasitic, or human SA.

Function and Regulation of Acid Resistance Antiporters.

Bacterial pathogens are a major cause of foodborne diseases and food poisoning. To cope with the acid conditions encountered in different environments such as in fermented food or in the gastric compartment, neutralophilic bacteria have developed several adaptive mechanisms. One of those mechanisms, the amino acid dependent system, consumes intracellular protons in biochemical reactions. It involves an antiporter that facilitates the exchange of external substrate amino acid for internal product and a cytoplasmic decarboxylase that catalyzes a proton-consuming decarboxylation of the substrate. So far, four acid resistance antiporters have been discovered, namely the glutamate-γ-aminobutyric acid antiporter GadC, the arginine-agmatine antiporter AdiC, the lysine-cadaverine antiporter CadB, and the ornithine-putrescine antiporter PotE. The 3D structures of AdiC and GadC, reveal an inverted-repeat fold of two times 5 transmembrane helices, typical of the amino acid-polyamine-organocation (APC) superfamily of transporters. This review summarizes our current knowledge on the transport mechanism, the pH regulation and the selectivity of these four acid resistance antiporters. It also highlights that AdiC is a paradigm for eukaryotic amino acid transporters of the APC superfamily as structural models of several of these transporters built using AdiC structures were exploited to unveil their mechanisms of amino acid recognition and translocation.

Monitoring glycosidase activity for clustered sugar substrates, a study on ß-glucuronidase.

Determination of glycosidase hydrolysis kinetics for a monovalent sugar substrate is relatively straightforward and classically achieved by monitoring the fluorescence signal released from the sugar-conjugated probe after enzymatic hydrolysis. Naturally occuring sugar epitopes are, however, often clustered on biopolymers or at biological surfaces, and previous reports have shown that glycosidase hydrolytic rates can differ greatly with multivalent presentation of the sugar epitopes. New probes are needed to make it easier to interpret the importance of substrate clustering towards a specific enzyme activity. In this work, we developed multivalent glucuronide substrates attached to fluorescent amino-coumarines through self-immolative linkers to enable real time-monitoring of the hydrolysing activity of E.coli ß-glucuronidases (GUS) towards clustered substrates. GUS are exoglycosidases of considerable therapeutic interest cleaving ß-d-glucuronides and are found in the lysosomes, in the tumoral microenvironment, and are expressed by gut microbiota. GUS showed a much lower catalytic efficiency in hydrolysing clustered glucuronides due to a significantly lower enzymatic velocity and affinity for the substrates. GUS was 52-fold less efficient in hydrolysing GlcA substrates presented on an octameric silsequioxane (COSS) compared with a monovalent GlcA of similar chemical structure. Thus, kinetic and thermodynamic data of GUS hydrolysis towards multivalent glucuronides were easily obtained with these new types of enzymatically-triggered probes. More generally, adapting the substrate nature and valency of these new probes, should improve understanding of the impact of multivalency for a specific enzyme.

Molecular mechanism of substrate selectivity of the arginine-agmatine Antiporter AdiC.

The arginine-agmatine antiporter (AdiC) is a component of an acid resistance system developed by enteric bacteria to resist gastric acidity. In order to avoid neutral proton antiport, the monovalent form of arginine, about as abundant as its divalent form under acidic conditions, should be selectively bound by AdiC for transport into the cytosol. In this study, we shed light on the mechanism through which AdiC distinguishes Arg+ from Arg2+ of arginine by investigating the binding of both forms in addition to that of divalent agmatine, using a combination of molecular dynamics simulations with molecular and quantum mechanics calculations. We show that AdiC indeed preferentially binds Arg+. The weaker binding of divalent compounds results mostly from their greater tendency to remain hydrated than Arg+. Our data suggests that the binding of Arg+ promotes the deprotonation of Glu208, a gating residue, which in turn reinforces its interactions with AdiC, leading to longer residence times of Arg+ in the binding site. Although the total electric charge of the ligand appears to be the determinant factor in the discrimination process, two local interactions formed with Trp293, another gating residue of the binding site, also contribute to the selection mechanism: a cation-π interaction with the guanidinium group of Arg+ and an anion-π interaction involving Glu208.

A Novel Integrated Way for Deciphering the Glycan Code for the FimH Lectin.

The fimbrial lectin FimH from uro- and enteropathogenic Escherichia coli binds with nanomolar affinity to oligomannose glycans exposing Manα1,3Man dimannosides at their non-reducing end, but only with micromolar affinities to Manα1,2Man dimannosides. These two dimannoses play a significantly distinct role in infection by E. coli. Manα1,2Man has been described early on as shielding the (Manα1,3Man) glycan that is more relevant to strong bacterial adhesion and invasion. We quantified the binding of the two dimannoses (Manα1,2Man and Manα1,3Man to FimH using ELLSA and isothermal microcalorimetry and calculated probabilities of binding modes using molecular dynamics simulations. Our experimentally and computationally determined binding energies confirm a higher affinity of FimH towards the dimannose Manα1,3Man. Manα1,2Man displays a much lower binding enthalpy combined with a high entropic gain. Most remarkably, our molecular dynamics simulations indicate that Manα1,2Man cannot easily take its major conformer from water into the FimH binding site and that FimH is interacting with two very different conformers of Manα1,2Man that occupy 42% and 28% respectively of conformational space. The finding that Manα1,2Man binding to FimH is unstable agrees with the earlier suggestion that E. coli may use the Manα1,2Man epitope for transient tethering along cell surfaces in order to enhance dispersion of the infection.

Targeting Dynamical Binding Processes in the Design of Non-Antibiotic Anti-Adhesives by Molecular Simulation-The Example of FimH.

Located at the tip of type I fimbria of Escherichia coli, the bacterial adhesin FimH is responsible for the attachment of the bacteria to the (human) host by specifically binding to highly-mannosylated glycoproteins located on the exterior of the host cell wall. Adhesion represents a necessary early step in bacterial infection and specific inhibition of this process represents a valuable alternative pathway to antibiotic treatments, as such anti-adhesive drugs are non-intrusive and are therefore unlikely to induce bacterial resistance. The currently available anti-adhesives with the highest affinities for FimH still feature affinities in the nanomolar range. A prerequisite to develop higher-affinity FimH inhibitors is a molecular understanding of the FimH-inhibitor complex formation. The latest insights in the formation process are achieved by combining several molecular simulation and traditional experimental techniques. This review summarizes how molecular simulation contributed to the current knowledge of the molecular function of FimH and the importance of dynamics in the inhibitor binding process, and highlights the importance of the incorporation of dynamical aspects in (future) drug-design studies.

Norovirus devours human milk oligosaccharides rich in α-fucose.

Human norovirus binding to histo-blood group antigens (HBGAs) is thought to direct their entry into host cells. However, the glycan epitopes characteristic of HBGAs are also present on oligosaccharides abundant in human milk. In this issue of JBC, Hanisch et al compared norovirus binding to human gastric mucins and human milk oligosaccharides, finding those bound most avidly are rich in α-fucose. Mimicry of these epitopes with α-fucose multivalently displayed on other carbohydrate scaffolds successfully scavenged this prevalent virus, providing new insights into norovirus biology and clues for future therapeutic development.

Lipidated apolipoprotein E4 structure and its receptor binding mechanism determined by a combined cross-linking coupled to mass spectrometry and molecular dynamics approach.

Apolipoprotein E (apoE) is a forefront actor in the transport of lipids and the maintenance of cholesterol homeostasis, and is also strongly implicated in Alzheimer's disease. Upon lipid-binding apoE adopts a conformational state that mediates the receptor-induced internalization of lipoproteins. Due to its inherent structural dynamics and the presence of lipids, the structure of the biologically active apoE remains so far poorly described. To address this issue, we developed an innovative hybrid method combining experimental data with molecular modeling and dynamics to generate comprehensive models of the lipidated apoE4 isoform. Chemical cross-linking combined with mass spectrometry provided distance restraints, characterizing the three-dimensional organization of apoE4 molecules at the surface of lipidic nanoparticles. The ensemble of spatial restraints was then rationalized in an original molecular modeling approach to generate monomeric models of apoE4 that advocated the existence of two alternative conformations. These two models point towards an activation mechanism of apoE4 relying on a regulation of the accessibility of its receptor binding region. Further, molecular dynamics simulations of the dimerized and lipidated apoE4 monomeric conformations revealed an elongation of the apoE N-terminal domain, whereby helix 4 is rearranged, together with Arg172, into a proper orientation essential for lipoprotein receptor association. Overall, our results show how apoE4 adapts its conformation for the recognition of the low density lipoprotein receptor and we propose a novel mechanism of activation for apoE4 that is based on accessibility and remodeling of the receptor binding region.

Transition of yeast Can1 transporter to the inward-facing state unveils an α-arrestin target sequence promoting its ubiquitylation and endocytosis.

Substrate-transport-elicited endocytosis is a common control mechanism of membrane transporters avoiding excess uptake of external compounds, though poorly understood at the molecular level. In yeast, endocytosis of transporters is triggered by their ubiquitylation mediated by the Rsp5 ubiquitin-ligase, recruited by α-arrestin-family adaptors. We here report that transport-elicited ubiquitylation of the arginine transporter Can1 is promoted by transition to an inward-facing state. This conformational change unveils a region of the N-terminal cytosolic tail targeted by the Art1 α-arrestin, which is activated via the TORC1 kinase complex upon arginine uptake. Can1 mutants altered in the arginine-binding site or a cytosolic tripeptide sequence permanently expose the α-arrestin-targeted region so that Art1 activation via TORC1 is sufficient to trigger their endocytosis. We also provide evidence that substrate-transport elicited endocytosis of other amino acid permeases similarly involves unmasking of a cytosolic Art1-target region coupled to activation of Art1 via TORC1. Our results unravel a mechanism likely involved in regulation of many other transporters by their own substrates. They also support the emerging view that transporter ubiquitylation relies on combinatorial interaction rules such that α-arrestins, stimulated via signaling cascades or in their basal state, recognize transporter regions permanently facing the cytosol or unveiled during transport.

Sites for Dynamic Protein-Carbohydrate Interactions of O- and C-Linked Mannosides on the E. coli FimH Adhesin.

Antagonists of the Escherichia coli type-1 fimbrial adhesin FimH are recognized as attractive alternatives for antibiotic therapies and prophylaxes against acute and recurrent bacterial infections. In this study α-d-mannopyranosides O- or C-linked with an alkyl, alkene, alkyne, thioalkyl, amide, or sulfonamide were investigated to fit a hydrophobic substituent with up to two aryl groups within the tyrosine gate emerging from the mannose-binding pocket of FimH. The results were summarized into a set of structure-activity relationships to be used in FimH-targeted inhibitor design: alkene linkers gave an improved affinity and inhibitory potential, because of their relative flexibility combined with a favourable interaction with isoleucine-52 located in the middle of the tyrosine gate. Of particular interest is a C-linked mannoside, alkene-linked to an ortho-substituted biphenyl that has an affinity similar to its O-mannosidic analog but superior to its para-substituted analog. Docking of its high-resolution NMR solution structure to the FimH adhesin indicated that its ultimate, ortho-placed phenyl ring is able to interact with isoleucine-13, located in the clamp loop that undergoes conformational changes under shear force exerted on the bacteria. Molecular dynamics simulations confirmed that a subpopulation of the C-mannoside conformers is able to interact in this secondary binding site of FimH.

The mitochondrial VDAC of bean seeds recruits phosphatidylethanolamine lipids for its proper functioning.

The voltage-dependent anion-selective channel (VDAC) is the main pathway for inorganic ions and metabolites through the mitochondrial outer membrane. Studies recently demonstrated that membrane lipids regulate its function. It remains, however, unclear how this regulation takes place. In this study, we show that phospholipids are key regulators of Phaseolus VDAC function and, furthermore, that the salt concentration modulates this regulation. Both selectivity and voltage dependence of Phaseolus VDAC are very sensitive to a change in the lipid polar head from PC to PE. Interestingly enough, this dependence is observed only at low salt concentration. Furthermore, significant changes in VDAC functional properties also occur with the gradual methylation of the PE group pointing to the role of subtle chemical variations in the lipid head group. The dependence of PcVDAC gating upon the introduction of a small mole fraction of PE in a PC bilayer has prompted us to propose the existence of a specific interaction site for PE on the outer surface of PcVDAC. Eventually, comparative modeling and molecular dynamics simulations suggest a potential mechanism to get insight into the anion selectivity enhancement of PcVDAC observed in PE relative to PC.

Mutation of Tyr137 of the universal Escherichia coli fimbrial adhesin FimH relaxes the tyrosine gate prior to mannose binding.

The most prevalent diseases manifested by Escherichia coli are acute and recurrent bladder infections and chronic inflammatory bowel diseases such as Crohn's disease. E. coli clinical isolates express the FimH adhesin, which consists of a mannose-specific lectin domain connected via a pilin domain to the tip of type 1 pili. Although the isolated FimH lectin domain has affinities in the nanomolar range for all high-mannosidic glycans, differentiation between these glycans is based on their capacity to form predominantly hydrophobic interactions within the tyrosine gate at the entrance to the binding pocket. In this study, novel crystal structures of tyrosine-gate mutants of FimH, ligand-free or in complex with heptyl α-d-O-mannopyranoside or 4-biphenyl α-d-O-mannopyranoside, are combined with quantum-mechanical calculations and molecular-dynamics simulations. In the Y48A FimH crystal structure, a large increase in the dynamics of the alkyl chain of heptyl α-d-O-mannopyranoside attempts to compensate for the absence of the aromatic ring; however, the highly energetic and stringent mannose-binding pocket of wild-type FimH is largely maintained. The Y137A mutation, on the other hand, is the most detrimental to FimH affinity and specificity: (i) in the absence of ligand the FimH C-terminal residue Thr158 intrudes into the mannose-binding pocket and (ii) ethylenediaminetetraacetic acid interacts strongly with Glu50, Thr53 and Asn136, in spite of multiple dialysis and purification steps. Upon mutation, pre-ligand-binding relaxation of the backbone dihedral angles at position 137 in the tyrosine gate and their coupling to Tyr48 via the interiorly located Ile52 form the basis of the loss of affinity of the FimH adhesin in the Y137A mutant.

Unveiling the Mechanism of Arginine Transport through AdiC with Molecular Dynamics Simulations: The Guiding Role of Aromatic Residues.

Commensal and pathogenic enteric bacteria have developed several systems to adapt to proton leakage into the cytoplasm resulting from extreme acidic conditions. One such system involves arginine uptake followed by export of the decarboxylated product agmatine, carried out by the arginine/agmatine antiporter (AdiC), which thus works as a virtual proton pump. Here, using classical and targeted molecular dynamics, we investigated at the atomic level the mechanism of arginine transport through AdiC of E. coli. Overall, our MD simulation data clearly demonstrate that global rearrangements of several transmembrane segments are necessary but not sufficient for achieving transitions between structural states along the arginine translocation pathway. In particular, local structural changes, namely rotameric conversions of two aromatic residues, are needed to regulate access to both the outward- and inward-facing states. Our simulations have also enabled identification of a few residues, overwhelmingly aromatic, which are essential to guiding arginine in the course of its translocation. Most of them belong to gating elements whose coordinated motions contribute to the alternating access mechanism. Their conservation in all known E. coli acid resistance antiporters suggests that the transport mechanisms of these systems share common features. Last but not least, knowledge of the functional properties of AdiC can advance our understanding of the members of the amino acid-carbocation-polyamine superfamily, notably in eukaryotic cells.

The Antiadhesive Strategy in Crohn's Disease: Orally Active Mannosides to Decolonize Pathogenic Escherichia coli from the Gut.

Blocking the adherence of bacteria to cells is an attractive complementary approach to current antibiotic treatments, which are faced with increasing resistance. This strategy has been particularly studied in the context of urinary tract infections (UTIs), in which the adhesion of pathogenic Escherichia coli strains to uroepithelial cells is prevented by blocking the FimH adhesin expressed at the tips of bacteria organelles called fimbriae. Recently, we extended the antiadhesive concept, showing that potent FimH antagonists can block the attachment of adherent-invasive E. coli (AIEC) colonizing the intestinal mucosa of patients with Crohn's disease (CD). In this work, we designed a small library of analogues of heptyl mannoside (HM), a previously identified nanomolar FimH inhibitor, but one that displays poor antiadhesive effects in vivo. The anomeric oxygen atom was replaced by a sulfur or a methylene group to prevent hydrolysis by intestinal glycosidases, and chemical groups were attached at the end of the alkyl tail. Importantly, a lead compound was shown to reduce AIEC levels in the feces and in the colonic and ileal mucosa after oral administration (10 mg kg(-1) ) in a transgenic mouse model of CD. The compound showed a low bioavailability, preferable in this instance, thus suggesting the possibility of setting up an innovative antiadhesive therapy, based on the water-soluble and non-cytotoxic FimH antagonists developed here, for the CD subpopulation in which AIEC plays a key role.

Function and Regulation of Fungal Amino Acid Transporters: Insights from Predicted Structure.

Amino acids constitute a major nutritional source for probably all fungi. Studies of model species such as the yeast Saccharomyces cerevisiae and the filamentous fungus Aspergillus nidulans have shown that they possess multiple amino acid transporters. These proteins belong to a limited number of superfamilies, now defined according to protein fold in addition to sequence criteria, and differ in subcellular location, substrate specificity range, and regulation. Structural models of several of these transporters have recently been built, and the detailed molecular mechanisms of amino acid recognition and translocation are now being unveiled. Furthermore, the particular conformations adopted by some of these transporters in response to amino acid binding appear crucial to promoting their ubiquitin-dependent endocytosis and/or to triggering signaling responses. We here summarize current knowledge, derived mainly from studies on S. cerevisiae and A. nidulans, about the transport activities, regulation, and sensing role of fungal amino acid transporters, in relation to predicted structure.

Dual mechanism of ion permeation through VDAC revealed with inorganic phosphate ions and phosphate metabolites.

In the exchange of metabolites and ions between the mitochondrion and the cytosol, the voltage-dependent anion channel (VDAC) is a key element, as it forms the major transport pathway for these compounds through the mitochondrial outer membrane. Numerous experimental studies have promoted the idea that VDAC acts as a regulator of essential mitochondrial functions. In this study, using a combination of molecular dynamics simulations, free-energy calculations, and electrophysiological measurements, we investigated the transport of ions through VDAC, with a focus on phosphate ions and metabolites. We showed that selectivity of VDAC towards small anions including monovalent phosphates arises from short-lived interactions with positively charged residues scattered throughout the pore. In dramatic contrast, permeation of divalent phosphate ions and phosphate metabolites (AMP and ATP) involves binding sites along a specific translocation pathway. This permeation mechanism offers an explanation for the decrease in VDAC conductance measured in the presence of ATP or AMP at physiological salt concentration. The binding sites occur at similar locations for the divalent phosphate ions, AMP and ATP, and contain identical basic residues. ATP features a marked affinity for a central region of the pore lined by two lysines and one arginine of the N-terminal helix. This cluster of residues together with a few other basic amino acids forms a "charged brush" which facilitates the passage of the anionic metabolites through the pore. All of this reveals that VDAC controls the transport of the inorganic phosphates and phosphate metabolites studied here through two different mechanisms.

Chemical cross-linking/mass spectrometry maps the amyloid ß peptide binding region on both apolipoprotein E domains.

Apolipoprotein E (apoE) binds the amyloid ß peptide (Aß), one of the major culprits in Alzheimer's disease development. The formation of apoE:Aß complexes is implicated in both Aß clearance and fibrillization. However, the binding interface between apoE and Aß is poorly defined despite substantial previous research efforts, and the exact role of apoE in the pathology of Alzheimer's disease remains largely elusive. Here, we compared the three main isoforms of apoE (E2, E3, and E4) for their interaction with Aß1-42 in an early stage of aggregation and at near physiological conditions. Using electron microscopy and Western blots, we showed that all three isoforms are able to prevent Aß fibrillization and form a noncovalent complex, with one molecule of Aß bound per apoE. Using chemical cross-linking coupled to mass spectrometry, we further examined the interface of interaction between apoE2/3/4 and Aß. Multiple high-confidence intermolecular apoE2/3/4:Aß cross-links confirmed that Lys16 is located in the region of Aß binding to apoE2/3/4. Further, we demonstrated that both N- and C-terminal domains of apoE2/3/4 are interacting with Aß. The cross-linked sites were mapped onto and evaluated in light of a recent structure of apoE. Our results support binding of the hydrophobic Aß at the apoE domain-domain interaction interface, which would explain how apoE is able to stabilize Aß and thereby prevent its subsequent aggregation.

Substrate-induced ubiquitylation and endocytosis of yeast amino acid permeases.

Many plasma membrane transporters are downregulated by ubiquitylation, endocytosis, and delivery to the lysosome in response to various stimuli. We report here that two amino acid transporters of Saccharomyces cerevisiae, the general amino acid permease (Gap1) and the arginine-specific permease (Can1), undergo ubiquitin-dependent downregulation in response to their substrates and that this downregulation is not due to intracellular accumulation of the transported amino acids but to transport catalysis itself. Following an approach based on permease structural modeling, mutagenesis, and kinetic parameter analysis, we obtained evidence that substrate-induced endocytosis requires transition of the permease to a conformational state preceding substrate release into the cell. Furthermore, this transient conformation must be stable enough, and thus sufficiently populated, for the permease to undergo efficient downregulation. Additional observations, including the constitutive downregulation of two active Gap1 mutants altered in cytosolic regions, support the model that the substrate-induced conformational transition inducing endocytosis involves remodeling of cytosolic regions of the permeases, thereby promoting their recognition by arrestin-like adaptors of the Rsp5 ubiquitin ligase. Similar mechanisms might control many other plasma membrane transporters according to the external concentrations of their substrates.

Origin of ion selectivity in Phaseolus coccineus mitochondrial VDAC.

The mitochondrial voltage-dependent a nion-selective channel (VDAC) is the major permeation pathway for small ions and metabolites. Although a wealth of electrophysiological data has been obtained on different VDAC species, the physical mechanisms of their ionic selectivity are still elusive. We addressed this issue using electrophysiological experiments performed on plant VDAC. A simple macroscopic electrodiffusion model accounting for ion diffusion and for an effective fixed charge of the channel describes well its selectivity. Brownian Dynamics simulations of ion permeation performed on plant and mammalian VDACs point to the role of specific charged residues located at about the middle of the pore.