Pore Structure and Synergy in Antimicrobial Peptides of the Magainin Family.

Magainin 2 and PGLa are among the best-studied cationic antimicrobial peptides. They bind preferentially to negatively charged membranes and apparently cause their disruption by the formation of transmembrane pores, whose detailed structure is still unclear. Here we report the results of 5-9 µs all-atom molecular dynamics simulations starting from tetrameric transmembrane helical bundles of these two peptides, as well as their stoichiometric mixture, and the analog MG-H2 in DMPC or 3:1 DMPC/DMPG membranes. The simulations produce pore structures that appear converged, although some effect of the starting peptide arrangement (parallel vs. antiparallel) is still observed on this timescale. The peptides remain mostly helical and adopt tilted orientations. The calculated tilt angles for PGLa are in excellent agreement with recent solid state NMR experiments. The antiparallel dimer structure in the magainin 2 simulations resembles previously determined NMR and crystal structures. More transmembrane orientations and a larger and more ordered pore are seen in the 1:1 heterotetramer with an antiparallel helix arrangement. Insights into the mechanism of synergy between these two peptides are obtained via implicit solvent modeling of homo- and heterodimers and analysis of interactions in the atomistic simulations. This analysis suggests stronger pairwise interactions in the heterodimer than in the two homodimers.

The structure of a melittin-stabilized pore.

Melittin has been reported to form toroidal pores under certain conditions, but the atomic-resolution structure of these pores is unknown. A 9-µs all-atom molecular-dynamics simulation starting from a closely packed transmembrane melittin tetramer in DMPC shows formation of a toroidal pore after 1 µs. The pore remains stable with a roughly constant radius for the rest of the simulation. Surprisingly, one or two melittin monomers frequently transition between transmembrane and surface states. All four peptides are largely helical. A simulation in a DMPC/DMPG membrane did not lead to a stable pore, consistent with the experimentally observed lower activity of melittin on anionic membranes. The picture that emerges from this work is rather close to the classical toroidal pore, but more dynamic with respect to the configuration of the peptides.

Implicit membrane treatment of buried charged groups: application to peptide translocation across lipid bilayers.

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

One-dimensional confinement of electric field and humidity dependent DNA conductivity.

The dependence of DNA assemblies conductance on relative humidity is investigated theoretically. Following earlier suggestions, we consider the ionic conductivity through the layers of water adsorbed by DNA molecules. The increase in humidity results in a growing water layer. The binding energy of ions depends on the thickness of the water layer due to change in water polarization. This dependence is very strong at smaller thicknesses of water layers due to the low-dimensional confinement of an electric field in water. We show that the associated change in ion concentration can explain the six orders of magnitude increase in conductivity, with relative humidity growing from 0.05 to 0.95.