Facilitated substrate transport through membrane proteins

Counter-Intuitive Features of Particle Dynamics in Nanopores

Using the framework of a continuous diffusion model based on the Smoluchowski equation, we analyze particle dynamics in the confinement of a transmembrane nanopore. We briefly review existing analytical results to highlight consequences of interactions between the channel nanopore and the translocating particles. These interactions are described within a minimalistic approach by lumping together multiple physical forces acting on the particle in the pore into a one-dimensional potential of mean force. Such radical simplification allows us to obtain transparent analytical results, often in a simple algebraic form. While most of our findings are quite intuitive, some of them may seem unexpected and even surprising at first glance. The focus is on five examples: (i) attractive interactions between the particles and the nanopore create a potential well and thus cause the particles to spend more time in the pore but, nevertheless, increase their net flux; (ii) if the potential well-describing particle-pore interaction occupies only a part of the pore length, the mean translocation time is a non-monotonic function of the well length, first increasing and then decreasing with the length; (iii) when a rectangular potential well occupies the entire nanopore, the mean particle residence time in the pore is independent of the particle diffusivity inside the pore and depends only on its diffusivity in the bulk; (iv) although in the presence of a potential bias applied to the nanopore the "downhill" particle flux is higher than the "uphill" one, the mean translocation times and their distributions are identical, i.e., independent of the translocation direction; and (v) fast spontaneous gating affects nanopore selectivity when its characteristic time is comparable to that of the particle transport through the pore.

Structural displacement model of chitooligosaccharide transport through chitoporin

VhChiP is a chitooligosaccharide-specific porin identified in the outer membrane of Vibrio campbellii type strain American Type Culture Collection BAA 1116. VhChiP contains three identical subunits, and in each subunit, the 19-amino acid N-terminal segment serves as a molecular plug (the "N-plug") that controls the closed/open dynamics of the neighboring pores. In this study, the crystal structures of VhChiP lacking the N-plug were determined in the absence and presence of chitohexaose. Binding studies of sugar-ligand interactions by single-channel recordings and isothermal microcalorimetry experiments suggested that the deletion of the N-plug peptide significantly weakened the sugar-binding affinity due to the loss of hydrogen bonds around the central affinity sites. Steered molecular dynamic simulations revealed that the movement of the sugar chain along the sugar passage triggered the ejection of the N-plug, while the H-bonds transiently formed between the reducing end GlcNAc units of the sugar chain with the N-plug peptide may help to facilitate sugar translocation. The findings enable us to propose the structural displacement model, which enables us to understand the molecular basis of chitooligosaccharide uptake by marine Vibrio bacteria.

Exaggerated trans-membrane charge of ammonium transporters in nutrient-poor marine environments

Transporter proteins are a vital interface between cells and their environment. In nutrient-limited environments, microbes with transporters that are effective at bringing substrates into their cells will gain a competitive advantage over variants with reduced transport function. Microbial ammonium transporters (Amt) bring ammonium into the cytoplasm from the surrounding periplasm space, but diagnosing Amt adaptations to low nutrient environments solely from sequence data has been elusive. Here, we report altered Amt sequence amino acid distribution from deep marine samples compared to variants sampled from shallow water in two important microbial lineages of the marine water column community-Marine Group I Archaea (Thermoproteota) and the uncultivated gammaproteobacterial lineage SAR86. This pattern indicates an evolutionary pressure towards an increasing dipole in Amt for these clades in deep ocean environments and is predicted to generate stronger electric fields facilitating ammonium acquisition. This pattern of increasing dipole charge with depth was not observed in lineages capable of accessing alternative nitrogen sources, including the abundant alphaproteobacterial clade SAR11. We speculate that competition for ammonium in the deep ocean drives transporter sequence evolution. The low concentration of ammonium in the deep ocean is therefore likely due to rapid uptake by Amts concurrent with decreasing nutrient flux.

Current noise of a protein-selective biological nanopore

1/f current noise is ubiquitous in protein pores, porins, and channels. We have previously shown that a protein-selective biological nanopore with an external protein receptor can function as a 1/f noise generator when a high-affinity protein ligand is reversibly captured by the receptor. Here, we demonstrate that the binding affinity and concentration of the ligand are key determinants for the nature of current noise. For example, 1/f was absent when a protein ligand was reversibly captured at a much lower concentration than its equilibrium dissociation constant against the receptor. Furthermore, we also analyzed the composite current noise that resulted from mixtures of low-affinity and high-affinity ligands against the same receptor. This study highlights the significance of protein recognition events in the current noise fluctuations across biological membranes.

Driven transport of soft Brownian particles through pore-like structures: Effective size method

Single-file transport in pore-like structures constitutes an important topic for both theory and experiment. For hardcore interacting particles, a good understanding of the collective dynamics has been achieved recently. Here, we study how softness in the particle interaction affects the emergent transport behavior. To this end, we investigate the driven Brownian motion of particles in a periodic potential. The particles interact via a repulsive softcore potential with a shape corresponding to a smoothed rectangular barrier. This shape allows us to elucidate effects of mutual particle penetration and particle crossing in a controlled manner. We find that even weak deviations from the hardcore case can have a strong impact on the particle current. Despite this fact, knowledge about the transport in a corresponding hardcore system is shown to be useful to describe and interpret our findings for the softcore case. This is achieved by assigning a thermodynamic effective size to the particles based on the equilibrium density functional of hard spheres.

