Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map

Structure and response to flow of the glycocalyx layer

The glycocalyx is a sugar-rich layer located at the luminal part of the endothelial cells. It is involved in key metabolic processes and its malfunction is related to several diseases. To understand the function of the glycocalyx, a molecular level characterization is necessary. In this article, we present large-scale molecular-dynamics simulations that provide a comprehensive description of the structure and dynamics of the glycocalyx. We introduce the most detailed, to-date, all-atom glycocalyx model, composed of lipid bilayer, proteoglycan dimers, and heparan sulfate chains with realistic sequences. Our results reveal the folding of proteoglycan ectodomain and the extended conformation of heparan sulfate chains. Furthermore, we study the glycocalyx response under shear flow and its role as a flypaper for binding fibroblast growth factors (FGFs), which are involved in diverse functions related to cellular differentiation, including angiogenesis, morphogenesis, and wound healing. The simulations show that the glycocalyx increases the effective concentration of FGFs, leading to FGF oligomerization, and acts as a lever to transfer mechanical stimulus into the cytoplasmic side of endothelial cells.

Geometrical membrane curvature as an allosteric regulator of membrane protein structure and function

Transmembrane proteins are embedded in cellular membranes of varied lipid composition and geometrical curvature. Here, we studied for the first time the allosteric effect of geometrical membrane curvature on transmembrane protein structure and function. We used single-channel optical analysis of the prototypic transmembrane β-barrel α-hemolysin (α-HL) reconstituted on immobilized single small unilamellar liposomes of different diameter and therefore curvature. Our data demonstrate that physiologically abundant geometrical membrane curvatures can enforce a dramatic allosteric regulation (1000-fold inhibition) of α-HL permeability. High membrane curvatures (1/diameter ~1/40 nm(-1)) compressed the effective pore diameter of α-HL from 14.2 ± 0.8 Å to 11.4 ± 0.6 Å. This reduction in effective pore area (~40%) when combined with the area compressibility of α-HL revealed an effective membrane tension of ~50 mN/m and a curvature-imposed protein deformation energy of ~7 kBT. Such substantial energies have been shown to conformationally activate, or unfold, β-barrel and α-helical transmembrane proteins, suggesting that membrane curvature could likely regulate allosterically the structure and function of transmembrane proteins in general.

Channel-forming bacterial toxins in biosensing and macromolecule delivery

To intoxicate cells, pore-forming bacterial toxins are evolved to allow for the transmembrane traffic of different substrates, ranging from small inorganic ions to cell-specific polypeptides. Recent developments in single-channel electrical recordings, X-ray crystallography, protein engineering, and computational methods have generated a large body of knowledge about the basic principles of channel-mediated molecular transport. These discoveries provide a robust framework for expansion of the described principles and methods toward use of biological nanopores in the growing field of nanobiotechnology. This article, written for a special volume on "Intracellular Traffic and Transport of Bacterial Protein Toxins", reviews the current state of applications of pore-forming bacterial toxins in small- and macromolecule-sensing, targeted cancer therapy, and drug delivery. We discuss the electrophysiological studies that explore molecular details of channel-facilitated protein and polymer transport across cellular membranes using both natural and foreign substrates. The review focuses on the structurally and functionally different bacterial toxins: gramicidin A of Bacillus brevis, α-hemolysin of Staphylococcus aureus, and binary toxin of Bacillus anthracis, which have found their "second life" in a variety of developing medical and technological applications.

Computational investigation of cholesterol binding sites on mitochondrial VDAC

The mitochondrial voltage-dependent anion channel (VDAC) allows passage of ions and metabolites across the mitochondrial outer membrane. Cholesterol binds mammalian VDAC, and we investigated the effects of binding to human VDAC1 with atomistic molecular dynamics simulations that totaled 1.4 μs. We docked cholesterol to specific sites on VDAC that were previously identified with NMR, and we tested the reliability of multiple docking results in each site with simulations. The most favorable binding modes were used to build a VDAC model with cholesterol occupying five unique sites, and during multiple 100 ns simulations, cholesterol stably and reproducibly remained bound to the protein. For comparison, VDAC was simulated in systems with identical components but with cholesterol initially unbound. The dynamics of loops that connect adjacent β-strands were most affected by bound cholesterol, with the averaged root-mean-square fluctuation (RMSF) of multiple residues altered by 20–30%. Cholesterol binding also stabilized charged residues inside the channel and localized the surrounding electrostatic potentials. Despite this, ion diffusion through the channel was not significantly affected by bound cholesterol, as evidenced by multi-ion potential of mean force measurements. Although we observed modest effects of cholesterol on the open channel, our model will be particularly useful in experiments that investigate how cholesterol affects VDAC function under applied electrochemical forces and also how other ligands and proteins interact with the channel.

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.

Adenine nucleotide translocase is acetylated in vivo in human muscle: Modeling predicts a decreased ADP affinity and altered control of oxidative phosphorylation

