Ion permeation through the alpha-hemolysin channel: theoretical studies based on Brownian dynamics and Poisson-Nernst-Plank electrodiffusion theory

Exploring ssDNA translocation through α-hemolysin using coarse-grained steered molecular dynamics

Protein nanopores have proven to be effective for single-molecule studies, particularly for single-stranded DNA (ssDNA) translocation. Previous experiments demonstrated their ability to distinguish differences in purine and pyrimidine bases and in the orientation of the ssDNA molecule inside nanopores. Unfortunately, the microscopic details of ssDNA translocation over experimental time scales, which are not accessible through all-atom molecular dynamics (MD), have yet to be examined. However, coarse-grained (CG) MD simulations enable systems to be simulated over longer characteristic times closer to experiments than all-atom MD. This paper studies ssDNA translocation through α-hemolysin nanopores exploiting steered MD using the MARTINI CG force field. The impacts of the sequence length, orientation inside the nanopore and DNA charges on translocation dynamics as well as the conformational dynamics of ssDNA during the translocation are explored. Our results highlight the efficacy of CG molecular dynamics in capturing the experimental properties of ssDNA translocation, including a wide distribution in translocation times per base. In particular, the phosphate charges of the DNA molecule are crucial in the translocation dynamics and impact the translocation rate. Additionally, the influence of the ssDNA molecule orientation on the translocation rate is explained by the conformational differences of ssDNA inside the nanopore during its translocation. Our study emphasizes the significance of obtaining sufficient statistics via CG MD, which can elucidate the great variety of translocation processes.

Electro-osmotic Flow Generation via a Sticky Ion Action

Selective transport of ions through nanometer-sized pores is fundamental to cell biology and central to many technological processes such as water desalination and electrical energy storage. Conventional methods for generating ion selectivity include placement of fixed electrical charges at the inner surface of a nanopore through either point mutations in a protein pore or chemical treatment of a solid-state nanopore surface, with each nanopore type requiring a custom approach. Here, we describe a general method for transforming a nanoscale pore into a highly selective, anion-conducting channel capable of generating a giant electro-osmotic effect. Our molecular dynamics simulations and reverse potential measurements show that exposure of a biological nanopore to high concentrations of guanidinium chloride renders the nanopore surface positively charged due to transient binding of guanidinium cations to the protein surface. A comparison of four biological nanopores reveals the relationship between ion selectivity, nanopore shape, composition of the nanopore surface, and electro-osmotic flow. Guanidinium ions are also found to produce anion selectivity and a giant electro-osmotic flow in solid-state nanopores via the same mechanism. Our sticky-ion approach to generate electro-osmotic flow can have numerous applications in controlling molecular transport at the nanoscale and for detection, identification, and sequencing of individual proteins.

Dynamics accelerate the kinetics of ion diffusion through channels: Continuous-time random walk models beyond the mean field approximation

Ion channels are proteins that play a significant role in physiological processes, including neuronal excitability and signal transduction. However, the precise mechanisms by which these proteins facilitate ion diffusion through cell membranes are not well understood. This is because experimental techniques to characterize ion channel activity operate on a time scale too large to understand the role of the various protein conformations on diffusion. Meanwhile, computational approaches operate on a time scale too short to rationalize the observed behavior at the microscopic scale. In this paper, we present a continuous-time random walk model that aims to bridge the scales between the atomistic models of ion channels and the experimental measurement of their conductance. We show how diffusion slows down in complex systems by using 3D lattices that map out the pore geometry of two channels: Nav1.7 and gramicidin. We also introduce spatial and dynamic site disorder to account for system heterogeneity beyond the mean field approximation. Computed diffusion coefficients show that an increase in spatial disorder slows down diffusion kinetics, while dynamic disorder has the opposite effect. Our results imply that microscopic or phenomenological models based on the potential of mean force data overlook the functional importance of protein dynamics on ion diffusion through channels.

Electro-Osmotic Flow Generation via a Sticky Ion Action

