The Location of the Antimicrobial Peptide Maculatin 1.1 in Model Bacterial Membranes

Maculatin 1.1 (Mac1) is an antimicrobial peptide (AMP) from the skin secretions of Australian tree frogs. In this work, the interaction of Mac1 with anionic phospholipid bilayers was investigated by NMR, circular dichroism (CD) spectroscopy, neutron reflectometry (NR) and molecular dynamics (MD). In buffer, the peptide is unstructured but in the presence of anionic (DPC/LMPG) micelles or (DMPC/DMPG/DHPC) bicelles adopts a helical structure. Addition of the soluble paramagnetic agent gadolinium (Gd-DTPA) into the Mac1-DPC/LMPG micelle solution showed that the N-terminus is more exposed to the hydrophilic Gd-DTPA than the C-terminus in micelles. 2H and 31P solid-state NMR showed that Mac1 had a greater effect on the anionic lipid (DMPG). A deuterium labeled Mac1 used in NR experiments indicated that the AMP spanned across anionic (PC/PG) bilayers, which was compatible with MD simulations. Simulations also showed that Mac1 orientation remained transmembrane in bilayers and wrapped on the surface of the micelles regardless of the lipid or detergent charge. Thus, the peptide orientation appears to be more susceptible to curvature than charged surface. These results support the formation of transmembrane pores by Mac1 in model bacterial membranes.


INTRODUCTION
Maculatin 1.1 (Mac1)s is an antimicrobial peptide (AMP) from the skin secretions of the Australian tree frog Litoria genimaculata (Rozek et al., 1998). The peptide forms part of the frog's innate immune system and is effective at killing a wide range of Gram-positive bacteria (Fernandez et al., 2009). This makes Mac1, along with a range of other AMPs found across nature, a possibility for development of alternative antibiotics (Lee et al., 2015). Alternatives are required due to the increasing prevalence of antibiotic resistance by bacteria to the most commonly used antibiotics currently available, making what were once treatable infections increasingly difficult to treat (Rice, 2009). Mechanistically, what makes most AMPs an attractive target is that they target disrupting the cell membrane rendering the cell unviable and making it improbable for resistance to occur, rather than targeting a metabolic process where bacteria can evolve to develop resistance (Lee et al., 2019).
Bacteria are broadly classified as either Gram-positive or Gram-negative based on the difference in the cell envelopes and their membranes consist of at least 15% anionic lipids (Epand and Epand, 2009). Membrane-active AMPs usually disrupt bacterial membranes in one of three ways: lysis of the membrane through an action known as the carpet mechanism, formation of a toroidal pore or formation of a barrel-stave pore (Sani and Separovic, 2016). The mode of action of AMPs is normally determined by the presence of positively charge residues, amphiphilicity, secondary structure, and lipid composition and charge of the target cell membrane (Koehbach and Craik, 2019). Previous investigations found that AMPs severely perturb anionic bacterial membranes but are less active against eukaryotic membranes (Balhara et al., 2013;Lee et al., 2013). The first driver of AMPs to approach a target cell membrane is the electrostatic interaction between the positively charged residues of the peptide and the negatively charged cell surface. The second driver is then the hydrophobic interactions between the amphipathic domains of the peptide and the acyl chains of the lipids that make up the membrane (Brogden, 2005). As peptides increase in length, net charge and/or hydrophobicity, the ability to disturb anionic lipid systems also increases (Jiang et al., 2008;Lee et al., 2014). Mac1 is 21 amino acids long and is disordered in aqueous solution but when bound to a membrane can span both leaflets of the membrane and adopt an alpha-helix structure (Chia et al., 2000;Fernandez et al., 2013a). At neutral pH Mac1 is cationic so has a strong affinity for anionic membranes such as bacterial membranes but will still bind to zwitterionic membranes (Mechler et al., 2007). The lipid tail composition as well as headgroup makeup will also influence Mac1 membrane binding with a preference for fluid phospholipid bilayers with alkyl chains that are 14 to 16 carbons in length (Sani et al., 2012;Lee et al., 2018). The final characteristic to note is that Mac1 has a proline residue that induces a kink in the structure of the peptide, and substitution of this residue with alanine or glycine results in reduced binding to membranes (Fernandez et al., 2013b). This is the result of the proline creating a wedge which allows the peptide to insert into the bilayer (Sani et al., 2015b).
Determining the structure and orientation of peptides bound within lipids is challenging due to the difficulty in having a single stable conformation of the peptide-lipid aggregate. A number of studies have investigated the location and orientation of Mac1 in lipid bilayers of various lipid compositions and a range of biomimetic membranes. Earlier work showed that Mac1 inserts into the hydrophobic core of the membrane at an angle of approximately 30 • to 40 • from the bilayer normal, depending on peptide concentration and lipid composition used, and with anionic lipids showing larger tilt angles (Chia et al., 2002;Marcotte et al., 2003). However, from this work, the location of Mac1 could not be ascertained. Previous molecular dynamics (MD) simulations concluded that Mac1 could have multiple insertion angles over a large range (0 • to 150 • ) but clusters of multiple peptides arranged in a parallel fashion in planar bilayers of dipalmitoylphosphatidylcholine, with the Nterminus first making contact with the outer leaflet to facilitate further insertion (Bond et al., 2008). Further MD simulations in an anionic POPE/POPG phospholipid system showed that maculatin spans the lipid bilayer in a parallel fashion but, although a higher density of the N-terminus was found in one leaflet, it was not sufficiently high enough to conclude a preferred orientation (Balatti et al., 2018). Other studies also found that Mac1 can adopt several different aggregate arrangements with dimers to pentamers dominant in phospholipids with short chains, such as DMPC, and hexamers or more dominating in lipids with longer chains (Wang et al., 2016). In the aggregates a mix of anti-parallel and parallel arrangements could be adopted but with the anti-parallel form dominating (Wang et al., 2016). In previous experimental and MD simulation work we have found that Mac1 adopts a preferred orientation in zwitterionic phosphatidylcholine (PC) membranes with the N-terminus in the outer leaflet (Sani et al., 2020). Mac1 will bind to a range of different model membranes that replicate various characteristics of the Gram-positive bacterial membrane (Sani et al., 2015a), and in this study we investigate the orientation of Mac1 in anionic phosphatidylglycerol (PG) membranes that better reflect the charge state of bacterial membranes which is the principal target of Mac1. We use neutron reflectometry and solid-state NMR spectroscopy to determine the location of deuterated Mac1 in PC/PG bilayers that better mimic the properties of membranes of Gram-positive bacteria.

