A Molecular Dynamics Study of the Short-Helical-Cytolytic Peptide

Jul 22, 2017 - However, the development of novel peptide drugs is still limited by the elusive molecular understanding of the membrane-lysis mechanism...
2 downloads 15 Views 3MB Size
Article pubs.acs.org/JPCC

A Molecular Dynamics Study of the Short-Helical-Cytolytic Peptide Assembling and Bioactive on Membrane Interface Fude Sun,† Xiufang Ding,† Lida Xu,† Jun F. Liang,‡ Long Chen,*,† and Shi-Zhong Luo*,† †

Beijing Key Laboratory of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China ‡ Department of Biomedical Engineering, Chemistry, and Biological Sciences, Stevens Institute of Technology, Castle Point on Hudson, Hoboken, New Jersey 07030, United States S Supporting Information *

ABSTRACT: Cytolytic peptides (CPs) have long been employed as broad-spectrum antibiotic agents to overcome multidrug resistance. However, the development of novel peptide drugs is still limited by the elusive molecular understanding of the membrane-lysis mechanism and modeling of CPs, especially of the short helical species. In this study, a known anticancer CP named PTP-7b (FLGALFKALSHLL) in disrupting membranes via selfassembling approach was studied by combining experiments and time-extended coarse-grained dynamic simulations. Effective membrane disintegration was induced by aggregation of the membrane-bound peptide individuals, rather than the preassembled peptide clusters. The disturbance level of lipid bilayers depended on the peptide concentrations, which was responsible for the long time-costing of PTP-7b in killing cells. On the basis of lines of simulations and energy-landscape calculations, the dynamics of membrane deformation evolving toward preliminary leakage resulted from the aggregated PTP-7b was demonstrated, which was subjected to the spatiotemporal cooperation of the membrane-inserted and the periplasmic peptides. The molecular mechanism incorporated the 11th histidine interaction coupled with the peptide amphiphilicity in accelerating phospholipid migration outward. This study revealed elaborate modeling and dynamics information about the short helical CPs in membrane lysis, which would be helpful to understand the underlying mechanisms and rational design of CPs for drug application.



INTRODUCTION

To our knowledge, the models of membrane disruption affected by AMPs/CPs can be classified into the transmembrane-pore and nonpore models. The former includes the “barrel-stave” model11 in which peptides interact laterally to form ion-channel-like transmembrane pores, and the “toroidal” model12 that peptides lead to local curvature of cell membranes and the lipid bilayer bends continuously through the pore neighboring with peptides. Due to the length requirement in formation of transmembrane pore, the involved peptide candidates are mostly long α-helical peptides, such as alamethicin,13 magainins,14 and mellitin.15 In addition, a“double-belt” pore structure was proposed as a novel model for the long amphiphilic peptides.16 However, for the majority of short CPs, it is difficult to parallel bundle together to form fixed channel-like peptide polymers. Depending on residue composition, the short CPs show alternative penetrating through or lying on the phospholipid bilayer to prepare for the further activities. A typical penetrating peptide is substrate

The increasing incidence of drug resistance, both in microbes and cancer cells, has become a serious public health problem worldwide.1,2 Cytolytic peptides (CPs), which belong to a major class of antimicrobial peptides (AMPs), have long been discussed as potential drugs to overcome multidrug resistance.3,4 Known as the broad-spectrum resistance against various guest pathogens, some peptide analogues show medical selectivity in anticancer treatment,5,6 providing an alternative approach with mild and sustainable advantages. Among the AMPs/CPs in APD2 database, the α-helical structure acts as the most general class occupying ∼14% of the total.7 In most cases, the helicity property strengthens when they approach to and attach on cell membranes, which plays as the crucial phase in measuring their antibiotic efficiency.8 On account of the diversity of CPs in length, sequence, and structure, a number of membrane disrupting models have been proposed to illustrate their action principles.9,10 However, there are still many limitations to obtaining a consensus in deciphering the modeling and mechanism of the short helical CPs in membrane lysis. © 2017 American Chemical Society

Received: May 7, 2017 Revised: July 19, 2017 Published: July 22, 2017 17263

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

study will deepen our comprehension into the membrane-lysis approach of short linear aggregated peptides and the future design.

of trans-activating transcriptor (TAT) which consists of 9−16 basic residues.17 A recent study of penetrating peptide R9 (nine arginines) shows it adopts alternative pore-forming or poreabsent model affected by lipid compositions.18 For the rest short CPs with helical structures, the “carpet” model is the most accepted style, in which the peptides parallel distribute on membranes to transport the cytoplasmic components via specific membrane perturbation.19 The “detergents” model is a “carpet” dependent phase when a peptide concentration threshold reaches, then the peptides form patches to damage bilayers.20,21 Self-assembling in solution and membrane has been observed for many CPs.22,23 Peptide aggregation help to stabilize and sustainable release of peptide drugs.6,24 In our previous work, a CP named PTP-7b with 13 residues is found possessing bioactive significance in killing lung cancer cells relying on aggregation.25 However, the membrane lysis by PTP-7b shows an evident time-costing feature compared with the rapid potency of its parent peptide PTP-7 (K11H forms PTP-7b).26 It is speculated that PTP-7b exerts its anticancer bioactivity based on peptide aggregation upon membrane association. Considering that the typical prototype of PTP-7b is amphipathic, αhelical, cationic, and short in sequence, revealing the molecular details underlying the membrane deformation triggered by peptide aggregation appears quite meaningful for the following improvement. However, to date it remains a challenge to describe the subtle membrane disruption by conventional experimental means. Recent years, relying on the advantages of time-scale expansion and accuracy elevation, the coarse-grained (CG) dynamic molecular simulations, especially the CG model based on the Martini force filed,27−30 have been extensively employed in exploring the molecular basis in various biological processes. Martini force field in principle integrates four heavy atoms into one CG bead that is classified into 18 subtypes to optimize their interaction. Extensive calibration of the nonbonded interactions between building blocks is established to match with the experimental data,31 especially the thermodynamic property such as oil/water partitioning energy. Various processes like peptide−membrane binding,32 lipid organization,33 and protein−protein interaction34,35 are substantially dependent on the preference of the constituent partition between polar and nonpolar environments. Based on good practicability of phospholipid bilayers and updating advances for proteins, it has been successfully used to model different aspects of helical bioactive peptides when interacting with membranes,36,37 including investigation of the membrane-pore modeling.38−40 In addition, CG models of β-amyloid peptide41 and cyclic peptide32,42 have been employed to present their effects on membranes by forming assembly on the surface or acyl core. Energy landscapes within the peptide-bilayer interplay generated from Martini CG model is used to establish the correlation between structural characteristics and membrane disruption mode.43,44 Herein, we investigated the membrane-lytic process of aggregated peptide PTP-7b by Martini CG model combining with experiment results to interpret the peptide−membrane interaction and the underlying driving forces. A stepwise “outstretch” cytolytic process was presented and analyzed, which incorporated the participations of the membrane-bound and the periplasmic peptides. The structural reorientation details and the crucial residues were explored to figure out the inherent principles. The detailed information revealed in the