Selective transport of ions through nanometer-sized pores is fundamental to cell biology and central to many technological processes such as water desalination and electrical energy storage. Conventional methods for generating ion selectivity include placement of fixed electrical charges at the inner surface of a nanopore through either point mutations in a protein pore or chemical treatment of a solid-state nanopore surface, with each nanopore type requiring a custom approach. Here, we describe a general method for transforming a nanoscale pore into a highly selective, anion-conducting channel capable of generating a giant electro-osmotic effect. Our molecular dynamics simulations and reverse potential measurements show that exposure of a biological nanopore to high concentrations of guanidinium chloride renders the nanopore surface positively charged due to transient binding of guanidinium cations to the protein surface. A comparison of four biological nanopores reveals the relationship between ion selectivity, nanopore shape, composition of the nanopore surface, and electro-osmotic flow. Remarkably, guanidinium ions are also found to produce anion selectivity and a giant electro-osmotic flow in solid-state nanopores via the same mechanism. Our sticky-ion approach to generate electro-osmotic flow can have numerous applications in controlling molecular transport at the nanoscale and for detection, identification, and sequencing of individual proteins.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

Translocation of Hydrophobic Polyelectrolytes under Electrical Field: Molecular Dynamics Study

We studied the translocation of polyelectrolyte (PE) chains driven by an electric field through a pore by means of molecular dynamics simulations of a coarse-grained HP model mimicking high salt conditions. Charged monomers were considered as polar (P) and neutral monomers as hydrophobic (H). We considered PE sequences that had equally spaced charges along the hydrophobic backbone. Hydrophobic PEs were in the globular form in which H-type and P-type monomers were partially segregated and they unfolded in order to translocate through the narrow channel under the electric field. We provided a quantitative comprehensive study of the interplay between translocation through a realistic pore and globule unraveling. By means of molecular dynamics simulations, incorporating realistic force fields inside the channel, we investigated the translocation dynamics of PEs at various solvent conditions. Starting from the captured conformations, we obtained distributions of waiting times and drift times at various solvent conditions. The shortest translocation time was observed for the slightly poor solvent. The minimum was rather shallow, and the translocation time was almost constant for medium hydrophobicity. The dynamics were controlled not only by the friction of the channel, but also by the internal friction related to the uncoiling of the heterogeneous globule. The latter can be rationalized by slow monomer relaxation in the dense phase. The results were compared with those from a simplified Fokker-Planck equation for the position of the head monomer.

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.

Current Rectification and Ionic Selectivity of α-Hemolysin: Coarse-Grained Molecular Dynamics Simulations

In order to understand the physical processes of nanopore experiments at the molecular level, microscopic information from molecular dynamics is greatly needed. Coarse-grained models are a good alternative to classical all-atom models since they allow longer and faster simulations. We performed coarse-grained molecular dynamics of the ionic transport through the α-hemolysin protein nanopore, inserted into a lipid bilayer surrounded by solvent and ions. For this purpose, we used the MARTINI coarse-grained force field and its polarizable water solvent (PW). Moreover, the electric potential difference applied experimentally was mimicked by the application of an electric field to the system. We present, in this study, the results of 1.5 μs long-molecular dynamics simulations of 12 different systems for which different charged amino acids were neutralized, each of them in the presence of nine different electric fields ranging between ±0.04 V/nm (a total of around 100 simulations). We were able to observe several specific features of this pore, current asymmetry and anion selectivity, in agreement with previous studies and experiments, and we identified the charged amino acids responsible for these current behaviors, therefore validating our coarse-grain approach to study ionic transport through nanopores. We also propose a microscopic explanation of these ionic current features using ionic density maps.

Translocation, Rejection and Trapping of Polyampholytes

Polyampholytes (PA) are a special class of polymers comprising both positive and negative monomers along their sequence. Most proteins have positive and negative residues and are PAs. Proteins have a well-defined sequence while synthetic PAs have a random charge sequence. We investigated the translocation behavior of random polyampholyte chains through a pore under the action of an electric field by means of Monte Carlo simulations. The simulations incorporated a realistic translocation potential profile along an extended asymmetric pore and translocation was studied for both directions of engagement. The study was conducted from the perspective of statistics for disordered systems. The translocation behavior (translocation vs. rejection) was recorded for all 220 sequences comprised of N = 20 charged monomers. The results were compared with those for 107 random sequences of N = 40 to better demonstrate asymptotic laws. At early times, rejection was mainly controlled by the charge sequence of the head part, but late translocation/rejection was governed by the escape from a trapped state over an antagonistic barrier built up along the sequence. The probability distribution of translocation times from all successful attempts revealed a power-law tail. At finite times, there was a population of trapped sequences that relaxed very slowly (logarithmically) with time. If a subensemble of sequences with prescribed net charge was considered the power-law decay was steeper for a more favorable net charge. Our findings were rationalized by theoretical arguments developed for long chains. We also provided operational criteria for the translocation behavior of a sequence, explaining the selection by the translocation process. From the perspective of protein translocation, our findings can help rationalize the behavior of intrinsically disordered proteins (IDPs), which can be modeled as polyampholytes. Most IDP sequences have a strong net charge favoring translocation. Even for sequences with those large net charges, the translocation times remained very dispersed and the translocation was highly sequence-selective.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Protein Transport through Nanopores Illuminated by Long-Time-Scale Simulations

