Understanding ion conductance on a molecular level: an all-atom modeling of the bacterial porin OmpF

Gating of β-Barrel Protein Pores, Porins, and Channels: An Old Problem with New Facets

β barrels are ubiquitous proteins in the outer membranes of mitochondria, chloroplasts, and Gram-negative bacteria. These transmembrane proteins (TMPs) execute a wide variety of tasks. For example, they can serve as transporters, receptors, membrane-bound enzymes, as well as adhesion, structural, and signaling elements. In addition, multimeric β barrels are common structural scaffolds among many pore-forming toxins. Significant progress has been made in understanding the functional, structural, biochemical, and biophysical features of these robust and versatile proteins. One frequently encountered fundamental trait of all β barrels is their voltage-dependent gating. This process consists of reversible or permanent conformational transitions between a large-conductance, highly permeable open state and a low-conductance, solute-restrictive closed state. Several intrinsic molecular mechanisms and environmental factors modulate this universal property of β barrels. This review article outlines the typical signatures of voltage-dependent gating. Moreover, we discuss recent developments leading to a better qualitative understanding of the closure dynamics of these TMPs.

Conformational Dynamics of Loop L3 in OmpF: Implications toward Antibiotic Translocation and Voltage Gating

In the present work, we delineate the molecular mechanism of a bulky antibiotic permeating through a bacterial channel and uncover the role of conformational dynamics of the constriction loop in this process. Using the temperature accelerated sliced sampling approach, we shed light onto the dynamics of the L3 loop, in particular the F118 to S125 segment, at the constriction regions of the OmpF porin. We complement the findings with single channel electrophysiology experiments and applied-field simulations, and we demonstrate the role of hydrogen-bond stabilization in the conformational dynamics of the L3 loop. A molecular mechanism of permeation is put forward wherein charged antibiotics perturb the network of stabilizing hydrogen-bond interactions and induce conformational changes in the L3 segment, thereby aiding the accommodation and permeation of bulky antibiotic molecules across the constriction region. We complement the findings with single channel electrophysiology experiments and demonstrate the importance of the hydrogen-bond stabilization in the conformational dynamics of the L3 loop. The generality of the present observations and experimental results regarding the L3 dynamics enables us to identify this L3 segment as the source of gating. We propose a mechanism of OmpF gating that is in agreement with previous experimental data that showed the noninfluence of cysteine double mutants that tethered the L3 tip to the barrel wall on the OmpF gating behavior. The presence of similar loop stabilization networks in porins of other clinically relevant pathogens suggests that the conformational dynamics of the constriction loop is possibly of general importance in the context of antibiotic permeation through porins.

Structure and Computational Electrophysiology of Ac-LS3, a Synthetic Ion Channel

Computer simulations are reported on Ac-LS3, a synthetic ion channel, containing 21 residues with a Leu-Ser-Ser-Leu-Leu-Ser-Leu heptad repeat, which forms ions channels upon application of voltage. A hexameric, coiled-coil bundle initially positioned perpendicular to the membrane settled into a stable, tilted structure after 1.5 μs, most likely to improve contacts between the non-polar exterior of the channel and the hydrophobic core of the membrane. Once tilted, the bundle remained in this state during subsequent simulations of nearly 10 μs at voltages ranging from 200 to -100 mV. In contrast, attempts to identify a stable pentameric structure failed, thus supporting the hypothesis that the channel is a hexamer. Results at 100 mV were used to reconstruct the free energy profiles for K⁺ and Cl^- in the channel. This was done by way of several methods in which results of molecular dynamics (MD) simulations were combined with the electrodiffusion model. Two of them developed recently do not require knowledge of the diffusivity. Instead, they utilize one-sided density profiles and committor probabilities. The consistency between different methods is very good, supporting the utility of the newly developed methods for reconstructing free energies of ions in channels. The flux of K⁺, which accounts for most of the current through the channel, calculated directly from MD matches well the total measured current. However, the current of Cl^- is somewhat overestimated, possibly due to a slightly unbalanced force field involving chloride. The current-voltage dependence was also reconstructed by way of a recently developed, efficient method that requires simulations only at a single voltage, yielding good agreement with the experiment. Taken together, the results demonstrate that computational electrophysiology has become a reliable tool for studying how channels mediate ion transport through membranes.

Rationalizing the generation of broad spectrum antibiotics with the addition of a positive charge

Antibiotic resistance of Gram-negative bacteria is largely attributed to the low permeability of their outer membrane (OM). Recently, we disclosed the eNTRy rules, a key lesson of which is that the introduction of a primary amine enhances OM permeation in certain contexts. To understand the molecular basis for this finding, we perform an extensive set of molecular dynamics (MD) simulations and free energy calculations comparing the permeation of aminated and amine-free antibiotic derivatives through the most abundant OM porin of E. coli, OmpF. To improve sampling of conformationally flexible drugs in MD simulations, we developed a novel, Monte Carlo and graph theory based algorithm to probe more efficiently the rotational and translational degrees of freedom visited during the permeation of the antibiotic molecule through OmpF. The resulting pathways were then used for free-energy calculations, revealing a lower barrier against the permeation of the aminated compound, substantiating its greater OM permeability. Further analysis revealed that the amine facilitates permeation by enabling the antibiotic to align its dipole to the luminal electric field of the porin and form favorable electrostatic interactions with specific, highly-conserved charged residues. The importance of these interactions in permeation was further validated with experimental mutagenesis and whole cell accumulation assays. Overall, this study provides insights on the importance of the primary amine for antibiotic permeation into Gram-negative pathogens that could help the design of future antibiotics. We also offer a new computational approach for calculating free-energy of processes where relevant molecular conformations cannot be efficiently captured.