METHODS Experimental Section. Atomic Force Microscopy. The nonassembled peptide PTP-7b (40 μM) was prepared freshly from the stock solution (10 mM). The self-assembled peptide was prepared by incubation at 37 °C for 24 h. In brief, 0.1 mL of PTP-7b solution was dropped on the top of a freshly cleaved mica sheet glued to a steel AFM sample plate. Samples were incubated under ambient conditions for 10 min, washed to remove unbound peptides, and then air-dried. The nonassembled and self-assembled PTP-7b samples were examined using an atomic force microscope (AFM) (Veeco Instruments Inc., USA). Measurements were performed using a silicon nitride probe mounted on a NanoScope V Controller (Veeco Instruments Inc., USA) and in ScanAsyst mode (scan size of 2 × 2 μm, scan rate of 1.0 Hz). MTT Assay. The cytotoxicity of the peptides was determined using the MTT assay as described previously. Briefly, cells in a complete medium were added into 96-well plates (5 × 103 cells/well) and cultured at 37 °C for 14−16 h. After the cells were washed, they were fed with serum-free medium containing various concentrations of nonassembled and self-assembled PTP-7b and incubated at 37 °C for 2 h. After that, 10 μL of MTT (5 mg/mL) were added into each well. Cell viability was determined after 4 h of incubation by dissolving crystallized MTT with 10% SDS solution containing 5% isopropanol and 0.1% HCl and measuring the absorbance at 570 nm. Peptide Stability Determination. The peptide stock (10 mM) was diluted by deionized water to the final peptide concentration of 100 μM. The peptide solution was incubated at 37 °C and was tested at time points of 0 and 24 h, respectively. The peptide stability was then analyzed by reversephase high-pressure liquid chromatography, and the intact peptides were identified by ESI-MS (Thermo LTQ Orbitrap XL equipped with an electrospray ionization source). Confocal Microscopy. Freshly trypsinized human cervical cancer Hela cell line were seeded in a collagen-coated eight-well glass chamber (2 × 104 cells/well) and cultured in DMEM medium containing 10% FBS at 37 °C (5% CO2) overnight. Before the assay, cell membranes were stained with Di-8ANEPPS (2.5 μM) for 10 min. After being washed three times with PBS to remove the excess Di-8-ANEPPS, cells were exposed to fluorescein isothiocyanate (FITC)-labeled PTP-7b. Peptide accumulation and self-assembly on cell surfaces were investigated using a Leica TCS SP8 confocal microscope (Leica Microsystems, Germany) by setting the excitation wavelength at 488 nm and emission wavelength at 505−530 nm (for peptides) and 560−600 nm (for cell membranes), respectively. CG Simulations. Model Setup. All the simulations were performed on the Gromacs-4.6.6 package.45 The Martini CG force field, in which four heavy atoms and associated hydrogens are considered as one bead,29−31 was employed to investigate the molecular basis of the membrane-peptide system. The input parameters were obtained from the Martini homepage (www. cgmartini.nl). A simulation box size 10 × 10 × 15 nm3 containing phosphatidylcholine (POPC)/phosphatidylserine (POPS) with a 3:1 ratio was constructed using the insane.py46,47 to imitate the cancer lipid bilayers.48,49 The small simulation box was used to avoid the unnecessary spontaneous curvature of the mixed bilayers happening in large 17264

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 1. (A) AFM images of fresh (nonassembled) and overnight-incubated (assembled) PTP-7b dissolved in pure water. (B) PTP-7b aggregation on cell membranes induced lysis of Hela cells. The cell membranes were labeled by di-8-ANEPPS (red) and the peptide was labeled by FITC (green). Peptide concentration was fixed at 40 μM. (C) The cytotoxicity of nonassembled and assembled PTP-7b measured on Hela cells. Three independent measurements were conducted. (D) Peptide-induced surface pressure changes of lipid monolayer with a composition of DPPCcholesterol-DPPS (50:10:2.5 in molar ratio) mimicking mammalian cell membranes. The initial surface pressure of the lipid monolayer was set at about 30 mN/m. The final concentration of loaded peptides was 10 μM.

Simulation Details. The systems were first energy minimized for 5000 steps by the steepest descent method. Then NPT equilibration for 50 ns was conducted with the bilayer lipids, peptides, solvent, and ions coupled separately in semi-isotropic style. The temperature adopted the weakcoupling scheme of Berendsen thermostat51 at 310 K with constant τT of 1 ps. The pressure was controlled with the Berendsen barostat, with a pressure reference of 1.0 bar, a compressibility of 4.5 × 10−5 bar−1, and a constant τP of 5 ps. Standard cutoff schemes were used for the nonbonded potentials: LJ interaction was shifted to zero in the range 0.9−1.2 nm, and the electrostatic interactions was in the range 0.0−1.2 nm. Periodic boundary condition was employed. The step time of the simulation was set as 20 fs. The neighbor list

scale. The atomistic structure of PTP-7b and its mutant were first constructed by Pymol software,50 and then the peptide CG model together with the topology restrains were generated by Martinize script from the atomistic structures. The secondary structure of the peptides studied here was restrained as α-helical that was confirmed in our previous study.6 Peptide individuals were randomly filled into an oblate box with a size of 10 × 10 × 2 nm3, and then the peptide groups were integrated into the upside of the membrane model with a distance of 2 nm between peptides down-edge and membrane surface. Alternatively, the distance was set with −1 nm to insert the peptide in membrane. The lipid bilayers were solvated with standard CG water.6 Sodium ions were added to keep the system neutral. 17265

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

94.1% after the overnight incubation, which excluded the interference of peptide degradation for estimations (Figure S1 and Table S1). On the basis of the results, the time-costing cell lysis caused by PTP-7b was mostly or initially attributed to the membrane-bound peptide individuals, which packed to aggregates subsequently. Molecular Characterization by CG Simulations. The membrane-lytic process occurs at a molecular level in long time-scale, which remains a challenge to be captured by the atomistic simulation and current experiment methods. In order to obtain the molecular perspective of the dynamic membrane deformation triggered by peptide aggregation, the Martini CG force field suitable for long-term biological process was employed in the study. The membrane model was composed of 436 POPC/POPS with a ratio of 3/1 to mimic the cancer cell membrane with anionic surface (all simulations have been listed in Table 1). Based on the above discussion, we first

for the pairwise nonbonded interactions was updated every 10 steps. The pair-list cutoff was set as 1.4 nm, which ensured that the nonbonded forces between nearly all atoms within the cutoff distance were calculated as the pair-list updated. Each simulation ran for an effective time of 6.0 μs. All the simulation representations were performed by VMD software.52 Potential of Mean Force (PMF). To compute the PMF between two packed peptides on the membrane surface, the same system setup in the simulation that contained the peptide couple, bilayer, and solvent was used. The selected structures were first energy minimized by the steepest descent method. The PMF was calculated using the umbrella sampling methods.53 For each system, more than 20 windows were selected to simulate according to an average distance interval of 0.2 nm between the two peptide helixes. The initial configurations for each window were generated by pulling the two peptides from their associated state to their window location using the umbrella potentials with a force constant of 1500 kJ·mol−1·nm−2 and pulling rate of 20 nm/ns for 250 ps. The fixed chain was constrained by force constant of 1000 kJ· mol−1·nm−2. Finally, we obtained the successive peptide-bilayer initial configurations with peptides remaining on the membrane surface showing successive separation distances. Each window was then equilibrated for 100 ns with the force constant of 1000 kJ·mol−1·nm−2, followed by a 1.0 μs production simulation. The WHAM method54 was used to unbias the umbrella potentials.

Table 1. Overview of Peptide Variants at Different Concentrations in Associating with the Membrane

peptide analogues



RESULTS Membrane Disintegration Caused by Peptide SelfAssembling. PTP-7b (FLGALFKALSHLL) formed evident peptide aggregates after overnight incubation compared with the fresh soluble peptide counterpart (Figure 1A). The results were consistent with the SEM views of our previous study, in which PTP-7b was proved showing comparable anticancer activity to its parent peptide PTP-7 after 60 min.25 Relying on confocal microscope, no peptide attachment on cell surface occurred (within the first half an hour) until 60 min when obvious peptide aggregates emerged on membrane localization, followed by a rapid cell lysis (Figure 1B). The biphasic cell lysis caused by PTP-7b was therefore interrelated with its selfassembling upon interacting with membranes. However, it was questioned that whether it was the preassembled PTP-7b or peptide aggregates formed by the membrane-inserted individuals that actually executed the subsequent cell lysis, since PTP7b self-assembling property was tested here in absence of the membrane environment. In order to evaluate the bioactivity of the membrane-inserted and preassembled PTP-7b, the cytotoxicity of the nonassembled (fresh) and preassembled (overnight incubation) peptide was respectively tested using MTT assays. The IC50 value of the preassembled PTP-7b was nearly 2-fold that of the nonassembled peptide (Figure 1C), suggesting the cytotoxicity of PTP-7b decreased severely when it preassembled to associate with cell surface. As a further inspection, the different cases of PTP-7b were estimated by the pressure shifts of a model lipid monolayer subjected by peptides (Figure 1D). The incubation of fresh PTP-7b induced a biphasic change of surface pressure: a rapid increase of 27.4 mN/m and a subsequent decrease of 2.7 mN/m, indicating the evident membrane binding of peptides followed by events of membrane-insertion. In contrast, the preassembled PTP-7b could not produce an equivalent membrane attachment (δ3 = 16.9 mN/m) either the membrane perturbation. PTP-7b could maintain intact in