The transport of molecules through nanoscale confined space is relevant in biology, biosensing, and industrial filtration. Microscopically modeling transport through nanopores is required for a fundamental understanding and guiding engineering, but the short duration and low replica number of existing simulation approaches limit statistically relevant insight. Here we explore protein transport in nanopores with a high-throughput computational method that realistically simulates hundreds of up to seconds-long protein trajectories by combining Brownian dynamics and continuum simulation and integrating both driving forces of electroosmosis and electrophoresis. Ionic current traces are computed to enable experimental comparison. By examining three biological and synthetic nanopores, our study answers questions about the kinetics and mechanism of protein transport and additionally reveals insight that is inaccessible from experiments yet relevant for pore design. The discovery of extremely frequent unhindered passage can guide the improvement of biosensor pores to enhance desired biomolecular recognition by pore-tethered receptors. Similarly, experimentally invisible nontarget adsorption to pore walls highlights how to improve recently developed DNA nanopores. Our work can be expanded to pressure-driven flow to model industrial nanofiltration processes.

Constructing ion channels from water-soluble α-helical barrels

The design of peptides that assemble in membranes to form functional ion channels is challenging. Specifically, hydrophobic interactions must be designed between the peptides and at the peptide-lipid interfaces simultaneously. Here, we take a multi-step approach towards this problem. First, we use rational de novo design to generate water-soluble α-helical barrels with polar interiors, and confirm their structures using high-resolution X-ray crystallography. These α-helical barrels have water-filled lumens like those of transmembrane channels. Then, we modify the sequences to facilitate their insertion into lipid bilayers. Single-channel electrical recordings and fluorescent imaging of the peptides in membranes show monodisperse, cation-selective channels of unitary conductance. Surprisingly, however, an X-ray structure solved from lipidic cubic phase for one peptide reveals an alternative state with tightly packed helices and a constricted channel. To reconcile these observations, we perform computational analyses to compare the properties of possible different states of the peptide.

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.

Laser-based temperature control to study the roles of entropy and enthalpy in polymer-nanopore interactions

Single-molecule approaches for probing the free energy of confinement for polymers in a nanopore environment are critical for the development of nanopore biosensors. We developed a laser-based nanopore heating approach to monitor the free energy profiles of such a single-molecule sensor. Using this approach, we measure the free energy profiles of two distinct polymers, polyethylene glycol and water-soluble peptides, as they interact with the nanopore sensor. Polyethylene glycol demonstrates a retention mechanism dominated by entropy with little sign of interaction with the pore, while peptides show an enthalpic mechanism, which can be attributed to physisorption to the nanopore (e.g., hydrogen bonding). To manipulate the energetics, we introduced thiolate-capped gold clusters [Au₂₅(SG)₁₈] into the pore, which increases the charge and leads to additional electrostatic interactions that help dissect the contribution that enthalpy and entropy make in this modified environment. These observations provide a benchmark for optimization of single-molecule nanopore sensors.

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.

Accurate modeling of a biological nanopore with an extended continuum framework

Despite the broad success of biological nanopores as powerful instruments for the analysis of proteins and nucleic acids at the single-molecule level, a fast simulation methodology to accurately model their nanofluidic properties is currently unavailable. This limits the rational engineering of nanopore traits and makes the unambiguous interpretation of experimental results challenging. Here, we present a continuum approach that can faithfully reproduce the experimentally measured ionic conductance of the biological nanopore Cytolysin A (ClyA) over a wide range of ionic strengths and bias potentials. Our model consists of the extended Poisson-Nernst-Planck and Navier-Stokes (ePNP-NS) equations and a computationally efficient 2D-axisymmetric representation for the geometry and charge distribution of the nanopore. Importantly, the ePNP-NS equations achieve this accuracy by self-consistently considering the finite size of the ions and the influence of both the ionic strength and the nanoscopic scale of the pore on the local properties of the electrolyte. These comprise the mobility and diffusivity of the ions, and the density, viscosity and relative permittivity of the solvent. Crucially, by applying our methodology to ClyA, a biological nanopore used for single-molecule enzymology studies, we could directly quantify several nanofluidic characteristics difficult to determine experimentally. These include the ion selectivity, the ion concentration distributions, the electrostatic potential landscape, the magnitude of the electro-osmotic flow field, and the internal pressure distribution. Hence, this work provides a means to obtain fundamental new insights into the nanofluidic properties of biological nanopores and paves the way towards their rational engineering.