Proteomics techniques have revealed that lysine acetylation is abundant in mitochondrial proteins. This study was undertaken (1) to determine the relationship between mitochondrial protein acetylation and insulin sensitivity in human skeletal muscle, identifying key acetylated proteins, and (2) to use molecular modeling techniques to understand the functional consequences of acetylation of adenine nucleotide translocase 1 (ANT1), which we found to be abundantly acetylated. Eight lean and eight obese nondiabetic subjects had euglycemic clamps and muscle biopsies for isolation of mitochondrial proteins and proteomics analysis. A number of acetylated mitochondrial proteins were identified in muscle biopsies. Overall, acetylation of mitochondrial proteins was correlated with insulin action (r = 0.60; P < 0.05). Of the acetylated proteins, ANT1, which catalyzes ADP-ATP exchange across the inner mitochondrial membrane, was acetylated at lysines 10, 23, and 92. The extent of acetylation of lysine 23 decreased following exercise, depending on insulin sensitivity. Molecular dynamics modeling and ensemble docking simulations predicted the ADP binding site of ANT1 to be a pocket of positively charged residues, including lysine 23. Calculated ADP-ANT1 binding affinities were physiologically relevant and predicted substantial reductions in affinity upon acetylation of lysine 23. Insertion of these derived binding affinities as parameters into a complete mathematical description of ANT1 kinetics predicted marked reductions in adenine nucleotide flux resulting from acetylation of lysine 23. Therefore, acetylation of ANT1 could have dramatic physiological effects on ADP-ATP exchange. Dysregulation of acetylation of mitochondrial proteins such as ANT1 therefore could be related to changes in mitochondrial function that are associated with insulin resistance.

An Overview of Electrostatic Free Energy Computations for Solutions and Proteins

Free energy simulations for electrostatic and charging processes in complex molecular systems encounter specific difficulties owing to the long-range, 1/r Coulomb interaction. To calculate the solvation free energy of a simple ion, it is essential to take into account the polarization of nearby solvent but also the electrostatic potential drop across the liquid-gas boundary, however distant. The latter does not exist in a simulation model based on periodic boundary conditions because there is no physical boundary to the system. An important consequence is that the reference value of the electrostatic potential is not an ion in a vacuum. Also, in an infinite system, the electrostatic potential felt by a perturbing charge is conditionally convergent and dependent on the choice of computational conventions. Furthermore, with Ewald lattice summation and tinfoil conducting boundary conditions, the charges experience a spurious shift in the potential that depends on the details of the simulation system such as the volume fraction occupied by the solvent. All these issues can be handled with established computational protocols, as reviewed here and illustrated for several small ions and three solvated proteins.

Rectification of Ion Current in Nanopores Depends on the Type of Monovalent Cations: Experiments and Modeling

Rectifying nanopores feature ion currents that are higher for voltages of one polarity compared to the currents recorded for corresponding voltages of the opposite polarity. Rectification of nanopores has been found to depend on the pore opening diameter and distribution of surface charges on the pore walls as well as pore geometry. Very little is known, however, on the dependence of ionic rectification on the type of transported ions of the same charge. We performed experiments with single conically shaped nanopores in a polymer film and recorded current-voltage curves in three electrolytes: LiCl, NaCl, and KCl. Rectification degrees of the pores, quantified as the ratio of currents recorded for voltages of opposite polarities, were the highest for KCl and the lowest for LiCl. The experimental observations could not be explained by a continuum modeling based on the Poisson-Nernst-Planck equations. All-atom molecular dynamics simulations revealed differential binding between Li⁺, Na⁺, and K⁺ ions and carboxyl groups on the pore walls, resulting in changes to both the effective surface charge of the nanopore and cation mobility within the pore.

Comparative analysis of nucleotide translocation through protein nanopores using steered molecular dynamics and an adaptive biasing force

The translocation of nucleotide molecules across biological and synthetic nanopores has attracted attention as a next generation technique for sequencing DNA. Computer simulations have the ability to provide atomistic-level insight into important states and processes, delivering a means to develop a fundamental understanding of the translocation event, for example, by extracting the free energy of the process. Even with current supercomputing facilities, the simulation of many-atom systems in fine detail is limited to shorter timescales than the real events they attempt to recreate. This imposes the need for enhanced simulation techniques that expand the scope of investigation in a given timeframe. There are numerous free energy calculation and translocation methodologies available, and it is by no means clear which method is best applied to a particular problem. This article explores the use of two popular free energy calculation methodologies in a nucleotide-nanopore translocation system, using the α-hemolysin nanopore. The first uses constant velocity-steered molecular dynamics (cv-SMD) in conjunction with Jarzynski's equality. The second applies an adaptive biasing force (ABF), which has not previously been applied to the nucleotide-nanpore system. The purpose of this study is to provide a comprehensive comparison of these methodologies, allowing for a detailed comparative assessment of the scientific merits, the computational cost, and the statistical quality of the data obtained from each technique. We find that the ABF method produces results that are closer to experimental measurements than those from cv-SMD, whereas the net errors are smaller for the same computational cost.

Structural determinants of water permeation through the sodium-galactose transporter vSGLT

