Flipped C-Terminal Ends of APOA1 Promote ABCA1-Dependent Cholesterol Efflux by Small HDLs.

BACKGROUND: Cholesterol efflux capacity (CEC) predicts cardiovascular disease independently of high-density lipoprotein (HDL) cholesterol levels. Isolated small HDL particles are potent promoters of macrophage CEC by the ABCA1 (ATP-binding cassette transporter A1) pathway, but the underlying mechanisms are unclear. METHODS: We used model system studies of reconstituted HDL and plasma from control and lecithin-cholesterol acyltransferase (LCAT)-deficient subjects to investigate the relationships among the sizes of HDL particles, the structure of APOA1 (apolipoprotein A1) in the different particles, and the CECs of plasma and isolated HDLs. RESULTS: We quantified macrophage and ABCA1 CEC of 4 distinct sizes of reconstituted HDL. CEC increased as particle size decreased. Tandem mass spectrometric analysis of chemically cross-linked peptides and molecular dynamics simulations of APOA1, the major protein of HDL, indicated that the mobility of C-terminus of that protein was markedly higher and flipped off the surface in the smallest particles. To explore the physiological relevance of the model system studies, we isolated HDL from LCAT-deficient subjects, whose small HDLs (like reconstituted HDLs) are discoidal and composed of APOA1, cholesterol, and phospholipid. Despite their very low plasma levels of HDL particles, these subjects had normal CEC. In both the LCAT-deficient subjects and control subjects, the CEC of isolated extra-small HDL (a mixture of extra-small and small HDL by calibrated ion mobility analysis) was 3- to 5-fold greater than that of the larger sizes of isolated HDL. Incubating LCAT-deficient plasma and control plasma with human LCAT converted extra-small and small HDL particles into larger particles, and it markedly inhibited CEC. CONCLUSIONS: We present a mechanism for the enhanced CEC of small HDLs. In smaller particles, the C-termini of the 2 antiparallel molecules of APOA1 are "flipped" off the lipid surface of HDL. This extended conformation allows them to engage with ABCA1. In contrast, the C-termini of larger HDLs are unable to interact productively with ABCA1 because they form a helical bundle that strongly adheres to the lipid on the particle. Enhanced CEC, as seen with the smaller particles, predicts decreased cardiovascular disease risk. Thus, extra-small and small HDLs may be key mediators and indicators of the cardioprotective effects of HDL.

Design principles of PI(4,5)P2 clustering under protein-free conditions: Specific cation effects and calcium-potassium synergy.

Phosphatidylinositol 4,5-bisphosphate (PIP2) clustering is a key component in cell signaling, yet little is known about the atomic-level features of this phenomenon. Network-theoretic analysis of multimicrosecond atomistic simulations of PIP2 containing asymmetric bilayers under protein-free conditions, presented here, reveals how design principles of PIP2 clustering are determined by the specific cation effects. Ca2+ generates large clusters (6% are pentamer or larger) by adding existing PIP2 dimers formed by strong O‒Ca2+‒O bridging interactions of unprotonated P4/P5 phosphates. In contrast, monovalent cations (Na+ and K+) form smaller and less-stable clusters by preferentially adding PIP2 monomers. Despite having the same net charge, the affinity to P4/P5 is higher for Na+, while affinity toward glycerol P1 is higher for K+. Consequently, a mixture of K+ and Ca2+ (as would be produced by Ca2+ influx) synergistically yields larger and more stable clusters than Ca2+ alone due to the different binding preferences of these cations.

Spontaneous curvature generation by peptides in asymmetric bilayers.

Peptides and proteins play crucial roles in membrane remodeling by inducing spontaneous curvature. However, extracting spontaneous curvatures from simulations of asymmetric bilayers is challenging because differential stress (i.e., the difference of the leaflet surface tensions) arising from leaflet area strains can vary substantially among initial conditions. This study investigates peptide-induced spontaneous curvature δc 0 p in asymmetric bilayers consisting of a single lipid type and a peptide confined to one leaflet; δc 0 p is calculated from the Helfrich equation using the first moment of the lateral pressure tensor and an alternative expression using the differential stress. It is shown that differential stress introduced during initial system generation is effectively relaxed by equilibrating using P21 periodic boundary conditions, which allows lipids to switch leaflets across cell boundaries and equalize their chemical potentials across leaflets. This procedure leads to robust estimates of δc 0 p for the systems simulated, and is recommended when equality of chemical potentials between the leaflets is a primary consideration.