P/La (peptide/ lipids)

initial positionb (with upper leaflet)

peptide assembling

membranec distortion level

× √ × √ √

− + − ++ +++

× × ×

− − −

PTP-7b (FLGALFKALSHLL)

15/168 20/168 20/168 40/168 40/168

PTP-7 (FLGALFKALSKLL)

15/168 20/168 20/168

inserted inserted upside inserted 20-inserted 20-upside inserted upside inserted

20/168

inserted





20/168

inserted

×



PTP-7b-H11A (FLGALFKALSALL) PTP-7b-H11Y (FLGALFKALSYLL) PTP-7b-H11 (+)d (FLGALFKALSHLL)

The lipid number in the “P/L” indicated the lipids in the upper leaflet. bThe initial position of “inserted” and “upside” means the distance between the lower edge of peptide group and the upper membranes surface of 2.0 nm and −1.0 nm, respectively. cThe “−” symbol mean the membrane was normal; the “+” mean membrane deformation caused by peptides, and multiply of it represented the membrane deformation degree. dProtonation of the 11th histidine residue. a

investigated the temporal-spatial relevance of PTP-7b in membrane disruption by peptide assembling inserted and outside the membrane, respectively. To this end, 20 separated peptides were placed upon the membrane with respective distances to the upper leaflet of 2 nm and −1 nm to build two simulation systems (the left images of Figure 2A,D). The peptide concentration in simulations did not exceed the saturated concentration calculated in our previous experiments6 and at a comparable magnitude order. After 6.0 μs simulation, the peptides in both systems were able to form aggregates. Peptides outside the membranes selfpacked rapidly before binding to the lipid bilayer (Figure S2). However, it showed an insignificant influence on the membrane structure (Figure 2A). In contrast, evident bilayer curvature was induced by aggregation of the membrane-inserted peptides (Figure 2B). Meanwhile, the lateral displacement of the lipid acyl chains of the upper leaflet was much closer to zero, 17266

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 2. Presentations of membrane deformations affected by association of PTP-7b in different situations. Peptides were, respectively, placed outside (A) and inserted into (B,D,E) the membrane. Only backbones of peptides were shown with blue chains. Tails of POPC and POPS were distinguished in green and purple, and the head groups were presented in red. For clarification, the water solvents were omitted. The upper images in A, B, D, and E represented the final states of the four lipid-bilayer systems at different concentrations. The lower images showed the mass density profiles55 of the upper lipid leaflet incorporating peptides within the last 1.0 μs, which were generated in the same view to the upper snapshots. (C) Disorder parameters (Sz) of the acyl carbon bonds of the lipids in the upper leaflet of the four systems.

reflecting the dramatically disorder of associated lipid and growing defects of the membrane structure (Figure S4). Investigation of Molecular Details. To reveal the molecular mechanism of PTP-7b in cytolysis, its parent peptide PTP-7 (FLGALFKALSKLL) with a weak self-assembling ability was subjected to simulations. Consistent with the mode in Figure 2D, simulations of PTP-7 at P/L of 15/168 and 20/168 were built to examine their dynamic behaviors. Similar with PTP-7b (Figure 3A), PTP-7 at the lower P/L ratio did not aggregate, and the situation maintained at the higher P/L ratio (Figure S3). The simulation results were consistent with the previous data that showed the distinct self-assembling capabilities between PTP-7 and PTP-7b.25 However, the peptide orientation of PTP-7 and PTP-7b at the lower concentration exhibited a common mode: the N-portion was concealed in the bilayer, retaining the C-portion outward (Figure 3B and S5). The specific peptide arrangement is interrelated with the fact that the N-portion of the peptides is relative hydrophobic (Figure 3C,D), which is favorable for minimizing the system energy. The likewise structural character

showing more distinct membrane distortion (Figure 2C) impacted by the membrane-inserted peptide block. Therefore, the simulation results matched well to the above experiments, indicating a reduced membrane-break ability of the preassembled PTP-7b aggregate. Plenty of evidence supports that the membrane disturbance by antimicrobial peptides requires a high peptide concentration.56,57 For this purpose, three systems with peptide concentration gradient of (P/L) 15/168, 20/168, and 40/168 were constructed. At the lower P/L of 15/168, the bilayer structure remained normal (Figure 2D), and self-assembling of PTP-7b was not observed (Figure S3). In comparison, notable aggregates of peptides formed when the P/L ratio elevated to 20/168. The upper leaflet showed an apparent curvature with the peptide aggregate protruding out to solvent phase (Figure 2B). As expected, increased lipid instability (Figure 2C) and a more profound bilayer curvature happened when 40 peptide chains self-assembled on the membranes, and even some lipids inclined to break away (Figure 2E). Correspondingly, destabilized by the increasing activity of aggregated peptides, the lipid tail gradually extended toward the membrane surface, 17267

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 3. Peptide arrangement details relative with membrane surface of PTP-7b. (A) Top view of the dispersive distribution of the peptide individuals and (B) lateral view of one random peptide inserting in the membrane. The peptide backbone was presented in blue. Phe, Lys and His were presented in red, green and olive, respectively. The solvent molecules were presented by iceblue beads. (C) Density distributions of the lipid head, tail, Phe and Lys&His residues along the membrane normal. (D) Respective distance involvements of the C-portion and N-portion to the upper membrane surface. The blue dotted line was used to connect the horizontal level of the upper leaflet surface in B and D.

asymmetrical contact distribution of the peptide hydrophobic interface with the lipid bilayer showed that only the 1-F, 5-L, 9L, and 12-L predominantly residing in the hydrophobic bilayer (Figure 5A). In contrast, the 6-F exhibited weaker affinity with the lipid bilayer. The abnormal residue alignment with bilayer not only favored for the peptide hydrophilic surface distributing in the polar solvent, it also allowed for the interaction occurrence between the histidine residues (Figure 5B). To quantitatively describe the association energy of PTP-7b, PMF of PTP-7b and the two mutants along with the separation distance between the paired peptides were estimated, respectively (Figure 5C). The association energy of PTP-7b showed two minimum regions with energy of −23.0 kJ/mol (0.9 nm) and −12.9 kJ/mol (∼1.4 nm). In comparison, the binding energy of H11A reached to only one lowest valley at the distance of ∼1.4 nm, and it showed a shallow energy valley of −6.8 kJ/mol. Despite that the association energy profile of H11Y presented a similar shape with PTP-7b, the typical well depth was −8.5 kJ/mol which was similar to that of H11A. The value was much lower than that of PTP-7b, implying the easier dissociation of the two histidine mutations. Therefore, the results confirmed the gift of the 11th histidine in stabilizing PTP-7b self-assembling on membrane interface. Moreover, the Coulomb potential evolvements of the different residues, ascalculated from the nonbonded potentials