Sodium-glucose transporters (SGLTs) facilitate the movement of water across the cell membrane, playing a central role in cellular homeostasis. Here, we present a detailed analysis of the mechanism of water permeation through the inward-facing state of vSGLT based on nearly 10 μs of molecular dynamics simulations. These simulations reveal the transient formation of a continuous water channel through the transporter that permits water to permeate the protein. Trajectories in which spontaneous release of galactose is observed, as well as those in which galactose remains in the binding site, show that the permeation rate, although modulated by substrate occupancy, is not tightly coupled to substrate release. Using a, to our knowledge, novel channel-detection algorithm, we identify the key residues that control water flow through the transporter and show that solvent gating is regulated by side-chain motions in a small number of residues on the extracellular face. A sequence alignment reveals the presence of two insertion sites in mammalian SGLTs that flank these outer-gate residues. We hypothesize that the absence of these sites in vSGLT may account for the high water permeability values for vSGLT determined via simulation compared to the lower experimental estimates for mammalian SGLT1.

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.

ATP transport through VDAC and the VDAC-tubulin complex probed by equilibrium and nonequilibrium MD simulations

Voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane, serves as a principal pathway for ATP, ADP, and other respiratory substrates across this membrane. Using umbrella-sampling simulations, we established the thermodynamic and kinetic components governing ATP transport across the VDAC1 channel. We found that there are several low-affinity binding sites for ATP along the translocation pathway and that the main barrier for ATP transport is located around the center of the channel and is formed predominantly by residues in the N-terminus. The binding affinity of ATP to an open channel was found to be in the millimolar to micromolar range. However, we show that this weak binding increases the ATP translocation probability by about 10-fold compared with the VDAC pore in which attractive interactions were artificially removed. Recently, it was found that free dimeric tubulin induces a highly efficient, reversible blockage of VDAC reconstituted into planar lipid membranes. It was proposed that by blocking VDAC permeability for ATP/ADP and other mitochondrial respiratory substrates tubulin controls mitochondrial respiration. Using the Rosetta protein-protein docking algorithm, we established a tentative structure of the VDAC-tubulin complex. An extensive set of equilibrium and nonequilibrium (under applied electric field) molecular dynamics (MD) simulations was used to establish the conductance of the open and blocked channel. It was found that the presence of the unstructured C-terminal tail of tubulin in the VDAC pore decreases its conductance by more than 40% and switches its selectivity from anionic to cationic. The subsequent 1D potential of mean force (PMF) computations for the VDAC-tubulin complex show that the state renders ATP transport virtually impossible. A number of residues pivotal for tubulin binding to the channel were identified that help to clarify the molecular details of VDAC-tubulin interaction and to provide new insight into the mechanism of the control of mitochondria respiration by VDAC.

Slow DNA transport through nanopores in hafnium oxide membranes

We present a study of double- and single-stranded DNA transport through nanopores fabricated in ultrathin (2-7 nm thick) freestanding hafnium oxide (HfO2) membranes. The high chemical stability of ultrathin HfO2 enables long-lived experiments with <2 nm diameter pores that last several hours, in which we observe >50 000 DNA translocations with no detectable pore expansion. Mean DNA velocities are slower than velocities through comparable silicon nitride pores, providing evidence that HfO2 nanopores have favorable physicochemical interactions with nucleic acids that can be leveraged to slow down DNA in a nanopore.

Bioinspired graphene nanopores with voltage-tunable ion selectivity for Na(+) and K(+)

Biological protein channels have many remarkable properties such as gating, high permeability, and selectivity, which have motivated researchers to mimic their functions for practical applications. Herein, using molecular dynamics simulations, we design bioinspired nanopores in graphene sheets that can discriminate between Na(+) and K(+), two ions with very similar properties. The simulation results show that, under transmembrane voltage bias, a nanopore containing four carbonyl groups to mimic the selectivity filter of the KcsA K(+) channel preferentially conducts K(+) over Na(+). A nanopore functionalized by four negatively charged carboxylate groups to mimic the selectivity filter of the NavAb Na(+) channel selectively binds Na(+) but transports K(+) over Na(+). Surprisingly, the ion selectivity of the smaller diameter pore containing three carboxylate groups can be tuned by changing the magnitude of the applied voltage bias. Under lower voltage bias, it transports ions in a single-file manner and exhibits Na(+) selectivity, dictated by the knock-on ion conduction and selective blockage by Na(+). Under higher voltage bias, the nanopore is K(+)-selective, as the blockage by Na(+) is destabilized and the stronger affinity for carboxylate groups slows the passage of Na(+) compared with K(+). The computational design of biomimetic ion-selective nanopores helps to understand the mechanisms of selectivity in biological ion channels and may also lead to a wide range of potential applications such as sensitive ion sensors, nanofiltration membranes for Na(+)/K(+) separation, and voltage-tunable nanofluidic devices.

Electric migration of α-hemolysin in supported n-bilayers: a model for transmembrane protein microelectrophoresis

Proteome analysis involves separating proteins as a preliminary step toward their characterization. This paper reports on the translational migration of a model transmembrane protein (α-hemolysin) in supported n-bilayers (n, the number of bilayers, varies from 1 to around 500 bilayers) when an electric field parallel to the membrane plane is applied. The migration changes in direction as the charge on the protein changes its sign. Its electrophoretic mobility is shown to depend on size and charge. The electrophoretic mobility varies as 1/R(2), with R the equivalent geometric radius of the embedded part of the protein. Measuring mobilities at differing pH in our system enables us to determine the pI and the charge of the protein. Establishing all these variations points to the feasibility of electrophoretic transport of a charged object in this medium and is a first step toward electrophoretic separation of membrane proteins in n-bilayer systems.