Dynamics of ionic flows via Poisson-Nernst-Planck systems with local hard-sphere potentials: Competition between cations

We study a quasi-one-dimensional steady-state Poisson-Nernst-Planck type model for ionic flows through a membrane channel with three ion species, two positively charged with the same valence and one negatively charged. Bikerman's local hard-sphere potential is included in the model to account for ion sizes. The problem is treated as a boundary value problem of a singularly perturbed differential system. Under the framework of a geometric singular perturbation theory, together with specific structures of this concrete model, the existence and uniqueness of solutions to the boundary value problem for small ion sizes is established. Furthermore, treating the ion sizes as small parameters, we derive an approximation of individual fluxes, from which one can further study the qualitative properties of ionic flows and extract concrete information directly related to biological measurements. Of particular interest is the competition between two cations due to the nonlinear interplay between finite ion sizes, diffusion coefficients and boundary conditions, which is closely related to selectivity phenomena of open ion channels with given protein structures. Furthermore, we are able to characterize the distinct effects of the nonlinear interplays between these physical parameters. Numerical simulations are performed to identify some critical potentials which play critical roles in examining properties of ionic flows in our analysis.

Modeling single-molecule stochastic transport for DNA exo-sequencing in nanopore sensors

We present a simulation framework for computing the probability that a single molecule reaches the recognition element in a nanopore sensor. The model consists of the Langevin equation for the diffusive motion of small particles driven by external forces and the Poisson-Nernst-Planck-Stokes equations to compute these forces. The model is applied to examine DNA exo-sequencing in α-hemolysin, whose practicability depends on whether isolated DNA monomers reliably migrate into the channel in their correct order. We find that, at moderate voltage, migration fails in the majority of trials if the exonuclease which releases monomers is located farther than 1 nm above the pore entry. However, by tuning the pore to have a higher surface charge, applying a high voltage of 1 V and ensuring the exonuclease stays close to the channel, success rates of over 95% can be achieved.

Polynucleotide differentiation using hybrid solid-state nanopore functionalizing with α-hemolysin

We report results from full atomistic molecular dynamics simulations on the properties of biomimetic nanopores. This latter result was obtained through the direct insertion of an α-hemolysin protein inside a hydrophobic solid-state nanopore. Upon translocation of different DNA strands, we demonstrate here that the theoretical system presents the same discrimination properties as the experimental one obtained previously. This opens an interesting way to promote the stability of a specific protein inside a solid nanopore to develop further biomimetic applications for DNA or protein sequencing.

Nanopore Fabrication and Application as Biosensors in Neurodegenerative Diseases

Since its conception as an applied biomedical technology nearly 30 years ago, nanopore is emerging as a promising, high-throughput, biomarker-targeted diagnostic tool for clinicians. The attraction of a nanopore-based detection system is its simple, inexpensive, robust, user-friendly, high-throughput blueprint with minimal sample preparation needed prior to analysis. The goal of clinical-based nanopore biosensing is to go from sample acquisition to a meaningful readout quickly. The most extensive work in nanopore applications has been targeted at DNA, RNA, and peptide identification. Although, biosensing of pathological biomarkers, which is covered in this review, is on the rise. This review is broken into two major sections: (i) the current state of existing biological, solid state, and hybrid nanopore systems and (ii) the applications of nanopore biosensors toward detecting neurodegenerative biomarkers.

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.

Selectivity of the KcsA potassium channel: Analysis and computation

Ion channels regulate the flux of ions through cell membranes and play significant roles in many physiological functions. Most of the existing literature focuses on computational approaches based on molecular dynamics simulation or numerical solution of the modified Poisson-Nernst-Planck (PNP) system. In this paper, we present an analytical and computational study of a mathematical model of the KcsA potassium channel, including the effects of ion size (Bikerman model) and solvation energy (Born model). Under equilibrium conditions, we obtain an analytical solution of our modified PNP system, which is used to explain selectivity of KcsA of various ions (K^{+}, Na^{+}, Cl^{-}, Ca^{2+}, and Ba^{2+}) due to negative permanent charges inside the filter region and the effect of ion sizes. Our results show that K^{+} is always selected over Na^{+}, as smaller Na^{+} ions have larger solvation energy. As the amount of negative charges in the filter exceeds a critical value, divalent ions (Ca^{2+} and Ba^{2+}) can enter the filter region and block the KcsA channel. For the nonequilibrium cases, due to difficulties associated with a pure analytical or numerical approach, we use a hybrid analytical-numerical method to solve the modified PNP system. Our predictions of selectivity of KcsA channels and saturation phenomenon of the current-voltage (I-V) curve agree with experimental observations.