Initiation and evolution of pores formed by influenza fusion peptides probed by lysolipid inclusion.

The fusion peptide (FP) domain is necessary for the fusogenic activity of spike proteins in a variety of enveloped viruses, allowing the virus to infect the host cell, and is the only part of the protein that interacts directly with the target membrane lipid tails during fusion. There are consistent findings of poration by this domain in experimental model membrane systems, and, in certain conditions, the isolated FPs can generate pores. Here, we use molecular dynamics simulations to investigate the specifics of how these FP-induced pores form in membranes with different compositions of lysolipid and POPC. The simulations show that pores form spontaneously at high lysolipid concentrations via hybrid intermediates, where FP aggregates in the cis leaflet tilt to form a funnel-like structure that spans the leaflet and locally reduces the hydrophobic thickness that must be traversed by water to form a pore. By restraining a single FP within an FP aggregate to this tilted conformation, pores can be formed in lower-lysolipid-content membranes, including pure POPC, on the 100-ns timescale, much more rapidly than in unbiased simulations in bilayers with the same composition. The pore formation pathway is similar to the spontaneous formation in high lysolipid concentrations. Depending on the membrane composition, the pores can be metastable (as seen in POPC) or lead to membrane rupture.

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations.

Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments alter their lipid composition in order to maintain a specific viscosity, and membrane viscosity has been linked to the rate of cellular respiration. To understand the factors that determine the viscosity of a membrane, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with the stress-tensor autocorrelation function modeled by a stretched exponential function. By simulating a series of lipids at different temperatures, we establish the dependence of viscosity on several aspects of lipid chemistry, including hydrocarbon chain length, unsaturation, and backbone structure. Sphingomyelin is found to have a remarkably high viscosity, roughly 20 times that of DPPC. Furthermore, we find that inclusion of the entire range of the dispersion interaction increases viscosity by up to 140%. The simulated viscosities are similar to experimental values obtained from the rotational dynamics of small chromophores and from the diffusion of integral membrane proteins but significantly lower than recent measurements based on the deformation of giant vesicles.

USP7 small-molecule inhibitors interfere with ubiquitin binding.

The ubiquitin system regulates essential cellular processes in eukaryotes. Ubiquitin is ligated to substrate proteins as monomers or chains and the topology of ubiquitin modifications regulates substrate interactions with specific proteins. Thus ubiquitination directs a variety of substrate fates including proteasomal degradation. Deubiquitinase enzymes cleave ubiquitin from substrates and are implicated in disease; for example, ubiquitin-specific protease-7 (USP7) regulates stability of the p53 tumour suppressor and other proteins critical for tumour cell survival. However, developing selective deubiquitinase inhibitors has been challenging and no co-crystal structures have been solved with small-molecule inhibitors. Here, using nuclear magnetic resonance-based screening and structure-based design, we describe the development of selective USP7 inhibitors GNE-6640 and GNE-6776. These compounds induce tumour cell death and enhance cytotoxicity with chemotherapeutic agents and targeted compounds, including PIM kinase inhibitors. Structural studies reveal that GNE-6640 and GNE-6776 non-covalently target USP7 12 Å distant from the catalytic cysteine. The compounds attenuate ubiquitin binding and thus inhibit USP7 deubiquitinase activity. GNE-6640 and GNE-6776 interact with acidic residues that mediate hydrogen-bond interactions with the ubiquitin Lys48 side chain, suggesting that USP7 preferentially interacts with and cleaves ubiquitin moieties that have free Lys48 side chains. We investigated this idea by engineering di-ubiquitin chains containing differential proximal and distal isotopic labels and measuring USP7 binding by nuclear magnetic resonance. This preferential binding protracted the depolymerization kinetics of Lys48-linked ubiquitin chains relative to Lys63-linked chains. In summary, engineering compounds that inhibit USP7 activity by attenuating ubiquitin binding suggests opportunities for developing other deubiquitinase inhibitors and may be a strategy more broadly applicable to inhibiting proteins that require ubiquitin binding for full functional activity.