Solution NMR Experiments
The solution NMR samples were made of: 1 mM multiple (x 11) 15 N-labeled Mac1 with 98% purity, dissolved in phosphate buffer (10 mM, pH 7.4) containing 50 mM KCl, 1 mM EDTA and 10% (v/v) D 2 O with 120 mM of d 38 -DPC/LMPG (9:1 mol/mol) and 1 mM unlabelled Mac1 in phosphate buffer containing 150 mM d 25 -SDS micelles. 1 H-15 N HSQC spectra of 15 N labeled Mac1 were acquired at 35 • C using a Bruker 600 MHz NMR spectrometer for experiments performed with d 38 -DPC/LMPG micelles. For each time increment, 16 transients were taken with 4k and 512 points in the 1 H and 15 N dimensions, respectively, and a 1.5 s recycle delay. Spectral width of 13 ppm centered at 4.7 ppm for 1 H and 40 ppm centered at 118 ppm for 15 N were used. 1 H-15 N-1 H HSQC-NOESY were performed with same parameters as HSQC experiments and with 4k, 64 and 128 points in the direct and indirect dimensions. For gadolinium titration, a concentrated stock solution of gadolinium diethylenetriamine-pentaacetic acid (Gd-DTPA) in MilliQ water was made and aliquots added directly into the NMR tube to reach 1 mM, 4 mM, 7 mM and 15 mM Gd 3+ . The total added volume of Gd 3+ was ca. 5% of the total sample volume. The NMR tube was allowed to equilibrate for 30 min at 35 • C prior to each measurement. The volume of each resonance after Gd-DTPA addition was obtained by integration, corrected for effect of dilution and normalized relative to the Gd 3+ -free resonance volume.
The unlabelled Mac1 in d 25 -SDS micelles sample was investigated at 37 • C on an 800 MHz Bruker Advance II spectrometer. 1 H homonuclear TOCSY (mixing time τ mix = 80 ms) and NOESY (τ mix = 150 and 300 ms) were acquired with 512 points and 1k points in the F1 dimension, respectively, and 4k points in the F2 dimension. Between 16 and 32 transients were accumulated with a 1.5 s recycle delay. The data were multiplied with a squared sine bell function shifted by 90 • . The 1 H spectral window was set to 9,600 Hz. 13 C-1 H HSQC experiments were performed with 256 points in the F1 dimension and 4k points in the F2 dimensions. 64 transients were accumulated with a 2 s recycle delay. The 13 C spectral window was set to 33,200 Hz. Non-uniform sampling 15 N-1 H HSQC experiments were performed with 37.5% of 128 points in the F1 dimension and 4k points in the F2 dimension. 1,024 Transients were accumulated with a 1.5s recycle delay. The 15 N spectral window was set to 3,240 Hz.
All data dimensions were zero-filled to twice the respective FID size. 1 H chemical shifts were referenced to DSS (sodium trimethylsilylpropanesulfonate) at 0 ppm. Data were processed in Topspin (Bruker) and analyzed using the CCPNmr Analysis program (Vranken et al., 2005). Backbone and side chains were assigned using all experiments.

Structure Calculations
The NOESY cross-peak assignments were subsequently used to generate distance restraints for the structure determination. The nOe distance restraints were supplemented with dihedral angle restraints predicted with DANGLE from H α, H N, N H, C α , C β chemical shifts (Cheung et al., 2010). A standard CNS 1.1-based protocol was employed using the ARIA 2.2 interface (Rieping et al., 2006). The 10 lowest energy structures were refined in a water shell and evaluated with MolProbity (Chen et al., 2010).