Dangerous liaisons between detergents and membrane proteins. The case of mitochondrial uncoupling protein 2

The extraction of membrane proteins from their native environment by detergents is central to their biophysical characterization. Recent studies have emphasized that detergents may perturb the structure locally and modify the dynamics of membrane proteins. However, it remains challenging to determine whether these perturbations are negligible or could be responsible for misfolded conformations, altering the protein's function. In this work, we propose an original strategy combining functional studies and molecular simulations to address the physiological relevance of membrane protein structures obtained in the presence of detergents. We apply our strategy to a structure of isoform 2 of an uncoupling protein (UCP2) binding an inhibitor recently obtained in dodecylphosphocholine detergent micelles. Although this structure shares common traits with the ADP/ATP carrier, a member of the same protein family, its functional and biological significance remains to be addressed. In the present investigation, we demonstrate how dodecylphosphocholine severely alters the structure as well as the function of UCPs. The proposed original strategy opens new vistas for probing the physiological relevance of three-dimensional structures of membrane proteins obtained in non-native environments.

Optimization of 3D Poisson-Nernst-Planck model for fast evaluation of diverse protein channels

We show the accuracy and applicability of our fast algorithmic implementation of a three-dimensional Poisson-Nernst-Planck (3D-PNP) flow model for characterizing different protein channels. Due to its high computational efficiency, our model can predict the full current-voltage characteristics of a channel within minutes, based on the experimental 3D structure of the channel or its computational model structure. Compared with other methods, such as Brownian dynamics, which currently needs a few weeks of the computational time, or even much more demanding molecular dynamics modeling, 3D-PNP is the only available method for a function-based evaluation of very numerous tentative structural channel models. Flow model tests of our algorithm and its optimal parametrization are provided for five native channels whose experimental structures are available in the protein data bank (PDB) in an open conductive state, and whose experimental current-voltage characteristics have been published. The channels represent very different geometric and structural properties, which makes it the widest test to date of the accuracy of 3D-PNP on real channels. We test whether the channel conductance, rectification, and charge selectivity obtained from the flow model, could be sufficiently sensitive to single-point mutations, related to unsignificant changes in the channel structure. Our results show that the classical 3D-PNP model, under proper parametrization, is able to achieve a qualitative agreement with experimental data for a majority of the tested characteristics and channels, including channels with narrow and irregular conductivity pores. We propose that although the standard PNP model cannot provide insight into complex physical phenomena due to its intrinsic limitations, its semiquantitative agreement is achievable for rectification and selectivity at a level sufficient for the bioinformatical purpose of selecting the best structural models with a great advantage of a very short computational time.

Disentangling steric and electrostatic factors in nanoscale transport through confined space

The voltage-driven passage of biological polymers through nanoscale pores is an analytically, technologically, and biologically relevant process. Despite various studies on homopolymer translocation there are still several open questions on the fundamental aspects of pore transport. One of the most important unresolved issues revolves around the passage of biopolymers which vary in charge and volume along their sequence. Here we exploit an experimentally tunable system to disentangle and quantify electrostatic and steric factors. This new, fundamental framework facilitates the understanding of how complex biopolymers are transported through confined space and indicates how their translocation can be slowed down to enable future sensing methods.

Role of the central arginine R133 toward the ion selectivity of the phosphate specific channel OprP: effects of charge and solvation

The outer membrane porin OprP of Pseudomonas aeruginosa forms a highly specific phosphate selective channel. This channel is responsible for the high-affinity uptake of phosphate ions into the periplasmic space of the bacteria. A detailed investigation of the structure-function relationship of OprP is inevitable to decipher the anion and phosphate selectivity of this porin in particular and to broaden the present understanding of the ion selectivity of different channels. To this end we investigated the role of the central arginine of OprP, R133, in terms of its effects in selectivity and ion transport properties of the pore. Electrophysiological bilayer measurements and free-energy molecular dynamics simulations were carried out to probe the transport of different ions through various R133 mutants. For these mutants, the change in phosphate binding specificity, ion conduction, and anion selectivity was determined and compared to previous molecular dynamic calculations and electrophysiological measurements with wild-type OprP. Molecular analysis revealed a rather particular role of arginine 133 and its charge, while at the same time this residue together with the network of other residues, namely, D94 and Y114, has the ability to dehydrate the permeating ion. These very specific features govern the ion selectivity of OprP.

Collective diffusion model for ion conduction through microscopic channels