The Whole Is Bigger than the Sum of Its Parts: Drug Transport in the Context of Two Membranes with Active Efflux

Cell envelope plays a dual role in the life of bacteria by simultaneously protecting it from a hostile environment and facilitating access to beneficial molecules. At the heart of this ability lie the restrictive properties of the cellular membrane augmented by efflux transporters, which preclude intracellular penetration of most molecules except with the help of specialized uptake mediators. Recently, kinetic properties of the cell envelope came into focus driven on one hand by the urgent need in new antibiotics and, on the other hand, by experimental and theoretical advances in studies of transmembrane transport. A notable result from these studies is the development of a kinetic formalism that integrates the Michaelis-Menten behavior of individual transporters with transmembrane diffusion and offers a quantitative basis for the analysis of intracellular penetration of bioactive compounds. This review surveys key experimental and computational approaches to the investigation of transport by individual translocators and in whole cells, summarizes key findings from these studies and outlines implications for antibiotic discovery. Special emphasis is placed on Gram-negative bacteria, whose envelope contains two separate membranes. This feature sets these organisms apart from Gram-positive bacteria and eukaryotic cells by providing them with full benefits of the synergy between slow transmembrane diffusion and active efflux.

How to Enter a Bacterium: Bacterial Porins and the Permeation of Antibiotics

Despite tremendous successes in the field of antibiotic discovery seen in the previous century, infectious diseases have remained a leading cause of death. More specifically, pathogenic Gram-negative bacteria have become a global threat due to their extraordinary ability to acquire resistance against any clinically available antibiotic, thus urging for the discovery of novel antibacterial agents. One major challenge is to design new antibiotics molecules able to rapidly penetrate Gram-negative bacteria in order to achieve a lethal intracellular drug accumulation. Protein channels in the outer membrane are known to form an entry route for many antibiotics into bacterial cells. Up until today, there has been a lack of simple experimental techniques to measure the antibiotic uptake and the local concentration in subcellular compartments. Hence, rules for translocation directly into the various Gram-negative bacteria via the outer membrane or via channels have remained elusive, hindering the design of new or the improvement of existing antibiotics. In this review, we will discuss the recent progress, both experimentally as well as computationally, in understanding the structure-function relationship of outer-membrane channels of Gram-negative pathogens, mainly focusing on the transport of antibiotics.

Fake It 'Till You Make It-The Pursuit of Suitable Membrane Mimetics for Membrane Protein Biophysics

Membrane proteins evolved to reside in the hydrophobic lipid bilayers of cellular membranes. Therefore, membrane proteins bridge the different aqueous compartments separated by the membrane, and furthermore, dynamically interact with their surrounding lipid environment. The latter not only stabilizes membrane proteins, but directly impacts their folding, structure and function. In order to be characterized with biophysical and structural biological methods, membrane proteins are typically extracted and subsequently purified from their native lipid environment. This approach requires that lipid membranes are replaced by suitable surrogates, which ideally closely mimic the native bilayer, in order to maintain the membrane proteins structural and functional integrity. In this review, we survey the currently available membrane mimetic environments ranging from detergent micelles to bicelles, nanodiscs, lipidic-cubic phase (LCP), liposomes, and polymersomes. We discuss their respective advantages and disadvantages as well as their suitability for downstream biophysical and structural characterization. Finally, we take a look at ongoing methodological developments, which aim for direct in-situ characterization of membrane proteins within native membranes instead of relying on membrane mimetics.

Protein Ligand-Induced Amplification in the 1/f Noise of a Protein-Selective Nanopore

Previous studies of transmembrane protein channels have employed noise analysis to examine their statistical current fluctuations. In general, these explorations determined a substrate-induced amplification in the Gaussian white noise of these systems at a low-frequency regime. This outcome implies a lack of slowly appearing fluctuations in the number and local mobility of diffusing charges in the presence of channel substrates. Such parameters are among the key factors in generating a low-frequency 1/f noise. Here, we show that a protein-selective biological nanopore exhibits a substrate-induced amplification in the 1/f noise. The modular composition of this biological nanopore includes a hydrophilic transmembrane protein pore fused to a water-soluble binding protein on its extramembranous side. In addition, this protein nanopore shows an open substate populated by a high-frequency current noise because of the flickering of an engineered polypeptide adaptor at the tip of the pore. However, the physical association of the protein ligand with the binding domain reversibly switches the protein nanopore from a high-frequency noise substate into a quiet substate. In the absence of the protein ligand, our nanopore shows a low-frequency white noise. Remarkably, in the presence of the protein ligand, an amplified low-frequency 1/f noise was detected in a ligand concentration-dependent fashion. This finding suggests slowly occurring equilibrium fluctuations in the density and local mobility of charge carriers under these conditions. Furthermore, we report that the excess in 1/f noise is generated by reversible switches between the noisy ligand-released substate and the quiet ligand-captured substate. Finally, quantitative aspects of the low-frequency 1/f noise are in accord with theoretical predictions of the current noise analysis of protein channel-ligand interactions.