Inward-facing glycine residues create sharp turns in β-barrel membrane proteins

The transmembrane region of outer-membrane proteins (OMPs) of Gram-negative bacteria are almost exclusively β-barrels composed of between 8 and 26 β-strands. To explore the relationship between β-barrel size and shape, we modeled and simulated engineered variants of the Escherichia coli protein OmpX with 8, 10, 12, 14, and 16 β-strands. We found that while smaller barrels maintained a roughly circular shape, the 16-stranded variant developed a flattened cross section. This flat cross section impeded its ability to conduct ions, in agreement with previous experimental observations. Flattening was determined to arise from the presence of inward-facing glycines at sharp turns in the β-barrel. An analysis of all simulations revealed that glycines, on average, make significantly smaller angles with residues on neighboring strands than all other amino acids, including alanine, and create sharp turns in β-barrel cross sections. This observation was generalized to 119 unique structurally resolved OMPs. We also found that the fraction of glycines in β-barrels decreases as the strand number increases, suggesting an evolutionary role for the addition or removal of glycine in OMP sequences.

Electrophysiological Properties from Computations at a Single Voltage: Testing Theory with Stochastic Simulations

We use stochastic simulations to investigate the performance of two recently developed methods for calculating the free energy profiles of ion channels and their electrophysiological properties, such as current–voltage dependence and reversal potential, from molecular dynamics simulations at a single applied voltage. These methods require neither knowledge of the diffusivity nor simulations at multiple voltages, which greatly reduces the computational effort required to probe the electrophysiological properties of ion channels. They can be used to determine the free energy profiles from either forward or backward one-sided properties of ions in the channel, such as ion fluxes, density profiles, committor probabilities, or from their two-sided combination. By generating large sets of stochastic trajectories, which are individually designed to mimic the molecular dynamics crossing statistics of models of channels of trichotoxin, p7 from hepatitis C and a bacterial homolog of the pentameric ligand-gated ion channel, GLIC, we find that the free energy profiles obtained from stochastic simulations corresponding to molecular dynamics simulations of even a modest length are burdened with statistical errors of only 0.3 kcal/mol. Even with many crossing events, applying two-sided formulas substantially reduces statistical errors compared to one-sided formulas. With a properly chosen reference voltage, the current–voltage curves can be reproduced with good accuracy from simulations at a single voltage in a range extending for over 200 mV. If possible, the reference voltages should be chosen not simply to drive a large current in one direction, but to observe crossing events in both directions.

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.

Millisecond-Long Simulations of Antibiotics Transport through Outer Membrane Channels

To reach their target site inside Gram-negative bacteria, almost all antibiotics need to cross the outer membrane. Computational modeling of such processes can be numerically demanding due to the size of the systems and especially due to the timescales involved. Recently, a hybrid Brownian and molecular dynamics approach, i.e., Brownian dynamics including explicit atoms (BRODEA), has been developed and evaluated for studying the transport of monoatomic ions through membrane channels. Later on, this numerically efficient scheme has been applied to determine the free energy surfaces of the ciprofloxacin and enrofloxacin translocation through the porin OmpC using temperature-accelerated simulations. To improve the usability and accuracy of the approach, schemes to approximate the position-dependent diffusion constant of the molecule while traversing the pore had to be established. To this end, we have studied the translocation of the charged phosphonic acid antibiotic fosfomycin through the porin OmpF from Escherichia coli devising and benchmarking several diffusion models. To test the efficiency and sensitivity of these models, the effect of OmpF mutations on the permeation of fosfomycin was analyzed. Permeation events have been recorded over millisecond-long biased and unbiased simulations, from which thermodynamics and kinetics quantities of the translocation processes were determined. As a result, the use of the BRODEA approach, together with the appropriate diffusion model, was seen to accurately reproduce the findings observed in electrophysiology experiments and all-atom molecular dynamics simulations. These results suggest that the BRODEA approach can become a valuable tool for screening numerous compounds to evaluate their outer membrane permeability, a property important in the development of new antibiotics.

Antibiotic uptake through porins located in the outer membrane of Gram-negative bacteria

