Computational prediction of the optimal oligomeric state for membrane-inserted β-barrels of protegrin-1 and related mutants

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Putative Pore Structures of Amyloid β 25-35 in Lipid Bilayers

The amyloid β peptide aggregates to form extracellular plaques in the brains of Alzheimer's disease patients. Certain of its fragments have been found to have similar properties to those of the full-length peptide. The best-studied of these is 25-35, which aggregates into fibrils, is toxic to neurons, and forms ion channels in synthetic lipid bilayers. Here, we investigate possible pore-forming structures of oligomers of this peptide in a POPC/POPG membrane. We consider octameric and decameric β-barrels of different topology, strand orientation, and shear, evaluate their stability in an implicit membrane model, and subject the best models to multimicrosecond all-atom molecular dynamics simulations. We find two decameric structures that are kinetically stable in membranes on this time scale: an imperfectly closed antiparallel β-barrel with K28 in the pore lumen and a short parallel β-barrel with K28 toward the membrane interface. Both structures exhibit dehydrated gaps in the pore lumen, which are larger for the antiparallel barrel. Based on these results, the experimental cation selectivity, the dependence of ion channel activity on voltage direction, and certain mutation data, the parallel model seems more compatible with experimental data.

Computational Methods and Tools in Antimicrobial Peptide Research

The evolution of antibiotic-resistant bacteria is an ongoing and troubling development that has increased the number of diseases and infections that risk going untreated. There is an urgent need to develop alternative strategies and treatments to address this issue. One class of molecules that is attracting significant interest is that of antimicrobial peptides (AMPs). Their design and development has been aided considerably by the applications of molecular models, and we review these here. These methods include the use of tools to explore the relationships between their structures, dynamics, and functions and the increasing application of machine learning and molecular dynamics simulations. This review compiles resources such as AMP databases, AMP-related web servers, and commonly used techniques, together aimed at aiding researchers in the area toward complementing experimental studies with computational approaches.

Experimental and Computational Characterization of Oxidized and Reduced Protegrin Pores in Lipid Bilayers

Protegrin-1 (PG-1), an 18-residue β-hairpin stabilized by two disulfide bonds, is a member of a family of powerful antimicrobial peptides which are believed to act through membrane permeabilization. Here we used a combination of experimental and computational approaches to characterize possible structural arrangements of PG-1 in lipid bilayers mimicking bacterial membranes. We have measured the dose-response function of the PG-1-induced leakage of markers of various sizes from vesicles and found it to be consistent with the formation of pores of two different sizes. The first one allows the release of small dyes and occurs at peptide:lipid ratios < 0.006. Above this ratio, larger pores are observed through which the smallest of dextrans FD4 can be released. In parallel with pore formation, we observe a general large-scale destabilization of vesicles which is probably related to complete rupture of some vesicles. The population of vesicles that are completely ruptured depends linearly on PG-1:lipid ratio. Neither pore size, nor vesicle rupture are influenced by the formation of disulfide bonds. Previous computational work on oxidized protegrin is complemented here by all-atom MD simulations of PG-1 with reduced disulfide bonds both in solution (monomer) and in a bilayer (dimer and octamer). The simulations provide molecular insights into the influence of disulfide bonds on peptide conformation, aggregation, and oligomeric structure.

The unusual conformation of cross-strand disulfide bonds is critical to the stability of β-hairpin peptides

The cross-strand disulfides (CSDs) found in β-hairpin antimicrobial peptides (β-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of β-AMPs. To address this, we solved the nuclear magnetic resonance (NMR) structure of the β-AMP gomesin (Gm) at 278, 298, and 310 K, examined the disulfide bond geometry of over 800 disulfide-rich peptides, and carried out extensive molecular dynamics (MD) simulation of the peptides Gm and protegrin. The NMR data suggests Cα-Cα distances characteristic for CSDs are independent of temperature. Analysis of disulfide-rich peptides from the Protein Data Bank revealed that right-handed and left-handed rotamers are equally likely in CSDs. The previously reported preference for right-handed rotamers was likely biased by restricting the analysis to peptides and proteins solved using X-ray crystallography. Furthermore, data from MD simulations showed that the short Cα-Cα distance is critical for the stability of these peptides. The unique disulfide geometry of CSDs poses a challenge to biomolecular force fields and to retain the stability of β-hairpin fold over long simulation times, restraints on the torsion angles might be required.

Insights into Membrane Translocation of Protegrin Antimicrobial Peptides by Multistep Molecular Dynamics Simulations

Protegrin-1 (PG-1) is a cationic arginine-rich antimicrobial peptide. It is widely accepted that PG-1 induces membrane disruption by forming pores that lead to cell death. However, the insertion mechanism for these highly cationic peptides into the hydrophobic membrane environment is still poorly understood at the molecular scale. It has previously been determined that the association of arginine guanidinium and lipid phosphate groups results in strong bidentate bonds that stabilize peptide-lipid complexes. It has also been suggested that arginine residues are able to drag phosphate groups as they insert inside the membrane to form a toroidal pore. However, whether bidentate bonds play a significant role in inducing a pore formation remains unclear. To investigate the role of bidentate complexes in PG-1 translocation, we conducted molecular dynamics simulations. Two computational electroporation methods were implemented to examine the translocation process. We found that PG-1 could insert into the membrane, provided the external electric potential is large enough to first induce a water column or a pore within the lipid bilayer membrane. We also found that the highly charged PG-1 is capable in itself of inducing molecular electroporation. Substitution of arginines with charge-equivalent lysines showed a markedly reduced tendency for insertion. This indicates that the guanidinium group likely facilitates PG-1 translocation. Potential of mean force calculations suggests that peptide insertion inside the hydrophobic environment of the membrane core is not favored. We found that formation of a water column or a pore might be a prerequisite for PG-1 translocation. We also found that PG-1 can stabilize the pore after insertion. We suggest that PG-1 could be a pore inducer and stabilizer. This work sheds some light on PG-1 translocation mechanisms at the molecular level. Methods presented in this study may be extended to other arginine-rich antimicrobial and cell-penetrating peptides.

Transmembrane Pore Structures of β-Hairpin Antimicrobial Peptides by All-Atom Simulations

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, whereas 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 six 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.