Effects of H-bonds on sugar binding to chitoporin from Vibrio harveyi

BACKGROUND

VhChiP is a sugar-specific-porin present in the outer membrane of the marine bacterium Vibrio harveyi and responsible for chitin uptake, with a high selectivity for chitohexaose.

METHODS

VhChiP and its mutants were expressed and purified from BL21 (DE3) Omp8 Rosetta strain. After reconstitution into planar lipid bilayers, the ion current fluctuations caused by chitohexaose entering the channel were measured in deuterium oxide and in water.

RESULTS

The role of hydrogen-bonding in sugar binding was investigated by comparing channel occlusion by chitohexaose in buffers containing H₂O and D₂O. The BLM results revealed the significant contribution of hydrogen bonding to the binding of chitohexaose in the constriction zone of VhChiP. Replacing H₂O as solvent by D₂O significantly decreased the on- and off-rates of sugar penetration into the channel. The importance of hydrogen bonding inside the channel was more noticeable when the hydrophobicity of the constriction zone was diminished by replacing Trp¹³⁶ with the charged residues Asp or Arg. The on- and off-rates decreased up to 2.5-fold and 4-fold when Trp¹³⁶ was replaced by Arg, or 5-fold and 3-fold for Trp¹³⁶ replacement by Asp, respectively. Measuring the on-rate at different temperatures and for different channel mutants revealed the activation energy for chitohexaose entrance into VhChiP channel.

CONCLUSIONS

Hydrogen-bonds contribute to sugar permeation.

Chitosugar translocation by an unexpressed monomeric protein channel

The outer membrane protein channel EcChiP, associated with a silent gene in E. coli, is a monomeric chitoporin. In a glucose-deficient environment, E. coli can express the ChiP gene to exploit chitin degradation products. Single-channel small ion current measurements, which reveal the dynamics of single sugar molecules trapped in channel, are used here to study the exotic transport of chitosugars by E. coli. Molecules escape from the channel on multiple timescales. Voltage-dependent trapping rates observed for charged chitosan molecules, as well as model calculations, indicate that the rapid escape processes are those in which the molecule escapes back to the side of the membrane from which it originated. The probability that a sugar molecule is translocated through the membrane is thus estimated from the current data and the dependence of this translocation probability on the length of the chitosugar molecule and the applied voltage analyzed. The described method for obtaining the translocation probability and related molecular translocation current is applicable to other transport channels.

Functional analysis of an unusual porin-like channel that imports chitin for alternative carbon metabolism in Escherichia coli

Escherichia coli have the genetic potential to use chitin as a carbon source in the absence of glucose, importing it via the chitin-uptake channel EcChiP for processing by the glucosamine catabolic pathway. The chip gene is usually not expressed when E. coli are grown on glucose-enriched nutrients, providing a general regulatory mechanism for the pathway. EcChiP is unusual in that it is homologous to porins and monomeric instead of trimeric, the typical form of sugar-specific channels, making it unclear how this channel operates. We recently reported that EcChiP could form a stable channel in lipid membranes and that the channel is specific for chitooligosaccharides. This report describes the biophysical nature of sugar-channel interactions and the kinetics of sugar association and dissociation. Titrating EcChiP with chitohexaose resulted in protein fluorescence enhancement in a concentration-dependent manner, yielding a binding constant of 2.9 × 10⁵ m^-1, consistent with the value of 2.5 × 10⁵ m^-1 obtained from isothermal titration microcalorimetry. Analysis of the integrated heat change suggested that the binding process was endothermic and driven by entropy. Single-channel recordings confirmed the voltage dependence of the penetration of chitohexaose molecules into and their release from EcChiP. Once inside the pore, the sugar release rate (k_off) from the affinity site increased with elevated voltage, regardless of the side of sugar addition. Our findings revealed distinct thermodynamic and kinetic features of the activity of sugar-specific EcChiP and advance our knowledge of the physiological possibility of chitin utilization by non-chitinolytic bacteria.

Correlated trapping of sugar molecules by the trimeric protein channel chitoporin

The protein channel chitoporin (ChiP), which is used by marine bacteria to translocate selected sugar molecules through the outer cell membrane, is studied via single channel current measurements in water and heavy water sugar solutions. The dynamic trapping and escape probabilities of sugar molecules from different monomers in the trimeric channel are characterized, including their dependence on channel orientation and sensitivity to a deuterium isotope effect. A detailed analysis of stochastic current fluctuations reveals that the trapping properties of chitoporin exhibit memory effects: the rate of trapping transitions depends on the previous sequence of transitions; and intermonomer correlations: the average trapping rate of an unblocked monomer is larger when its neighboring monomers are blocked. The latter, likely resulting from rapid re-trapping of recently escaped sugar molecules, is considered as a possible design strategy to enhance sugar transport.