Introduction: Making selective inhibitors of novel Gram-negative targets is not a substantial challenge - getting them into Gram-negative bacteria to reach their lethal target is the bottleneck. Poor permeability of the antibiotic requires high concentration causing off target activity. The lack of simple experimental techniques to measure antibiotic uptake as well as the local concentration at the target site creates a particular bottleneck in understanding and in improving the antibiotic activity.Areas covered: Here we recall current approaches to quantify the uptake. For a few antibiotics with known evidence for channel-limited permeation, the flux across a single OmpF or OmpC channel has been measured. For a typical concentration gradient of 1 µM of antibiotics the uptake varies between one up to few hundred molecules per second and per channel.Expert opinion: The current research effort is on quantifying the flux for a larger list of compounds on a cellular (mass spectra, fluorescence) or at single channel level (electrophysiology). A larger dataset of single channel permeabilities under various condition will be a powerful tool for understanding and improving the activity of antibiotics.

Ionic transport through a protein nanopore: a Coarse-Grained Molecular Dynamics Study

The MARTINI coarse-grained (CG) force field is used to test the ability of CG models to simulate ionic transport through protein nanopores. The ionic conductivity of CG ions in solution was computed and compared with experimental results. Next, we studied the electrostatic behavior of a solvated CG lipid bilayer in salt solution under an external electric field. We showed this approach correctly describes the experimental conditions under a potential bias. Finally, we performed CG molecular dynamics simulations of the ionic transport through a protein nanopore (α-hemolysin) inserted in a lipid bilayer, under different electric fields, for 2-3 microseconds. The resulting I - V curve is qualitatively consistent with experiments, although the computed current is one order of magnitude smaller. Current saturation was observed for potential biases over ±350 mV. We also discuss the time to reach a stationary regime and the role of the protein flexibility in our CG simulations.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.

Fosfomycin Permeation through the Outer Membrane Porin OmpF

Fosfomycin is a frequently prescribed drug in the treatment of acute urinary tract infections. It enters the bacterial cytoplasm and inhibits the biosynthesis of peptidoglycans by targeting the MurA enzyme. Despite extensive pharmacological studies and clinical use, the permeability of fosfomycin across the bacterial outer membrane is largely unexplored. Here, we investigate the fosfomycin permeability across the outer membrane of Gram-negative bacteria by electrophysiology experiments as well as by all-atom molecular dynamics simulations including free-energy and applied-field techniques. Notably, in an electrophysiological zero-current assay as well as in the molecular simulations, we found that fosfomycin can rapidly permeate the abundant Escherichia coli porin OmpF. Furthermore, two triple mutants in the constriction region of the porin have been investigated. The permeation rates through these mutants are slightly lower than that of the wild type but fosfomycin can still permeate. Altogether, this work unravels molecular details of fosfomycin permeation through the outer membrane porin OmpF of E. coli and moreover provides hints for understanding the translocation of phosphonic acid antibiotics through other outer membrane pores.

Outer Membrane Porins

The transport of small molecules across membranes is essential for the import of nutrients and other energy sources into the cell and, for the export of waste and other potentially harmful byproducts out of the cell. While hydrophobic molecules are permeable to membranes, ions and other small polar molecules require transport via specialized membrane transport proteins . The two major classes of membrane transport proteins are transporters and channels. With our focus here on porins-major class of non-specific diffusion channel proteins , we will highlight some recent structural biology reports and functional assays that have substantially contributed to our understanding of the mechanism that mediates uptake of small molecules, including antibiotics, across the outer membrane of Enterobacteriaceae . We will also review advances in the regulation of porin expression and porin biogenesis and discuss these pathways as new therapeutic targets.

Lipid binding attenuates channel closure of the outer membrane protein OmpF

Outer-membrane porins are often considered as passive conduits of small molecules across lipid bilayers. Using native mass spectrometry experiments we identify a pH-sensitive lipid-binding mechanism of outer membrane porin F, which enables increased threading of a colicin-derived peptide through open channels. Supported by molecular dynamics simulations and channel recording experiments, we posit that this mechanism attenuates channel opening in response to changes in environmental conditions, specifically pH. These findings have important consequences for mass spectrometry experiments, wherein the role of charge is often overlooked, and they also could help provide understanding of antibiotics that gain access to Gram-negative bacteria through porin-mediated pathways.

Strong interactions between lipids and proteins occur primarily through association of charged headgroups and amino acid side chains, rendering the protonation status of both partners important. Here we use native mass spectrometry to explore lipid binding as a function of charge of the outer membrane porin F (OmpF). We find that binding of anionic phosphatidylglycerol (POPG) or zwitterionic phosphatidylcholine (POPC) to OmpF is sensitive to electrospray polarity while the effects of charge are less pronounced for other proteins in outer or mitochondrial membranes: the ferripyoverdine receptor (FpvA) or the voltage-dependent anion channel (VDAC). Only marginal charge-induced differences were observed for inner membrane proteins: the ammonia channel (AmtB) or the mechanosensitive channel. To understand these different sensitivities, we performed an extensive bioinformatics analysis of membrane protein structures and found that OmpF, and to a lesser extent FpvA and VDAC, have atypically high local densities of basic and acidic residues in their lipid headgroup-binding regions. Coarse-grained molecular dynamics simulations, in mixed lipid bilayers, further implicate changes in charge by demonstrating preferential binding of anionic POPG over zwitterionic POPC to protonated OmpF, an effect not observed to the same extent for AmtB. Moreover, electrophysiology and mass-spectrometry–based ligand-binding experiments, at low pH, show that POPG can maintain OmpF channels in open conformations for extended time periods. Since the outer membrane is composed almost entirely of anionic lipopolysaccharide, with similar headgroup properties to POPG, such anionic lipid binding could prevent closure of OmpF channels, thereby increasing access of antibiotics that use porin-mediated pathways.