Finite-Size Effects in Simulations of Peptide/Lipid Assembly.

Molecular dynamics simulations are an attractive tool for understanding lipid/peptide self-assembly but can be plagued by inaccuracies when the system sizes are too small. The general guidance from self-assembly simulations of homogeneous micelles is that the total number of surfactants should be three to five times greater than the equilibrium aggregate number of surfactants per micelle. Herein, the heuristic is tested on the more complicated self-assembly of lipids and amphipathic peptides using the Cooke and Martini 3 coarse-grained models. Cooke model simulations with 50 to 1000 lipids and no peptide are dominated by finite-size effects, with usually one aggregate (micelle or nanodisc) containing most of the lipids forming at each system size. Approximately 200 systems of different peptide/lipid (P/L) ratios and sizes of up to 1000 lipids yield a "finite-size phase diagram" for peptide driven self-assembly, including a coexistence region of micelles and discs. Insights from the Cooke model are applied to the assembly of dimyristoylphosphatidylcholine and the ELK-neutral peptide using the Martini 3 model. Systems of 150, 450, and 900 lipids with P/L = 1/6.25 form mixtures of lipid-rich discs that agree in size with experiment and peptide-rich micelles. Only the 150-lipid system shows finite-size effects, which arise from the long-tailed distribution of aggregate sizes. The general rule of three to five times the equilibrium aggregate size remains a practical heuristic for the Cooke and Martini 3 systems investigated here.

Developing initial conditions for simulations of asymmetric membranes: a practical recommendation.

It has been proposed that the surface tension difference between leaflets (or differential stress) in asymmetric bilayers is generally nonvanishing. This implies that there is no unique approach to generate initial conditions for simulations of asymmetric bilayers in the absence of experimentally derived constraints. Current generation methods include individual area per lipid (APL) based, leaflet surface area (SA) matching, and zero leaflet tension based (0-DS). This work adds a bilayer-based approach that aims for achieving partial chemical equilibrium by interleaflet switching of selected lipids via P21 periodic boundary conditions. Based on a recently proposed theoretical framework, we obtained expressions for tensions in asymmetric bilayers from both the bending and area strains. We also developed a quantitative measure for the energetic penalty from the differential stress. The impacts of APL-, SA-, and 0-DS-based approaches on mechanical properties are assessed for two different asymmetric bilayers. The lateral pressure profile and its moments differ significantly for each method, whereas the area compressibility modulus is relatively insensitive. Application of P21 periodic boundary conditions (APL/P21, SA/P21, and 0-DS/P21) results in better agreement in mechanical properties between asymmetric bilayers generated by APL-, SA-, and 0-DS-based approaches, in which changes are the smallest for bilayers from the SA-based method. The estimated differential stress from the theory shows good agreement with that from the simulations. These simulation results and the good agreement between the predicted and observed differential stress further support the theoretical framework in which bilayer mechanical properties are outcomes of the interplay between intrinsic bending and asymmetric lipid packing. Based on the simulation results and theoretical predictions, the SA/P21-based, or at least the SA-based (when the differential stress is small), approach is recommended as a practical method for developing initial conditions for asymmetric bilayer simulations.

Molecular Dynamics Simulations of Membrane Permeability.

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

Tertiary structure of apolipoprotein A-I in nascent high-density lipoproteins.

