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.

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.

Sensing and generation of lipid membrane curvature, mediated by the binding of specific proteins onto the membrane surface, play crucial roles in cell biology. A number of mechanisms have been proposed, but the molecular understanding of these processes is incomplete. All-atom molecular dynamics simulations have offered valuable insights but are extremely demanding computationally. Implicit membrane simulations could provide a viable alternative, but current models apply only to planar membranes. In this work, the implicit membrane model 1 is extended to spherical and tubular membranes. The geometric change from planar to curved shapes is straightforward but insufficient for capturing the full curvature effect, which includes changes in lipid packing. Here, these packing effects are taken into account via the lateral pressure profile. The extended implicit membrane model 1 is tested on the wild-types and mutants of the antimicrobial peptide magainin, the ALPS motif of arfgap1, α-synuclein, and an ENTH domain. In these systems, the model is in qualitative agreement with experiments. We confirm that favorable electrostatic interactions tend to weaken curvature sensitivity in the presence of strong hydrophobic interactions but may actually have a positive effect when those are weak. We also find that binding to vesicles is more favorable than binding to tubes of the same diameter and that the long helix of α-synuclein tends to orient along the axis of tubes, whereas shorter helices tend to orient perpendicular to it. Adoption of a specific orientation could provide a mechanism for coupling protein oligomerization to tubule formation.

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.

No abstract