Understanding pH-dependent selectivity of alamethicin K18 channels by computer simulation

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8Å pore and the octamer in an 11Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted "barrel-stave" model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.

Lysine methylation strategies for characterizing protein conformations by NMR

In the presence of formaldehyde and a mild reducing agent, reductive methylation is known to achieve near complete dimethylation of protein amino groups under non-denaturing conditions, in aqueous media. Amino methylation of proteins is employed in mass spectrometric, crystallographic, and NMR studies. Where biosynthetic labeling is prohibitive, amino (13)C-methylation provides an attractive option for monitoring folding, kinetics, protein-protein and protein-DNA interactions by NMR. Here, we demonstrate two improvements over traditional (13)C-reductive methylation schemes: (1) By judicious choice of stoichiometry and pH, ε-aminos can be preferentially monomethylated. Monomethyl tags are less perturbing and generally exhibit improved resolution over dimethyllysines, and (2) By use of deuterated reducing agents and (13)C-formaldehyde, amino groups can be labeled with (13)CH(2)D tags. Use of deutero-(13)C-formaldehyde affords either (13)CHD(2), or (13)CD(3) probes depending on choice of reducing agent. Making use of (13)C-(2)H scalar couplings, we demonstrate a filtering scheme that eliminates natural abundance (13)C signal.

Ligand-induced extramembrane conformation switch controlling alamethicin assembly and the channel current

In this review, we describe our approach to creating artificial receptor-channel proteins or sensor systems, using an extramembrane segment conformationally switchable by external stimuli. Alamethicin is known to self-assemble in membranes to form ion channels with various open states. Employment of an alpha-helical leucine-zipper segment resulted in the effective modulation of the association states of alamethicin to produce a single predominant channel-open state. A decrease in the helical content of the extramembrane segments was found to induce a channel-current increase. Therefore, conformational changes in the extramembrane segments induced by the interaction with ligands can be reflected in the current levels.

Channel-Forming Activity of Alamethicin: Effects of Covalent Tethering

In this short review article, the effects of covalent tethering of alamethicin molecules on channel-forming behavior are described. Broadly speaking, these chemical modifications have provided insight into all three aspects of channel behavior: the structure of the conducting state, the ion-selectivity and ion-permeation properties, and the voltage dependence. Each of these aspects are discussed in turn.

Implicit solvent simulations of peptide interactions with anionic lipid membranes

A recently developed implicit membrane model (IMM1) is supplemented with a Gouy-Chapman term describing counterion-screened electrostatic interactions of a solute with negatively charged membrane lipids. The new model is tested on peptides that bind to anionic membranes. Pentalysine binds just outside the plane of negative charge, whereas Lys-Phe peptides insert their aromatic rings into the hydrophobic core. Melittin and magainin 2 bind more strongly to anionic than to neutral membranes and in both cases insert their hydrophobic residues into the hydrocarbon core. The third domain of Antennapedia homeodomain (penetratin) binds as an alpha-helix in the headgroup region. Cardiotoxin II binds strongly to anionic membranes but marginally to neutral ones. In all cases, the location and configuration of the peptides are consistent with experimental data, and the effective energy changes upon binding compare favorably with experimental binding free energies. The model opens the way to exploring the effect of membrane charge on the location, conformation, and dynamics of a large variety of biologically active peptides on membranes.

Voltage-dependent energetics of alamethicin monomers in the membrane

The implicit membrane model IMM1 is extended to include the effect of transmembrane potential and used to investigate the optimal membrane binding configurations and energies for alamethicin helices. In the absence of voltage, the lowest energy configuration is on the membrane surface with a tilt allowing the N terminus to be fully buried. Slightly higher in energy is an also tilted configuration with the N terminus deeper in the membrane and almost crossing the membrane. In 26A membranes and in the presence of 0.1V voltage, the TM orientation becomes lower in energy. This is consistent with the assumption that voltage induces a transition from the interfacial to the inserted (TM) orientation. This effect of voltage is smaller in thicker membranes. The results are compared to previous experimental and theoretical studies and the findings are discussed in relation to the mechanism of channel formation by alamethicin.

Total synthesis in solution of alamethicin F50/5 by an easily tunable segment condensation approach