Understanding the function of high-density lipoprotein (HDL) requires detailed knowledge of the structure of its primary protein, apolipoprotein A-I (APOA1). However, APOA1 flexibility and HDL heterogeneity have confounded decades of efforts to determine high-resolution structures and consistent models. Here, molecular dynamics simulations totaling 30 µs on two nascent HDLs, each with 2 APOA1 and either 160 phospholipids and 24 cholesterols or 200 phospholipids and 20 cholesterols, show that residues 1-21 of the N-terminal domains of APOA1 interact via strong salt bridges. Residues 26-43 of one APOA1 in the smaller particle form a hinge on the disc edge, which displaces the C-terminal domain of the other APOA1 to the phospholipid surface. The proposed structures are supported by chemical cross-linking, Rosetta modeling of the N-terminal domain, and analysis of the lipid-free ∆185APOA1 crystal structure. These structures provide a framework for understanding HDL maturation and revise all previous models of nascent HDL.

Incorporation of α-methylated amino acids into Apolipoprotein A-I mimetic peptides improves their helicity and cholesterol efflux potential.

Apolipoprotein A-I (ApoA-I) mimetic peptides are potential therapeutic agents for promoting the efflux of excess cellular cholesterol, which is dependent upon the presence of an amphipathic helix. Since α-methylated Ala enhances peptide helicity, we hypothesized that incorporating other types of α-methylated amino acids into ApoA-I mimetic peptides may also increase their helicity and cholesterol efflux potential. The last helix of apoA-I, peptide 'A' (VLESFKVSFLSALEEYTKKLNT), was used to design peptides containing a single type of α-methylated amino acid substitution (Ala/Aα, Glu/Dα, Lys/Kα, Leu/Lα), as well as a peptide containing both α-methylated Lys and Leu (6α). Depending on the specific residue, the α-helical content as measured by CD-spectroscopy and calculated hydrophobic moments were sometimes higher for peptides containing other types of α-methylated amino acids than those with α-methylated Ala. In ABCA1-transfected cells, cholesterol efflux to the peptides showed the following order of potency: 6α>Kα≈Lα≈Aα≫Dα≈A. In general, α-methylated peptides were resistant to proteolysis, but this varied depending on the type of protease and specific amino acid substitution. In summary, increased helicity and amphilicity due to α-methylated amino acid substitutions in ApoA-I mimetic peptides resulted in improved cholesterol efflux capacity and resistance to proteolysis, indicating that this modification may be useful in the future design of therapeutic ApoA-I mimetic peptides.

Optimization of globomycin analogs as novel gram-negative antibiotics.

Discovery of novel classes of Gram-negative antibiotics with activity against multi-drug resistant infections is a critical unmet need. As an essential member of the lipoprotein biosynthetic pathway, lipoprotein signal peptidase II (LspA) is an attractive target for antibacterial drug discovery, with the natural product inhibitor globomycin offering a modestly-active starting point. Informed by structure-based design, the globomycin depsipeptide was optimized to improve activity against E. coli. Backbone modifications, together with adjustment of physicochemical properties, afforded potent compounds with good in vivo pharmacokinetic profiles. Optimized compounds such as 51 (E. coli MIC 3.1 µM) and 61 (E. coli MIC 0.78 µM) demonstrate broad spectrum activity against gram-negative pathogens and may provide opportunities for future antibiotic discovery.

Membrane permeability of small molecules from unbiased molecular dynamics simulations.

Permeation of many small molecules through lipid bilayers can be directly observed in molecular dynamics simulations on the nano- and microsecond timescale. While unbiased simulations provide an unobstructed view of the permeation process, their feasibility for computing permeability coefficients depends on various factors that differ for each permeant. The present work studies three small molecules for which unbiased simulations of permeation are feasible within less than a microsecond, one hydrophobic (oxygen), one hydrophilic (water), and one amphiphilic (ethanol). Permeabilities are computed using two approaches: counting methods and a maximum-likelihood estimation for the inhomogeneous solubility diffusion (ISD) model. Counting methods yield nearly model-free estimates of the permeability for all three permeants. While the ISD-based approach is reasonable for oxygen, it lacks precision for water due to insufficient sampling and results in misleading estimates for ethanol due to invalid model assumptions. It is also demonstrated that simulations using a Langevin thermostat with collision frequencies of 1/ps and 5/ps yield oxygen permeabilities and diffusion constants that are lower than those using Nosé-Hoover by statistically significant margins. In contrast, permeabilities from trajectories generated with Nosé-Hoover and the microcanonical ensemble do not show statistically significant differences. As molecular simulations become more affordable and accurate, calculation of permeability for an expanding range of molecules will be feasible using unbiased simulations. The present work summarizes theoretical underpinnings, identifies pitfalls, and develops best practices for such simulations.

