Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4-5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations-except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.

The voltage-dependent anion channel (VDAC) in the outer membrane of mitochondria serves an essential role in the transport of metabolites and electrolytes between the cell matrix and mitochondria. To examine its structure, dynamics, and the mechanisms underlying its electrophysiological properties, we performed a total of 1.77 μs molecular dynamics simulations of human VDAC isoform 1 in DOPE/DOPC mixed bilayers in 1 M KCl solution with transmembrane potentials of 0, ±25, ±50, ±75, and ±100 mV. The calculated conductance and ion selectivity are in good agreement with the experimental measurements. In addition, ion density distributions inside the channel reveal possible pathways for different ion species. Based on these observations, a mechanism underlying the anion selectivity is proposed; both ion species are transported across the channel, but the rate for K(+) is smaller than that for Cl(-) because of the attractive interactions between K(+) and residues on the channel wall. This difference leads to the anion selectivity of VDAC.

It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.

Brownian dynamics (BD) in a suitably constructed potential of mean force is an efficient and accurate method for simulating ion transport through wide ion channels. Here, a web-based graphical user interface (GUI) is presented for grand canonical Monte Carlo (GCMC) BD simulations of channel proteins: http://www.charmm-gui.org/input/gcmcbd. The webserver is designed to help users avoid most of the technical difficulties and issues encountered in setting up and simulating complex pore systems. GCMC/BD simulation results for three proteins, the voltage dependent anion channel (VDAC), α-Hemolysin, and the protective antigen pore of the anthrax toxin (PA), are presented to illustrate system setup, input preparation, and typical output (conductance, ion density profile, ion selectivity, and ion asymmetry). Two models for the input diffusion constants for potassium and chloride ions in the pore are compared: scaling of the bulk diffusion constants by 0.5, as deduced from previous all-atom molecular dynamics simulations of VDAC; and a hydrodynamics based model (HD) of diffusion through a tube. The HD model yields excellent agreement with experimental conductances for VDAC and α-Hemolysin, while scaling bulk diffusion constants by 0.5 leads to underestimates of 10–20%. For PA, simulated ion conduction values overestimate experimental values by a factor of 1.5 to 7 (depending on His protonation state and the transmembrane potential), implying that the currently available computational model of this protein requires further structural refinement.

The voltage dependent anion channel (VDAC) in the outer membrane of mitochondria serves an essential role in transport of metabolites and electrolytes between the cell matrix and mitochondria. To examine its structure, dynamics, and mechanisms underlying its electrophysiological properties, we have performed a total of 1.77 ms molecular dynamics simulations of human VDAC isoform 1 in DOPE/DOPC mixed bilayers in 1M KCl solution with transmembrane potentials of 0, 525, 550, 575, and 5100 mV. The calculated conductance and ion selectivity are in good agreement with the experimental measurements. In addition, ion density distributions inside the channel reveal possible pathways for different ion species. Based on these observations, a mechanism underlying the anion selectivity is proposed; both ion species are transported across the channel, but the rate for Kþ is smaller than that of Cl-because of the attractive interactions between Kþ and residues on the channel wall. This difference leads to the anion selectivity of VDAC.