Web interface for Brownian dynamics simulation of ion transport and its applications to beta-barrel pores

A Deep Dive into VDAC1 Conformational Diversity Using All-Atom Simulations Provides New Insights into the Structural Origin of the Closed States

The voltage-dependent anion channel 1 (VDAC1) is a crucial mitochondrial transporter that controls the flow of ions and respiratory metabolites entering or exiting mitochondria. As a voltage-gated channel, VDAC1 can switch between a high-conducting “open” state and a low-conducting “closed” state emerging at high transmembrane (TM) potentials. Although cell homeostasis depends on channel gating to regulate the transport of ions and metabolites, structural hallmarks characterizing the closed states remain unknown. Here, we performed microsecond accelerated molecular dynamics to highlight a vast region of VDAC1 conformational landscape accessible at typical voltages known to promote closure. Conformers exhibiting durable subconducting properties inherent to closed states were identified. In all cases, the low conductance was due to the particular positioning of an unfolded part of the N-terminus, which obstructed the channel pore. While the N-terminal tail was found to be sensitive to voltage orientation, our models suggest that stable low-conducting states of VDAC1 predominantly take place from disordered events and do not result from the displacement of a voltage sensor or a significant change in the pore. In addition, our results were consistent with conductance jumps observed experimentally and corroborated a recent study describing entropy as a key factor for VDAC gating.

Anthrax toxin channel: What we know based on over 30 years of research

Protective antigen channel is the central component of the deadly anthrax exotoxin responsible for binding and delivery of the toxin's enzymatic lethal and edema factor components into the cytosol. The channel, which is more than three times longer than the lipid bilayer membrane thickness and has a 6-Å limiting diameter, is believed to provide a sophisticated unfoldase and translocase machinery for the foreign protein transport into the host cell cytosol. The tripartite toxin can be reengineered, one component at a time or collectively, to adapt it for the targeted cancer therapeutic treatments. In this review, we focus on the biophysical studies of the protective antigen channel-forming activity, small ion transport properties, enzymatic factor translocation, and blockage comparing it with the related clostridial binary toxin channels. We address issues linked to the anthrax toxin channel structural dynamics and lipid dependence, which are yet to become generally recognized as parts of the toxin translocation machinery.

Hydrophobic Gating and 1/f Noise of the Anthrax Toxin Channel

"Pink" or 1/f noise is a natural phenomenon omnipresent in physics, economics, astrophysics, biology, and even music and languages. In electrophysiology, the stochastic activity of a number of biological ion channels and artificial nanopores could be characterized by current noise with a 1/f power spectral density. In the anthrax toxin channel (PA₆₃), it appears as fast voltage-independent current interruptions between conducting and nonconducting states. This behavior hampers potential development of PA₆₃ as an ion-channel biosensor. On the bright side, the PA₆₃ flickering represents a mesmerizing phenomenon to investigate. Notably, similar 1/f fluctuations are observed in the channel-forming components of clostridial binary C2 and iota toxins, which share functional and structural similarities with the anthrax toxin channel. Similar to PA₆₃, they are evolved to translocate the enzymatic components of the toxins into the cytosol. Here, using high-resolution single-channel lipid bilayer experiments and all-atom molecular dynamic simulations, we suggest that the 1/f noise in PA₆₃ occurs as a result of "hydrophobic gating" at the ϕ-clamp region, the phenomenon earlier observed in several water-filled channels "fastened" inside by the hydrophobic belts. The ϕ-clamp is a narrow "hydrophobic ring" in the PA₆₃ lumen formed by seven or eight phenylalanine residues at position 427, conserved in the C2 and iota toxin channels, which catalyzes protein translocation. Notably, the 1/f noise remains undetected in the F427A PA₆₃ mutant. This finding can elucidate the functional purpose of 1/f noise and its possible role in the transport of the enzymatic components of binary toxins.

The intrinsically disordered N-terminus of the voltage-dependent anion channel