Chitoporin from the Marine Bacterium Vibrio harveyi: PROBING THE ESSENTIAL ROLES OF TRP136 AT THE SURFACE OF THE CONSTRICTION ZONE

VhChiP is a sugar-specific porin present in the outer membrane of the marine bacterium Vibrio harveyi. VhChiP is responsible for the uptake of chitin oligosaccharides, with particular selectivity for chitohexaose. In this study, we employed electrophysiological and biochemical approaches to demonstrate that Trp(136), located at the mouth of the VhChiP pore, plays an essential role in controlling the channel's ion conductivity, chitin affinity, and permeability. Kinetic analysis of sugar translocation obtained from single channel recordings indicated that the Trp(136) mutations W136A, W136D, W136R, and W136F considerably reduce the binding affinity of the protein channel for its best substrate, chitohexaose. Liposome swelling assays confirmed that the Trp(136) mutations decreased the rate of bulk chitohexaose permeation through the VhChiP channel. Notably, all of the mutants show increases in the off-rate for chitohexaose of up to 20-fold compared with that of the native channel. Furthermore, the cation/anion permeability ratio Pc/Pa is decreased in the W136R mutant and increased in the W136D mutant. This demonstrates that the negatively charged surface at the interior of the protein lumen preferentially attracts cationic species, leading to the cation selectivity of this trimeric channel.

Excess white noise to probe transport mechanisms in a membrane channel

Current fluctuation analysis has been successfully used over the years to investigate the physical properties of different systems. Here, we perform single-channel time-resolved current experiments in a protein channel to evaluate the different transport mechanisms governing the channel function. Using different salts of monovalent and divalent cations in a wide range of concentrations and applied potentials, we analyze current fluctuations focusing on the voltage dependence of the additional white noise that appears in the low-frequency range of the spectra. We demonstrate that the channel displays two characteristic transport regimes: at low salt concentrations (10 mM to 1 M) ion permeation is controlled by the protein fixed charges that induce accumulation or exclusion of ions to preserve local electroneutrality. At high salt concentrations (>1M) adsorption processes associated to the binding of cations to the channel charges regulate the transport properties.

Analysis of fast channel blockage: revealing substrate binding in the microsecond range

For an antibiotic to be effective, it needs to cross the outer membrane barrier and reach the target inside the cell. Hydrophilic antibiotics, e.g.β-lactams, use porin channels to cross the outer membrane and accumulate in the periplasm. Experimental determination of antibiotic interactions with porin is performed by using electrophysiology on a single channel level by noise analysis or single event analysis methods. We report a novel framework for analyzing the ion-current noise, taking into account the corrections due to the analogous filter and the sampling procedure, with the goal of extending the time resolution to a range previously inaccessible by event analysis or by conventional noise analysis. The new method allows one to analyse fast binding events and/or the case when the single channel is not completely blocked by the substrate. We demonstrate the power of this approach by using as an example the interactions of meropenem, an antibiotic of the carbapenem family, with the OmpF porin that is considered to be one of the main pathways for antibiotics to enter Escherichia coli. The presence of meropenem in OmpF is detected by ion current blockages, and the on and off rates are estimated from the concentration dependence of the average ion current and of its power spectral density. The obtained average residence time of the antibiotic inside the channel is in the range of a few microseconds, i.e. more than 50 times smaller than the inverse cut-off frequency of the analogous filter.

Role of the potential landscape on the single-file diffusion through channels

Transport of colloid particles through narrow channels is ubiquitous in cell biology as well as becoming increasingly important for microfluidic applications or targeted drug delivery. Membrane channels in cells are useful models for artificial designs because of their high efficiency, selectivity, and robustness to external fluctuations. Here, we model the passive channels that let cargo simply diffuse through them, affected by a potential profile along the way. Passive transporters achieve high levels of efficiency and specificity from binding interactions with the cargo inside the channel. This however leads to a paradox: why should channels which are so narrow that they are blocked by their cargo evolve to have binding regions for their cargo if that will effectively block them? Using Brownian dynamics simulations, we show that different potentials, notably symmetric, increase the flux through narrow passive channels - and investigate how shape and depth of potentials influence the flux. We find that there exist optimal depths for certain potential shapes and that it is most efficient to apply a small force over an extended region of the channel. On the other hand, having several spatially discrete binding pockets will not alter the flux significantly. We also explore the role of many-particle effects arising from pairwise particle interactions with their neighbours and demonstrate that the relative changes in flux can be accounted for by the kinetics of the absorption reaction at the end of the channel.

Channel-facilitated diffusion boosted by particle binding at the channel entrance

We investigate single-file diffusion of Brownian particles in arrays of closely confining microchannels permeated by a variety of attractive optical potentials and connecting two baths with equal particle concentration. We simultaneously test free diffusion in the channel, diffusion in optical traps coupled in the center of the channel, and diffusion in traps extending into the baths. We found that both classes of attractive optical potentials enhance the translocation rate through the channel with respect to free diffusion. Surprisingly, for the latter class of potentials we measure a 40-fold enhancement in the translocation rate with respect to free diffusion and find a sublinear power law dependence of the translocation rate on the average number of particles in the channel. Our results reveal the function of particle binding at the channel entrances for diffusive transport and open the way to a better understanding of membrane transport and design of synthetic membranes with enhanced diffusion rate.