Rapid and Accurate Determination of Nanopore Ionic Current Using a Steric Exclusion Model

Nanopore sensing has emerged as a versatile approach to detection and identification of biomolecules. Presently, researchers rely on experience and intuition for choosing or modifying the nanopores to detect a target analyte. The field would greatly benefit from a computational method that could relate the atomic-scale geometry of the nanopores and analytes to the blockade nanopore currents they produce. Existing computational methods are either computationally too expensive to be used routinely in experimental laboratories or not sensitive enough to account for the atomic structure of the pore and the analytes. Here, we demonstrate a robust and inexpensive computational approach-the steric exclusion model (SEM) of nanopore conductance-that is orders of magnitude more efficient than all-atom MD and yet is sensitive enough to account for the atomic structure of the nanopore and the analyte. The method combines the computational efficiency of a finite element solver with the atomic precision of a nanopore conductance map to yield unprecedented speed and accuracy of ionic current prediction. We validate our SEM approach through comparison with the current blockades computed using the all-atom molecular dynamics method for a range of proteins confined to a solid-state nanopore, biological channels embedded in a lipid bilayer membranes, and blockade currents produced by DNA homopolymers in MspA. We illustrate potential applications of SEM by computing blockade currents produced by nucleosome proteins in a solid-state nanopore, individual amino acids in MspA, and by testing the effect of point mutations on amino acid distinguishability. We expect our SEM approach to become an integral part of future development of the nanopore sensing field.

Orientational and Folding Thermodynamics via Electric Dipole Moment Restraining

The projection of molecular processes onto a small set of relevant descriptors, the so-called reaction coordinates or collective variables (CVs), is a technique nowadays routinely employed by the biomolecular simulation community. In this work, we implemented two CVs to manipulate the orientation (i.e., angle) (μ⃗_a) and magnitude (|μ⃗|) of the electric dipole moment. In doing so, we studied the thermodynamics of water orientation under the application of external voltages and the folding of two polypeptides at zero-field conditions. The projection of the free-energy [potential of mean force (PMF)] along water orientation defined an upper limit of around 0.3 V for irrelevant thermodynamic effects. On the other hand, sufficiently strong μ⃗_a restraints applied on 12-alanine (Ala₁₂) triggered structural effects because of the alignment of local dipoles; for lower restraints, a full-body rotation is achieved. The manipulation of |μ⃗| produced strong perturbations on the secondary structure of Ala₁₂, promoting an enhanced sampling to its configurational space. Rigorous free-energy calculations in the form of 2-D PMFs for deca-alanine showed the utility of |μ⃗| as a reaction coordinate to study folding in small α helices. As a whole, we propose that the manipulation of both components of the dipole moment, μ⃗_a and |μ⃗|, provides thermodynamics insights into the structural conformation and stability of biomolecules. These new CVs are implemented in the Colvars module, available for NAMD and LAMMPS.

Alpha-hemolysin nanopore allows discrimination of the microcystins variants

Microcystins (MCs) are a class of cyclic heptapeptides with more than 100 variants produced by cyanobacteria present in surface waters. MCs are potent hepatotoxic agents responsible for fatal poisoning in animals and humans. Several techniques are employed in the detection of MCs, however, there is a shortage of methods capable of discriminating variants of MCs. In this work we demonstrate that the α-hemolysin (αHL) nanopore can detect and discriminate the variants (LR, YR and RR) of MCs in aqueous solution. The discrimination process is based on the analysis of the residence times of each variant of MCs within the unitary nanopore, as well as, on the amplitudes of the blockages in the ionic current flowing through it. Simulations of molecular dynamics and calculation of the electrostatic potential revealed that the variants of MCs present different charge distribution and correlated with the three patterns on the amplitudes of the blockages in the ionic current. Additionally, molecular docking analysis indicates different patterns of interaction of the variants of MCs with two specific regions of the nanopore. We conclude that αHL nanopore can discriminate variants of microcystins by a mechanism based mainly on electrostatic interaction. Finally, we propose the use of nanopore-based technology as a promising method for analyzing microcystins in aqueous solutions.