The voltage-dependent anion channel (VDAC) is a critical β-barrel membrane protein of the mitochondrial outer membrane, which regulates the transport of ions and ATP between mitochondria and the cytoplasm. In addition, VDAC plays a central role in the control of apoptosis and is therefore of great interest in both cancer and neurodegenerative diseases. Although not fully understood, it is presumed that the gating mechanism of VDAC is governed by its N-terminal region which, in the open state of the channel, exhibits an α-helical structure positioned midway inside the pore and strongly interacting with the β-barrel wall. In the present work, we performed molecular simulations with a recently developed force field for disordered systems to shed new light on known experimental results, showing that the N-terminus of VDAC is an intrinsically disordered region (IDR). First, simulation of the N-terminal segment as a free peptide highlighted its disordered nature and the importance of using an IDR-specific force field to properly sample its conformational landscape. Secondly, accelerated dynamics simulation of a double cysteine VDAC mutant under applied voltage revealed metastable low conducting states of the channel representative of closed states observed experimentally. Related structures were characterized by partial unfolding and rearrangement of the N-terminal tail, that led to steric hindrance of the pore. Our results indicate that the disordered properties of the N-terminus are crucial to properly account for the gating mechanism of VDAC.

The voltage-dependent anion channel (VDAC) is a membrane protein playing a pivotal role in the transport of ions or ATP across the mitochondrial outer membrane as well as in the induction of apoptosis. At high enough membrane potential, VDAC is known to transition from an open state to multiple closed states, reducing the flow of ions through the channel and blocking the passage of large metabolites. While the structure of the open state was resolved more than a decade ago, a molecular description of the gating mechanism of the channel is still missing. Here we show that the N-terminus of VDAC is an intrinsically disordered region and that such a property has a profound impact on its dynamics either as a free peptide or as part of the channel. By taking disordered properties of the N-terminus into account, we managed to generate long-lived closed conformations of the channel at experimental values of the membrane potential. Our results provide new insights into the molecular mechanism driving the gating of VDAC.

Single-molecule sensing of peptides and nucleic acids by engineered aerolysin nanopores

Nanopore sensing is a powerful single-molecule approach for the detection of biomolecules. Recent studies have demonstrated that aerolysin is a promising candidate to improve the accuracy of DNA sequencing and to develop novel single-molecule proteomic strategies. However, the structure–function relationship between the aerolysin nanopore and its molecular sensing properties remains insufficiently explored. Herein, a set of mutated pores were rationally designed and evaluated in silico by molecular simulations and in vitro by single-channel recording and molecular translocation experiments to study the pore structural variation, ion selectivity, ionic conductance and capabilities for sensing several biomolecules. Our results show that the ion selectivity and sensing ability of aerolysin are mostly controlled by electrostatics and the narrow diameter of the double β-barrel cap. By engineering single-site mutants, a more accurate molecular detection of nucleic acids and peptides has been achieved. These findings open avenues for developing aerolysin nanopores into powerful sensing devices.

Aerolysin pores have potential to improve the accuracy of DNA sequencing and single-molecule proteomics. Here, the authors rationally design a set of mutated pores to achieve a more accurate detection of peptides and nucleic acids.

Multivalent Inhibitors of Channel-Forming Bacterial Toxins

Rational design of multivalent molecules represents a remarkable modern tool to transform weak non-covalent interactions into strong binding by creating multiple finely-tuned points of contact between multivalent ligands and their supposed multivalent targets. Here, we describe several prominent examples where the multivalent blockers were investigated for their ability to directly obstruct oligomeric channel-forming bacterial exotoxins, such as the pore-forming bacterial toxins and B component of the binary bacterial toxins. We address problems related to the blocker/target symmetry match and nature of the functional groups, as well as chemistry and length of the linkers connecting the functional groups to their multivalent scaffolds. Using the anthrax toxin and AB5 toxin case studies, we briefly review how the oligomeric toxin components can be successfully disabled by the multivalent non-channel-blocking inhibitors, which are based on a variety of multivalent scaffolds.

A Crucial Role for Side-Chain Conformation in the Versatile Charge Selectivity of Cys-Loop Receptors