Solid-State NMR Experiments
The solid-state NMR sample consisted of the unlabelled Mac1 dissolved in anisotropic bicelles at a lipid to peptide molar ratio of 50:1. The anisotropic bicelles were composed of either d 54 -DMPC and DMPG, or DMPC and d 54 -DMPG, at a molar ratio of 4:1 and DHPC lipids mixed at a molar ratio (DMPC+DMPG)/DHPC (q) of 3.6. The lipid concentration (C L ) was 20% (w/v) in imidazole buffer. The sample was then packed into a 4 mm Bruker MAS rotor. 31 P NMR Experiments 31 P NMR experiments were performed on a 400 MHz Bruker Avance III NMR spectrometer at a frequency of 161.5 MHz. A 4 mm triple resonance probe was used in a double resonance mode. 62.5 kHz direct excitation pulse experiments were used under 31.25 kHz 1 H SPINAL64 decoupling scheme and a recycle delay of 3 s. Typically, 1k scans were acquired and processed with 8k zero-filling and linebroadening from 20 Hz to 100 Hz were used.

H NMR Experiments
The static 2 H solid-state NMR experiments were performed on a 400 MHz Bruker Avance III NMR spectrometer at a frequency of 61.5 MHz. A 4 mm triple resonance probe was used in a double resonance mode. The solid echo pulse sequence was used with 45.5 kHz 2 H excitation, an echo delay of 26 µs and a recycle delay of 0.5 s. The spectra were recorded with a 500 kHz spectral window and typically 128 k transients were accumulated. The FIDs were processed using 4 k zero filling and a 150 Hz line broadening.

Circular Dichroism Experiments
The CD samples were made with phosphate buffer since imidazole prevents signal acquisition below 210 nm. Similar lipid to peptide molar ratios were used as for NMR. CD spectra were acquired on a Chirascan spectropolarimeter (Applied Photophysics Ltd, UK) between 180 and 260 nm using a 0.1 mm path length cylindrical quartz cell (Starna, Hainault, UK). Spectra were acquired at 1 nm intervals, 1 s integration time and 3 scans accumulated. CD signal was recorded in milli-degree units at 35 • C and reported as mean residue ellipticity (MRE) using the conversion formula: where CD is the signal in mdegree, [Mac1] the concentration of peptide in mol.L −1 , L the cell path length in cm and N r the number of residues.
Neutron Reflectometry D 2 O and H 2 O solutions for NR experiments were buffered with 10 mM MOPS pH/D 7.0 and 150 mM NaCl. Samples for NR were prepared on round silicon wafers that were 100 mm in diameter and 10 mm thick. Before use the silicon wafers were cleaned in a Jelight UV-ozone cleaner for 20 min and then extensively washed with 2% (v/v) Hellmanex solution, followed by ultrapure (>18 M ) water, and then analytical grade ethanol. The wafers were dried under a stream of nitrogen. The polished side of the wafer was held against a roughened backing silicon wafer that had inlet and outlet holes that were connected to a HPLC pump for solution exchange. There was a 100 µm thick PTFE gasket that separated the sample and backing wafer creating a 283 µL volume for sample/solution injection. The wafers were held together in aluminum clamping plates that were also connected to a Julabo water bath for temperature control. Solid-supported membranes were created through vesicle deposition as described previously (Fernandez et al., 2013a) by injecting a 0.1 mg mL −1 vesicle solution in the reservoir, incubating for 1 h at 30 • C and then rising excess vesicles away with 5 mL of buffer. After characterization of the membranes 2 mL of d-Mac1 at 10 µM was injected and incubated for 1 h before excess peptide was removed with a 5 mL wash of buffer.
Neutron reflectometry measurements were conducted on the blank silicon wafer, the formed solid-supported membrane, and after Mac1 incubation on the Platypus time-of-flight neutron reflectometer at the 20 MW OPAL Research Reactor (Sydney, Australia) (James et al., 2011). The instrument views the coldneutron source and utilizes a useable wavelength bandwidth of 2.5 to 20 Å using a disc chopper system that was set to a λ/λ ∼ 8.4% and 24 Hz. The instrument has horizontal sample geometry and collimation slits of 0.71 and 2.91 mm were used for an angle of incidence of 0.85 • and 3.5 • . Neutrons were detected on a 2D 3 He detector and the raw data was reduced using in-house software that stitches the two angles together at the appropriate overlap region, re-bins the data to instrument resolution, correctly scales the data so that any critical edge is unity and corrects for detector efficiency (Nelson, 2010). The final data is then presented as reflectivity vs. momentum transfer, Q, which is given according the equation below: where θ is the angle of incidence and λ the wavelength. Data was analyzed using an Abele's matrix method in the MOTOFIT software package (Nelson, 2006). The membrane system is divided into a series of layers and each layer is defined by its thickness (in Å), its scattering length density (SLD) (in Å −2 ) and interfacial roughness (in Å). An initial model is created and then a reflectivity profile calculated and compared against the experimental data. Parameters are then adjusted until the calculated profile adequately matches the experimental data using a least-squares regression. A genetic algorithm is used so that reasonable limits can be placed on each varied parameter and to avoid becoming trapped in global minima. A scattering length density is analogous to a neutron's refractive index and can be calculated according to: where V m is the molecular volume and b i is the scattering length for each isotope within the molecule. The theoretical SLD and molecular volumes of the materials used can be found in Supplementary Table S1. As scattering length varies across each different isotope (e.g., for hydrogen b H = −3.741 ×10 −5 Å and for its stable isotope deuterium b D = +6.667 ×10 −5 Å), different components of each layer can be determined through contrast variation and the total SLD of a layer is determined through the sum of the SLD for each component according to: where Ø is the volume fraction. The phase of the lipid bilayer at the different temperatures was determined through calculating the area per lipid, A lipid , using the fitted properties of the lipid tails as follows: where τ is the thickness of the lipid tails. Error values for each fitted parameter were determined using Monte Carlo resampling as described previously Holt et al., 2009). Briefly, 1008 fits of the data were completed and the distribution of fitted values plotted for each parameter varied. The error is then the 95% confidence interval of the distribution.