Microcystins (MCs) are a class of cyclic heptapeptides with more than 100 variants produced by cyanobacteria present in surface waters.

Comparative biosensing of glycosaminoglycan hyaluronic acid oligo- and polysaccharides using aerolysin and [Formula: see text]-hemolysin nanopores⋆

Seeking new tools for the analysis of glycosaminoglycans, we have compared the translocation of anionic oligosaccharides from hyaluronic acid using aerolysin and [Formula: see text]-hemolysin nanopores. We show that pores of similar channel length and diameter lead to distinct translocation behavior of the same macromolecules, due to different structural properties of the nanopores. When passing from the vestibule side of the nanopores, short hyaluronic acid oligosaccharides could be detected during their translocation across an aerolysin nanopore but not across an [Formula: see text]-hemolysin nanopore. We were however able to detect longer oligosaccharide fragments, resulting from the in situ enzymatic depolymerization of hyaluronic acid polysaccharides, with both nanopores, meaning that short oligosaccharides were crossing the [Formula: see text]-hemolysin nanopore with a speed too high to be detected. The translocation speed was an order of magnitude higher across [Formula: see text]-hemolysin compared to aerolysin. These results show that the choice of a nanopore to be used for resistive pulse sensing experiments should not rely only on the diameter of the channel but also on other parameters such as the charge repartition within the pore lumen.

The Role of Lipid Interactions in Simulations of the α-Hemolysin Ion-Channel-Forming Toxin

Molecular dynamics simulations were performed to describe the function of the ion-channel-forming toxin α-hemolysin (αHL) in lipid membranes that were composed of either 1,2-diphytanoyl-sn-glycero-3-phospho-choline or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline. The simulations highlight the importance of lipid type in maintaining αHL structure and function, enabling direct comparison to experiments for biosensing applications. We determined that although the two lipids studied are similar in structure, 1,2-diphytanoyl-sn-glycero-3-phospho-choline membranes better match the hydrophobic thickness of αHL compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline membranes. This hydrophobic match is essential to maintaining proper alignment of β-sheet loops at the trans entrance of αHL, which, when disrupted, creates an additional constriction to ion flow that decreases the channel current below experimental values and creates greater variability in channel conductance. Agreement with experiments was further improved with sufficient lipid membrane equilibration and allowed the discrimination of subtle αHL conduction states with lipid type. Finally, we explore the effects of truncating the extramembrane cap of αHL and its role in maintaining proper alignment of αHL in the membrane and channel conductance. Our results demonstrate the essential role of lipid type and lipid-protein interactions in simulations of αHL and will considerably improve the interpretation of experimental data.

Boundary Layer Effects on Ionic Flows Via Classical Poisson-Nernst-Planck Systems