Whether synaptic transmission is excitatory or inhibitory depends, to a large extent, on whether the ion channels that open upon binding the released neurotransmitter conduct cations or anions. The mechanistic basis of the opposite charge selectivities of Cys-loop receptors has only recently begun to emerge. It is now clear that ionized side chains-whether pore-facing or buried-in the first α-helical turn of the second transmembrane segments underlie this phenomenon and that the electrostatics of backbone atoms are not critically involved. Moreover, on the basis of electrophysiological observations, it has recently been suggested that not only the sign of charged side chains but also their conformation are crucial determinants of cation-anion selectivity. To challenge these ideas with the chemical and structural rigor that electrophysiological observations naturally lack, we performed molecular dynamics, Brownian dynamics, and electrostatics calculations of ion permeation. To this end, we used structural models of the open-channel conformation of the α1 glutamate-gated Cl^- channel and the α1 glycine receptor. Our results provided full support to the notion that the conformation of charged sides chains matters for charge selectivity. Indeed, whereas some rotamers of the buried arginines at position 0' conferred high selectivity for anions, others supported the permeation of cations and anions at similar rates or even allowed the faster permeation of cations. Furthermore, we found that modeling glutamates at position -1' of the anion-selective α1 glycine receptor open-state structure-instead of the five native alanines-switches charge selectivity also in a conformation-dependent manner, with some glutamate rotamers being much more effective at conferring selectivity for cations than others. Regarding pore size, we found that the mere expansion of the pore has only a minimal impact on cation-anion selectivity. Overall, these results bring to light the previously unappreciated impact of side-chain conformation on charge selectivity in Cys-loop receptors.

Lipid composition and salt concentration as regulatory factors of the anion selectivity of VDAC studied by coarse-grained molecular dynamics simulations

The voltage-dependent anion channel (VDAC) is a mitochondrial outer membrane protein whose fundamental function is to facilitate and regulate the flow of metabolites between the cytosol and the mitochondrial intermembrane space. Using coarse-grained molecular dynamics simulations, we investigated the dependence of VDAC selectivity towards small inorganic anions on two factors: the ionic strength and the lipid composition. In agreement with experimental data we found that VDAC becomes less anion selective with increasing salt concentration due to the screening of a few basic residues that point into the pore lumen. The molecular dynamics simulations provide insight into the regulation mechanism of VDAC selectivity by the composition in the lipid membrane and suggest that the ion distribution is differently modulated by POPE compared to the POPC bilayer. This occurs through the more persistent interactions of acidic residues located at both rims of the β-barrel with head groups of POPE which in turn impact the electrostatic potential and thereby the selectivity of the pore. This mechanism occurs not only in POPE single component membranes but also in a mixed POPE/POPC bilayer by an enrichment of POPE over POPC lipids on the surface of VDAC. Thus we show here that computationally-inexpensive coarse-grained simulations are able to capture, in a semi-quantitative way, essential features of VDAC anion selectivity and could pave the way toward a molecular level understanding of metabolite transport in natural membranes.

Effect of late endosomal DOBMP lipid and traditional model lipids of electrophysiology on the anthrax toxin channel activity

Anthrax toxin action requires triggering of natural endocytic transport mechanisms whereby the binding component of the toxin forms channels (PA₆₃) within endosomal limiting and intraluminal vesicle membranes to deliver the toxin's enzymatic components into the cytosol. Membrane lipid composition varies at different stages of anthrax toxin internalization, with intraluminal vesicle membranes containing ~70% of anionic bis(monoacylglycero)phosphate lipid. Using model bilayer measurements, we show that membrane lipids can have a strong effect on the anthrax toxin channel properties, including the channel-forming activity, voltage-gating, conductance, selectivity, and enzymatic factor binding. Interestingly, the highest PA₆₃ insertion rate was observed in bis(monoacylglycero)phosphate membranes. The molecular dynamics simulation data show that the conformational properties of the channel are different in bis(monoacylglycero)phosphate compared to PC, PE, and PS lipids. The anthrax toxin protein/lipid bilayer system can be advanced as a novel robust model to directly investigate lipid influence on membrane protein properties and protein/protein interactions.

