Exploration of the structural features defining the conduction properties of a synthetic ion channel

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

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

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

De novo protein design, a retrospective

Proteins are molecular machines whose function depends on their ability to achieve complex folds with precisely defined structural and dynamic properties. The rational design of proteins from first-principles, or de novo, was once considered to be impossible, but today proteins with a variety of folds and functions have been realized. We review the evolution of the field from its earliest days, placing particular emphasis on how this endeavor has illuminated our understanding of the principles underlying the folding and function of natural proteins, and is informing the design of macromolecules with unprecedented structures and properties. An initial set of milestones in de novo protein design focused on the construction of sequences that folded in water and membranes to adopt folded conformations. The first proteins were designed from first-principles using very simple physical models. As computers became more powerful, the use of the rotamer approximation allowed one to discover amino acid sequences that stabilize the desired fold. As the crystallographic database of protein structures expanded in subsequent years, it became possible to construct proteins by assembling short backbone fragments that frequently recur in Nature. The second set of milestones in de novo design involves the discovery of complex functions. Proteins have been designed to bind a variety of metals, porphyrins, and other cofactors. The design of proteins that catalyze hydrolysis and oxygen-dependent reactions has progressed significantly. However, de novo design of catalysts for energetically demanding reactions, or even proteins that bind with high affinity and specificity to highly functionalized complex polar molecules remains an importnant challenge that is now being achieved. Finally, the protein design contributed significantly to our understanding of membrane protein folding and transport of ions across membranes. The area of membrane protein design, or more generally of biomimetic polymers that function in mixed or non-aqueous environments, is now becoming increasingly possible.

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

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