Structure and Function in Antimicrobial Piscidins: Histidine Position, Directionality of Membrane Insertion, and pH-Dependent Permeabilization.

Piscidins are histidine-enriched antimicrobial peptides that interact with lipid bilayers as amphipathic α-helices. Their activity at acidic and basic pH in vivo makes them promising templates for biomedical applications. This study focuses on p1 and p3, both 22-residue-long piscidins with 68% sequence identity. They share three histidines (H3, H4, and H11), but p1, which is significantly more permeabilizing, has a fourth histidine (H17). This study investigates how variations in amphipathic character associated with histidines affect the permeabilization properties of p1 and p3. First, we show that the permeabilization ability of p3, but not p1, is strongly inhibited at pH 6.0 when the conserved histidines are partially charged and H17 is predominantly neutral. Second, our neutron diffraction measurements performed at low water content and neutral pH indicate that the average conformation of p1 is highly tilted, with its C-terminus extending into the opposite leaflet. In contrast, p3 is surface bound with its N-terminal end tilted toward the bilayer interior. The deeper membrane insertion of p1 correlates with its behavior at full hydration: an enhanced ability to tilt, bury its histidines and C-terminus, induce membrane thinning and defects, and alter membrane conductance and viscoelastic properties. Furthermore, its pH-resiliency relates to the neutral state favored by H17. Overall, these results provide mechanistic insights into how differences in the histidine content and amphipathicity of peptides can elicit different directionality of membrane insertion and pH-dependent permeabilization. This work features complementary methods, including dye leakage assays, NMR-monitored titrations, X-ray and neutron diffraction, oriented CD, molecular dynamics, electrochemical impedance spectroscopy, surface plasmon resonance, and quartz crystal microbalance with dissipation.

Gramicidin A Channel Formation Induces Local Lipid Redistribution II: A 3D Continuum Elastic Model.

To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.

Gramicidin A Channel Formation Induces Local Lipid Redistribution I: Experiment and Simulation.

Integral membrane protein function can be modulated by the host bilayer. Because biological membranes are diverse and nonuniform, we explore the consequences of lipid diversity using gramicidin A channels embedded in phosphatidylcholine (PC) bilayers composed of equimolar mixtures of di-oleoyl-PC and di-erucoyl-PC (dC18:1+dC22:1, respectively), di-palmitoleoyl-PC and di-nervonoyl-PC (dC16:1+dC24:1, respectively), and di-eicosenoyl-PC (pure dC20:1), all of which have the same average bilayer chain length. Single-channel lifetime experiments, molecular dynamics simulations, and a simple lipid compression model are used in tandem to gain insight into lipid redistribution around the channel, which partially alleviates the bilayer deformation energy associated with channel formation. The average single-channel lifetimes in the two-component bilayers (95 ± 10 ms for dC18:1+dC22:1 and 195 ± 20 ms for dC16:1+dC24:1) were increased relative to the single-component dC20:1 control bilayer (65 ± 10 ms), implying lipid redistribution. Using a theoretical treatment of thickness-dependent changes in channel lifetimes, the effective local enrichment of lipids around the channel was estimated to be 58 ± 4% dC18:1 and 66 ± 2% dC16:1 in the dC18:1+dC22:1 and dC16:1+dC24:1 bilayers, respectively. 3.5-µs molecular dynamics simulations show 66 ± 2% dC16:1 in the first lipid shell around the channel in the dC16:1+dC24:1 bilayer, but no significant redistribution (50 ± 4% dC18:1) in the dC18:1+dC22:1 bilayer; these simulated values are within the 95% confidence intervals of the experimental averages. The strong preference for the better matching lipid (dC16:1) near the channel in the dC16:1+dC24:1 mixture and lesser redistribution in the dC18:1+dC22:1 mixture can be explained by the energetic cost associated with compressing the lipids to match the channel's hydrophobic length.