Ion conduction through microscopic channels is of central importance in both biology and nanotechnology. To better understand the current-voltage (I-V) dependence of ion channels, here we describe and prove a collective diffusion model that quantitatively relates the spontaneous ion permeation at equilibrium to the stationary ionic fluxes driven by small voltages. The model makes it possible to determine the channel conductance in the linear I-V range from equilibrium simulations without the application of a voltage. To validate the theory, we perform molecular-dynamics simulations on two channels-a conical-shaped nanopore and the transmembrane pore of an α-hemolysin-under both equilibrium and nonequilibrium conditions. The simulations reveal substantial couplings between the motions of cations and anions, which are effectively captured by the collective coordinate in the model. Although the two channels exhibit very different linear ranges in the I-V curves, in both cases the channel conductance at small voltages is in reasonable agreement with the prediction from the equilibrium simulation. The simulations also suggest that channel charges, rather than geometric asymmetry, play a more prominent role in current rectification.

The effects of diffusion on an exonuclease/nanopore-based DNA sequencing engine

Over 15 years ago, the ability to electrically detect and characterize individual polynucleotides as they are driven through a single protein ion channel was suggested as a potential method for rapidly sequencing DNA, base-by-base, in a ticker tape-like fashion. More recently, a variation of this method was proposed in which a nanopore would instead detect single nucleotides cleaved sequentially by an exonuclease enzyme in close proximity to one pore entrance. We analyze the exonuclease/nanopore-based DNA sequencing engine using analytical theory and computer simulations that describe nucleotide transport. The available data and analytical results suggest that the proposed method will be limited to reading <80 bases, imposed, in part, by the short lifetime each nucleotide spends in the vicinity of the detection element within the pore and the ability to accurately discriminate between the four mononucleotides.

Designing a polycationic probe for simultaneous enrichment and detection of microRNAs in a nanopore

The nanopore sensor can detect cancer-derived nucleic acid biomarkers such as microRNAs (miRNAs), providing a noninvasive tool potentially useful in medical diagnostics. However, the nanopore-based detection of these biomarkers remains confounded by the presence of numerous other nucleic acid species found in biofluid extracts. Their nonspecific interactions with the nanopore inevitably contaminate the target signals, reducing the detection accuracy. Here we report a novel method that utilizes a polycationic peptide-PNA probe as the carrier for selective miRNA detection in the nucleic acid mixture. The cationic probe hybridized with microRNA forms a dipole complex, which can be captured by the pore using a voltage polarity that is opposite the polarity used to capture negatively charged nucleic acids. As a result, nontarget species are driven away from the pore opening, and the target miRNA can be detected accurately without interference. In addition, we demonstrate that the PNA probe enables accurate discrimination of miRNAs with single-nucleotide difference. This highly sensitive and selective nanodielectrophoresis approach can be applied to the detection of clinically relevant nucleic acid fragments in complex samples.

Theory of polymer-nanopore interactions refined using molecular dynamics simulations

Molecular dynamics simulations were used to refine a theoretical model that describes the interaction of single polyethylene glycol (PEG) molecules with α-hemolysin (αHL) nanopores. The simulations support the underlying assumptions of the model, that PEG decreases the pore conductance by binding cations (which reduces the number of mobile ions in the pore) and by volume exclusion, and provide bounds for fits to new experimental data. Estimation of cation binding indicates that four monomers coordinate a single K(+) in a crown-ether-like structure, with, on average, 1.5 cations bound to a PEG 29-mer at a bulk electrolyte concentration of 4 M KCl. Additionally, PEG is more cylindrical and has a larger cross-section area in the pore than in solution, although its volume is similar. Two key experimental quantities of PEG are described by the model: the ratio of single channel current in the presence of PEG to that in the polymer's absence (blockade depth) and the mean residence time of PEG in the pore. The refined theoretical model is simultaneously fit to the experimentally determined current blockade depth and the mean residence times for PEGs with 15 to 45 monomers, at applied transmembrane potentials of -40 to -80 mV and for three electrolyte concentrations. The model estimates the free energy of the PEG-cation complexes to be -5.3 kBT. Finally the entropic penalty of confining PEG to the pore is found to be inversely proportional to the electrolyte concentration.

Oroxylin A inhibits hemolysis via hindering the self-assembly of α-hemolysin heptameric transmembrane pore

Alpha-hemolysin (α-HL) is a self-assembling, channel-forming toxin produced by most Staphylococcus aureus strains as a 33.2-kDa soluble monomer. Upon binding to a susceptible cell membrane, the monomer self-assembles to form a 232.4-kDa heptamer that ultimately causes host cell lysis and death. Consequently, α-HL plays a significant role in the pathogenesis of S. aureus infections, such as pneumonia, mastitis, keratitis and arthritis. In this paper, experimental studies show that oroxylin A (ORO), a natural compound without anti-S. aureus activity, can inhibit the hemolytic activity of α-HL. Molecular dynamics simulations, free energy calculations, and mutagenesis assays were performed to understand the formation of the α-HL-ORO complex. This combined approach revealed that the catalytic mechanism of inhibition involves the direct binding of ORO to α-HL, which blocks the conformational transition of the critical “Loop” region of the α-HL protein thereby inhibiting its hemolytic activity. This mechanism was confirmed by experimental data obtained from a deoxycholate-induced oligomerization assay. It was also found that, in a co-culture system with S. aureus and human alveolar epithelial (A549) cells, ORO could protect against α-HL-mediated injury. These findings indicate that ORO hinders the lytic activity of α-HL through a novel mechanism, which should facilitate the design of new and more effective antibacterial agents against S. aureus.