Enrofloxacin Permeation Pathways across the Porin OmpC

In Gram-negative bacteria, the lack or quenching of antibiotic translocation across the outer membrane is one of the main factors for acquiring antibiotic resistance. An atomic-level comprehension of the key features governing the transport of drugs by outer-membrane protein channels would be very helpful in developing the next generation of antibiotics. In a previous study [ J. D. Prajapati et al. J. Chem. Theory Comput. 2017 , 13 , 4553 ], we characterized the diffusion pathway of a ciprofloxacin molecule through the outer membrane porin OmpC of Escherichia coli by combining metadynamics and a zero-temperature string method. Here, we evaluate the diffusion route through the OmpC porin for a similar fluoroquinolone, that is, the enrofloxacin molecule, using the previously developed protocol. As a result, it was found that the lowest-energy pathway was similar to that for ciprofloxacin; namely, a reorientation was required on the extracellular side with the carboxyl group ahead before enrofloxacin reached the constriction region. In turn, the free-energy basins for both antibiotics are located at similar positions in the space defined by selected reaction coordinates, and their affinity sites share a wide number of porin residues. However, there are some important deviations due to the chemical differences of these two drugs. On the one hand, a slower diffusion process is expected for enrofloxacin, as the permeation pathway exhibits higher overall energy barriers, mainly in the constriction region. On the other hand, enrofloxacin needs to replace some polar interactions in its affinity sites with nonpolar ones. This study demonstrates how minor chemical modifications can qualitatively affect the translocation mechanism of an antibiotic molecule.

Label-Free Single-Molecule Thermoscopy Using a Laser-Heated Nanopore

When light is used to excite electronic transitions in a material, nonradiative energy during relaxation is often released in the form of heat. In this work, we show that photoexcitation of a silicon nitride nanopore using a focused visible laser results in efficient localized photothermal heating, which reduces the nearby electrolyte viscosity and increases the ionic conductance. In addition, a strong localized thermal gradient in the pore vicinity is produced, evidenced by finite-element simulations and experimental observation of both ion and DNA thermophoresis. After correcting for thermophoresis, the nanopore current can be used as a nanoscale thermometer, enabling rapid force thermoscopy. We utilize this to probe thermal melting transitions in synthetic and native biomolecules that are heated at the nanopore. Our results on single molecules are validated by correspondence to bulk measurements, which paves the way to various biophysical experiments that require rapid temperature and force control on individual molecules.

Conversion of OprO into an OprP-like Channel by Exchanging Key Residues in the Channel Constriction

Under phosphate-limiting conditions, the channels OprP and OprO are induced and expressed in the outer membrane of Pseudomonas aeruginosa. Despite their large homology, the phosphate-specific OprP and the diphosphate-specific OprO pores show structural differences in their binding sites situated in the constriction region. Previously, it was shown that the mutation of amino acids in OprP (Y62F and Y114D) led to an exchange in substrate specificity similar to OprO. To support the role of these key amino acids in the substrate sorting of these specific channels, the reverse mutants for OprO (F62Y, D114Y, and F62Y/D114Y) were created in this study. The phosphate and diphosphate binding of the generated channels was studied in planar lipid bilayers. Our results show that mutations of key residues indeed reverse the substrate specificity of OprO to OprP and support the view that just a few strategically positioned amino acids are mainly responsible for its substrate specificity.

Exploring bacterial outer membrane barrier to combat bad bugs

One of the main fundamental mechanisms of antibiotic resistance in Gram-negative bacteria comprises an effective change in the membrane permeability to antibiotics. The Gram-negative bacterial complex cell envelope comprises an outer membrane that delimits the periplasm from the exterior environment. The outer membrane contains numerous protein channels, termed as porins or nanopores, which are mainly involved in the influx of hydrophilic compounds, including antibiotics. Bacterial adaptation to reduce influx through these outer membrane proteins (Omps) is one of the crucial mechanisms behind antibiotic resistance. Thus to interpret the molecular basis of the outer membrane permeability is the current challenge. This review attempts to develop a state of knowledge pertinent to Omps and their effective role in antibiotic influx. Further, it aims to study the bacterial response to antibiotic membrane permeability and hopefully provoke a discussion toward understanding and further exploration of prospects to improve our knowledge on physicochemical parameters that direct the translocation of antibiotics through the bacterial membrane protein channels.

Role of Electroosmosis in the Permeation of Neutral Molecules: CymA and Cyclodextrin as an Example