Anisotropic diffusion of spherical particles in closely confining microchannels

We present here the measurement of the diffusivity of spherical particles closely confined by narrow microchannels. Our experiments yield a two-dimensional map of the position-dependent diffusion coefficients parallel and perpendicular to the channel axis with a resolution down to 129 nm. The diffusivity was measured simultaneously in the channel interior, the bulk reservoirs, as well as the channel entrance region. In the channel interior we found strongly anisotropic diffusion. While the perpendicular diffusion coefficient close to the confining walls decreased down to approximately 25% of the value on the channel axis, the parallel diffusion coefficient remained constant throughout the entire channel width. In addition to the experiment, we performed finite element simulations for the diffusivity in the channel interior and found good agreement with the measurements. Our results reveal the distinctive influence of strong confinement on Brownian motion, which is of significance to microfluidics as well as quantitative models of facilitated membrane transport.

Single-molecule trapping dynamics of sugar-uptake channels in marine bacteria

Stochastic fluctuations of ion current through one chitoporin (ChiP) channel within a bilayer lipid membrane in sugar solution are analyzed. These reflect single-molecule dynamics, which indicate that ChiP has multiple binding sites for sugar and exploits interactions between bound molecules to direct sugar passage. Since ChiP is used by marine bacteria, this is likely an adaptive strategy to enhance sugar translocation from rough water.

Chitoporin from Vibrio harveyi, a channel with exceptional sugar specificity

Chitoporin (VhChiP) is a sugar-specific channel responsible for the transport of chitooligosaccharides through the outer membrane of the marine bacterium Vibrio harveyi. Single channel reconstitution into black lipid membrane allowed single chitosugar binding events in the channel to be resolved. VhChiP has an exceptionally high substrate affinity, with a binding constant of K = 5.0 × 10(6) M(-1) for its best substrate (chitohexaose). The on-rates of chitosugars depend on applied voltages, as well as the side of the sugar addition, clearly indicating the inherent asymmetry of the VhChiP lumen. The binding affinity of VhChiP for chitohexaose is 1-5 orders of magnitude larger than that of other known sugar-specific porins for their preferred substrates. Thus, VhChiP is the most potent sugar-specific channel reported to date, with its high efficiency presumably reflecting the need for the bacterium to take up chitin-containing nutrients promptly under turbulent aquatic conditions to exploit them efficiently as its sole source of energy.

Molecular uptake of chitooligosaccharides through chitoporin from the marine bacterium Vibrio harveyi

BACKGROUND

Chitin is the most abundant biopolymer in marine ecosystems. However, there is no accumulation of chitin in the ocean-floor sediments, since marine bacteria Vibrios are mainly responsible for a rapid turnover of chitin biomaterials. The catabolic pathway of chitin by Vibrios is a multi-step process that involves chitin attachment and degradation, followed by chitooligosaccharide uptake across the bacterial membranes, and catabolism of the transport products to fructose-6-phosphate, acetate and NH(3).

PRINCIPAL FINDINGS

This study reports the isolation of the gene corresponding to an outer membrane chitoporin from the genome of Vibrio harveyi. This porin, expressed in E. coli, (so called VhChiP) was found to be a SDS-resistant, heat-sensitive trimer. Immunoblotting using anti-ChiP polyclonal antibody confirmed the expression of the recombinant ChiP, as well as endogenous expression of the native protein in the V. harveyi cells. The specific function of VhChiP was investigated using planar lipid membrane reconstitution technique. VhChiP nicely inserted into artificial membranes and formed stable, trimeric channels with average single conductance of 1.8±0.13 nS. Single channel recordings at microsecond-time resolution resolved translocation of chitooligosaccharides, with the greatest rate being observed for chitohexaose. Liposome swelling assays showed no permeation of other oligosaccharides, including maltose, sucrose, maltopentaose, maltohexaose and raffinose, indicating that VhChiP is a highly-specific channel for chitooligosaccharides.

CONCLUSION/SIGNIFICANCE

We provide the first evidence that chitoporin from V. harveyi is a chitooligosaccharide specific channel. The results obtained from this study help to establish the fundamental role of VhChiP in the chitin catabolic cascade as the molecular gateway that Vibrios employ for chitooligosaccharide uptake for energy production.

OccK channels from Pseudomonas aeruginosa exhibit diverse single-channel electrical signatures but conserved anion selectivity