The mechanism controlling protein-ligand interactions is one of the most important processes in rational drug design. X-ray crystallography is a traditional tool used to investigate the interaction of ligands and proteins in a complex. However, protein crystallography is inefficient, and the development of crystal technology and research remains unequally distributed. Thus, it seems impractical to explore the structure of the α-hemolysin-ORO monomer complex by crystallography. Therefore, we used molecular dynamics simulations to investigate the receptor-ligand interaction in the α-HL-ORO monomer complex. In this study, we found that oroxylin A (ORO), a natural compound with little anti-S. aureus activity, can inhibit the hemolytic activity of α-HL at low concentrations. Through molecular docking and molecular dynamics simulations, we determined the potential binding mode of the protein-ligand interaction. The data revealed that ORO directly binds to α-HL, an interaction that blacks the conformational transition of the critical “Loop” region in α-HL and thus prevents the formation of the α-HL heptameric transmembrane pore, which ultimately inhibits the hemolytic activity of α-HL. This mechanism was confirmed by experimental data. Furthermore, we demonstrated that ORO could protect against α-HL-mediated injury in human alveolar epithelial (A549) cells.

Molecular insight into the inhibition mechanism of cyrtominetin to α-hemolysin by molecular dynamics simulation

The protein α-hemolysin (α-HL) is a self-assembling exotoxin that binds to the membrane of a susceptible host cell. In this paper, experimental studies show that cyrtominetin (CTM) can inhibit the hemolytic activity of α-HL. To understand how CTM can affect hemolytic activity, molecular dynamics simulations were carried out for α-HL-CTM complex and these results were compared with the crystal structure of monomeric α-HL. With this approach, the analysis revealed that the inhibition of CTM involves CTM directly binding to α-HL. Due to the binding of CTM, the conformation of the critical "Loop" region was restrained. This mechanism was confirmed by the experimental data. These findings indicate that CTM hinders the lysis activity of α-HL through a novel mechanism.

Lipid nanobilayers to host biological nanopores for DNA translocations

We characterize a recently introduced novel nanobilayer technique [Gornall, J. L., Mahendran, K. R., Pambos, O. J., Steinbock, L. J., Otto, O., Chimerel, C., Winterhalter, M., and Keyser, U. F. Simple reconstitution of protein pores in nano lipid bilayers. Nano Lett. 2011, 11 (8), 3334-3340] and its practical aspects for incorporating the biological nanopore α-hemolysin from Staphylococcus aureus and subsequent studies on the translocation of biomolecules under various conditions. This technique provides advantages over classical bilayer methods, especially the quick formation and extended stability of a bilayer. We have also developed a methodology to prepare a uniform quality of giant unilamellar vesicles (GUVs) in a reproducible way for producing nanobilayers. The process and the characteristics of the reconstitution of α-hemolysin in nanobilayers were examined by exploiting various important parameters, including pH, applied voltage, salt concentration, and number of nanopores. Protonation of α-hemolysin residues in the low pH region affects the translocation durations, which, in turn, changes the statistics of event types as a result of electrostatics and potentially the structural changes in DNA. When the pH and applied voltage were varied, it was possible to investigate and partly control the capture rates and type of translocation events through α-hemolysin nanopores. This study could be helpful to use the nanobilayer technique for further explorations, particularly owing to its advantages and technical ease compared to existing bilayer methods.

Pulling peptides across nanochannels: resolving peptide binding and translocation through the hetero-oligomeric channel from Nocardia farcinica

We investigated translocation of cationic peptides through nanochannels derived from the Gram-positive bacterium Nocardia farcinica at the single-molecule level. The two subunits NfpA and NfpB form a hetero-oligomeric cation selective channel. On the basis of amino acid comparison we performed homology modeling and obtained a channel structurally related to MspA of Mycobacterium smegmatis. The quantitative single-molecule measurements provide an insight into transport processes of solutes through nanochannels. High-resolution ion conductance measurements in the presence of peptides of different charge and length revealed the kinetics of peptide binding. The observed asymmetry in peptide binding kinetics indicated a unidirectional channel insertion in the lipid bilayer. In the case of cationic peptides, the external voltage acts as a driving force that promotes the interaction of the peptide with the channel surface. At low voltage, the peptide just binds to the channel, whereas at higher voltage, the force is strong enough to pull the peptide across the channel. This allows distinguishing quantitatively between peptide binding and translocation through the channel.

Protein nanopore-based, single-molecule exploration of copper binding to an antimicrobial-derived, histidine-containing chimera peptide