Molecular dynamics simulations of homo-oligomeric bundles embedded within a lipid bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)3-NH2 (LS2) or Ac-(LSSLLSL)3-NH2 (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine-serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a "dimer of dimers." In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.

Interacting ions in biophysics: real is not ideal

Ions in water are important throughout biology, from molecules to organs. Classically, ions in water were treated as ideal noninteracting particles in a perfect gas. Excess free energy of each ion was zero. Mathematics was not available to deal consistently with flows, or interactions with other ions or boundaries. Nonclassical approaches are needed because ions in biological conditions flow and interact. The concentration gradient of one ion can drive the flow of another, even in a bulk solution. A variational multiscale approach is needed to deal with interactions and flow. The recently developed energetic variational approach to dissipative systems allows mathematically consistent treatment of the bio-ions Na(+), K(+), Ca(2+), and Cl(-) as they interact and flow. Interactions produce large excess free energy that dominate the properties of the high concentration of ions in and near protein active sites, ion channels, and nucleic acids: the number density of ions is often >10 M. Ions in such crowded quarters interact strongly with each other as well as with the surrounding protein. Nonideal behavior found in many experiments has classically been ascribed to allosteric interactions mediated by the protein and its conformation changes. The ion-ion interactions present in crowded solutions-independent of conformation changes of the protein-are likely to change the interpretation of many allosteric phenomena. Computation of all atoms is a popular alternative to the multiscale approach. Such computations involve formidable challenges. Biological systems exist on very different scales from atomic motion. Biological systems exist in ionic mixtures (like extracellular and intracellular solutions), and usually involve flow and trace concentrations of messenger ions (e.g., 10(-7) M Ca(2+)). Energetic variational methods can deal with these characteristic properties of biological systems as we await the maturation and calibration of all-atom simulations of ionic mixtures and divalents.

Rectification of the channelrhodopsin early conductance

We analyzed the nonlinear current-voltage relationships of the early conducting state of channelrhodopsin-2 expressed in Xenopus oocytes and human embryonic kidney cells with respect to changes of the electrochemical gradients of H(+), Na(+)/K(+), and Ca(2+)/Mg(2+). Several models were tested for wild-type ChR2 and mutations at positions E90, E123, H134, and T159. Voltage-gating was excluded as cause for the nonlinearity. However, a general enzyme kinetic model with one predominant binding site yielded good fits throughout. The empty site with an apparent charge number of about -0.3 and strong external cation binding causes some inward rectification of the uniport function. Additional inward rectification is due to asymmetric competition from outside between the transported ion species. Significant improvement of the fits was achieved by introducing an elastic voltage-divider formed by the voltage-sensitive barriers.

Soft wall ion channel in continuum representation with application to modeling ion currents in α-hemolysin

A soft repulsion (SR) model of short-range interactions between mobile ions and protein atoms is introduced in the framework of continuum representation of the protein and solvent. The Poisson-Nernst-Plank (PNP) theory of ion transport through biological channels is modified to incorporate this soft wall protein model. Two sets of SR parameters are introduced. The first is parametrized for all essential amino acid residues using all atom molecular dynamic simulations; the second is a truncated Lennard-Jones potential. We have further designed an energy-based algorithm for the determination of the ion accessible volume, which is appropriate for a particular system discretization. The effects of these models of short-range interactions were tested by computing current-voltage characteristics of the α-hemolysin channel. The introduced SR potentials significantly improve prediction of channel selectivity. In addition, we studied the effect of the choice of some space-dependent diffusion coefficient distributions on the predicted current-voltage properties. We conclude that the diffusion coefficient distributions largely affect total currents and have little effect on rectifications, selectivity, or reversal potential. The PNP-SR algorithm is implemented in a new efficient parallel Poisson, Poisson-Boltzmann, and PNP equation solver, also incorporated in a graphical molecular modeling package HARLEM.

Model generation of viral channel forming 2B protein bundles from polio and coxsackie viruses

2B is a 99 amino acid membrane protein encoded by enteroviruses such as polio and coxsackie viruses with two transmembrane domains. The protein is found to make membranes of infected cells permeable. Using a computational approach which positions the models and assesses stability by molecular dynamics (MD) simulations a putative tetrameric bundle model of 2B is generated. The bundles show a pore lining motif of three lysines followed by a serine. The bundle is discussed in terms of different possible orientations of the helices in the membrane and the consequences this has on the in vivo activity of 2B.

Viral channel-forming proteins

Channel-forming proteins are found in a number of viral genomes. In some cases, their role in the viral life cycle is well understood, in some cases it needs still to be elucidated. A common theme is that their mode of action involves a change of electrochemical or proton gradient across the lipid membrane which modulates the viral or cellular activity. Blocking these proteins can be a suitable therapeutic strategy as for some viruses this may be "lethal." Besides the many biological relevant questions still to be answered, there are also many open questions concerning the biophysical side as well as structural information and the mechanism of function on a molecular level. The immanent biophysical issues are addressed and the work in the field is summarized.

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.

Self-energy-limited ion transport in subnanometer channels

The current-voltage characteristics of the alpha-Hemolysin protein pore during the passage of single-stranded DNA under varying ionic strength C are studied experimentally. We observe strong blockage of the current, weak superlinear growth of the current as a function of voltage, and a minimum of the current as a function of C. These observations are interpreted as the result of the ion electrostatic self-energy barrier originating from the large difference in the dielectric constants of water and the lipid bilayer. The dependence of DNA capture rate on C also agrees with our model.

Ion exchange phase transitions in water-filled channels with charged walls

Ion transport through narrow water-filled channels is impeded by a high electrostatic barrier. The latter originates from the large ratio of the dielectric constants of the water and the surrounding media. We show that "doping," i.e., immobile charges attached to the walls of the channel, substantially reduces the barrier. This explains why most of the biological ion channels are "doped." We show that at rather generic conditions the channels may undergo ion exchange phase transitions (typically of the first order). Upon such a transition a finite latent concentration of ions may either enter or leave the channel, or be exchanged between the ions of different valences. We discuss possible implications of these transitions for the Ca-vs-Na selectivity of biological Ca channels. We also show that transport of divalent Ca ions is assisted by their fractionalization into two separate excitations.

Topological Equilibria of Ion Channel Peptides in Oriented Lipid Bilayers Revealed by 15N Solid-State NMR Spectroscopy

Ion channel peptides have been prepared by solid-phase peptide synthesis, labeled with 15N at selected sites, and reconstituted into oriented lipid bilayers. The (Leu-Ser-Ser-Leu-Leu-Ser-Leu)3-CONH2 peptide has previously been shown to exhibit well-defined and discrete ionic conductances when investigated by single-channel measurements [Lear, J. D., et al. (1988) Science 240, 1177]. Proton-decoupled 15N solid-state NMR spectroscopy indicates that (Leu-Ser-Ser-Leu-Leu-Ser-Leu)3-CONH2 preferentially aligns parallel to the membrane surface in excellent agreement with its amphipathic helical structure. However, by carefully choosing the conditions of the membrane environment, significant contributions that are indicative of transmembrane alignments become obvious in the 15N chemical shift solid-state NMR spectra. The data thereby provide experimental evidence for an equilibrium between in-plane and transmembrane-oriented helix configurations where the transmembrane and surface-oriented peptide fractions are in slow exchange. Similar topological equilibria are observed when the N-terminus of the LS21 peptide is acetylated. These observations provide experimental support for previous models, suggesting that the channels observed in single-channel conductance measurements are indeed formed by hexameric transmembrane helical bundles. In contrast, the shorter peptide (Leu-Ser-Ser-Leu-Leu-Ser-Leu)2-CONH2 is oriented parallel to the membrane surface under all conditions tested. This peptide exhibits erratic conductance changes when investigated by electrophysiological methods, probably because it is too short to span the lipid bilayer.

Computing numerically the access resistance of a pore

The access resistance (AR) of a channel is an important component of the conductance of ion channels, particularly in wide and short channels, where it accounts for a substantial fraction of the total resistance to the movement of ions. The AR is usually calculated by using a classical and simple expression derived by Hall from electrostatics (J.E. Hall 1975 J. Gen. Phys. 66:531-532), though other expressions, both analytical and numerical, have been proposed. Here we report some numerical results for the AR of a channel obtained by solving the Poisson-Nernst-Planck equations at the entrance of a circular pore. Agreement is found between numerical calculations and analytical results from Hall's equation for uncharged pores in neutral membranes. However, for channels embedded in charged membranes, Hall's expression overestimates the AR, which is much lower and can even be neglected in some cases. The weak dependence of AR on the pore radius for charged membranes at low salt concentration can be exploited to separate the channel and the access contributions to the measured conductance.

Poisson–Nernst–Planck Theory Approach to the Calculation of Current Through Biological Ion Channels

The Poisson-Nernst-Planck (PNP) theory of electro-diffusion is reviewed. Techniques for numerical solution of the three-dimensional PNP equations are summarized, and several illustrative applications to ion transport through protein channels are presented. Strengths and weaknesses of the theory are discussed, as well as attempts to improve it via increasingly realistic evaluation of the force acting on each ion due to the protein/membrane environment.

Structure and dynamics of model pore insertion into a membrane

A cylindrical transmembrane molecule is constructed by linking hydrophobic sites selected from a coarse grain model. The resulting hollow tube assembly serves as a representation of a transmembrane channel, pore, or a carbon nanotube. The interactions of a coarse grain di-myristoyl-phosphatidyl-choline hydrated bilayer with both a purely hydrophobic tube and a tube with hydrophilic caps are studied. The hydrophobic tube rotates in the membrane and becomes blocked by lipid tails after a few tens of nanoseconds. The hydrophilic sites of the capped tube stabilize it by anchoring the tube in the lipid headgroup/water interfacial region of each membrane leaflet. The capped tube remains free of lipid tails. The capped tube spontaneously conducts coarse grain water sites; the free-energy profile of this process is calculated using three different methods and is compared to the barrier for water permeation through the lipid bilayer. Spontaneous tube insertion into an undisturbed lipid bilayer is also studied, which we reported briefly in a previous publication. The hydrophobic tube submerges into the membrane core in a carpetlike manner. The capped tube laterally fuses with the closest leaflet, and then, after plunging into the membrane interior, rotates to assume a transbilayer orientation. Two lipids become trapped at the end of the tube as it penetrates the membrane. The hydrophilic headgroups of these lipids associate with the lower tube cap and assist the tube in crossing the interior of the membrane. When the rotation is complete these lipids detach from the tube caps and fuse with the lower leaflet lipids.

Gramicidin-based channel systems for the detection of protein-ligand interaction

To detect protein-ligand interaction a gramicidin-based sensor was developed. Biotin was tagged to the C-terminus of gramicidin (Gram-bio 1). The biotin-moiety, which faces the electrolyte, gave little effect on single-channel conductance. Streptavidin added to the electrolyte was detected by Gram-bio 1 through the monitoring channel current using the planar bilayer system. The suppression of macroscopic currents and the acceleration of their decaying time course were observed in a concentration dependent manner. In the single-channel level, however, no significant effect on the single-channel conductance and the open dwell time was observed upon addition of streptavidin. Therefore, streptavidin neither blocked the open channel nor changed the stability of the conducting dimer. Insertion of a linker between gramicidin and biotin did not change the streptavidin-sensitivity of the current reduction. We conclude that the binding of streptavidin to the Gram-bio 1 shifted the distribution of the complex from the membrane to the electrolyte and, thus, reduced the formation of conducting dimer of Gram-bio 1 in the membrane. Interaction of biotin with an anti-biotin antibody was also observed using this system, indicating that this system is applicable for the detection of protein-ligand interaction having a binding constant of approximately 10(8-9) M(-1) or more. Both the adamantane-tagged gramicidin for detection of beta-cyclodextrin and the Strep Tag-II-tagged gramicidin for detection of streptavidin (binding constant: approximately 10(5) M(-1) or less) failed to respond. Thus, high-affinity ligands upon tagging to gramicidin render the gramicidin-based sensor able to execute as a real-time monitoring system for protein-ligand interaction.

Understanding nature's design for a nanosyringe

Synthetic and natural peptide assemblies can possess transport or conductance activity across biomembranes through the formation of nanopores. The fundamental mechanisms of membrane insertion necessary for antimicrobial or synthetic pore formation are poorly understood. We observe a lipid-assisted mechanism for passive insertion into a model membrane from molecular dynamics simulations. The assembly used in the study, a generic nanotube functionalized with hydrophilic termini, is assisted in crossing the membrane core by transleaflet lipid flips. Lipid tails occlude a purely hydrophobic nanotube. The observed insertion mechanism requirements for hydrophobic-hydrophilic matching have implications for the design of synthetic channels and antibiotics.

Design, synthesis, and characterization of peptide nanostructures having ion channel activity

We report the synthesis and the functional studies of multiple crown alpha-helical peptides designed to form artificial ion channels. The approach combines the versatility of solid phase peptide synthesis, the conformational predictability of peptidic molecules, and the solution synthesis of crown ethers with engineerable ion-binding abilities. Several biophysical methods were employed to characterize the activity and the mode of action of these crown peptide nanostructures. The 21 residue peptides bearing six 21-EC-7 turned out to facilitate the translocation of ions in a similar fashion to natural ion channels.

Electrodiffusion near an ion channel and the effect of mobile buffer

A numerical method was set up to calculate the dynamic concentration behavior of charged particles in the vicinity of an ion channel. It takes into account the electric potential due to the charge of the transported ions. Additionally, the finite on- and off-kinetics of mobile ion-buffers such as EGTA can be added to the simulations. The calculations were carried out using a modified Crank-Nicolson algorithm to solve the partial differential equations describing the problem. It was found that the electrostatic effect on the concentration of permeating ions is negligible in the presence of physiological salt concentrations. Nevertheless, there are electrostatic effects on other ion species near the channel mouth. Studies on the effect of a Ba2+ -current through a Ca2+ -channel onto the Ca2+ -concentration in the bath, and on the amplification of the Ca2+ -effect on the BK-channel due to the K+ -flux are presented. Additionally, the effect of mobile buffers was simulated and the numerical results are compared with some common analytical approximations.

Improved 3D continuum calculations of ion flux through membrane channels

A continuum model, based on the Poisson-Nernst-Planck (PNP) theory, is applied to simulate steady-state ion flux through protein channels. The PNP equations are modified to explicitly account (1) for the desolvation of mobile ions in the membrane pore and (2) for effects related to ion sizes. The proposed algorithm for a three-dimensional self-consistent solution of PNP equations, in which final results are refined by a focusing technique, is shown to be suitable for arbitrary channel geometry and arbitrary protein charge distribution. The role of the pore shape and protein charge distribution in formation of basic electrodiffusion properties, such as channel conductivity and selectivity, as well as concentration distributions of mobile ions in the pore region, are illustrated by simulations on model channels. The influence of the ionic strength in the bulk solution and of the externally applied electric field on channel properties are also discussed.

Dielectric boundary force and its crucial role in gramicidin

In an electrostatic problem with nonuniform geometry, a charge Q in one region induces surface charges [called dielectric boundary charges (DBC)] at boundaries between different dielectrics. These induced surface charges, in return, exert a force [called dielectric boundary force (DBF)] on the charge Q that induced them. The DBF is often overlooked. It is not present in standard continuum theories of (point) ions in or near membranes and proteins, such as Gouy-Chapman, Debye-Huckel, Poisson-Boltzmann or Poisson-Nernst- Planck. The DBF is important when a charge Q is near dielectric interfaces, for example, when ions permeate through protein channels embedded in biological membranes. In this paper, we define the DBF and calculate it explicitly for a planar dielectric wall and for a tunnel geometry resembling the ionic channel gramicidin. In general, we formulate the DBF in a form useful for continuum theories, namely, as a solution of a partial differential equation with boundary conditions. The DBF plays a crucial role in the permeation of ions through the gramicidin channel. A positive ion in the channel produces a DBF of opposite sign to that of the fixed charge force (FCF) produced by the permanent charge of the gramicidin polypeptide, and so the net force on the positive ion is reduced. A negative ion creates a DBF of the same sign as the FCF and so the net (repulsive) force on the negative ion is increased. Thus, a positive ion can permeate the channel, while a negative ion is excluded from it. In gramicidin, it is this balance between the FCF and DBF that allows only singly charged positive ions to move into and through the channel. The DBF is not directly responsible, however, for selectivity between the alkali metal ions (e.g., Li+, Na+, K+): we prove that the DBF on a mobile spherical ion is independent of the ion's radius.

Multistep mechanism of chloride translocation in a strongly anion-selective porin channel

The strongly anion-selective porin channel Omp32 from the bacterium Delftia acidovorans differs from other unspecific porins by its pronounced selectivity for anions and its particularly small channel cross-section. Multinanosecond molecular dynamics simulations of chloride ion movement in this pore protein suggest that translocated anions interact intimately with the charges of a "basic ladder", whose dynamics lead the anions in a stepwise manner through the constriction zone of the channel. The ladder-steps comprise the central clustered arginine groups and flanking basic residues at its exoplasmic and periplasmic sides. The computed free energy profile of ion movement in and around the constriction zone shows a corresponding succession of free energy minima and barriers. A number of polar atoms from other amino acids contribute to the coordination of Cl(-) at certain sites and to its temporary immobilization in the channel. A special binding site occurs at the transition of the constriction zone to the periplasmic funnel, binding the chloride ion over significant lengths of time. The results from our MD study offer a possible explanation for the nonlinear conductance properties and unusual salt-dependent characteristics of Omp32 observed earlier in experimental measurements.

The role of the dielectric barrier in narrow biological channels: a novel composite approach to modeling single-channel currents

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.

Controlled assembly of carbon nanotubes by designed amphiphilic Peptide helices

Carbon nanotubes have properties potentially useful in diverse electrical and mechanical nanoscale devices and for making strong, light materials. However, carbon nanotubes are difficult to solubilize and organize into architectures necessary for many applications. In the present paper, we describe an amphiphilic alpha-helical peptide specifically designed not only to coat and solubilize carbon nanotubes, but also to control the assembly of the peptide-coated nanotubes into macromolecular structures through peptide-peptide interactions between adjacent peptide-wrapped nanotubes. The data presented herein show that the peptide folds into an amphiphilic alpha-helix in the presence of carbon nanotubes and disperses them in aqueous solution by noncovalent interactions with the nanotube surface. Electron microscopy and polarized Raman studies reveal that the peptide-coated nanotubes assemble into fibers with the nanotubes aligned along the fiber axis. Most importantly, the size and morphology of the fibers can be controlled by manipulating solution conditions that affect peptide-peptide interactions.

Analysis and evaluation of channel models: simulations of alamethicin

Alamethicin is an antimicrobial peptide that forms stable channels with well-defined conductance levels. We have used extended molecular dynamics simulations of alamethicin bundles consisting of 4, 5, 6, 7, and 8 helices in a palmitoyl-oleolyl-phosphatidylcholine bilayer to evaluate and analyze channel models and to link the models to the experimentally measured conductance levels. Our results suggest that four helices do not form a stable water-filled channel and might not even form a stable intermediate. The lowest measurable conductance level is likely to correspond to the pentamer. At higher aggregation numbers the bundles become less symmetrical. Water properties inside the different-sized bundles are similar. The hexamer is the most stable model with a stability comparable with simulations based on crystal structures. The simulation was extended from 4 to 20 ns or several times the mean passage time of an ion. Essential dynamics analyses were used to test the hypothesis that correlated motions of the helical bundles account for high-frequency noise observed in open channel measurements. In a 20-ns simulation of a hexameric alamethicin bundle, the main motions are those of individual helices, not of the bundle as a whole. A detailed comparison of simulations using different methods to treat long-range electrostatic interactions (a twin range cutoff, Particle Mesh Ewald, and a twin range cutoff combined with a reaction field correction) shows that water orientation inside the alamethicin channels is sensitive to the algorithms used. In all cases, water ordering due to the protein structure is strong, although the exact profile changes somewhat. Adding an extra 4-nm layer of water only changes the water ordering slightly in the case of particle mesh Ewald, suggesting that periodicity artifacts for this system are not serious.

Detection of protein-ligand interaction on the membranes using C-terminus biotin-tagged alamethicin

C-terminal biotin-tagged alamethicin, which has several alpha-aminoisobutyric acid (Aib) residues in its sequence, was synthesized by the preparation of the protected peptide segment using the 2-chlorotrityl resin, followed by conjugation with biotin hydrazide. Suppression of the channel current of the biotin-tagged alamethicin by the addition of streptavidin to the electrolyte was monitorable in real time using the planar lipid-bilayer method. The system was also applicable to the detection of interaction of the biotin-tagged alamethicin with the anti-biotin antibody.

Modifications of alamethicin ion channels by substitution of Glu-7 for Gln-7

To evaluate the role of charged residues facing a pore lumen in stability of channel structure and ion permeation, we studied electrical properties of ion channels formed by synthesized native alamethicins (Rf50 (alm-Q7Q18) and Rf30 (alm-Q7E18)) and their analogs with Glu-7 (alm-E7Q18 and alm-E7E18). The single-channel currents were measured over a pH range of 3.5 to 8.7 using planar bilayers of diphytanoyl PC. The peptides all showed multi-level current fluctuations in this pH range. At pH 3.5 the channels formed by the four peptides were similar to each other irrespective of the side chain differences at positions 7 and 18. The ionization of Glu-7 (E7) and Glu-18 (E18) above neutral pH reduced the relative probabilities of low-conductance states (levels 1 and 2) and increased those of high-conductance states (levels 4-6). The channel conductance of the peptides with E7 and/or E18, which was distinct from that of alm-Q7Q18, showed a marked pH-dependence, especially for low-conductance states. The ionization of E7 further reduced the stability of channel structure, altered the current-voltage curve from a superlinear relation to a sublinear one, and enhanced cation selectivity. These results indicate that ionized E7 strongly influences the channel structure and the ion permeation, in contrast to ionized E18.

Alamethicin-leucine zipper hybrid peptide: a prototype for the design of artificial receptors and ion channels

In this report, we describe a novel concept of extramembrane control of channel peptide assembly and the eventual channel current modulation. Alamethicin is a peptide antibiotic, which usually forms ion channels in various association states. By introducing an extramembrane leucine zipper segment (Alm-LeuZ), the association number of alamethicin was effectively controlled to produce a single predominant channel open state. The assembly was estimated to be a tetramer, by comparison of the channel conductance with that of the template-assembled Alm-LeuZ tetramer, which was prepared by the conjugation of a maleimide-functionalized peptide template with cysteine-derivatized Alm-LeuZ segments. Employment of an extramembrane segment of a random conformation provided higher levels of channel conductance. The result exemplified the possibility of channel current control by a conformational switch of the extramembrane segments.

Conformational and orientation studies of artificial ion channels incorporated into lipid bilayers

The conformational and orientation studies in lipid bilayers of 21 amino acid peptides bearing six crown ethers are reported. The compounds were designed to form artificial ion channels by stacking the crown rings, and were shown to be functional in bilayer membranes. We used Fourier transform infrared spectroscopy and CD spectropolarimetry to study the conformation of the peptides in solution and in lipid bilayers. These studies revealed that hexacrown peptides retain their alpha-helical conformation when incorporated in a lipid bilayer environment. Attenuated total reflectance spectroscopy was used to investigate the orientation of the peptides in a lipid bilayer. Results demonstrated that the peptides are not oriented at a fixed angle in membrane, but rather are in incorporation equilibrium between an active state parallel to the lipid chain and an inactive state adsorbed at the surface of the bilayer. From these results, we propose a model for the channel activity and the gating mechanism of these hexacrown peptides in bilayer membranes.

Framework model for single proton conduction through gramicidin

This paper describes a framework model for proton conduction through gramicidin; a model designed to incorporate information from molecular dynamics and use this to predict conductance properties. The state diagram describes both motion of an excess proton within the pore as well as the reorientation of waters within the pore in the absence of an excess proton. The model is constructed as the diffusion limit of a random walk, allowing control over the boundary behavior of trajectories. Simple assumptions about the boundary behavior are made, which allow an analytical solution for the proton current and conductance. This is compared with corresponding expressions from statistical mechanics. The random walk construction allows diffusing trajectories underlying the model to be simulated in a simple way. Details of the numerical algorithm are described.

Simulations of ion channels--watching ions and water move

Ion channels mediate electrical excitability in neurons and muscle. Three-dimensional structures for model peptide channels and for a potassium (K+) channel have been combined with computer simulations to permit rigorous exploration of structure-function relations of channels. Water molecules and ions within transbilayer pores tend to diffuse more slowly than in bulk solutions. In the narrow selectivity filter of the bacterial K+ channel (i.e. the region of the channel that discriminates between different species of ions) a column of water molecules and K+ ions moves in a concerted fashion. By combining atomistic simulations (in which all atoms of the channel molecule, water and ions are treated explicitly) with continuum methods (in which the description of the channel system is considerably simplified) it is possible to simulate some of the physiological properties of channels.

Three-dimensional Poisson-Nernst-Planck theory studies: influence of membrane electrostatics on gramicidin A channel conductance

A recently introduced real-space lattice methodology for solving the three-dimensional Poisson-Nernst-Planck equations is used to compute current-voltage relations for ion permeation through the gramicidin A ion channel embedded in membranes characterized by surface dipoles and/or surface charge. Comparisons to a variety of experimental results, presented herein, have proven largely successful. Strengths and weaknesses of the method are discussed.

Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics

Continuum theories of electrolytes are widely used to describe physical processes in various biological systems. Although these are well-established theories in macroscopic situations, it is not clear from the outset that they should work in small systems whose dimensions are comparable to or smaller than the Debye length. Here, we test the validity of the mean-field approximation in Poisson-Boltzmann theory by comparing its predictions with those of Brownian dynamics simulations. For this purpose we use spherical and cylindrical boundaries and a catenary shape similar to that of the acetylcholine receptor channel. The interior region filled with electrolyte is assumed to have a high dielectric constant, and the exterior region representing protein a low one. Comparisons of the force on a test ion obtained with the two methods show that the shielding effect due to counterions is overestimated in Poisson-Boltzmann theory when the ion is within a Debye length of the boundary. As the ion gets closer to the boundary, the discrepancy in force grows rapidly. The implication for membrane channels, whose radii are typically smaller than the Debye length, is that Poisson-Boltzmann theory cannot be used to obtain reliable estimates of the electrostatic potential energy and force on an ion in the channel environment.