Pseudomonas aeruginosa is a Gram-negative bacterium that utilizes substrate-specific outer membrane (OM) proteins for the uptake of small, water-soluble nutrients employed in the growth and function of the cell. In this paper, we present for the first time a comprehensive single-channel examination of seven members of the OM carboxylate channel K (OccK) subfamily. Recent biochemical, functional, and structural characterization of the OccK proteins revealed their common features, such as a closely related, monomeric, 18-stranded β-barrel conformation with a kidney-shaped transmembrane pore and the presence of a basic ladder within the channel lumen. Here, we report that the OccK proteins exhibited fairly distinct unitary conductance values, in a much broader range than previously expected, which includes low (~40-100 pS) and medium (~100-380 pS) conductance. These proteins showed diverse single-channel dynamics of current gating transitions, revealing one-open substate (OccK3), two-open substate (OccK4-OccK6), and three-open substate (OccK1, OccK2, and OccK7) kinetics with functionally distinct conformations. Interestingly, we discovered that anion selectivity is a conserved trait among the members of the OccK subfamily, confirming the presence of a net pool of positively charged residues within their central constriction. Moreover, these results are in accord with an increased specificity and selectivity of these protein channels for negatively charged, carboxylate-containing substrates. Our findings might ignite future functional examinations and full atomistic computational studies for unraveling a mechanistic understanding of the passage of small molecules across the lumen of substrate-specific, β-barrel OM proteins.

Effect of interactions on molecular fluxes and fluctuations in the transport across membrane channels

Transport of molecules across membrane channels is investigated theoretically using exactly solvable one-dimensional discrete-state stochastic models. An interaction between molecules and membrane pores is modeled via a set of binding sites with different energies. It is shown that the interaction potential strongly influences the particle currents as well as fluctuations in the number of translocated molecules. For small concentration gradients, the attractive sites lead to largest currents and fluctuations, while the repulsive interactions yield the largest fluxes and dispersions for large concentration gradients. Interaction energies that lead to maximal currents and maximal fluctuations are the same only for locally symmetric potentials, where transition states are equally distant from the neighboring binding sites, while they differ for the locally asymmetric potentials. The conditions for the most optimal translocation transport with maximal current and minimal dispersion are discussed. It is argued that, in this case, the interaction strength is independent of local symmetry of the potential of mean forces. In addition, the effect of the global asymmetry of the interaction potential is investigated, and it is shown that it also strongly affects the particle translocation dynamics. These phenomena can be explained by analyzing the details of the particle entering and leaving the binding sites in the channel.

Sugar transport across lactose permease probed by steered molecular dynamics

Escherichia coli lactose permease (LacY) transports sugar across the inner membrane of the bacterium using the proton motive force to accumulate sugar in the cytosol. We have probed lactose conduction across LacY using steered molecular dynamics, permitting us to follow molecular and energetic details of lactose interaction with the lumen of LacY during its permeation. Lactose induces a widening of the narrowest parts of the channel during permeation, the widening being largest within the periplasmic half-channel. During permeation, the water-filled lumen of LacY only partially hydrates lactose, forcing it to interact with channel lining residues. Lactose forms a multitude of direct sugar-channel hydrogen bonds, predominantly with residues of the flexible N-domain, which is known to contribute a major part of LacY's affinity for lactose. In the periplasmic half-channel lactose predominantly interacts with hydrophobic channel lining residues, whereas in the cytoplasmic half-channel key protein-substrate interactions are mediated by ionic residues. A major energy barrier against transport is found within a tight segment of the periplasmic half-channel where sugar hydration is minimal and protein-sugar interaction maximal. Upon unbinding from the binding pocket, lactose undergoes a rotation to permeate either half-channel with its long axis aligned parallel to the channel axis. The results hint at the possibility of a transport mechanism, in which lactose permeates LacY through a narrow periplasmic half-channel and a wide cytoplasmic half-channel, the opening of which is controlled by changes in protonation states of key protein side groups.

Channel-facilitated molecular transport across membranes: attraction, repulsion, and asymmetry

Transport of molecules across membrane channels is investigated theoretically using exactly solvable discrete stochastic site-binding models. It is shown that the interaction potential between molecules and the channel has a strong effect on translocation dynamics. The presence of attractive binding sites in the pore accelerates the particle current for small concentrations outside the membrane, while for large concentrations, surprisingly, repulsive binding sites yield the most optimal transport. In addition, the asymmetry of the interaction potential also strongly influences the channel transport. The mechanism underlying these phenomena is discussed using the details of particle dynamics at the binding sites.

Molecular transport through channels and pores: effects of in-channel interactions and blocking

Facilitated translocation of molecules through channels and pores is of fundamental importance for transmembrane transport in biological systems. Several such systems have specific binding sites inside the channel, but a clear understanding of how the interaction between channel and molecules affects the flow is still missing. We present a generic analytical treatment of the problem that relates molecular flow to the first passage time across and the number of particles inside the channel. Both quantities depend in different ways on the channel properties. For the idealized case of noninteracting molecules, we find an increased flow whenever there is a binding site in the channel, despite an increased first passage time. In the more realistic case that molecules may block the channel, we find an increase of flow only up to a certain threshold value of the binding strength and a dependence on the sign of the concentration gradient, i.e., asymmetric transport. The optimal binding strength in that case is analyzed. In all cases the reason for transport facilitation is an increased occupation probability of a particle inside the channel that overcomes any increase in the first passage time because of binding.

Probing the orientation of reconstituted maltoporin channels at the single-protein level