To quantify the flow of small uncharged molecules into and across nanopores, one often uses ion currents. The respective ion-current fluctuations caused by the presence of the analyte make it possible to draw some conclusions about the direction and magnitude of the analyte flow. However, often this flow appears to be asymmetric with respect to the applied voltage. As a possible reason for this asymmetry, we identified the electroosmotic flow (EOF), which is the water transport associated with ions driven by the external transmembrane voltage. As an example, we quantify the contribution of the EOF through a nanopore by investigating the permeation of α-cyclodextrin through CymA, a cyclodextrin-specific channel from Klebsiella oxytoca. To understand the results from electrophysiology on a molecular level, all-atom molecular dynamics simulations are used to detail the effect of the EOF on substrate entry to and exit from a CymA channel in which the N-terminus has been deleted. The combined experimental and computational results strongly suggest that one needs to account for the significant contribution of the EOF when analyzing the penetration of cyclodextrins through the CymA pore. This example study at the same time points to the more general finding that the EOF needs to be considered in translocation studies of neutral molecules and, at least in many cases, should be able to help in discriminating between translocation and binding events.

Dynamics and Interactions of OmpF and LPS: Influence on Pore Accessibility and Ion Permeability

The asymmetric outer membrane of Gram-negative bacteria is formed of the inner leaflet with phospholipids and the outer leaflet with lipopolysaccharides (LPS). Outer membrane protein F (OmpF) is a trimeric porin responsible for the passive transport of small molecules across the outer membrane of Escherichia coli. Here, we report the impact of different levels of heterogeneity in LPS environments on the structure and dynamics of OmpF using all-atom molecular dynamics simulations. The simulations provide insight into the flexibility and dynamics of LPS components that are highly dependent on local environments, with lipid A being the most rigid and O-antigen being the most flexible. Increased flexibility of O-antigen polysaccharides is observed in heterogeneous LPS systems, where the adjacent O-antigen repeating units are weakly interacting and thus more dynamic, compared to homogeneous LPS systems in which LPS interacts strongly with each other with limited overall flexibility due to dense packing. The model systems were validated by comparing molecular-level details of interactions between OmpF surface residues and LPS core sugars with experimental data, establishing the importance of LPS core oligosaccharides in shielding OmpF surface epitopes recognized by monoclonal antibodies. There are LPS environmental influences on the movement of bulk ions (K(+) and Cl(-)), but the ion selectivity of OmpF is mainly affected by bulk ion concentration.

Validity of the Electrodiffusion Model for Calculating Conductance of Simple Ion Channels

We examine the validity and utility of the electrodiffusion (ED) equation, i.e., the generalized Nernst-Planck equation, to characterize, in combination with molecular dynamics, the electrophysiological behavior of simple ion channels. As models, we consider three systems-two naturally occurring channels formed by α-helical bundles of peptaibols, trichotoxin, and alamethicin, and a synthetic, hexameric channel, formed by a peptide that contains only leucine and serine. All these channels mediate transport of potassium and chloride ions. Starting with equilibrium properties, such as the potential of mean force experienced by an ion traversing the channel and diffusivity, obtained from molecular dynamics simulations, the ED equation can be used to determine the full current-voltage dependence with modest or no additional effort. The potential of mean force can be obtained not only from equilibrium simulations, but also, with comparable accuracy, from nonequilibrium simulations at a single voltage. The main assumptions underlying the ED equation appear to hold well for the channels and voltages studied here. To expand the utility of the ED equation, we examine what are the necessary and sufficient conditions for Ohmic and nonrectifying behavior and relate deviations from this behavior to the shape of the ionic potential of mean force.

Simulating Current-Voltage Relationships for a Narrow Ion Channel Using the Weighted Ensemble Method

Ion channels are responsible for a myriad of fundamental biological processes via their role in controlling the flow of ions through water-filled membrane-spanning pores in response to environmental cues. Molecular simulation has played an important role in elucidating the mechanism of ion conduction, but connecting atomistically detailed structural models of the protein to electrophysiological measurements remains a broad challenge due to the computational cost of reaching the necessary time scales. Here, we introduce an enhanced sampling method for simulating the conduction properties of narrow ion channels using the Weighted ensemble (WE) sampling approach. We demonstrate the application of this method to calculate the current–voltage relationship as well as the nonequilibrium ion distribution at steady-state of a simple model ion channel. By direct comparisons with long brute force simulations, we show that the WE simulations rigorously reproduce the correct long-time scale kinetics of the system and are capable of determining these quantities using significantly less aggregate simulation time under conditions where permeation events are rare.

What controls open-pore and residual currents in the first sensing zone of alpha-hemolysin nanopore? Combined experimental and theoretical study