is also found by previous experimental data, which suggests that PTP-7 adopts the “carpet” like model in killing cancer cells.26 In order to explain the different routines of PTP-7b and PTP-7 behavior following on membranes, the influence of the unique 11th histidine was taken into consideration. It was interesting that the histidine residues situated on the top of peptide aggregates, arranged linearly and exposed to the water phase (Figure 4A). By calculating the residue mean smallest distance between the neighboring peptides, the histidine residues presented remarkable on the interacting matrix (Figure S6). As a better comparison, the distance between the 11th histidine backbones on the neighboring peptide pair was calculated, and it was obviously lower than the distance between the neighboring peptides (Figure 4B). The result therefore revealed the importance of the histidine in peptide association. Furthermore, the 11th histidine was substituted by alanine (H11A) to remain the same charge amount and by tyrosine (H11Y) to keep a comparable size of the side chain. The peptide variants were subjected to simulation at the ratio of 20/168. As a result, the clustering of H11A exhibited low ability as weak as PTP-7 (Figure 4C). Despite that H11Y showed a propensity to clustering, the lipid disorder of the upper leaflet affected by the two mutants were less than that caused by PTP-7b (Figure 4D). PTP-7b is an amphipathic peptide. However, an 17268

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 4. (A) Arrangement of the histidine residues in the peptide self-assembly. The histidine residues were presented by yellow beads, and the presentation methods of the bilayer were consistent with Figure 2. Top and lateral views were shown, respectively. (B) The distance evolvements of the two interacting peptides and the two histidine residues. Only the backbone beads were considered in the calculations. (C) Number of peptide clusters (NC) evolvements of PTP-7, PTP-7b, H11A, and H11Y as a function of time. (D) Lipid disorder parameters of the upper leaflets subjected by PTP-7, PTP-7b, H11A, and H11Y.

Figure 5. (A) Contact intensity between the respective residues on PTP-7b and the upper bilayer lipids. The data was generated from the system with P/L ratio of 20/168. All the peptides were taken into account and then averaged. (B) Presentation of one interacting peptide pair on the membrane surface. The assembling-involved histidine residues were denoted by red dotted circle, and the bilayer was shown with a red square. (C) PMFs of PTP-7b-WT, H11A, and H11Y calculated by umbrella sampling calculations at a sampling generation of 1.0 μs per window (configure samplings were shown in Figure S7). The association energy was set to zero at the largest pulling distance.

the histidine residues. The results revealed that the high affinity between PTP-7b peptides did not merely depend on the π−π stacking of the histidine aromatic side chains.58 A previous work shows that the nitrogen atoms of the imidazole ring can act as a hydrogen bond acceptor and/or donor in the neutral state.59 It is speculated that interaction of the N−H···N type formed between histidine imidazole and the lysine side chain of the paired peptide. Despite that the explicit hydrogen bond is not available in Martini CG model, the interaction between

of PTP-7b and the histidine mutations, were similar and changed little (Figure 6A). The Lennard-Jones (L-J) potential of H11A also remained unchanged, corresponding to the dispersive peptide distribution on the membrane (Figure S3). H11Y showed an evident L-J potential decrease, which acted as the motivation for the obvious aggregation of H11Y (Figure 4C). In comparison, wild type (WT) of PTP-7b generated the most dramatic decline of the L-J potential (Figure 6B), reflecting the remarkable molecular interaction derived from 17269

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 6. (A) Coulomb potentials and (B) L-J potentials of PTP-7b, H11A, H11Y as a function of time. (C) NC evolvements of the three peptides on membranes. (D) Disorder level of acyl carbon bonds of upper lipids impacted by three different peptides.

Figure 7. (A) Membrane-outside and -inserted peptides were distinguished with yellow and blue without the side chains. Acyl tails of POPC and POPS were colored as green and purple, respectively. The images were captured at the final frame in lateral and top view. The initial state can be seen in Figure 9. (B) Peptide self-assembling evolvement (black line) and breaking away of the peptide-associated lipids from the membrane (blue line).

lysine side chain and histidine imidazole is at the “attractive”, ranking the second strongest level of the L-J interaction (cgmartini.nl). In fact, obvious crossed interaction between the 11-histidine and 7′-lysine was observed in the residue contact distribution (Figure S6). Therefore, the histidine enabled PTP7b to tightly bind together relying on the π−π stacking and additive hydrogen bonds with the lysine on the paired peptide. Because PTP-7b adopted α-helical structure, the lysine and histidine locate at the same surface of the peptide helical projection, which is favorable for the crossed interaction between them (Figure 5B).

The histidine residue can be protonated when the pH is below 5.5. The cell surface of cancer especially the carcinoma is more acidic than the normal cells. Our previous study shows that the human lung carcinoma cell A549 has a more acidic surface (pH = 7.15) than that of the normal human lung cell CCD-13Lu (pH = 7.4) under physiological condition.60 However, the pH level of A549 did not reach the pH threshold to protonate the histidine residue of PTP-7b.6 When the histidine residue was protonated, PTP-7b could not behave the biphasic cell lysis because of the loss of peptide self-assembling, resulting the PTP-7b in an identical pathway to that of PTP17270

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 8. (A) Contacts between peptide and liberated lipids (red line) and solvent-accessible surface area (SASA, green line) of the liberated lipids along with the time. (B) Distances of the membrane-inserted peptides from the bilayer center with and without the periplasmic peptides, respectively.

7.25 The histidine-protonated PTP-7b at the P/L ratio of 20/ 168 was conducted in simulation. In consistence with the experiments, the protonated PTP-7b could not aggregate on the membrane surface (Figure 6C), and the membrane structure remained normal (Figure 6D). Evidence shows that the acidic periplasmic space of tumor tissues was mainly generated from the glycolytic metabolism.61 Additionally, the pH of cancer cells changes in different periods, and the protons transport from the cytosol and anionic lipids can be obstructed by the membrane-anchored glycocalyx. Subsequently, they form an acid layer outside of the bilayer.62 The bilayer interface thus forms a negative net potential that is responsible for attracting the cationic peptides.48,63 Because PTP-7b destabilizes the bilayer primarily relying on the assembling of the membrane-inserted peptides, the protonation state is hard to maintain because of the negative net potential of the membrane interface. However, as the protons of anionic lipids such as the phosphatidylserines transfer toward the periplasmic space to form the acidic layer, PTP-7b is possibly protonated during its approaching to the lipid membrane while passing the acidic layer. The protonated PTP-7b shows a weak self-assembling in periplasmic space (Figure S8), which is favorable for its selective binding to cancer cells. Modeling Details of PTP-7b Aggregates To Induce Bilayer Leakage. Enlighted by the dynamics evidence of some lipids extending up in the system at the highest P/L ratio of 40/ 168 (Figure 2E), further cell lysis is believed to occur by the continued peptide−lipid interplay. Considering the successive adhesion of the periplasmic peptides toward the target cells, which resulted in a concentration gradient of peptides distributing in vicinity of membranes,64 the peptide-bilayer system with the P/L ratio of 40/168 was performed again; however, we divided the peptide pool into two classes, which were positioned outside and inserted into the model membrane, respectively. As a result, accompanying the clustering of peptides, some individual lipids were completely abstracted from the upper leaflet (Figure 7A,B-black line), which was considered as the directive phenomenon of cell leakage.25,65 By calculating the distance of the translocated lipids from the membrane, the separated lipids reached to a high level of 2.5 nm from the bilayer center (Figure 7B, blue line). Moreover, the peptides from the membrane outside encircled the liberated lipids by a neatly parallel arrangement, which was consistent with that on the membrane surface shown in Figure 4A. The arrangement pattern is structurally favorable