Effect of endosomal acidification on small ion transport through the anthrax toxin PA63 channel

Tight regulation of pH is critical for the structure and function of cells and organelles. The pH environment changes dramatically along the endocytic pathway, an internalization transport process that is 'hijacked' by many intracellularly active bacterial exotoxins, including the anthrax toxin. Here, we investigate the role of pH on single-channel properties of the anthrax toxin protective antigen (PA₆₃ ). Using conductance and current noise analysis, blocker binding, ion selectivity, and poly(ethylene glycol) partitioning measurements, we show that the channel exists in two different open states ('maximum' and 'main') at pH ≥ 5.5, while only a maximum conductance state is detected at pH < 5.5. We describe two substantially distinct patterns of PA₆₃ conductance dependence on KCl concentration uncovered at pH 6.5 and 4.5.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

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

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

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.

Concerted motions networking pores and distant ferroxidase centers enable bacterioferritin function and iron traffic

X-ray crystallography, molecular dynamics (MD) simulations, and biochemistry were utilized to investigate the effect of introducing hydrophobic interactions in the 4-fold (N148L and Q151L) and B-pores (D34F) of Pseudomonas aeruginosa bacterioferritin B (BfrB) on BfrB function. The structures show only local structural perturbations and confirm the anticipated hydrophobic interactions. Surprisingly, structures obtained after soaking crystals in Fe2+-containing crystallization solution revealed that although iron loads into the ferroxidase centers of the mutants, the side chains of ferroxidase ligands E51 and H130 do not reorganize to bind the iron ions, as is seen in the wt BfrB structures. Similar experiments with a double mutant (C89S/K96C) prepared to introduce changes outside the pores show competent ferroxidase centers that function akin to those in wt BfrB. MD simulations comparing wt BfrB with the D34F and N148L mutants show that the mutants exhibit significantly reduced flexibility and reveal a network of concerted motions linking ferroxidase centers and 4-fold and B-pores, which are important for imparting ferroxidase centers in BfrB with the required flexibility to function efficiently. In agreement, the efficiency of Fe2+ oxidation and uptake of the 4-fold and B-pore mutants in solution is significantly compromised relative to wt or C89S/K96C BfrB. Finally, our structures show a large number of previously unknown iron binding sites in the interior cavity and B-pores of BfrB, which reveal in unprecedented detail conduits followed by iron and phosphate ions across the BfrB shell, as well as paths in the interior cavity that may facilitate nucleation of the iron phosphate mineral.

BROMOC suite: Monte Carlo/Brownian dynamics suite for studies of ion permeation and DNA transport in biological and artificial pores with effective potentials

The transport of ions and solutes by biological pores is central for cellular processes and has a variety of applications in modern biotechnology. The time scale involved in the polymer transport across a nanopore is beyond the accessibility of conventional MD simulations. Moreover, experimental studies lack sufficient resolution to provide details on the molecular underpinning of the transport mechanisms. BROMOC, the code presented herein, performs Brownian dynamics simulations, both serial and parallel, up to several milliseconds long. BROMOC can be used to model large biological systems. IMC-MACRO software allows for the development of effective potentials for solute-ion interactions based on radial distribution function from all-atom MD. BROMOC Suite also provides a versatile set of tools to do a wide variety of preprocessing and postsimulation analysis. We illustrate a potential application with ion and ssDNA transport in MspA nanopore.

Microsecond simulations of DNA and ion transport in nanopores with novel ion-ion and ion-nucleotides effective potentials