Metal ions binding exert a crucial influence upon the aggregation properties and stability of peptides, and the propensity of folding in various substates. Herein, we demonstrate the use of the α-HL protein as a powerful nanoscopic tool to probe Cu(2+)-triggered physicochemical changes of a 20 aminoacids long, antimicrobial-derived chimera peptide with a His residue as metal-binding site, and simultaneously dissect the kinetics of the free- and Cu(2+)-bound peptide interaction to the α-HL pore. Combining single-molecule electrophysiology on reconstituted lipid membranes and fluorescence spectroscopy, we show that the association rate constant between the α-HL pore and a Cu(2+)-free peptide is higher than that of a Cu(2+)-complexed peptide. We posit that mainly due to conformational changes induced by the bound Cu(2+) on the peptide, the resulting complex encounters a higher energy barrier toward its association with the protein pore, stemming most likely from an extra entropy cost needed to fit the Cu(2+)-complexed peptide within the α-HL lumen region. The lower dissociation rate constant of the Cu(2+)-complexed peptide from α-HL pore, as compared to that of Cu(2+)-free peptide, supports the existence of a deeper free energy well for the protein interaction with a Cu(2+)-complexed peptide, which may be indicative of specific Cu(2+)-mediated contributions to the binding of the Cu(2+)-complexed peptide within the pore lumen.

Modeling the Ion Selectivity of the Phosphate Specific Channel OprP

Ion selectivity of transport systems is an essential property of membranes from living organisms. These entities are used to regulate multifarious biological processes by virtue of selective participation of specific ions in transport processes. To understand this process, we studied the phosphate selectivity of the OprP porin from Pseudomonas aeruginosa using all-atom free-energy molecular dynamics simulations. These calculations were performed to define the energetics of phosphate, sulfate, chloride, and potassium ion transport through OprP. Atomic-level analysis revealed that the overall electrostatic environment of the channel was responsible for the anion selectivity of the channel, whereas the particular balance of interactions between the permeating ions and water as well as channel residues drove the selectivity between different anions. The selectivity of OprP is discussed in light of well-studied ion channels that are highly selective for potassium or chloride.

Toward detection of DNA-bound proteins using solid-state nanopores: insights from computer simulations

Through all-atom molecular dynamics simulations, we explore the use of nanopores in thin synthetic membranes for detection and identification of DNA binding proteins. Reproducing the setup of a typical experiment, we simulate electric field driven transport of DNA-bound proteins through nanopores smaller in diameter than the proteins. As model systems, we use restriction enzymes EcoRI and BamHI specifically and nonspecifically bound to a fragment of dsDNA, and streptavidin and NeutrAvidin proteins bound to dsDNA and ssDNA via a biotin linker. Our simulations elucidate the molecular mechanics of nanopore-induced rupture of a protein-DNA complex, the effective force applied to the DNA-protein bond by the electrophoretic force in a nanopore, and the role of DNA-surface interactions in the rupture process. We evaluate the ability of the nanopore ionic current and the local electrostatic potential measured by an embedded electrode to report capture of DNA, capture of a DNA-bound protein, and rupture of the DNA-protein bond. We find that changes in the strain on dsDNA can reveal the rupture of a protein-DNA complex by altering both the nanopore ionic current and the potential of the embedded electrode. Based on the results of our simulations, we suggest a new method for detection of DNA binding proteins that utilizes peeling of a nicked double strand under the electrophoretic force in a nanopore.

Synthetic lipid membrane channels formed by designed DNA nanostructures

We created nanometer-scale transmembrane channels in lipid bilayers by means of self-assembled DNA-based nanostructures. Scaffolded DNA origami was used to create a stem that penetrated and spanned a lipid membrane, as well as a barrel-shaped cap that adhered to the membrane, in part via 26 cholesterol moieties. In single-channel electrophysiological measurements, we found similarities to the response of natural ion channels, such as conductances on the order of 1 nanosiemens and channel gating. More pronounced gating was seen for mutations in which a single DNA strand of the stem protruded into the channel. Single-molecule translocation experiments show that the synthetic channels can be used to discriminate single DNA molecules.

Directly observing the motion of DNA molecules near solid-state nanopores

We investigate the diffusion and the drift motion of λ DNA molecules near solid-state nanopores prior to their translocation through the nanopores using fluorescence microscopy. The radial dependence of the electric field near a nanopore generated by an applied voltage in ionic solution can be estimated quantitatively in 3D by analyzing the motion of negatively charged DNA molecules. We find that the electric field is approximately spherically symmetric around the nanopore under the conditions investigated. In addition, DNA clogging at the nanopore was directly observed. Surprisingly, the probability of the clogging event increases with increasing external bias voltage. We also find that DNA molecules clogging the nanopore reduce the electric field amplitude at the nanopore membrane surface. To better understand these experimental results, analytical method with Ohm's law and computer simulation with Poisson and Nernst-Planck (PNP) equations are used to calculate the electric field near the nanopore. These results are of great interest in both experimental and theoretical considerations of the motion of DNA molecules near voltage-biased nanopores. These findings will also contribute to the development of solid-state nanopore-based DNA sensing devices.

Inspection of the engineered FhuA ΔC/Δ4L protein nanopore by polymer exclusion