A total synthesis in solution of the 19-mer peptide component F50/5 of alamethicin, the most extensively investigated among the channel-former peptaibol antibiotics, is reported. Three peptide segments (A, B, C) were prepared and assembled, followed by incorporation of the acetylated N-terminal amino acid. The synthetic modules B and C are characterized by three Glu(OMe) residues (at positions 7, 18, and 19) that, after completion of the synthesis, were reacted with ammonia to provide alamethicin F50/5. By use of this general strategy, we also prepared the [Gln7, Glu(OMe)18,19] alamethicin F50/5 analogue. The purity and conformation of the final products were assessed by chromatographic, spectrometric, and spectroscopic techniques. This tunable segment condensation approach will pave the way for an easy synthesis of alamethicin analogues bearing amino acid residues with desired side-chain probes even at the N-terminus and in internal positions of the sequence.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Self-assembly and properties of diblock copolymers by coarse-grain molecular dynamics

Block-copolymer amphiphiles have been observed to assemble into vesicles and other morphologies long known for lipids but with remarkably different properties. Coarse-grain molecular dynamics (CG-MD) is used herein to elaborate the structures and properties of diblock copolymer assemblies in water. By varying the hydrophilic/hydrophobic ratio of the copolymer in line with experiment, bilayer, cylindrical and spherical micelle morphologies spontaneously assemble. Varying the molecular weight (MW) with hydrophilic/hydrophobic ratio appropriate to a bilayer yields a hydrophobic core thickness that scales for large MW as a random coil polymer, in agreement with experiment. The extent of hydrophobic-segment overlap in the core increases nonlinearly with MW, indicative of chain entanglements and consistent with the dramatic decrease reported for lateral mobility in polymer vesicles. Calculated trends with MW as well as hydrophilic/hydrophobic ratio thus agree with experiment, demonstrating that CG-MD simulations provide a rational design tool for diblock copolymer assemblies.

Electrostatic potentials and fields in the vicinity of engineered nanostructures

We have developed a method for calculating the electrostatic potentials and fields in the vicinity of geometrically complex engineered nanostructures composed of varying materials in electrolytes of arbitrary pH and ionic strength. The method involves direct summation of charged Debye-Hückel spheres composing the nanostructural surfaces and, by including charge redistribution on the surface of conducting materials held at constant potential, is applicable to mixed boundary conditions. The method is validated by comparison to analytical solutions for an infinite plane (Gouy-Chapman), an infinite cylinder (Bessel functions), and an infinite plane which contains a hole and which is held at constant potential. Excellent agreement between the potentials obtained by our numerical method and the closed form solutions is found for these conditions. The method is applied to the calculation of the electric field enhancement in the vicinity of a nanomembrane whose pore wall is held at constant charge and whose membrane surfaces are held at constant potential. The electric field is found to be enhanced by the charge buildup in the rim of the hole of the nanomembrane; the buildup results from the potential being held constant in the conducting region. Ion concentrations are also calculated. Positive ion rejection is found to be enhanced by this charge buildup in the region of the rim when a constant positive potential is applied.

Engineering charge selectivity in model ion channels

Most ion channel proteins exhibit some degree of charge selectivity, that is, an ability to conduct ions of one charge more efficiently than ions of the opposite charge. The structural origins of charge selectivity remain incompletely understood despite recent advances in the determination of cation-selective and anion-selective channel protein structures. Helix bundle channels formed via self-assembly of the peptide alamethicin provide a tractable model system for exploring the structural basis of charge selectivity. We synthesized covalently-linked alamethicin dimers, with amino acid substitutions at position 18 [lysine (Lys), arginine (Arg), glutamine (Gln), 2,3-diaminopropionic acid (Dpr)] in each helix, to assess the role of this position as a charge-selectivity determinant in alamethicin channels. Of the position 18 substitutions investigated, the Lys derivative exhibited the greatest degree of anion selectivity. Arg-containing channels were slightly less anion-selective than Lys. Interestingly, Dpr channels showed cation selectivity nearly equivalent to that exhibited by the neutral Gln derivative. We suggest that this result is due to a wider pore diameter that permits a greater number of counter-ions leading to enhanced charge screening and a lower effective side-chain positive charge.