since it is able to tie the neighboring peptides in the most efficient way to block the lipids. In addition, it was favorable for the stability of the system by coordinating the peptide amphipathicity and histidine interaction. Meanwhile, segregated subdomain of POPS was observed to be enriched below the peptide assembly (Figure S9), which has been considered as a benchmark of membrane degradation.66,67 Note that the liberated lipids were mainly POPC. It was speculated that the leakage of the POPC lipids increased the positional possibility of POPS lipids for concentrative translocation. Moreover, the enrichment of POPS could be additionally guided by the electrostatic attraction derived from the growing positive PTP-7b aggregate. The continuous peptide aggregation produced the lipid− peptide mixed cluster that extended to the water phase. In order to define the structure character of the mixture, the compositions on the shell and core of the mixed cluster were analyzed. An obvious increase of the contact intensity between the liberated lipids and the peptides was observed during the peptide-bilayer interplay (Figure 8A, red line), which implied the enhanced integration of peptide−lipid mixture. Meanwhile, the solvent accessible surface area of the liberated lipids declined evidently during the dynamic evolution (Figure 8A, green line). The results indicated that though the liberated lipids moved away from the bilayer, they were kept from the solvent environment. In other words, during the out-leakage of the bilayer lipids, they were sealed inside the peptide shell that performed as the “cage-like” structure. Comparing the distance evolvements of the membrane-inserted peptides from bilayer center with and without the periplasmic peptides, the periplasmic peptides could assist the membrane-inserted peptides to extend increasingly out to the solvent environment (Figure 8B). Taken together, the results indicated that the bilayer destruction induced by aggregated peptides also relied on the additive contribution of the periplasmic peptides.



DISCUSSION Peptide self-assembling has been extensively studied for years in drug delivery68 and sustained release.25 Antibiotic assessment of natural and artificial peptides by forming ordered aggregates have been extensively explored recently.69,70 As the majority class of AMPs, CPs in assembling action attracts increasing attentions and the interrelated cytolytic mechanism is different from the existing pore-forming models. It is important to establish molecular knowledge of the modeling details, 17271

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

Figure 9. Evolvement of the bilayer lysis triggered by the spatiotemporal aggregation of PTP-7b. The periplasmic and the membrane-inserted peptides were respectively rendered with orange and blue to clarify their initial positions and the subsequent cooperation. The POPC and POPS lipids were represented with green and purple beads, respectively.

peptides initially bind to the anionic membrane as individuals relying on the electrostatic interaction. As continued accretion of peptides attach on the localized membrane, the membrane suffers unfavorable stress, and the membrane integrity starts to worsen. Furthermore, the membrane distortion becomes more dramatic by addition of the peripheral peptides. Relying on the traction of the volume-increased peptide aggregate toward the water phase, the associated lipids gradually migrate from the hydrophobic membrane and were concealed in peptides (Movie S1). The resulting mixture of lipids and peptides might act as the preliminary micelle that cracks away from the membrane, followed by the membrane disruption and cell lysis. The phenomenon is likewise with the report of another CP,80 which highlights the importance of peptide behaving at the water−lipid interface. During the process, the completed lipid migration caused by aggregated peptides incorporates the membrane-inserted as well as the periplasmic peptides, indicating the significance of the periplasmic peptides in membrane disruption. The same phenomenon is also observed in the formation of “pre-pore” by melittin, which requires a cooperative response of periplasmic melittin molecules to lower the free energy in creating small transmembrane-pore.81 Note that in the delivery and uptake of peptide drugs, the situations of soluble and assembled peptides usually coexist. Moreover, affected by environment factors like ion buffering and changeable local pH, the conversion between the soluble and assembled peptides therefore changes occasionally, especially on the drug targets of membrane surface. Based on the above results, we suggest some propositions underlying the short CP in associating with lipid bilayer. First, the cationic PTP-7b shows propensity to binding the anionic cancer cell membranes, which enables its anticancer selectivity. To this point, PTP-7 possesses a high cytolysis efficacy probably relying on the enhanced electrostatic interaction with cancer cells.5 However, PTP-7 loses the ability to form aggregates because of electrostatic repulsion. In contrast, a neutral peptide variant can form ordered nanofibers in pure water, as shown in a previous study.24 Note that aggregation of PTP-7b occurred at the membrane interface. On the one hand, it ensures stability for the peptide, depending on the hydrophobic residues of 9th-Phe, 1st-Leu, and 12th-Leu anchoring on membrane. In fact, the amphipathic distribution of polar and hydrophobic residues acts as a benchmark for membrane-active peptides in disturbing bilayer integrity, either by creation of lipid disruption or pore formation.4 On the other hand, the residue arrangement is responsible for the interaction between the 11th histidine residue, which stabilizes the peptide assembly in the membrane.82,83 Moreover, the histidine

structure relevance, and dynamics underlying the membrane distortion followed by ultimate leakage. In this study, the molecular details of membrane disruption induced by aggregated CP of PTP-7b were revealed combining experiments and extended CG dynamic simulations, which is hopeful to provide referential information in the development of new CP drugs. The membrane disruption process shares several hallmarks similar to previously proposed models. For example, when peptides bind on membrane surface, CPs adopt an orientation with the hydrophobic terminal inserted into the membrane and distribute as dispersive individuals at low concentration. This phase is generally considered as the preparation for the following membrane insertion or pore formation.71 PTP-7b in assembled and nonassembled cases do not present “pore” style that is suitable for Magainin-2 and Melittin studied in CG modeling.40 It is suspected that PTP-7b in the initial phase of binding membrane surface conforms to the description of the “carpet” model, which has been used to explain the mechanism of the cationic lytic peptides,72 such as PTP-7.73 The membrane coated by “carpet” peptides is often followed by membrane disruption alternatively realizing in a “detergents” manner, or forming transient pores to disintegrate membranes. One review describes that the “detergents” model derives from the “carpet” model and gradually generates micelle-like mixtures of lipids and peptides.74 Likewise, the active processing of aggregated PTP-7b then exceeds the critical concentration. However, a controversial report describes that the onset of “detergents” model starts from the attachment of peptide clusters concomitant with pore formation.4,75 In this regard, more efforts are required to depict the principles of the “detergents” style. In this study, the preassembled PTP-7b is proved showing lower efficiency to disintegrate membrane without formation of transmembrane pore, and PTP-7b mostly adopts another variant of “carpet” model to extend out the localized membrane by aggregation. In a previous study exploring the modeling of maculatin, a multitude of disputed evidence is reported that maculatin adopts a pore-like model76,77 or carpet-like model,78,79 depending on giant unilamellar vesicles, NMR spectroscopy, and so on. It indicates that the modeling presentation is affected by the methods and conditions. The CG simulations were hence executed and membrane insertion of maculatin aggregates was observed at a higher P/L ratio,39 which compensated the molecular details and explanation for experiments. To explore the cytolysis of PTP-7b, we built lines of simulations based on known experimental data to explore the substantial molecular details. As shown in Figure 9, 17272

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C

aggregation of the membrane-inserted peptides, rather than the preassembled peptides. (2) The self-assembling and membrane-lytic abilities of PTP-7b on membranes depend on the amphipathic property and peptide interaction based on the histidine residues. (3) The processing of lipid leakage triggered by peptide aggregation requires for a cooperation of the membrane-inserted and the periplasmic peptides. Taken together, the molecular perspective of this study is favorable for the understanding about action modeling of short helical CPs and development of novel clinical peptides.