The electrophoretic transport of single-stranded DNA through biological nanopores such as alpha-hemolysin (αHL) is a promising and cost-effective technology with the potential to revolutionize genomics. The rational design of pores with the controlled polymer translocation rates and high contrast between different nucleotides could improve significantly nanopore sequencing applications. Here, we apply a combination of theoretical and experimental methods in an attempt to elucidate several selective modifications in the pore which were proposed to be central for the effective discrimination between purines and pyrimidines. Our nanopore test set includes the wild type αHL and six mutants (E111N/M113X/K147N) in which the cross-section and chemical functionality of the first constriction zone of the pore are modified. Electrophysiological recordings were combined with all-atom Molecular Dynamics simulations (MD) and a recently developed Brownian Dynamics (BROMOC) protocol to investigate residual ion currents and pore-DNA interactions for two homo-polymers e.g. poly(dA)40 or poly(dC)40 blocking the pore. The calculated residual currents and contrast in the poly(dA)40/poly(dC)40 blocked pore are in qualitative agreement with the experimental recordings. We showed that a simple structural metric allows rationalization of key elements in the emergent contrast between purines and pyrimidines in the modified αHL mutants. The shape of the pore and its capacity for hydrogen bonding to a translocated polynucleotide are two essential parameters for contrast optimization. To further probe the impact of these two factors in the ssDNA sensing, we eliminated the effect of the primary constriction using serine substitutions (i.e. E111S/M113S/T145S/K147S) and increased the hydrophobic volume of the central residue in the secondary constriction (L135I). This pore modification sharply increased the contrast between Adenine (A) and Cytosine (C).

In Silico Structure and Sequence Analysis of Bacterial Porins and Specific Diffusion Channels for Hydrophilic Molecules: Conservation, Multimericity and Multifunctionality

Diffusion channels are involved in the selective uptake of nutrients and form the largest outer membrane protein (OMP) family in Gram-negative bacteria. Differences in pore size and amino acid composition contribute to the specificity. Structure-based multiple sequence alignments shed light on the structure-function relations for all eight subclasses. Entropy-variability analysis results are correlated to known structural and functional aspects, such as structural integrity, multimericity, specificity and biological niche adaptation. The high mutation rate in their surface-exposed loops is likely an important mechanism for host immune system evasion. Multiple sequence alignments for each subclass revealed conserved residue positions that are involved in substrate recognition and specificity. An analysis of monomeric protein channels revealed particular sequence patterns of amino acids that were observed in other classes at multimeric interfaces. This adds to the emerging evidence that all members of the family exist in a multimeric state. Our findings are important for understanding the role of members of this family in a wide range of bacterial processes, including bacterial food uptake, survival and adaptation mechanisms.

Characteristics of Sucrose Transport through the Sucrose-Specific Porin ScrY Studied by Molecular Dynamics Simulations

Sucrose-specific porin (ScrY) is a transmembrane protein that allows for the uptake of sucrose under growth-limiting conditions. The crystal structure of ScrY was resolved before by X-ray crystallography, both in its uncomplexed form and with bound sucrose. However, little is known about the molecular characteristics of the transport mechanism of ScrY. To date, there has not yet been any clear demonstration for sucrose transport through the ScrY. Here, the dynamics of the ScrY trimer embedded in a phospholipid bilayer as well as the characteristics of sucrose translocation were investigated by means of atomistic molecular dynamics (MD) simulations. The potential of mean force (PMF) for sucrose translocation through the pore showed two main energy barriers within the constriction region of ScrY. Energy decomposition allowed to pinpoint three aspartic acids as key residues opposing the passage of sucrose, all located within the L3 loop. Mutation of two aspartic acids to uncharged residues resulted in an accordingly modified electrostatics and decreased PMF barrier. The chosen methodology and results will aid in the design of porins with modified transport specificities.

A computational study of ion current modulation in hVDAC3 induced by disulfide bonds

The human VDAC channel exists in three isoforms characterized by high sequence homology and structural similarity. Yet the function and mode of action of hVDAC3 are still elusive. The presence of six surface cysteines exposed to the oxidizing environment of the mitochondrial inter-membrane space suggests the possible establishment of intramolecular disulfide bonds. Two natural candidates for disulfide bridge formation are Cys2 and Cys8 that, located on the flexible N-terminal domain, can easily come in contact. A third potentially important residue is Cys122 that is close to Cys2 in the homology model of VDAC3. Here we analyzed the impact of SS bonds through molecular dynamics simulations of derivatives of hVDAC3 (dubbed SS-2-8, SS-2-122, SS-8-122) including a single disulfide bond. Simulations showed that in SS-8-122, the fragment 1-7 crosses the top part of the barrel partially occluding the pore and causing a 20% drop of conductance. In order to identify other potential channel-occluding disulfide bonds, we used a set of neural networks and structural bioinformatics algorithms, after filtering with the steric constraints imposed by the 3D-structure. We identified other three species, namely SS-8-65, SS-2-36 and SS-8-36. While the conductance of SS-8-65 and SS-2-36 is about 30% lower than that of the species without disulfide bonds, the conductance of SS-8-36 was 40-50% lower. The results show how VDAC3 is able to modulate its pore size and current by exploiting the mobility of the N-terminal and forming, upon external stimuli, disulfide bridges with cysteine residues located on the barrel and exposed to the inter-membrane space.

Simulations of outer membrane channels and their permeability

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

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.

Physical methods to quantify small antibiotic molecules uptake into Gram-negative bacteria