We developed a novel scheme based on the grand-canonical Monte Carlo/Brownian dynamics simulations and have extended it to studies of ion currents across three nanopores with the potential for single-stranded DNA (ssDNA) sequencing: solid-state nanopore Si₃N₄, α-hemolysin, and E111N/M113Y/K147N mutant. To describe nucleotide-specific ion dynamics compatible with ssDNA coarse-grained model, we used the inverse Monte Carlo protocol, which maps the relevant ion-nucleotide distribution functions from all-atom molecular dynamics (MD) simulations. Combined with the previously developed simulation platform for Brownian dynamics simulations of ion transport, it allows for microsecond- and millisecond-long simulations of ssDNA dynamics in the nanopore with a conductance computation accuracy that equals or exceeds that of all-atom MD simulations. In spite of the simplifications, the protocol produces results that agree with the results of previous studies on ion conductance across open channels and provide direct correlations with experimentally measured blockade currents and ion conductances that have been estimated from all-atom MD simulations.

CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residues

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface to prepare molecular simulation systems and input files to facilitate the usage of common and advanced simulation techniques. Since it is originally developed in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to setup a broad range of simulations including free energy calculation and large-scale coarse-grained representation. Here, we describe functionalities that have recently been integrated into CHARMM-GUI PDB Manipulator, such as ligand force field generation, incorporation of methanethiosulfonate spin labels and chemical modifiers, and substitution of amino acids with unnatural amino acids. These new features are expected to be useful in advanced biomolecular modeling and simulation of proteins.

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.

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.

Side-chain conformation at the selectivity filter shapes the permeation free-energy landscape of an ion channel

On the basis of single-channel currents recorded from the muscle nicotinic acetylcholine receptor (AChR), we have recently hypothesized that the conformation adopted by the glutamate side chains at the first turn of the pore-lining α-helices is a key determinant of the rate of ion permeation. In this paper, we set out to test these ideas within a framework of atomic detail and stereochemical rigor by conducting all-atom molecular dynamics and Brownian dynamics simulations on an extensively validated model of the open-channel muscle AChR. Our simulations provided ample support to the notion that the different rotamers of these glutamates partition into two classes that differ markedly in their ability to catalyze ion conduction, and that the conformations of the four wild-type glutamates are such that two of them "fall" in each rotamer class. Moreover, the simulations allowed us to identify the mm (χ(1) ≅ -60°; χ(2) ≅ -60°) and tp (χ(1) ≅ 180°; χ(2) ≅ +60°) rotamers as the likely conduction-catalyzing conformations of the AChR's selectivity-filter glutamates. More generally, our work shows an example of how experimental benchmarks can guide molecular simulations into providing a type of structural and mechanistic insight that seems otherwise unattainable.

Cationic PAMAM dendrimers as pore-blocking binary toxin inhibitors

Dendrimers are unique highly branched macromolecules with numerous groundbreaking biomedical applications under development. Here we identified poly(amido amine) (PAMAM) dendrimers as novel blockers for the pore-forming B components of the binary anthrax toxin (PA63) and Clostridium botulinum C2 toxin (C2IIa). These pores are essential for delivery of the enzymatic A components of the internalized toxins from endosomes into the cytosol of target cells. We demonstrate that at low μM concentrations cationic PAMAM dendrimers block PA63 and C2IIa to inhibit channel-mediated transport of the A components, thereby protecting HeLa and Vero cells from intoxication. By channel reconstitution and high-resolution current recording, we show that the PAMAM dendrimers obstruct transmembrane PA63 and C2IIa pores in planar lipid bilayers at nM concentrations. These findings suggest a new potential role for the PAMAM dendrimers as effective polyvalent channel-blocking inhibitors, which can protect human target cells from intoxication with binary toxins from pathogenic bacteria.

Origin of ion selectivity in Phaseolus coccineus mitochondrial VDAC

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

Theoretical frameworks for multiscale modeling and simulation

Biomolecular systems have been modeled at a variety of scales, ranging from explicit treatment of electrons and nuclei to continuum description of bulk deformation or velocity. Many challenges of interfacing between scales have been overcome. Multiple models at different scales have been used to study the same system or calculate the same property (e.g., channel conductance). Accurate modeling of biochemical processes under in vivo conditions and the bridging of molecular and subcellular scales will likely soon become reality.

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.

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.

Molecular origin of VDAC selectivity towards inorganic ions: a combined molecular and Brownian dynamics study

