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.

Melittin is a short cationic peptide that exerts cytolytic effects on bacterial and eukaryotic cells. Experiments suggest that in zwitterionic membranes, melittin forms transmembrane toroidal pores supported by four to eight peptides. A recently constructed melittin variant with a reduced cationic charge, MelP5, is active at 10-fold lower concentrations. In previous work, we performed molecular dynamics simulations on the microsecond timescale to examine the supramolecular pore structure of a melittin tetramer in zwitterionic and partially anionic membranes. We now extend that study to include the effects of peptide charge, initial orientation, and number of monomers on the pore formation and stabilization processes. Our results show that parallel transmembrane orientations of melittin and MelP5 are more consistent with experimental data. Whereas a MelP5 parallel hexamer forms a large stable pore during the 5-μs simulation time, a melittin hexamer and an octamer are not fully stable, with several monomers dissociating during the simulation time. Interaction-energy analysis shows that this difference in behavior between melittin and MelP5 is not due to stronger electrostatic repulsion between neighboring melittin peptides but to peptide-lipid interactions that disfavor the isolated MelP5 transmembrane monomer. The ability of melittin monomers to diffuse freely in the 1,2-dimyristoyl-SN-glycero-3-phosphocholine membrane leads to dynamic pores with varying molecularity.

Protegrin-1 is an 18-residue β-hairpin antimicrobial peptide (AMP) that has been suggested to form transmembrane β-barrels in biological membranes. However, alternative structures have also been proposed. Here, we performed multimicrosecond, all-atom molecular dynamics simulations of various protegrin-1 oligomers on the membrane surface and in transmembrane topologies. The membrane surface simulations indicated that protegrin dimers are stable, while trimers and tetramers break down. Tetrameric arcs remained stably inserted in lipid membranes, but the pore water was displaced by lipid molecules. Unsheared protegrin β-barrels opened into β-sheets that surrounded stable aqueous pores, whereas tilted barrels with sheared hydrogen bonding patterns were stable in most topologies. A third type of observed pore consisted of multiple small oligomers surrounding a small, partially lipidic pore. We also considered the β-hairpin AMP tachyplesin, which showed less tendency to oligomerize than protegrin: the octameric bundle resulted in small pores surrounded by 6 peptides as monomers and dimers, with some peptides returning to the membrane surface. The results imply that multiple configurations of protegrin oligomers may produce aqueous pores and illustrate the relationship between topology and putative steps in protegrin-1’s pore formation. However, the long-term stability of these structures needs to be assessed further.

For a long time the structural and molecular features of mammalian histidine decarboxylase (EC 4.1.1.22), the enzyme that produces histamine, have evaded characterization. We overcome the experimental problems for the study of this enzyme by using a computer-based modelling and simulation approach, and have now the conditions to use histidine decarboxylase as a target in histamine pharmacology. In this review, we present the recent (last 5 years) advances in the structure-function relationship of histidine decarboxylase and the strategy for the discovery of new drugs.