Recently we have shown that maltoporin channels reconstituted into black lipid membranes have pronounced asymmetric properties in both ion conduction and sugar binding. This asymmetry revealed also that maltoporin insertion is directional. However, the orientation in the lipid bilayer remained an open question. To elucidate the orientation, we performed point mutations at each side of the channel and analyzed the ion current fluctuation caused by an asymmetric maltohexaose addition. In a second series we used a chemically modified maltohexaose sugar molecule with inhibited entry possibility from the periplasmic side. In contrast to the natural outer cell wall of bacteria, we found that the maltoporin inserts in artificial lipid bilayer in such a way that the long extracellular loops are exposed to the same side of the membrane than protein addition. Based on this orientation, the directional properties of sugar binding were correlated to physiological conditions. We found that nature has optimized maltoporin channels by lowering the activation barriers at each extremity of the pore to trap sugar molecules from the external medium and eject them most efficiently to the periplasmic side.

On translocation through a membrane channel via an internal binding site: kinetics and voltage dependence

Here we present a model for maltodextrin translocation through maltoporin channels. In a first step, our theoretical analysis does consider the case of a single binding site for a given substrate in a structurally unaffected channel with a possibly different entrance barrier on either side. It is shown how by means of conventional electrical conductance measurements (including current noise analysis) the basic equilibrium and rate constants can be determined as functions of the applied voltage. Then also the net translocation rate of the substrate becomes accessible quantitatively. This most simple model mechanism has been extended to include a voltage-dependent fast conformational change of the channel that prevents the binding process. The so developed approach has been tested with experimental data for a single maltoporin trimer being reconstituted in black lipid membranes when studied in the presence of maltohexaose as the substrate. The experimental results turned out to be clearly incompatible with binding alone. They are, however, very satisfactorily fitted by pertinent theoretical curves if also inhibition of binding by a conformational transition is taken into account. Accordingly, quantitative evaluations of the underlying parameters and eventually of the translocation rate have been carried out successfully. Our analysis reveals a set of parameters necessary for an optimal translocation that nicely corresponds to natural conditions.

Site-directed mutagenesis of tyrosine 118 within the central constriction site of the LamB (maltoporin) channel of Escherichia coli. II. Effect on maltose and maltooligosaccharide binding kinetics

The 3-D structure of the maltooligosaccharide-specific LamB channel of Escherichia coli (also called maltoporin) is known from x-ray crystallography. The central constriction of the channel formed by the external loop 3 is controlled by tyrosine 118. Y118 was replaced by site-directed mutagenesis by 10 other amino acids (alanine (A), isoleucine (I), asparagine (N), serine (S), cysteine (C), aspartic acid (D), arginine (R), histidine (H), phenylalanine (F), and tryptophan (W)) including neutral ones, negatively and positively charged amino acids to study the effect of their size, their hydrophobicity index, and their charge on maltose and maltooligosaccharide binding to LamB. The mutants were reconstituted into lipid bilayer membranes and the stability constants for binding of maltose, maltotriose, maltopentaose, and maltoheptaose to the channel were measured using titration experiments. The mutation of Y118 to any other non-aromatic amino acid led to a substantial decrease of the stability constant of binding by factors between about two and six. The highest effect was observed for the mutant Y118A. Replacement of Y118 by the two other aromatic amino acids, phenylalanine (F) and tryptophan (W), resulted in a substantial increase of the stability constant maximally by a factor of almost 400 for the Y118W mutant. The carbohydrate-induced block of the channel function was used for the study of current noise through the different mutant LamB channels. The analysis of the power density spectra allowed the evaluation of the on- and off-rate constants (k(1) and k(-1)) of sugar binding. The results suggest that both rate constants were affected by the mutations. For most mutants, with the exception of Y118F and Y118W, k(1) decreased and k(-1) increased, whereas the opposite was found for the aromatic amino acid mutants. The results suggest that tyrosine 118 has a crucial effect on carbohydrate transport through LamB.

Sugar transport through maltoporin of Escherichia coli: role of the greasy slide

The lining of the maltodextrin-specific maltoporin (LamB) channel exhibits a string of aromatic residues, the greasy slide, part of which has been shown previously by crystallography to be involved in substrate binding. To probe the functional role of the greasy slide, alanine scanning mutagenesis has been performed on the six greasy slide residues and Y118 at the channel constriction. The mutants were characterized by an in vivo uptake assay and sugar-induced-current-noise analysis. Crystallographic analysis of the W74A mutant showed no perturbation of the structure. All mutants showed considerably decreased maltose uptake rates in vivo (<10% of the wild-type value), indicating the functional importance of the investigated residues. Substitutions at the channel center revealed appreciably increased (up to 100-fold) in vitro half-saturation concentrations for maltotriose and maltohexaose binding to the channel. Sugar association rates, however, were significantly affected also by the mutations at either end of the slide (W74A, W358A, and F227A), an effect which became most apparent upon nonsymmetrical sugar addition. The kinetic data are discussed on the basis of an asymmetric one-site two-barrier model, which suggests that, at low substrate concentrations, as are found under physiological conditions, only the heights of the extracellular and periplasmic barriers, which are reduced by the presence of the greasy slide, determine the efficiency of this facilitated diffusion channel.