The development of antibiotics against Gram-negative bacteria is a challenge: any active compound must cross the outer cell envelope composed of a hydrophilic highly charged lipopolysaccharide layer followed by a tight hydrophobic layer containing water filled gates called porins to reach the hydrophilic periplasmic space and depending on the target with the further need to cross the hydrophobic inner membrane. In addition to a possible rapid enzymatic deactivation efflux pumps shuffle compounds back outside. The resulting low permeability of cell envelope requires high dose and leads therefore to toxicity problems. Despite its relevance the permeability barrier in Gram-negative bacteria is not well understood partially caused by the lack of appropriate direct assays. Here we give a brief introduction on current available techniques to quantify passive diffusion of small hydrophilic molecules into Gram-negative bacteria.

Understanding Voltage Gating of Providencia stuartii Porins at Atomic Level

Bacterial porins are water-filled β-barrel channels that allow translocation of solutes across the outer membrane. They feature a constriction zone, contributed by the plunging of extracellular loop 3 (L3) into the channel lumen. Porins are generally in the open state, but undergo gating in response to external voltages. To date the underlying mechanism is unclear. Here we report results from molecular dynamics simulations on the two porins of Providenica stuartii, Omp-Pst1 and Omp-Pst2, which display distinct voltage sensitivities. Voltage gating was observed in Omp-Pst2, where the binding of cations in-between L3 and the barrel wall results in exposing a conserved aromatic residue in the channel lumen, thereby halting ion permeation. Comparison of Omp-Pst1 and Omp-Pst2 structures and trajectories suggests that their sensitivity to voltage is encoded in the hydrogen-bonding network anchoring L3 onto the barrel wall, as we observed that it is the strength of this network that governs the probability of cations binding behind L3. That Omp-Pst2 gating is observed only when ions flow against the electrostatic potential gradient of the channel furthermore suggests a possible role for this porin in the regulation of charge distribution across the outer membrane and bacterial homeostasis.

Porins are the main conduits for hydrophilic nutrients and ions uptake into the periplasm of Gram-negative bacteria. Their translocation permeability is determined by the amino-acid distribution on their extracellular loop L3. Bacterial porin channels have long been known to undergo step-wise gating, under the application of a transmembrane potential. Yet the exact molecular mechanism by which gating is achieved and the exact relevance of this evolved characteristic remain elusive. In the present study, we report on electrophysiology experiments and molecular dynamics simulations on the two general-diffusion porins of Providencia stuartii, Omp-Pst1 and Omp-Pst2. Our results show that gating in Omp-Pst2 occurs as the result of L3 displacement, which follows from the binding of cations in acidic niches between L3 and the barrel wall and effects in exposing the side chain of a highly conserved aromatic residue at the tip of L3 in the channel lumen. That Omp-Pst2 displays asymmetric voltage sensitivity and that the likelihood of gating is increased when cations transit from the extracellular to the intracellular side suggests voltage-gating underlies a regulatory role in bacterial homeostasis. Rational antibiotic-design strategies based on the maximization of antibiotic penetration and accumulation at their target sites, should take this role into account.

Breaking the hydrophobicity of the MscL pore: insights into a charge-induced gating mechanism

The mechanosensitive channel of large conductance (MscL) is a protein that responds to membrane tension by opening a transient pore during osmotic downshock. Due to its large pore size and functional reconstitution into lipid membranes, MscL has been proposed as a promising artificial nanovalve suitable for biotechnological applications. For example, site-specific mutations and tailored chemical modifications have shown how MscL channel gating can be triggered in the absence of tension by introducing charged residues at the hydrophobic pore level. Recently, engineered MscL proteins responsive to stimuli like pH or light have been reported. Inspired by experiments, we present a thorough computational study aiming at describing, with atomistic detail, the artificial gating mechanism and the molecular transport properties of a light-actuated bacterial MscL channel, in which a charge-induced gating mechanism has been enabled through the selective cleavage of photo-sensitive alkylating agents. Properties such as structural transitions, pore dimension, ion flux and selectivity have been carefully analyzed. Besides, the effects of charge on alternative sites of the channel with respect to those already reported have been addressed. Overall, our results provide useful molecular insights into the structural events accompanying the engineered MscL channel gating and the interplay of electrostatic effects, channel opening and permeation properties. In addition, we describe how the experimentally observed ionic current in a single-subunit charged MscL mutant is obtained through a hydrophobicity breaking mechanism involving an asymmetric inter-subunit motion.

Theoretical analysis of ion conductance and gating transitions in the OpdK (OccK1) channel

Molecular simulations have been performed on the pore OpdK elucidating molecular details of ion conductance and a possible gating mechanism.

Entropy-enthalpy compensation at the single protein level: pH sensing in the bacterial channel OmpF

The pH sensing mechanism of the OmpF channel operates via ligand modification: increasing acidity induces the replacement of cations with protons in critical binding sites decreasing the channel conductance. Aside from the change in enthalpy associated with the binding, there is also a change in the microscopic arrangements of ligands, receptors and the surrounding solvent. We show that the pH-modulation of the single channel conduction involves small free energy changes because large enthalpic and entropic contributions change in opposite ways, demonstrating an approximate enthalpy-entropy compensation for different salts and concentrations.

VCA1008: An Anion-Selective Porin of Vibrio Cholerae