Identification of a novel lipid binding motif in apolipoprotein B by the analysis of hydrophobic cluster domains.

Apolipoprotein B (apoB) is a large amphipathic protein that is the structural scaffold for the formation of several classes of lipoproteins involved in lipid transport throughout the body. The goal of the present study was to identify specific domains in the apoB sequence that contribute to its lipid binding properties. A sequence analysis algorithm was developed to identify stretches of hydrophobic amino acids devoid of charged amino acids, which are referred to as hydrophobic cluster domains (HCDs). This analysis identified 78 HCDs in apoB with hydrophobic stretches ranging from 6 to 26 residues. Each HCD was analyzed in silico for secondary structure and lipid binding properties, and a subset was synthesized for experimental evaluation. One HCD peptide, B38, showed high affinity binding to both isolated HDL and LDL, and could exchange between lipoproteins. All-atom molecular dynamics simulations indicate that B38 inserts 3.7Å below the phosphate plane of the bilayer. B38 forms an unusual α-helix with a broad hydrophobic face and polar serine and threonine residues on the opposite face. Based on this structure, we hypothesized that B38 could efflux cholesterol from cells. B38 showed a 12-fold greater activity than the 5A peptide, a bihelical Class A amphipathic helix (EC50 of 0.2658 vs. 3.188µM; p<0.0001), in promoting cholesterol efflux from ABCA1 expressing BHK-1 cells. In conclusion, we have identified novel domains within apoB that contribute to its lipid biding properties. Additionally, we have discovered a unique amphipathic helix design for efficient ABCA1-specific cholesterol efflux.

Simulations of Membrane-Disrupting Peptides I: Alamethicin Pore Stability and Spontaneous Insertion.

An all-atom molecular dynamics simulation of the archetype barrel-stave alamethicin (alm) pore in a 1,2-dioleoyl-sn-glycero-3-phosphocholine bilayer at 313 K indicates that ∼7 µs is required for equilibration of a preformed 6-peptide pore; the pore remains stable for the duration of the remaining 7 µs of the trajectory, and the structure factors agree well with experiment. A 5 µs simulation of 10 surface-bound alm peptides shows significant peptide unfolding and some unbinding, but no insertion. Simulations at 363 and 413 K with a -0.2 V electric field yield peptide insertion in 1 µs. Insertion is initiated by the folding of residues 3-11 into an α-helix, and mediated by membrane water or by previously inserted peptides. The stability of five alm pore peptides at 413 K with a -0.2 V electric field demonstrates a significant preference for a transmembrane orientation. Hence, and in contrast to the cationic antimicrobial peptide described in the following article, alm shows a strong preference for the inserted over the surface-bound state.

Simulations of Membrane-Disrupting Peptides II: AMP Piscidin 1 Favors Surface Defects over Pores.

Antimicrobial peptides (AMPs) that disrupt bacterial membranes are promising therapeutics against the growing number of antibiotic-resistant bacteria. The mechanism of membrane disruption by the AMP piscidin 1 was examined with multimicrosecond all-atom molecular dynamics simulations and solid-state NMR spectroscopy. The primary simulation was initialized with 20 peptides in four barrel-stave pores in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol bilayer. The four pores relaxed to toroidal by 200 ns, only one porelike structure containing two transmembrane helices remained at 26 µs, and none of the 18 peptides released to the surface reinserted to form pores. The simulation was repeated at 413 K with an applied electric field and all peptides were surface-bound by 200 ns. Trajectories of surface-bound piscidin with and without applied fields at 313 and 413 K and totaling 6 µs show transient distortions of the bilayer/water interface (consistent with (31)P NMR), but no insertion to transmembrane or pore states. (15)N chemical shifts confirm a fully surface-bound conformation. Taken together, the simulation and experimental results imply that transient defects rather than stable pores are responsible for membrane disruption by piscidin 1, and likely other AMPs.