Molecular Dynamics
The starting conformation of Mac1 peptide was generated from NMR data. The CHARMM-GUI membrane builder (Jo et al., 2007(Jo et al., , 2008(Jo et al., , 2009Wu et al., 2014) was then used to prepare the peptide-bilayer system with 40 DMPC and 10 DMPG lipids per leaflet and the peptide-micelle complex with 58 DPC and 7 LMPG lipid molecules randomly distributed. Each system was generated using a rectangular box containing 50 mM KCl and a 12.5 Å layer of water. All histidine residues were singly protonated to model the ionization state expected at pH 7.4 and the peptide C-terminus was amidated to match the experimental conditions. The simulations were performed using the CHARMM-36m force fields. The minimization, equilibration and production runs were performed with the NAMD package on a desktop machine fitting a GPU GeForce GTX 1080 titanium and a CPU with 12 cores.
Each system was first minimized for 2500 steps using the steepest descent method followed by 2500 steps of the conjugate gradient method with a 12 Å non-bonded interaction cut-off. The peptides and lipids were restrained with a 10 and 2.5 kCal.mol −1 potential, respectively. Then, 25 ps equilibration MD simulations were run at 308 K for each system, using a 10, 5, 2.5 kCal/mol restraint to maintain the peptide backbone and 2.5, 0 and 0 kCal/mol restraint on the lipids atom positions with 0.001 ps time step. Next, the peptide positional restraints were reduced to 1, 0.5, and 0.1 kCal/mol for 100 ps with 0.002 ps time step. All covalent bonds involving hydrogen atoms were constrained using the SHAKE algorithm (Ryckaert et al., 1977) and the rigid internal geometry for TIP3P water molecules was constrained with the SETTLE algorithm (Miyamoto and Kollman, 1992). The system temperature was maintained at 308 K using a Langevin thermostat (Izaguirre et al., 2001) with a 1 ps −1 collision frequency. The system pressure was controlled at 1 bar using a Langevin barostat for the DPC/LMPG micelle and a semi-isotropic Berendsen barostat with a xy surface tension for the DMPC/DMPG bilayer. Finally, the systems were run for 1 ns prior the production runs. 100 ns long simulations were Frontiers in Chemistry | www.frontiersin.org Assignments were made using TOCSY and NOESY experiments and previously published data (Sani et al., 2020). The NMR experiments were performed at 35 • C in imidazole buffer (pH 7.4). (B) CD spectra of Mac1in aqueous buffer (black squares), in the presence of zwitterionic DPC micelles (blue triangles) or anionic SDS micelles (red circles). The CD experiments were performed at 37 • C in phosphate buffer (pH 5). (C) Resonance peak intensity of 15 N labeled Mac1 residues as a function of the Gd-DTPA concentration in the presence of d 38 -DPC/LMPG (9:1) micelles. The NMR experiments were performed at 35 • C in imidazole buffer (pH 7.4).
(D) 10 lowest energy structures of Mac1 calculated from the solution NMR experiments (Supplementary Figure S1 and Supplementary Table S2) performed in SDS micelles (pH 5) at 37 • C. (E) Electrostatic surface representation of Mac1 averaged structure with the hydrophobic (red), hydrophilic (blue) and neutral (white) residues according to Kyte and Doolittle (1982).
performed with a random restart every 1 ns, the trajectories were then concatenated together, and the analysis was performed on the full 100 ns simulation.
The MD trajectories were visualized and analyzed using VMD (Humphrey et al., 1996) with custom scripts and the CPPTRAJ (Roe and Cheatham, 2013) software and fitting procedures and plots were created in Gnuplot.

Effect of Negatively Charged Membranes on Mac1 Helical Conformation and Preferential C-Terminus Insertion
The effect of the membrane curvature on the secondary structure and insertion of Mac1 has been recently characterized in neutral membrane mimetics using solution NMR and CD techniques (Sani et al., 2020). Here, similar investigations were performed in anionic membrane mimetics to probe the role of electrostatic interactions in AMP mode of action. The addition of 10% mol of the anionic LMPG in DPC micelles induced slightly different chemical shift perturbations along the peptide sequence. Comparing the 15 N HSQC obtained in DPC micelles and DPC/LMPG (9:1) showed that the N-terminus of Mac1 exhibited small chemical shift changes at Val 5 , Leu 6 , and Ala 7 while the middle and C-terminus sections exhibited more pronounced chemical shift perturbations, the highest observed for Ala16 ( Figure 1A). Interestingly, the shielding from the hydrophilic gadolinium moiety was greater for the middle section and the C-terminus while the N-terminus was severely exposed to paramagnetic effects as seen in the signal intensity loss upon Gd 3+ titration ( Figure 1C). The circular dichroism spectra of Mac1 in the presence of DPC and SDS micelles showed small differences: the lineshapes exhibited the typical Hα features of two minima at 222 nm and 209 nm and a maximum at about 195 nm (Figure 1B). The 10 lowest energy structures of Mac1 ( Figure 1D) were determined using solution NMR restraints obtained from the NOESY experiments performed in SDS micelles (Supplementary Figure S1). Mac1 showed a continuous Hα stretch with a bend located near Ala16, adjacent to the Pro15 residue. The electrostatic surface of Mac1 was computed using the averaged structure and exhibited an amphipathic distribution. Altogether, these results are consistent with the previously reported structure of Mac1 in DPC micelles and indicate that Mac1 retained a preferential insertion through its C-terminus, despite a slightly greater hydrophilic surface.

Preferential Interaction of Mac1 With Anionic Lipids
Highly curved surfaces are not optimal to understand the molecular mechanism of interaction between Mac1 and cell membranes. Thus, lipid bilayers of DMPC and DMPG (4:1 mol/mol) were used and anisotropic bicelles were formed using a short chain DHPC lipids (q = 3.6). This lipid system increases the spectral resolution due to spontaneous alignment in the magnetic field. By either adding d 54 -DMPC or d 54 -DMPG, specific modulations were observed for neutral or anionic lipids, respectively, and Mac1. The 2 H spectra displayed in Figure 2A clearly showed that Mac1 decreased strongly the dynamics (or increased the order) of the d 54 -DMPG lipid chains as seen in the increase of the 2 H quadrupolar splitting. The 2 H spectrum of the d 54 -DMPC showed little change in the quadrupolar splitting but significant line broadening was observed.
A similar broadening of the overall 31 P line shape was observed in the presence of Mac1, and interestingly, the 31 P chemical shift of DMPG and DHPC was significantly shifted while DMPC 31 P chemical shifts remained mainly unchanged but broadened. This also indicates that PG headgroup is in more curved regions near DHPC; and supports that the cationic Mac1 peptide had a stronger impact on anionic lipids than on the neutral lipids, yet, without heavily disturbing the lipid packing at this lipid to peptide molar ratio.

Mac1 Bound Within Solid-Supported Phospholipid Membranes
Neutron reflectometry is a technique able to probe the structure of surfaces and interfaces at nanoscale dimensions. The information obtained is a one-dimensional description of the structures through the different layers of the system (Penfold and Thomas, 2014). This information is useful for membrane systems as many surface sensitive and microscopy techniques can only view the surface topology whereas NR can probe through the entire bilayer (Wacklin, 2010;Lakey, 2019). The first part of the measurements is creation of a solid-supported phospholipid membrane on a silicon wafer. The membranes, consisting of d 54 -DMPC/d 54 -DMPG or h-DMPC/h-DMPG with a PC to PG mole ratio of 3:1, were deposited using the widely used vesicle deposition technique onto silicon wafers which were found to have a 8.3 ± 3.3 Å and 10.7 ± 2.1 Å oxide layer with a 3.2 Å and 2.6 Å roughness for the deuterated and hydrogenous bilayer, respectively. All reflectometry experiments were characterized in a buffered D 2 O and H 2 O contrast which were fitted simultaneously (Figure 3). For the modeling the bilayer was divided into three separate layers consisting of inner headgroups (those closest to the silicon oxide surface), tails and outer headgroups (those closest to the bulk solvent). The bilayers were characterized at 30 • C and were found to have a good coverage with a volume fraction of 0.894 ± 0.133 and 0.965 ± 0.030, and total thicknesses of 52.0 ± 1.7 Å and 59.4 ± 3.6 Å, which corresponds to an area per lipid of 57.4 ± 8.6 Å 2 and 49.0 ± 1.8 Å 2 each for the deuterated and hydrogenous bilayers, respectively (Supplementary Table S3). These results are consistent with previous observations of solidsupported membranes with similar lipid composition at 30 • C and correspond to being in the liquid crystalline (L α ) phase (Fernandez et al., 2013a). The bilayers were cooled to 15 • C and the bilayer re-measured. Upon cooling the fringe in the reflectivity profile shifts to lower Q-values ( Figure 3B) indicating the bilayer has become thicker. After complete data analysis the bilayer properties at 15 • C the area per lipid reduces to 45.4 ± 0.8 Å 2 and 45.0 ± 5.2 Å 2 for the deuterated and hydrogenous bilayers, respectively (Supplementary Table S3), indicating that the bilayer is now in the gel (L β ) phase, with the thickness for each bilayer also increasing with the volume fraction remaining largely unchanged (within error). It should be noted that bilayers deposited directly onto solid-supports through methods such a vesicle deposition do not display the ripple (P β ) phase (Naumann et al., 1992).
After characterization of the bilayers 2 mL of d-Mac1 prepared in buffered D 2 O at 10 µM was flowed into the sample cell and left to incubate on the bilayer for 1 h at 30 • C. After the incubation period the excess peptide was removed with a wash of buffered D 2 O the NR was measured. A version of Mac1 where six of the amino acids in the N-terminal half of the peptide are deuterated creates a substantially different SLD between the N-and C-terminal halves of the peptide (see Supplementary Table S1 and section Materials for details). This means that if Mac1 embeds within the membrane in a specific orientation this will be highlighted as the two leaflets of the membrane will have different SLD values. As can be seen in Figure 3A the reflectivity profiles before and after peptide look different suggesting that d-Mac1 has bound to the bilayer. Initially, the bilayers with peptide bound were fitted using the same three layer model above. However, a four layer model, where the tails were split into two separate layers of inner and outer tails with the SLD allowed to vary for each layer, was more suitable. After peptide addition it was observed that, overall the combined thickness of the tails layers increased in the presence of d-Mac1 (Figure 4 and Supplementary Table S4), which is consistent with previous observations of h-Mac1 binding to anionic lipid bilayers (Fernandez et al., 2013a). The SLD of the tail layer changed in both leaflets with an overall decrease for the deuterated bilayer ( Figure 4A) and increase for the hydrogenous bilayer (Figure 4B), which indicated that the d-Mac1 has inserted across both leaflets of the bilayer.
Examining more closely the SLD of each leaflet, a difference between the two sides of the bilayer (Figure 4 and Supplementary Table S5), which is particularly notable in the h-DMPC/h-DMPG (3:1) membrane (Figure 4B), can be seen. One interpretation is that the variation in SLD across the bilayer leaflet after peptide addition is due to an uneven distribution of d-Mac1 across the lipid leaflets. Upon calculating the volume fractions of d-Mac1 using this scenario, the volume fraction for the inner/outer leaflet are 0.172/0.209 and 0.118/0.182 for the deuterated and hydrogenous bilayer, respectively (Supplementary Table S4). However, as d-Mac1 is 21 amino acids and long enough to span a bilayer then, if the peptide is assembling in the lipid bilayer in a parallel fashion, the uneven SLD could be due to a preferred orientation. In most cases the SLD of the inner leaflet is slightly higher than that of the outer leaflet. This would suggest an orientation where the N-terminal half of the peptide is buried in the lower leaflet, which is curiously the opposite orientation of the d-Mac1 in a pure DMPC bilayer where the N-terminal half is found in the outer leaflet (Sani et al., 2020). The volume fraction of peptide in this scenario would be 0.256 for the deuterated bilayer and 0.176 for the hydrogenous bilayer. In either scenario, lipid is lost from the surface suggesting some lytic activity. It should be noted that for the anionic  membranes used in this study, the difference in SLD this time is within error of each other and there is also a contribution to the SLD of the layer from the volume fraction of solvent present, which would have an uneven distribution across the two bilayer leaflets.
The bilayers with d-Mac1 bound were cooled to 15 • C and remeasured. As for peptide-free bilayers the reflectivity shifts to lower Q-values (Figure 3), indicating that the membrane has become thicker. Looking at the values for thickness the bilayer does indeed become thicker with an increase of 2.2 Å (Supplementary Table S4). A difference in SLD across the two leaflets was also maintained showing that changes in the phase state of the membrane for a DMPC/DMPG (3:1) composition is not influenced by the presence of d-Mac1. Additionally, cooling the DMPC/DMPG (3:1) bilayers to the gel phase did not alter the orientation or distribution of d-Mac1 within the bilayer with no change observed in the orientation of the peptide (Figure 5). The difference in SLD in the inner tails increases when the temperature was changed from 30 • C to 15 • C in the d 54 -DMPC/d 54 -DMPG (3:1) case (Supplementary Table S5) with this difference corresponding to a loss of lipid and, therefore, an increase in the peptide volume fraction (Supplementary Table S4).
FIGURE 6 | Analysis of the trajectories obtained from 100 ns simulation of Mac1 in lipid environments. N-Cα-C backbone rmsd from the equilibrated structure in: (A) a DPC/LMPG (9:1) micelle, and (B) a DMPC/DMPG (4:1) bilayer using residues 1-21 (black line). The averaged rmsd was fitted with a linear function (red line). Distance fluctuation between the 1 Gly and the 21 Phe nitrogens of Mac1 inserted in: (C) a DPC/LMPG (9:1) micelle, and (D) a DMPC/DMPG (4:1) bilayer. The averaged distance was obtained by fitting the data to a linear function (red line).

Mac1 Exposure to Solvent Is Dependent on Membrane Curvature Rather Than Surface Charge
Full atom MD simulations were used to gain residue specific details of Mac1 interaction with anionic micelles and bilayers. The initial configurations were set as a single Mac1 peptide fully inserted into a micelle made up of 54 DPC and 6 LMPG or a bilayer made of 40 DMPC and 10 DMPG molecules per leaflet. Each system was equilibrated for ca. 1 ns at 35 • C as described in the Methods section. The 35 • C temperature was chosen as this is well above the phase transition temperature for both lipid types and are thus unlikely to form different phase domains (Lewis et al., 2005).

Effect of Anionic Lipids on Mac1 Topology and Length
Mac1 migrated to a peripheral location with a wrapped configuration around the micelle but remained in a transmembrane orientation in the bilayer, a similar outcome as for MD simulations performed in neutral lipid systems (Sani et al., 2020). The N-Cα-C backbone structure fluctuations were also similar for both anionic systems and neutral systems across the 100 ns simulations, indicating that Mac1 secondary structure was not greatly different in neutral and anionic lipid systems. Furthermore, the peptide length, measured as the head to tail distance between the nitrogen backbone atoms, (Figure 6) showed little variation in the micelle (about 3 nm) but some fluctuation in the bilayer system were observed with a significantly lower averaged peptide length (2.7 nm). As seen in Figure 7, the difference corresponded to a flexible bend initiated at Val 14 (before Pro 15 ), bringing the lowest peptide length to 2.3 nm, well below the optimal theoretical length for a 21 amino acid long linear peptide in an α-helical configuration of 3.15 nm.

Insertion Depth and Water Penetration
The depth of Mac1 penetration into the lipid core (Supplementary Figures S2, S3) and its exposure to water (Supplementary Figures S4, S5) were calculated over the 100 ns simulations. On average, Mac1 was positioned mainly just below the phosphates in the micelle (Figure 8A), the middle section of the peptide (from Leu 6 to Val 14 , excepting Lys 8 ) was slightly deeper within the hydrophobic core and the N-terminus was a little more exposed than the C-terminus ( Figure 8C) as previously observed with MD simulation performed in neutral DPC micelle (Sani et al., 2020). In the bilayer system, the peptide maintained a transmembrane orientation (Figure 8B), with residues 6-18 significantly shielded from water exposure. Interestingly, some water molecules were found to reside for a notable amount of time near Lys 8 (side chain), His 12 and Ala 16 ( Figure 8D). A snapshot showed Mac1 with a bent structure and 8 Lys extending to reach the phosphate-water interface and water flowing along the hydrophilic i,i+4 side of the peptide as shown in Figure 1E. Note, however, that the simulations performed herein were with an inserted peptide but earlier MD results (Wang et al., 2016) show that Mac1 inserts transmembrane into a phospholipid bilayer from the aqueous phase.
FIGURE 7 | Depth-cued snapshot of Mac1 in DMPC/DMPG (4:1) bilayers (sticks) at 57 ns showing a bent structure with water (dark blue surface) penetrating through the bilayer. Lys 8 is snorkeling out to reach the water interface and His 12 , Pro 15 /Ala 16 , and Asp 19 provide the hydrophilic surface for water molecules to hydrate the peptide within the hydrophobic core.

Peptide Pairing With Anionic vs. Neutral Lipids
Mac1 is a cationic peptide, but with a rather low overall charge of +1 at neutral pH. It has been shown to interact with many types of lipid membranes, inducing severe leakage in neutral and anionic membranes, but to a greater extent in the former (Fernandez et al., 2013a). However, electrostatic attractions have been shown to dominate the first binding step, as competition between neutral and anionic lipids have a larger preference to the latter surface. Nonetheless, once bound and inserted within the lipid hydrophobic core, MD simulations showed that the lipid headgroup has little influence on the peptide-lipid pairing as no segregation was observed within the simulation timeframe (Figure 9). Several possibilities are to be considered: (1) the simulation was not long enough to reach a pairing equilibrium; (2) Mac1 is mainly located below the phosphate headgroup and thus less prone to sense the headgroup difference between phosphocholine and phosphoglycerol moieties; and 3) the density of anionic lipids is not sufficient for statistically significant pairing to occur.

DISCUSSION
Previous studies investigating the location and orientation of Mac1 within phospholipid membranes focused on using zwitterionic PC headgroups (Sani et al., 2020), which generally model most eukaryotic cell membranes well, but do not replicate the anionic nature of bacterial membranes which are rich in FIGURE 8 | Averaged distance from the N th residue to the center of mass of (A) the DPC/LMPG (9:1) micelle, and (B) DMPC/DMPG (4:1) bilayer, as a proxy of Mac1 insertion depth into the hydrophobic core of the lipids. Averaged contact number between the n th residue nitrogen and water molecules within 5 Å in (C) DPC/LMPG (9:1) micelle, and (D) DMPC/DMPG (4:1) bilayer. The dashed lines correspond to the phosphates averaged position to the center of mass of the lipid systems.
PG and other anionic lipids. Whilst knowing how AMPs behave in zwitterionic membranes environments is important as it has to be understood how host cells are going to be affected by AMP treatment, ultimately it is the interactions with the bacterial pathogen that are of interest and this work has addressed the structure of the peptide-lipid aggregates that form. We find that the secondary structure of Mac1 remains unchanged whether in a zwitterionic lipid or an anionic environment (Figure 1). This suggests that headgroup composition plays little or no role in Mac1 structure once embedded within a membrane and that, as previous studies have shown, tail composition is the important factor in driving Mac1 conformation when embedded within a membrane (Sani et al., 2012). In addition, Mac1 still embeds deeply into the membrane in the presence of anionic lipids and assembles in a parallel fashion as observed for zwitterionic membranes (Sani et al., 2020). However, with zwitterionic membranes and the use of d-Mac1 it was observed that Mac1 adopts a preferred orientation with the N-terminus in the outer leaflet of the membrane (Sani et al., 2020). In this work, however, a preferred orientation within anionic membranes could not be conclusively determined in solid-supported membranes as the difference in SLD between the tails in each leaflet was with within error of each other (Figure 4 and Supplementary Table S5). This does not mean that a preferred orientation is ruled out as a preferred orientation in this case can be masked either due to the greater solvent content observed with d-Mac1 bound to anionic membranes, or some uneven distribution of the peptide across the bilayer leaflets. The location of d-Mac1 also was not influenced by the phase of the membrane (Figure 5).
Mac1 orientation was inconclusive in solid-supported anionic bilayers with the difference in SLD either being due to a uneven distribution of peptide across the different leaflets or parallel assembly with a preferential orientation. Other studies have examined the orientation of Mac1 within a phopholipid bilayer and show a variety of orientations based on different lipid compositions and model membranes (Chia et al., 2002;Bond et al., 2008;Wang et al., 2016), with a defined orientation being inconclusive for anionic bilayers (Balatti et al., 2018). However, in each case Mac1 spans the bilayer and, given that the volume fractions between each leaflet in the solid-supported bilayers are within error of each other (Supplementary Table S4) and the MD simulations suggest a similar orientation as for zwitterionic membranes with the N-terminus more exposed to water (Figure 8), it is likely that in this case d-Mac1 also spans the anionic bilayer. Where previous MD work has shown uneven distribution of Mac1 across the leaflets of a planar bilayer (Bond et al., 2008;Balatti et al., 2018), the difference in distribution is small (<15% difference) and these small difference would not greatly impact on the NR results, further supporting a membrane spanning model. Interestingly, the helical structure of Mac1 generates a hydrophilic surface, which in addition to Lys 8 snorkeling out to reach the water interface, allowed water to penetrate deeply into the bilayer hydrophobic core (Figures 7, 8). This may be the initial step for promoting peptide-peptide selfassembly that would form larger pore structures, as previously demonstrated experimentally (Sani et al., 2013) and in µs-long MD simulations (Wang et al., 2016).
It has been previously shown that membranes containing anionic lipids bind more peptide molecules per lipid (Mechler et al., 2007;Fernandez et al., 2013a), despite neutral membranes being more prone to lytic activity (Sani et al., 2014), and this work is no exception to that. However, what is shown here is that anionic lipids are important in their preferential binding to cationic Mac1. Figure 2 demonstrates that Mac1 has a preference for anionic lipids and has little impact on zwitterionic lipids in a mixed phospholipid bilayer. Possibly the presence of Mac1 causes the PC and PG headgroups to laterally separate (Fernandez et al., 2013b) but the techniques used in this work do not demonstrate this, which is the subject of future studies.
There is a wide range of AMPs that target anionic lipids either for entry into a bacterial cell to reach its final target or because the mode of action is membrane disruption (Epand and Epand, 2009). How different AMPs target and assemble in bacterial membranes varies greatly depending on the peptide and the lipids used. In the case of indolicidin, when bound to model bacterial membranes the distribution of lipids and the membrane thickness are unaffected (Nielsen et al., 2018). Whereas, the short 13 amino acid peptide aurein 1.2, which operates through the carpet mechanism causes lateral segregation of mixed lipid systems (Sharma and Qian, 2019). Our results show that for the case of Mac1 when presented to a phospholipid bilayer of mixed anionic and zwitterionic lipids, the AMP preferentially assembles in a transmembrane fashion and targets the anionic lipids.
In conclusion, Mac1 embeds itself within a model bacterial membrane in a fashion that spans the whole membrane. The presence of anionic lipids does not dramatically alter the overall mechanism of Mac1 compared to membranes composed entirely of zwitterionic lipids; however, the presence of anionic lipids does attract Mac1 to the membrane which allows the peptide to better target membranes from pathogenic bacteria over host membranes.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
AL, SZ, and M-AS conducted the experiments and the data analysis. FS conceived and lead the work. All authors wrote the paper.

FUNDING
The authors thank the Australian Research Council for funding (DP160100959 to FS and DE140101788 to AL) and ANSTO for the awarding of neutron beam time (P3890).