A putative porin function has been assigned to VCA1008 of Vibrio cholerae. Its coding gene, vca1008, is expressed upon colonization of the small intestine in infant mice and human volunteers, and is essential for infection. In vitro, vca1008 is expressed under inorganic phosphate limitation and, in this condition, VCA1008 is the major outer membrane protein of the bacterium. Here, we provide the first functional characterization of VCA1008 reconstituted into planar lipid bilayers. Our main findings were: 1) VCA1008 forms an ion channel that, at high voltage (~±100 mV), presents a voltage-dependent activity and displays closures typical of trimeric porins, with a conductance of 4.28±0.04 nS (n=164) in 1M KCl; 2) It has a preferred selectivity for anions over cations; 3) Its conductance saturates with increasing inorganic phosphate concentration, suggesting VCA1008 contains binding site(s) for this anion; 4) Its ion selectivity is controlled by both fixed charged residues within the channel and diffusion along the pore; 5) Partitioning of poly (ethylene glycol)s (PEGs) of different molecular mass suggests that VCA1008 channel has a pore exclusion limit of 0.9 nm.

A structural study of ion permeation in OmpF porin from anomalous X-ray diffraction and molecular dynamics simulations

OmpF, a multiionic porin from Escherichia coli, is a useful protypical model system for addressing general questions about electrostatic interactions in the confinement of an aqueous molecular pore. Here, favorable anion locations in the OmpF pore were mapped by anomalous X-ray scattering of Br(–) ions from four different crystal structures and compared with Mg(2+) sites and Rb(+) sites from a previous anomalous diffraction study to provide a complete picture of cation and anion transfer paths along the OmpF channel. By comparing structures with various crystallization conditions, we find that anions bind in discrete clusters along the entire length of the OmpF pore, whereas cations find conserved binding sites with the extracellular, surface-exposed loops. Results from molecular dynamics simulations are consistent with the experimental data and help highlight the critical residues that preferentially contact either cations or anions during permeation. Analysis of these results provides new insights into the molecular mechanisms that determine ion selectivity in OmpF porin.

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations

A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific β-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26–27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K⁺ channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K⁺ and Na⁺ concentration gradients and the electrostatic potential difference of living cells.

Stable closure of the cytoplasmic half-channel is required for efficient proton transport at physiological membrane potentials in the bacteriorhodopsin catalytic cycle

The bacteriorhodopsin (BR) Asp96Gly/Phe171Cys/Phe219Leu triple mutant has been shown to translocate protons 66% as efficiently as the wild-type protein. Light-dependent ATP synthesis in haloarchaeal cells expressing the triple mutant is 85% that of the wild-type BR expressing cells. Therefore, the functional activity of BR seems to be largely preserved in the triple mutant despite the observations that its ground-state structure resembles that of the wild-type M state (i.e., the so-called cytoplasmically open state) and that the mutant shows no significant structural changes during its photocycle, in sharp contrast to what occurs in the wild-type protein in which a large structural opening and closing occurs on the cytoplasmic side. To resolve the contradiction between the apparent functional robustness of the triple mutant and the presumed importance of the opening and closing that occurs in the wild-type protein, we conducted additional experiments to compare the behavior of wild-type and mutant proteins under different operational loads. Specifically, we characterized the ability of the two proteins to generate light-driven proton currents against a range of membrane potentials. The wild-type protein showed maximal conductance between −150 and −50 mV, whereas the mutant showed maximal conductance at membrane potentials >+50 mV. Molecular dynamics (MD) simulations of the triple mutant were also conducted to characterize structural changes in the protein and in solvent accessibility that might help to functionally contextualize the current–voltage data. These simulations revealed that the cytoplasmic half-channel of the triple mutant is constitutively open and dynamically exchanges water with the bulk. Collectively, the data and simulations help to explain why this mutant BR does not mediate photosynthetic growth of haloarchaeal cells, and they suggest that the structural closing observed in the wild-type protein likely plays a key role in minimizing substrate back flow in the face of electrochemical driving forces present at physiological membrane potentials.

A Stabilized Finite Element Method for Modified Poisson-Nernst-Planck Equations to Determine Ion Flow Through a Nanopore

The conventional Poisson-Nernst-Planck equations do not account for the finite size of ions explicitly. This leads to solutions featuring unrealistically high ionic concentrations in the regions subject to external potentials, in particular, near highly charged surfaces. A modified form of the Poisson-Nernst-Planck equations accounts for steric effects and results in solutions with finite ion concentrations. Here, we evaluate numerical methods for solving the modified Poisson-Nernst-Planck equations by modeling electric field-driven transport of ions through a nanopore. We describe a novel, robust finite element solver that combines the applications of the Newton's method to the nonlinear Galerkin form of the equations, augmented with stabilization terms to appropriately handle the drift-diffusion processes. To make direct comparison with particle-based simulations possible, our method is specifically designed to produce solutions under periodic boundary conditions and to conserve the number of ions in the solution domain. We test our finite element solver on a set of challenging numerical experiments that include calculations of the ion distribution in a volume confined between two charged plates, calculations of the ionic current though a nanopore subject to an external electric field, and modeling the effect of a DNA molecule on the ion concentration and nanopore current.