residues can form hydrogen bonds with the opposite lysine on the paired peptide. It therefore compensates the availability for further fusion between the peptide assemblies on membranes and the outside counterpart. The peptide arrangement makes the positive lysine outward to water instead of interacting to the lipid phosphates. The phenomenon has been considered as one of “carpet” model principles,84 which supports that PTP-7b acts as a derivative “carpet” model. Note that the peptide aggregation ability and bioactivity was sensitive to histidine position, since the other variant sequence of FLGALFHALSKL (exchange the lysine and histidine of PTP-7b) shows degraded aggregation and membrane-lysis.25 In a previous CG simulation study,42 an aggregated cyclic peptide induces bilayer lipids to form micelle-like mixtures. Note that the preassembled cyclic peptides show higher membrane destruction than the membrane-inserted peptides, which is different from the linear helical PTP-7b in this study. Given the difference on structure and sequence, they would adopt a distinct mechanism during their bilayer leakage processes. For example, for the cyclic peptide, no evidence shows the cooperation between the membrane-inserted and periplasmic peptides at high peptide concentration. However, the linear helical PTP-7b needs extra assistance of the periplasmic peptides for lipid further migrating outward. In addition, the generated micelle-like mixture formed by the cyclic peptide is disordered except for the limited peptide interstacking. In contrast, PTP-7b forms a “cage-like” structure where the peptides display a parallel arrangement to better pack together. Therefore, PTP-7b assembles on the membrane surface and induces cell lysis, suggesting it adopts a special pathway. It attributes to the specific residue alignments in balancing the amphipathic property and interaction based on the 11th histidine residues. Moreover, it is important for the anticancer activity of PTP-7b since the balance between charges and hydrophobicity is important for selectivity of amphipathic helical peptides.85,86 The phenomenon-based models, such as “carpet” model and “membrane-pore” model, give limited understanding of short helical CPs to predict the activity and design novel peptide drugs.9 The modeling of action described in this study represents a new type of aggregated CP, which comprises the peptide distribution between solvent and hydrophobic membranes by specific peptide interaction mode. Although this study has considered the increasing demand of antibiotic agents with improving properties and practicability, as well as molecular details involved the modeling and driving forces of the modified CPs with enhancing stability, both efficiency and simplification still require continued efforts.Despite the successfulness of CG simulations in presenting AMPs poreforming model,39,40 the deficiency of the CG simulation by a reduction of atom resolution and screening limitation of the electrostatic interaction should be carefully regarded. A more accurate definition is promising for the future research of the bioactive aggregated peptide upon membrane association. In addition to anticancer activity and sustainable release of the short helical CPs, the unknown clinical applications like drug delivery and drug modification deserve to be further exploited.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.7b04347. Respective self-assembling of the membrane-inserted and the periplasmic peptides, as well as the connection between them (MPG) Further association of the peptide blocks (MPG) Lipids were pulled outward by coordination of the peptide assembly (MPG) The liberated lipids were caged in a peptide capsule and remained steady (MPG) Additional description of analysis and supplemental data (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. Tel: 86-10-64438015. *E-mail: [email protected]. ORCID

Fude Sun: 0000-0003-2914-6356 Shi-Zhong Luo: 0000-0002-4880-5962 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the National Basic Research Program of China (973 program) (2013CB910700), the National Natural Science Foundation of China (21672019, 21372026, 21402006, 61304147), the Fundamental Research Funds for the Central Universities (Grant Nos. YS1407, 2050205, 20130801, PT1613-08, and PYBZ1711), and Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase).



REFERENCES

(1) Arias, C. A.; Murray, B. E. Antibiotic-Resistant Bugs in the 21st Centurya Clinical Super-Challenge. N. Engl. J. Med. 2009, 360, 439−443. (2) Wicki, A.; Mandalà, M.; Massi, D.; Taverna, D.; Tang, H.; Hemmings, B. A.; Xue, G. Acquired Resistance to Clinical Cancer Therapy: A Twist in Physiological Signaling. Physiol. Rev. 2016, 96, 805−829. (3) Zasloff, M. Antimicrobial Peptides of Multicellular Organisms. Nature 2002, 415, 389−395. (4) Bechinger, B.; Lohner, K. Detergent-Like Actions of Linear Amphipathic Cationic Antimicrobial Peptides. Biochim. Biophys. Acta, Biomembr. 2006, 1758, 1529−1539. (5) Mader, J. S.; Hoskin, D. W. Cationic Antimicrobial Peptides as Novel Cytotoxic Agents for Cancer Treatment. Expert Opin. Invest. Drugs 2006, 15, 933−946.



CONCLUSIONS In summary, we employed CG molecular dynamic simulations combined with experiments to explore the molecular basis of the cytolysis mechanism of the self-assembling peptide PTP-7b. (1) Efficient membrane perturbation is induced by lateral 17273

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C (6) Chen, L.; Tu, Z.; Voloshchuk, N.; Liang, J. F. Lytic Peptides with Improved Stability and Selectivity Designed for Cancer Treatment. J. Pharm. Sci. 2012, 101, 1508−1517. (7) Wang, G.; Li, X.; Wang, Z. Apd2: The Updated Antimicrobial Peptide Database and Its Application in Peptide Design. Nucleic Acids Res. 2009, 37, D933−D937. (8) Yeaman, M. R.; Yount, N. Y. Mechanisms of Antimicrobial Peptide Action and Resistance. Pharmacol. Rev. 2003, 55, 27−55. (9) Wimley, W. C.; Hristova, K. Antimicrobial Peptides: Successes, Challenges and Unanswered Questions. J. Membr. Biol. 2011, 239, 27− 34. (10) He, K.; Ludtke, S. J.; Huang, H. W.; Worcester, D. L. Antimicrobial Peptide Pores in Membranes Detected by Neutron inPlane Scattering. Biochemistry 1995, 34, 15614−15618. (11) Rapaport, D.; Shai, Y. Interaction of Fluorescently Labeled Pardaxin and Its Analogues with Lipid Bilayers. J. Biol. Chem. 1991, 266, 23769−23775. (12) Ludtke, S. J.; He, K.; Heller, W. T.; Harroun, T. A.; Yang, L.; Huang, H. W. Membrane Pores Induced by Magainin. Biochemistry 1996, 35, 13723−13728. (13) Yan, H.; Hancock, R. E. Synergistic Interactions between Mammalian Antimicrobial Defense Peptides. Antimicrob. Agents Chemother. 2001, 45, 1558−1560. (14) Berkowitz, B. A.; Bevins, C. L.; Zasloff, M. A. Magainins: A New Family of Membrane-Active Host Defense Peptides. Biochem. Pharmacol. 1990, 39, 625−629. (15) Zanetti, M. The Role of Cathelicidins in the Innate Host Defenses of Mammals. Curr. Issues. Mol. Biol. 2005, 7, 179−196. (16) Vácha, R.; Frenkel, D. Simulations Suggest Possible Novel Membrane Pore Structure. Langmuir 2014, 30, 1304−1310. (17) Rao, K. S.; Labhasetwar, V. Trans-Activating Transcriptional Activator (Tat) Peptide-Mediated Brain Drug Delivery. J. Biomed. Nanotechnol. 2006, 2, 173−185. (18) Sharmin, S.; Islam, M. Z.; Karal, M. A. S.; Alam Shibly, S. U.; Dohra, H.; Yamazaki, M. Effects of Lipid Composition on the Entry of Cell-Penetrating Peptide Oligoarginine into Single Vesicles. Biochemistry 2016, 55, 4154−4165. (19) Shai, Y. Mode of Action of Membrane Active Antimicrobial Peptides. Biopolymers 2002, 66, 236−248. (20) Strömstedt, A. A.; Ringstad, L.; Schmidtchen, A.; Malmsten, M. Interaction between Amphiphilic Peptides and Phospholipid Membranes. Curr. Opin. Colloid Interface Sci. 2010, 15, 467−478. (21) Wimley, W. C. Describing the Mechanism of Antimicrobial Peptide Action with the Interfacial Activity Model. ACS Chem. Biol. 2010, 5, 905−917. (22) Liu, L.; Xu, K.; Wang, H.; Tan, P. J.; Fan, W.; Venkatraman, S. S.; Li, L.; Yang, Y.-Y. Self-Assembled Cationic Peptide Nanoparticles as an Efficient Antimicrobial Agent. Nat. Nanotechnol. 2009, 4, 457− 463. (23) Dehsorkhi, A.; Castelletto, V.; Hamley, I. W. Self-Assembling Amphiphilic Peptides. J. Pept. Sci. 2014, 20, 453−467. (24) Tu, Z.; Hao, J.; Kharidia, R.; Meng, X. G.; Liang, J. F. Improved Stability and Selectivity of Lytic Peptides through Self-Assembly. Biochem. Biophys. Res. Commun. 2007, 361, 712−717. (25) Chen, L.; Patrone, N.; Liang, J. F. Peptide Self-Assembly on Cell Membranes to Induce Cell Lysis. Biomacromolecules 2012, 13, 3327− 3333. (26) Tu, Z. Design and Mechanism Study of Anti-Cancer Lytic Peptides. Ph.D. Thesis, University of Hong Kong, 2009. (27) Marrink, S. J.; Risselada, H. J.; Yefimov, S.; Tieleman, D. P.; De Vries, A. H. The Martini Force Field: Coarse Grained Model for Biomolecular Simulations. J. Phys. Chem. B 2007, 111, 7812−7824. (28) Monticelli, L.; Kandasamy, S. K.; Periole, X.; Larson, R. G.; Tieleman, D. P.; Marrink, S.-J. The Martini Coarse-Grained Force Field: Extension to Proteins. J. Chem. Theory Comput. 2008, 4, 819− 834. (29) Marrink, S. J.; Tieleman, D. P. Perspective on the Martini Model. Chem. Soc. Rev. 2013, 42, 6801−6822.