A quasi-one-dimensional steady-state Poisson-Nernst-Planck model of two oppositely charged ion species through a membrane channel is analyzed. The model problem is treated as a boundary value problem of a singularly perturbed differential system. Our analysis is based on the geometric singular perturbation theory but, most importantly, on specific structures of this concrete model. The existence and (local ) uniqueness of solutions to the boundary value problem is established. In particular, an approximation of both the individual flux and the I-V (current-voltage) relation are derived explicitly from the zeroth order approximation (in ") solutions, from which the boundary layer effects on ionic flows are studied in great details.

Conic shapes have higher sensitivity than cylindrical ones in nanopore DNA sequencing

Nanopores have emerged as helpful research tools for single molecule detection. Through continuum modeling, we investigated the effects of membrane thickness, nanopore size, and pore shape on current signal characteristics of DNA. The simulation results showed that, when reducing the pore diameter, the amplitudes of current signals of DNA increase. Moreover, we found that, compared to cylindrically shaped nanopores, conical-shaped nanopores produce greater signal amplitudes from biomolecules translocation. Finally, we demonstrated that continuum model simulations for the discrimination of DNA and RNA yield current characteristics approximately consistent with experimental measurements and that A-T and G-C base pairs can be distinguished using thin conical solid-state nanopores. Our study not only suggests that computational approaches in this work can be used to guide the designs of nanopore for single molecule detection, but it also provides several possible ways to improve the current amplitudes of nanopores for better resolution.

Single-Molecule Titration in a Protein Nanoreactor Reveals the Protonation/Deprotonation Mechanism of a C:C Mismatch in DNA

Measurement of single-molecule reactions can elucidate microscopic mechanisms that are often hidden from ensemble analysis. Herein, we report the acid-base titration of a single DNA duplex confined within the wild-type α-hemolysin (α-HL) nanopore for up to 3 h, while monitoring the ionic current through the nanopore. Modulation between two states in the current-time trace for duplexes containing the C:C mismatch in proximity to the latch constriction of α-HL is attributed to the base flipping of the C:C mismatch. As the pH is lowered, the rate for the C:C mismatch to flip from the intra-helical state to the extra-helical state ( k_intra-extra) decreases, while the rate for base flipping from the extra-helical state to the intra-helical state ( k_extra-intra) remains unchanged. Both k_intra-extra and k_extra-intra are on the order of 1 × 10^-2 s^-1 to 1 × 10^-1 s^-1 and remain stable over the time scale of the measurement (several hours). Analysis of the pH-dependent kinetics of base flipping using a hidden Markov kinetic model demonstrates that protonation/deprotonation occurs while the base pair is in the intra-helical state. We also demonstrate that the rate of protonation is limited by transport of H⁺ into the α-HL nanopore. Single-molecule kinetic isotope experiments exhibit a large kinetic isotope effect (KIE) for k_intra-extra ( k_H/ k_D ≈ 5) but a limited KIE for k_extra-intra ( k_H/ k_D ≈ 1.3), supporting our model. Our experiments correspond to the longest single-molecule measurements performed using a nanopore, and demonstrate its application in interrogating mechanisms of single-molecule reactions in confined geometries.

Membrane Position Dependency of the pKa and Conductivity of the Protein Ion Channel

The dependency of current-voltage characteristics of the α-hemolysin channel on the channel position within the membrane was studied using Poisson-Nernst-Planck theory of ion conductivity with soft repulsion between mobile ions and protein atoms (SP-PNP). The presence of the membrane environment also influences the protonation state of the residues at the boundary of the water-lipid interface. In this work, we predict that Asp and Lys residues at the protein rim change their protonation state upon penetration to the lipid environment. Free energies of protein insertion in the membrane for different penetration depths were estimated using the Poisson-Boltzmann/solvent-accessible surface area (PB/SASA) model. The results show that rectification and reversal potentials are very sensitive to the relative position of channel in the membrane, which in turn contributes to alternative protonation states of lipid-penetrating ionizable groups. The prediction of channel position based on the matching of calculated rectification with experimentally determined rectification is in good agreement with recent neutron reflection experiments. Based on the results, we conclude that α-hemolysin membrane position is determined by a combination of factors and not only by the pattern of the surface hydrophobicity as is typically assumed.

Multi-Dielectric Brownian Dynamics and Design-Space-Exploration Studies of Permeation in Ion Channels

This paper proposes a multi-dielectric Brownian dynamics simulation framework for design-space-exploration (DSE) studies of ion-channel permeation. The goal of such DSE studies is to estimate the channel modeling-parameters that minimize the mean-squared error between the simulated and expected "permeation characteristics." To address this computational challenge, we use a methodology based on statistical inference that utilizes the knowledge of channel structure to prune the design space. We demonstrate the proposed framework and DSE methodology using a case study based on the KcsA ion channel, in which the design space is successfully reduced from a 6-D space to a 2-D space. Our results show that the channel dielectric map computed using the framework matches with that computed directly using molecular dynamics with an error of 7%. Finally, the scalability and resolution of the model used are explored, and it is shown that the memory requirements needed for DSE remain constant as the number of parameters (degree of heterogeneity) increases.

Gating of Connexin Channels by transjunctional-voltage: Conformations and models of open and closed states

Voltage is an important physiologic regulator of channels formed by the connexin gene family. Connexins are unique among ion channels in that both plasma membrane inserted hemichannels (undocked hemichannels) and intercellular channels (aggregates of which form gap junctions) have important physiological roles. The hemichannel is the fundamental unit of gap junction voltage-gating. Each hemichannel displays two distinct voltage-gating mechanisms that are primarily sensitive to a voltage gradient formed along the length of the channel pore (the transjunctional voltage) rather than sensitivity to the absolute membrane potential (V_m or V_i-o). These transjunctional voltage dependent processes have been termed V_j- or fast-gating and loop- or slow-gating. Understanding the mechanism of voltage-gating, defined as the sequence of voltage-driven transitions that connect open and closed states, first and foremost requires atomic resolution models of the end states. Although ion channels formed by connexins were among the first to be characterized structurally by electron microscopy and x-ray diffraction in the early 1980's, subsequent progress has been slow. Much of the current understanding of the structure-function relations of connexin channels is based on two crystal structures of Cx26 gap junction channels. Refinement of crystal structure by all-atom molecular dynamics and incorporation of charge changing protein modifications has resulted in an atomic model of the open state that arguably corresponds to the physiologic open state. Obtaining validated atomic models of voltage-dependent closed states is more challenging, as there are currently no methods to solve protein structure while a stable voltage gradient is applied across the length of an oriented channel. It is widely believed that the best approach to solve the atomic structure of a voltage-gated closed ion channel is to apply different but complementary experimental and computational methods and to use the resulting information to derive a consensus atomic structure that is then subjected to rigorous validation. In this paper, we summarize our efforts to obtain and validate atomic models of the open and voltage-driven closed states of undocked connexin hemichannels. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

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.

Current state of theoretical and experimental studies of the voltage-dependent anion channel (VDAC)

Voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane provides a controlled pathway for respiratory metabolites in and out of the mitochondria. In spite of the wealth of experimental data from structural, biochemical, and biophysical investigations, the exact mechanisms governing selective ion and metabolite transport, especially the role of titratable charged residues and interactions with soluble cytosolic proteins, remain hotly debated in the field. The computational advances hold a promise to provide a much sought-after solution to many of the scientific disputes around solute and ion transport through VDAC and hence, across the mitochondrial outer membrane. In this review, we examine how Molecular Dynamics, Free Energy, and Brownian Dynamics simulations of the large β-barrel channel, VDAC, advanced our understanding. We will provide a short overview of non-conventional techniques and also discuss examples of how the modeling excursions into VDAC biophysics prospectively aid experimental efforts. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Molecular dynamics simulation of water permeation through the alpha-hemolysin channel

The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (kBT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.

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.

Diffusive Models of Membrane Permeation with Explicit Orientational Freedom

Accurate calculation of permeabilities from first-principles has been a long-standing challenge for computer simulations, notably in the context of drug discovery, as a route to predict the propensity of small, organic molecules to spontaneously translocate biological membranes. Of equal importance is the understanding of the permeation process and the pathway followed by the permeant from the aqueous medium to the interior of the lipid bilayer, and back out again. A convenient framework for the computation of permeabilities is provided by the solubility-diffusion model, which requires knowledge of the underlying free-energy and diffusivity landscapes. Here, we develop a formalism that includes an explicit description of the orientational motion of the solute as it diffuses across the membrane. Toward this end, we have generalized a recently proposed method that reconciles thermodynamics and kinetics in importance-sampling simulations by means of a Bayesian-inference scheme to reverse-solve the underlying Smoluchowski equation. Performance of the proposed formalism is examined in the model cases of a water and an ethanol molecule crossing a fully hydrated lipid bilayer. Our analysis reveals a conspicuous dependence of the free-energy and rotational diffusivity on the orientation of ethanol when it lies within the headgroup region of the bilayer. Specifically, orientations for which the hydroxyl group lies among the polar lipid head groups, while the ethyl group recedes toward the hydrophobic interior are associated with free-energy minima ∼2kBT deep, as well as significantly slower orientational kinetics compared to the bulk solution or the core of the bilayer. The conspicuous orientational anisotropy of ethanol at the aqueous interface is suggestive of a complete rotation of the permeant as it crosses the hydrophobic interior of the membrane.

Fast assessment of structural models of ion channels based on their predicted current-voltage characteristics

Computational prediction of protein structures is a difficult task, which involves fast and accurate evaluation of candidate model structures. We propose to enhance single-model quality assessment with a functionality evaluation phase for proteins whose quantitative functional characteristics are known. In particular, this idea can be applied to evaluation of structural models of ion channels, whose main function - conducting ions - can be quantitatively measured with the patch-clamp technique providing the current-voltage characteristics. The study was performed on a set of KcsA channel models obtained from complete and incomplete contact maps. A fast continuous electrodiffusion model was used for calculating the current-voltage characteristics of structural models. We found that the computed charge selectivity and total current were sensitive to structural and electrostatic quality of models. In practical terms, we show that evaluating predicted conductance values is an appropriate method to eliminate models with an occluded pore or with multiple erroneously created pores. Moreover, filtering models on the basis of their predicted charge selectivity results in a substantial enrichment of the candidate set in highly accurate models. Tests on three other ion channels indicate that, in addition to being a proof of the concept, our function-oriented single-model quality assessment method can be directly applied to evaluation of structural models of some classes of protein channels. Finally, our work raises an important question whether a computational validation of functionality should be included in the evaluation process of structural models, whenever possible.