The voltage-dependent anion channel (VDAC) serves as the major pore for metabolites and electrolytes in the outer mitochondrial membrane. To refine our understanding of ion permeation through this channel we performed an extensive Brownian (BD) and molecular dynamics (MD) study on the mouse VDAC isoform 1 wild-type and mutants (K20E, D30K, K61E, E158K and K252E). The selectivity and the conductance of the wild-type and of the variant channels computed from the BD trajectories are in agreement with experimental data. The calculated selectivity is shown to be very sensitive to slight conformational changes which may have some bearing on the variability of the selectivity values measured on the VDAC open state. The MD and BD free energy profiles of the ion permeation suggest that the pore region comprising the N-terminal helix and the barrel band encircling it predominantly controls the ion transport across the channel. The overall 12μs BD and 0.9μs MD trajectories of the mouse VDAC isoform 1 wild-type and mutants feature no distinct pathways for ion diffusion and no long-lived ion-protein interactions. The dependence of ion distribution in the wild-type channel with the salt concentration can be explained by an ionic screening of the permanent charges of the protein arising from the pore. Altogether these results bolster the role of electrostatic features of the pore as the main determinant of VDAC selectivity towards inorganic anions.

Protein dynamics and ion traffic in bacterioferritin

Bacterioferritin (Bfr) is a spherical protein composed of 24 subunits and 12 heme molecules. Bfrs contribute to regulate iron homeostasis in bacteria by capturing soluble but potentially toxic Fe(2+) and by compartmentalizing it in the form of a bioavailable ferric mineral inside the protein's hollow cavity. When iron is needed, Fe(3+) is reduced and mobilized into the cytosol as Fe(2+). Hence, key to the function of Bfr is its ability to permeate iron ions in and out of its interior cavity, which is likely imparted by a flexible protein shell. To examine the conformational flexibility of Bfrs in a native-like environment and the way in which the protein shell interacts with monovalent cations, we have performed molecular dynamics (MD) simulations of BfrB from Pseudomonas aeruginosa (Pa BfrB) in K(2)HPO(4) solutions at different ionic strengths. The results indicate the presence of coupled thermal fluctuations (dynamics) in the 4-fold pores and B-pores of the protein, which is key to allowing passage of monovalent cations through the protein shell using B-pores as conduits. The MD simulations also show that Pa BfrB ferroxidase centers are highly dynamic and permanently populated by transient cations exchanging with other cations in the interior cavity, as well as the solution bathing the protein. Taken together, the findings suggest that Fe(2+) passes across the Pa BfrB shell via B-pores and that the ferroxidase pores allow the capture and oxidation of Fe(2+), followed by translocation of Fe(3+) to the interior cavity, aided by the conformationally active H130.

BROMOC-D: Brownian Dynamics/Monte-Carlo Program Suite to Study Ion and DNA Permeation in Nanopores

A theoretical framework is presented to model ion and DNA translocation across a nanopore confinement under an applied electric field. A combined Grand Canonical Monte Carlo Brownian Dynamics (GCMC/BD) algorithm offers a general approach to study ion permeation through wide molecular pores with a direct account of ion-ion and ion-DNA correlations. This work extends previously developed theory by incorporating the recently developed coarse-grain polymer model of DNA by de Pablo and colleagues [Knotts, T. A.; Rathore, N.; Schwartz, D. C.; de Pablo, J. J. J. Chem. Phys. 2007, 126] with explicit ions for simulations of polymer dynamics. Atomistic MD simulations were used to guide model developments. The power of the developed scheme is illustrated with studies of single-stranded DNA (ss-DNA) oligomer translocation in two model cases: a cylindrical pore with a varying radius and a well-studied experimental system, the staphylococcal α-hemolysin channel. The developed model shows good agreement with experimental data for model studies of two homopolymers: ss-poly(dA)(n) and ss-poly(dC)(n). The developed protocol allows for direct evaluation of different factors (charge distribution and pore shape and size) controlling DNA translocation in a variety of nanopores.