(30) Marrink, S. J.; De Vries, A. H.; Mark, A. E. Coarse Grained Model for Semiquantitative Lipid Simulations. J. Phys. Chem. B 2004, 108, 750−760. (31) de Jong, D. H.; Singh, G.; Bennett, W. D.; Arnarez, C.; Wassenaar, T. A.; Schäfer, L. V.; Periole, X.; Tieleman, D. P.; Marrink, S. J. Improved Parameters for the Martini Coarse-Grained Protein Force Field. J. Chem. Theory Comput. 2013, 9, 687−697. (32) Khalfa, A.; Treptow, W.; Maigret, B.; Tarek, M. Self Assembly of Peptides near or within Membranes Using Coarse Grained Md Simulations. Chem. Phys. 2009, 358, 161−170. (33) Ingólfsson, H. I.; Melo, M. N.; van Eerden, F. J.; Arnarez, C. m.; Lopez, C. A.; Wassenaar, T. A.; Periole, X.; de Vries, A. H.; Tieleman, D. P.; Marrink, S. J. Lipid Organization of the Plasma Membrane. J. Am. Chem. Soc. 2014, 136, 14554−14559. (34) Sun, F.; Chen, L.; Wei, P.; Chai, M.; Ding, X.; Xu, L.; Luo, S.-Z. The Dimerization and Structural Stability of Amyloid Precursor Proteins Affected by the Membrane Microenvironments. J. Chem. Inf. Model. 2017, 57, 1375−1387. (35) Wei, P.; Sun, F.; Zuo, L.-m.; Qu, J.; Chen, P.; Xu, L.-d.; Luo, S.Z. Critical Residues and Motifs for Homodimerization of the First Transmembrane Domain of the Plasma Membrane Glycoprotein Cd36. J. Biol. Chem. 2017, 292, 8683−8693. (36) Parton, D. L.; Akhmatskaya, E. V.; Sansom, M. S. Multiscale Simulations of the Antimicrobial Peptide Maculatin 1.1: Water Permeation through Disordered Aggregates. J. Phys. Chem. B 2012, 116, 8485−8493. (37) Zhao, X.; Yu, H.; Yang, L.; Li, Q.; Huang, X. Simulating the Antimicrobial Mechanism of Human B-Defensin-3 with CoarseGrained Molecular Dynamics. J. Biomol. Struct. Dyn. 2015, 33, 2522−2529. (38) Rzepiela, A. J.; Sengupta, D.; Goga, N.; Marrink, S. J. Membrane Poration by Antimicrobial Peptides Combining Atomistic and CoarseGrained Descriptions. Faraday Discuss. 2010, 144, 431−443. (39) Bond, P. J.; Parton, D. L.; Clark, J. F.; Sansom, M. S. CoarseGrained Simulations of the Membrane-Active Antimicrobial Peptide Maculatin 1.1. Biophys. J. 2008, 95, 3802−3815. (40) Santo, K. P.; Berkowitz, M. L. Difference between Magainin-2 and Melittin Assemblies in Phosphatidylcholine Bilayers: Results from Coarse-Grained Simulations. J. Phys. Chem. B 2012, 116, 3021−3030. (41) Pannuzzo, M.; Milardi, D.; Raudino, A.; Karttunen, M.; La Rosa, C. Analytical Model and Multiscale Simulations of Aβ Peptide Aggregation in Lipid Membranes: Towards a Unifying Description of Conformational Transitions, Oligomerization and Membrane Damage. Phys. Chem. Chem. Phys. 2013, 15, 8940−8951. (42) Khalfa, A.; Tarek, M. On the Antibacterial Action of Cyclic Peptides: Insights from Coarse-Grained Md Simulations. J. Phys. Chem. B 2010, 114, 2676−2684. (43) Gkeka, P.; Sarkisov, L. Interactions of Phospholipid Bilayers with Several Classes of Amphiphilic A-Helical Peptides: Insights from Coarse-Grained Molecular Dynamics Simulations. J. Phys. Chem. B 2010, 114, 826−839. (44) Zhao, J.; Zhao, C.; Liang, G.; Zhang, M.; Zheng, J. Engineering Antimicrobial Peptides with Improved Antimicrobial and Hemolytic Activities. J. Chem. Inf. Model. 2013, 53, 3280−3296. (45) Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. Gromacs: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701−1718. (46) Wassenaar, T. A.; Ingólfsson, H. I.; Böckmann, R. A.; Tieleman, D. P.; Marrink, S. J. Computational Lipidomics with Insane: A Versatile Tool for Generating Custom Membranes for Molecular Simulations. J. Chem. Theory Comput. 2015, 11, 2144−2155. (47) Qi, Y.; Ingólfsson, H. I.; Cheng, X.; Lee, J.; Marrink, S. J.; Im, W. Charmm-Gui Martini Maker for Coarse-Grained Simulations with the Martini Force Field. J. Chem. Theory Comput. 2015, 11, 4486−4494. (48) Hu, J.; Chen, C.; Zhang, S.; Zhao, X.; Xu, H.; Zhao, X.; Lu, J. R. Designed Antimicrobial and Antitumor Peptides with High Selectivity. Biomacromolecules 2011, 12, 3839−3843. (49) Sani, M.-A.; Whitwell, T. C.; Separovic, F. Lipid Composition Regulates the Conformation and Insertion of the Antimicrobial 17274

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275

Article