Extensive engineering of protein nanopores for biotechnological applications using native scaffolds requires further inspection of their internal geometry and size. Recently, we redesigned ferric hydroxamate uptake component A (FhuA), a 22-β-stranded protein containing an N-terminal 160-residue cork domain (C). The cork domain and four large extracellular loops (4L) were deleted to obtain an unusually stiff engineered FhuA ΔC/Δ4L nanopore. We employed water-soluble poly(ethylene glycols) and dextran polymers to examine the interior of FhuA ΔC/Δ4L. When this nanopore was reconstituted into a synthetic planar lipid bilayer, addition of poly(ethylene glycols) produced modifications in the single-channel conductance, allowing for the evaluation of the nanopore diameter. Here, we report that FhuA ΔC/Δ4L features an approximate conical internal geometry with the cis entrance smaller than the trans entrance, in accord with the asymmetric nature of the crystal structure of the wild-type FhuA protein. Further experiments with impermeable dextran polymers indicated an average internal diameter of ~2.4 nm, a conclusion we arrived at based upon the polymer-induced alteration of the access resistance contribution to the nanopore's total resistance. Molecular insights inferred from this work represent a platform for future protein engineering of FhuA that will be employed for specific tasks in biotechnological applications.

The role of Lys147 in the interaction between MPSA-gold nanoparticles and the α-hemolysin nanopore

Single channel recordings were used to determine the effect of direct electrostatic interactions between sulfonate-coated gold nanoparticles and the constriction of the Staphylococcus aureus α-hemolysin protein channel on the ionic current amplitude. We provide evidence that Lys147 of α-hemolysin can interact with the sulfonate groups at the nanoparticle surface, and these interactions can reversibly block 100% of the residual ionic current. Lys147 is normally involved in a salt bridge with Glu111. The capture of a nanoparticle leads to a partial current block at neutral pH values, but protonation of Glu111 at pH 2.8 results in a full current block when the nanoparticle is captured. At pH 2.8, we suggest that Lys147 is free to engage in electrostatic interactions with sulfonates at the nanoparticle surface. To verify our results, we engineered a mutation in the α-hemolysin protein, where Glu111 is substituted by Ala (E111A), thus removing Glu111-Lys147 interactions and facilitating Lys147-sulfonate electrostatic interactions. This mutation leads to a 100% current block at pH 2.8 and a 92% block at pH 8.0, showing that electrostatic interactions are formed between the nanopore and the nanoparticle surface. Besides demonstrating the effect of electrostatic interactions on cross channel ionic current, this work offers a novel approach to controlling open and closed states of the α-hemolysin nanopore as a function of external gears.

Beyond Saffman-Delbruck approximation: a new regime for 2D diffusion of α-hemolysin complexes in supported lipid bilayer

Cell mechanisms are actively modulated by membrane dynamics. We studied the dynamics of a first-stage biomimetic system by Fluorescence Recovery After Patterned Photobleaching. Using this simple biomimetic system, constituted by α -hemolysin from Staphylococcus aureus inserted as single heptameric pore or complexes of pores in a glass-supported DMPC bilayer, we observed true diffusion behavior, with no immobile fraction. We find two situations: i) when incubation is shorter than 15 hours, the protein inserts as a heptameric pore and diffuses roughly three times more slowly than its host lipid bilayer; ii) incubation longer than 15 hours leads to the formation of larger complexes which diffuse more slowly. Our results indicate that, while the Saffman-Delbruck model adequately describes the diffusion coefficient D for small radii, D of the objects decreases as 1/R(2) for the size range explored in this study. Additionally, in the presence of inserted proteins, the gel-to-fluid transition of the supported bilayer as well as a temperature shift in the gel-to-fluid transition are observed.

Ion transport through a graphene nanopore

Molecular dynamics simulation is utilized to investigate the ionic transport of NaCl in solution through a graphene nanopore under an applied electric field. Results show the formation of concentration polarization layers in the vicinity of the graphene sheet. The nonuniformity of the ion distribution gives rise to an electric pressure which drives vortical motions in the fluid if the electric field is sufficiently strong to overcome the influence of viscosity and thermal fluctuations. The relative importance of hydrodynamic transport and thermal fluctuations in determining the pore conductivity is investigated. A second important effect that is observed is the mass transport of water through the nanopore, with an average velocity proportional to the applied voltage and independent of the pore diameter. The flux arises as a consequence of the asymmetry in the ion distribution which can be attributed to differing mobilities of the sodium and chlorine ions and to the polarity of water molecules. The accumulation of liquid molecules in the vicinity of the nanopore due to re-orientation of the water dipoles by the local electric field is seen to result in a local increase in the liquid density. Results confirm that the electric conductance is proportional to the nanopore diameter for the parameter regimes that we simulated. The occurrence of fluid vortices is found to result in an increase in the effective electrical conductance.

Mechanism of voltage gating in potassium channels

The mechanism of ion channel voltage gating-how channels open and close in response to voltage changes-has been debated since Hodgkin and Huxley's seminal discovery that the crux of nerve conduction is ion flow across cellular membranes. Using all-atom molecular dynamics simulations, we show how a voltage-gated potassium channel (KV) switches between activated and deactivated states. On deactivation, pore hydrophobic collapse rapidly halts ion flow. Subsequent voltage-sensing domain (VSD) relaxation, including inward, 15-angstrom S4-helix motion, completes the transition. On activation, outward S4 motion tightens the VSD-pore linker, perturbing linker-S6-helix packing. Fluctuations allow water, then potassium ions, to reenter the pore; linker-S6 repacking stabilizes the open pore. We propose a mechanistic model for the sodium/potassium/calcium voltage-gated ion channel superfamily that reconciles apparently conflicting experimental data.