Energetics of glycerol conduction through aquaglyceroporin GlpF

Aquaglyceroporin GlpF selectively conducts water and linear polyalcohols, such as glycerol, across the inner membrane of Escherichia coli. We report steered molecular dynamics simulations of glycerol conduction through GlpF, in which external forces accelerate the transchannel conduction in a manner that preserves the intrinsic conduction mechanism. The simulations reveal channel-glycerol hydrogen bonding interactions and the stereoselectivity of the channel. Employing Jarzynski's identity between free energy and irreversible work, we reconstruct the potential of mean force along the conduction pathway through a time series analysis of molecular dynamics trajectories. This potential locates binding sites and barriers inside the channel; it also reveals a low energy periplasmic vestibule suited for efficient uptake of glycerol from the environment.

Transport of maltodextrins through maltoporin: a single-channel study

Transport of sugars through maltoporin channels reconstituted into planar lipid membranes has traditionally been addressed using multichannel preparations. Here we show that single-channel experiments offer new possibilities to reveal molecular details of the interaction between the sugar and the channel. We analyze time-resolved transient interruptions in the maltoporin ionic current in the presence of differently sized maltodextrins. We find for all studied sugars, from maltotriose to maltoheptaose, that only one sugar molecule is required to completely block one of the pores in the maltoporin trimer. The probability of simultaneous blockage of different pores increases with sugar concentration in a manner that demonstrates their mutual independence. The maltoporin channel is asymmetric and, added from one side only, predominantly inserts in an oriented manner. The asymmetry of the channel structure manifests itself in two ways. First, it is seen as an asymmetrical response to applied voltage at otherwise symmetrical conditions; second, as asymmetrical rates of sugar entry into the channel with asymmetrical (one-sided) sugar addition. Importantly, we find that the sugar residence time in the pore does not depend on which side the sugar is added. This voltage-dependent time is the same for symmetrical, cis, or trans sugar addition. This observation suggests that once a sugar molecule is captured by the "greasy slide" of the channel, it spends enough time there to "forget" from what entrance it was captured. This also means that the blockage events studied here represent sugar translocation events, and not just binding at and release from the same entrance of the channel.

Dynamics of nucleotides in VDAC channels: structure-specific noise generation

Nucleotide penetration into the voltage-dependent mitochondrial ion channel (VDAC) reduces single-channel conductance and generates excess current noise through a fully open channel. VDAC channels were reconstituted into planar phospholipid membranes bathed in 1.0 M NaCl. At a given nucleotide concentration, the average decrease in small-ion channel conductance induced by mononucleotides ATP, ADP, AMP, and UTP and dinucleotides beta- and alpha-NADH, NAD, and NADPH are very close. However, the excess current noise is about seven times higher in the presence of NADPH than in the presence of ATP and is about 40 times higher than in the presence of UTP. The nucleotide-generated low-frequency noise obeys the following sequence: beta-NADPH > beta-NADH = alpha-NADH > ATP > ADP > beta-NAD > or = AMP > UTP. Measurements of bulk-phase diffusion coefficients and of the effective charge of the nucleotides in 1.0 M NaCl suggest that differences in size and charge cannot be the major factors responsible for the ability to generate current noise. Thus, although the ability of nucleotides to partition into the channel's pore, as assessed by the reduction in conductance, is very similar, the ability to generate current noise involves a detailed recognition of the three-dimensional structure of the nucleotide by the VDAC channel. A possible mechanism for this selectivity is two noise-generating processes operating in parallel.

The mechanism of glycerol conduction in aquaglyceroporins

BACKGROUND

The E. coli glycerol facilitator, GlpF, selectively conducts glycerol and water, excluding ions and charged solutes. The detailed mechanism of the glycerol conduction and its relationship to the characteristic secondary structure of aquaporins and to the NPA motifs in the center of the channel are unknown.

RESULTS

Molecular dynamics simulations of GlpF reveal spontaneous glycerol and water conduction driven, on a nanosecond timescale, by thermal fluctuations. The bidirectional conduction, guided and facilitated by the secondary structure, is characterized by breakage and formation of hydrogen bonds for which water and glycerol compete. The conduction involves only very minor changes in the protein structure, and cooperativity between the GlpF monomers is not evident. The two conserved NPA motifs are strictly linked together by several stable hydrogen bonds and their asparagine side chains form hydrogen bonds with the substrates passing the channel in single file.

CONCLUSIONS

A complete conduction of glycerol through the GlpF was deduced from molecular dynamics simulations, and key residues facilitating the conduction were identified. The nonhelical parts of the two half-membrane-spanning segments expose carbonyl groups towards the channel interior, establishing a curve-linear pathway. The conformational stability of the NPA motifs is important in the conduction and critical for selectivity. Water and glycerol compete in a random manner for hydrogen bonding sites in the protein, and their translocations in single file are correlated. The suggested conduction mechanism should apply to the whole family.