The Journal of Physical Chemistry C Peptide Maculatin 1.1. Biochim. Biophys. Acta, Biomembr. 2012, 1818, 205−211. (50) DeLano, W. L. Pymol Molecular Graphics System; Schrödinger, LLC, 2002 (51) Berendsen, H. J.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684−3690. (52) Humphrey, W.; Dalke, A.; Schulten, K. Vmd: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33−38. (53) Torrie, G. M.; Valleau, J. P. Nonphysical Sampling Distributions in Monte Carlo Free-Energy Estimation: Umbrella Sampling. J. Comput. Phys. 1977, 23, 187−199. (54) Kumar, S.; Rosenberg, J. M.; Bouzida, D.; Swendsen, R. H.; Kollman, P. A. The Weighted Histogram Analysis Method for FreeEnergy Calculations on Biomolecules. I. The Method. J. Comput. Chem. 1992, 13, 1011−1021. (55) Castillo, N.; Monticelli, L.; Barnoud, J.; Tieleman, D. P. Free Energy of Walp23 Dimer Association in Dmpc, Dppc, and Dopc Bilayers. Chem. Phys. Lipids 2013, 169, 95−105. (56) Rausch, J. M.; Marks, J. R.; Rathinakumar, R.; Wimley, W. C. BSheet Pore-Forming Peptides Selected from a Rational Combinatorial Library: Mechanism of Pore Formation in Lipid Vesicles and Activity in Biological Membranes. Biochemistry 2007, 46, 12124−12139. (57) Melo, M. N.; Ferre, R.; Castanho, M. A. Antimicrobial Peptides: Linking Partition, Activity and High Membrane-Bound Concentrations. Nat. Rev. Microbiol. 2009, 7, 245−250. (58) Wang, L.; Sun, N.; Terzyan, S.; Zhang, X.; Benson, D. R. A Histidine/Tryptophan Π-Stacking Interaction Stabilizes the HemeIndependent Folding Core of Microsomal Apocytochrome B5 Relative to That of Mitochondrial Apocytochrome B5. Biochemistry 2006, 45, 13750−13759. (59) Krishna Deepak, R.; Sankararamakrishnan, R. N−H··· N Hydrogen Bonds Involving Histidine Imidazole Nitrogen Atoms: A New Structural Role for Histidine Residues in Proteins. Biochemistry 2016, 55, 3774−3783. (60) Chen, L.; Liang, J. F. The Potential Roles of Cell Surface Phs in Bioactive Peptide Activation. Chem. Biol. Drug Des. 2015, 85, 208− 215. (61) Helmlinger, G.; Sckell, A.; Dellian, M.; Forbes, N. S.; Jain, R. K. Acid Production in Glycolysis-Impaired Tumors Provides New Insights into Tumor Metabolism. Clin. Cancer. Res. 2002, 8, 1284− 1291. (62) Du, J.; Yarema, K. J. Carbohydrate Engineered Cells for Regenerative Medicine. Adv. Drug Delivery Rev. 2010, 62, 671−682. (63) Dobrzyńska, I.; Szachowicz-Petelska, B.; Sulkowski, S.; Figaszewski, Z. Changes in Electric Charge and Phospholipids Composition in Human Colorectal Cancer Cells. Mol. Cell. Biochem. 2005, 276, 113−119. (64) Yang, L.; Weiss, T. M.; Lehrer, R. I.; Huang, H. W. Crystallization of Antimicrobial Pores in Membranes: Magainin and Protegrin. Biophys. J. 2000, 79, 2002−2009. (65) Sullivan, R.; Santarpia, P.; Lavender, S.; Gittins, E.; Liu, Z.; Anderson, M.; He, J.; Shi, W.; Eckert, R. Clinical Efficacy of a Specifically Targeted Antimicrobial Peptide Mouth Rinse: Targeted Elimination of Streptococcus Mutans and Prevention of Demineralization. Caries Res. 2011, 45, 415−428. (66) Polyansky, A. A.; Ramaswamy, R.; Volynsky, P. E.; Sbalzarini, I. F.; Marrink, S. J.; Efremov, R. G. Antimicrobial Peptides Induce Growth of Phosphatidylglycerol Domains in a Model Bacterial Membrane. J. Phys. Chem. Lett. 2010, 1, 3108−3111. (67) Oreopoulos, J.; Epand, R. F.; Epand, R. M.; Yip, C. M. PeptideInduced Domain Formation in Supported Lipid Bilayers: Direct Evidence by Combined Atomic Force and Polarized Total Internal Reflection Fluorescence Microscopy. Biophys. J. 2010, 98, 815−823. (68) Shah, R. N.; Shah, N. A.; Lim, M. M. D. R.; Hsieh, C.; Nuber, G.; Stupp, S. I. Supramolecular Design of Self-Assembling Nanofibers for Cartilage Regeneration. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 3293−3298.

(69) Liao, H.-S.; Lin, J.; Liu, Y.; Huang, P.; Jin, A.; Chen, X. SelfAssembly Mechanisms of Nanofibers from Peptide Amphiphiles in Solution and on Substrate Surfaces. Nanoscale 2016, 8, 14814−14820. (70) Tripathi, J. K.; Pal, S.; Awasthi, B.; Kumar, A.; Tandon, A.; Mitra, K.; Chattopadhyay, N.; Ghosh, J. K. Variants of Self-Assembling Peptide, Kld-12 That Show Both Rapid Fracture Healing and Antimicrobial Properties. Biomaterials 2015, 56, 92−103. (71) Huang, H. W. Action of Antimicrobial Peptides: Two-State Model. Biochemistry 2000, 39, 8347−8352. (72) Pouny, Y.; Rapaport, D.; Mor, A.; Nicolas, P.; Shai, Y. Interaction of Antimicrobial Dermaseptin and Its Fluorescently Labeled Analogues with Phospholipid Membranes. Biochemistry 1992, 31, 12416−12423. (73) Fjell, C. D.; Hiss, J. A.; Hancock, R. E.; Schneider, G. Designing Antimicrobial Peptides: Form Follows Function. Nat. Rev. Neurosci. 2012, 11, 37−51. (74) Zhang, M.; Zhao, J.; Zheng, J. Molecular Understanding of a Potential Functional Link between Antimicrobial and Amyloid Peptides. Soft Matter 2014, 10, 7425−7451. (75) Stella, L.; Mazzuca, C.; Venanzi, M.; Palleschi, A.; Didone, M.; Formaggio, F.; Toniolo, C.; Pispisa, B. Aggregation and WaterMembrane Partition as Major Determinants of the Activity of the Antibiotic Peptide Trichogin Ga Iv. Biophys. J. 2004, 86, 936−945. (76) Boland, M. P.; Separovic, F. Membrane Interactions of Antimicrobial Peptides from Australian Tree Frogs. Biochim. Biophys. Acta, Biomembr. 2006, 1758, 1178−1183. (77) Chia, C.; Torres, J.; Cooper, M.; Arkin, I.; Bowie, J. The Orientation of the Antibiotic Peptide Maculatin 1.1 in Dmpg and Dmpc Lipid Bilayers. Support for a Pore-Forming Mechanism. FEBS Lett. 2002, 512, 47−51. (78) Chia, B.; Carver, J. A.; Mulhern, T. D.; Bowie, J. H. Maculatin 1.1, an Anti-Microbial Peptide from the Australian Tree Frog, Litoria Genimaculata. Eur. J. Biochem. 2000, 267, 1894−1908. (79) Balla, M.; Bowie, J.; Separovic, F. Solid-State Nmr Study of Antimicrobial Peptides from Australian Frogs in Phospholipid Membranes. Eur. Biophys. J. 2004, 33, 109−116. (80) Fan, R.; Yuan, Y.; Zhang, Q.; Zhou, X.-R.; Jia, L.; Liu, Z.; Yu, C.; Luo, S.-Z.; Chen, L. Isoleucine/Leucine Residues at ″a″ and ″D″ Positions of a Heptad Repeat Sequence Are Crucial for the Cytolytic Activity of a Short Anticancer Lytic Peptide. Amino Acids 2017, 49, 193. (81) Sun, D.; Forsman, J.; Woodward, C. E. Multistep Molecular Dynamics Simulations Identify the Highly Cooperative Activity of Melittin in Recognizing and Stabilizing Membrane Pores. Langmuir 2015, 31, 9388−9401. (82) Sun, F.; Xu, L.; Chen, P.; Wei, P.; Qu, J.; Chen, J.; Luo, S.-Z. Insights into the Packing Switching of the Epha2 Transmembrane Domain by Molecular Dynamic Simulations. J. Phys. Chem. B 2015, 119, 7816−7824. (83) De Smet, K.; Contreras, R. Human Antimicrobial Peptides: Defensins, Cathelicidins and Histatins. Biotechnol. Lett. 2005, 27, 1337−1347. (84) Bechinger, B.; Zasloff, M.; Opella, S. J. Structure and Orientation of the Antibiotic Peptide Magainin in Membranes by Solid-State Nuclear Magnetic Resonance Spectroscopy. Protein Sci. 1993, 2, 2077−2084. (85) Dathe, M.; Schümann, M.; Wieprecht, T.; Winkler, A.; Beyermann, M.; Krause, E.; Matsuzaki, K.; Murase, O.; Bienert, M. Peptide Helicity and Membrane Surface Charge Modulate the Balance of Electrostatic and Hydrophobic Interactions with Lipid Bilayers and Biological Membranes. Biochemistry 1996, 35, 12612−12622. (86) Sharon, M.; Oren, Z.; Shai, Y.; Anglister, J. 2d-Nmr and Atr-Ftir Study of the Structure of a Cell-Selective Diastereomer of Melittin and Its Orientation in Phospholipids. Biochemistry 1999, 38, 15305−15316.

17275

DOI: 10.1021/acs.jpcc.7b04347 J. Phys. Chem. C 2017, 121, 17263−17275