Thomas G.
Meikle
ab,
Charlotte E.
Conn
c,
Frances
Separovic
a and
Calum J.
Drummond
*bc
aSchool of Chemistry, Bio21 Institute, University of Melbourne, VIC 3010, Australia
bCSIRO Manufacturing Flagship, Bag 10, Clayton, VIC 3169, Australia. E-mail: calum.drummond@rmit.edu.au
cSchool of Science, College of Science, Engineering and Health, RMIT University, Melbourne, VIC 3000, Australia
First published on 19th July 2016
Lipid based bicontinuous cubic mesophases provide a low-cost, robust membrane mimetic nanomaterial which allows for the incorporation of membrane peptides and proteins. However, the relationship between the mesostructure of the host lipidic bicontinuous mesophase, the chemical structure of its constituents, and the secondary structure of encapsulated biomolecules is complex, and at present largely unclear. Here we have examined the effects of adding gramicidin A′, an anti-microbial peptide, to lipidic bicontinuous cubic phases composed of a number of different lipids. We demonstrate, using a combination of synchrotron small angle X-ray scattering and circular dichroism, that fundamental physicochemical parameters of the lipid mesophase impact both the structural response to peptide addition, and the conformation of the encapsulated peptide. We have rationalised the results with reference to hydrophobic mismatch, the putative lateral pressure profile and the intrinsic surface curvature of each lipid system. Results should be of use for several applications of hybrid peptide-lipid materials including peptide based drug delivery and the design of in meso crystallization trials.
Bicontinuous cubic phases consist of a lipid bilayer, curved in three-dimensional space into the shape of a triply periodic minimal surface. Three bicontinuous cubic structures are experimentally observed in pure lipid systems, including the Schwarz diamond (QDII with crystallographic space group Pn3m), the primitive (QPII with space group Im3m) and the Schoen gyroid (QGII with space group Ia3d).6 The bilayer in each of these cubic mesophases features zero mean curvature and a high surface area to volume ratio, and results in a network of two interpenetrating but unconnected water channels (Fig. 1).
![]() | ||
Fig. 1 3D representation of the cubic Im3m mesophase, with the lipid bilayer shown in red, and the two interpenetrating water channels in blue and green. The bilayer is removed on the right side to clearly show the water channel structure. Modelled with K3Dsurf and Blender.7,8 |
Lipidic bicontinuous cubic mesophases offer a range of advantages for the incorporation of proteins and peptides. The properties of lipidic cubic mesophases are comparable to those existing in biological membranes, and thus encapsulated proteins and peptides can potentially retain their activity. These properties make the lipidic cubic phase a highly prospective material for drug delivery applications,5,9 including the encapsulation and controlled release of therapeutic proteins and peptides within the body. They are particularly useful for hydrophobic or amphiphilic molecules not suited to traditional administration routes. A number of proteins and peptides have demonstrated controlled release from lipidic cubic mesophases.10,11 Complex, triggered release materials are possible,12 whereby proteins encapsulated in the bilayer function to increase and decrease the rate of drug release via an opening and closing of pores or channels. The lipidic cubic phase can also be used to probe fundamental peptide–lipid interactions.
Guest biomolecules, when inserted into the bilayer, can influence the symmetry and structural unit dimensions of lipidic cubic mesophases.13 If we are to take full advantage of these materials we must understand the complex relationship between the chemical composition of the bilayer, the lipidic mesostructure, and the guest biomolecule. We have studied the effects of the antimicrobial peptide gA′ on cubic phases comprised of a range of commercially available lipids, tracking the changes over a period of two days. The evolution of the cubic structure as the concentration of peptide was increased was correlated with both the structure of the host lipid and the structure of gA′ within the bilayer.
Gramicidin A′ is a mixture of the linear pentadecapeptides gramicidin A, B and C; comprising ∼80%, ∼6% and ∼14%, respectively (Fig. 2). This antimicrobial peptide is produced by Bacillus brevis and exerts its effect by forming a dimeric channel spanning the bilayer and facilitating the uncontrolled movement of monovalent cations (Na+) between the intra- and extra-cellular environment, leading to cell death. Gramicidin A′ is well studied,14–16 cheap, stable, and structure grade crystals have been grown in meso,17 indicating that the peptide is readily reconstituted in the lipidic cubic phase. In addition, pore forming peptides have applications as drugs18 and can be used to further our understanding of a range of diseases where the formation of pores in the cell membrane are believed to play a role, including Alzheimer's disease19 and bacterial infections.20
Since its sequence elucidation in 1964,21 there has been extensive work to understand the function and structure of gA′. The peptide can exist as a variety of different conformations depending on its local environment22 and forms two major folding types: the double helix and the helical dimer (Fig. 3). The double helix, also termed the “pore” conformation is the result of two monomers forming a series of intermolecular hydrogen bonds along their length in either a parallel or anti-parallel β-sheet like motif, forming a double helix. This is the dominant motif present in solution in most organic solvents,23 and usually consists of four interconverting conformational variant species at equilibrium.
The helical dimer or “channel” conformation results from a series of intramolecular hydrogen bonds occurring along the length of the molecule which warps to form a single helical monomer. The dimer is formed when two of these monomers are joined anti-parallel, N-terminal to N-terminal through 6 intermolecular hydrogen bonds linking the formyl groups and the N-terminal amino acids with one another. This conformation is also referred to as a β6.3 structure,24 due to the 6.3 residues per turn in the helix.
Caffrey et al. has demonstrated that addition of gA′ to a series of cubic phase monoacylglycerol (MAG) lipids results in a change in symmetry, an effect closely correlated with hydrophobic mismatch.25 Herein we report the effect of peptide addition to a less studied class of cubic-phase forming lipids, the branched chain or isoprenoid-type lipids. These include phytantriol (PT), a saturated branched chain lipid frequently used as a cosmetic ingredient, and phytanoyl monoethanolamide (PE), a custom lipid synthesized in-house and designed to maintain cubic symmetry at low temperatures suitable for in meso crystallization.26 Similar lipids are found in archaebacteria and show enhanced stability compared to MAGs. The branched hydrocarbon chain architecture has been shown to result in increased local lateral pressure when compared to equivalent straight chain lipids.27 In addition, we expand upon previous studies on the effect of hydrophobic mismatch using MAG lipids by including another MAG, monovaccenin (MV), in addition to monoolein (MO) and monopalmitolein (MP). We have investigated a mixture of monoolein + 8 wt% cholesterol, an important biomolecule which, at this concentration, has been found to facilitate crystal growth in meso. These experiments were designed to explore the mesophase behaviour of a variety of commonly used cubic phase lipids subsequent to peptide addition, and to compare between two distinct classes of lipid: the monoacylglycerols and the branched chain isoprenoids. For direct relevance to applications such as in meso crystallization and drug delivery, we have varied the lipid between predominately commercially available lipids (MO, MV, MP, PT), and one custom synthesized lipid (PE).
Gramicidin A′, dissolved in either methanol or ethanol at 8 different concentrations (one for each final gramicidin A′ concentration), was dispensed into each well of a 96-well plate using a Mosquito robot (TTP Labtech Ltd, Melbourn, UK). Note that due to robotic protocols a set volume of each peptide stock was used for all lipids, such that the final mol% value differs slightly for each lipid. Lipids were dissolved in ethanol to a concentration of 200 mg mL−1 and 1.05 μL was deposited in each well over the peptide giving a lipid mass of 0.21 mg. The plate was then dried for 5 days under vacuum to remove all solvent. The dry plates were hydrated at various times prior to SAXS screening by dispensing 0.84 μL of water into each well, resulting in a solvent to lipid ratio of 80:
20. Lastly, 20 μL of water was added to the reservoir of each well and the plate sealed.
![]() | ||
Fig. 5 The chemical structures of each of the five bicontinuous cubic phase forming lipids used in our SAXS experiments. |
To assess the overall reproducibility of the data, duplicate samples were run at each data-point. The lattice parameters obtained from both data-sets are plotted against each other in the ESI Fig. S1.† The line y = x, representative of an entirely reproducible data-set is provided as a guide to the eye. Note that where coexisting phases were observed, the phase consistent between the two data-points was plotted. The data were highly reproducible for all plates tested, with a difference of 3 Å or less between top and bottom wells in 89% and 96% of wells in 2 and 48 h plates, respectively. Wells that proved irreproducible were discarded from further analysis.
The addition of gA′ at increasing concentrations produced similar phase behaviour in all three MAG lipid systems; an initial decrease in lattice parameter with increasing gramicidin concentration was observed with a decrease in the magnitude of the slope approaching 15 mol% peptide. The data points were fitted with a first order exponential decay (Fig. 6), with the fitting parameters listed in the ESI.† For each lipid, the relationship between lattice parameter and gA′ concentration for the Pn3m range was found to be statistically significant (in all cases p < 0.01, with the exception of MO + Chol. at 2 hours, for which p < 0.05). Despite a considerable variation in the initial unit cell size, the lattice parameter at approximately 15 mol% peptide was relatively similar for all three MAGs. For monoolein and monovaccenin, increasing peptide concentration resulted in a phase transition from the QDII phase, to a QDII phase coexisting with an inverse hexagonal (HII) mesophase (at 10 mol% and 3 mol% peptide, respectively). These HII phases coexisted with the QDII phase up to 15 mol% peptide, the highest peptide concentration studied. Additionally, at approximately 15 mol% peptide, both monoolein and monopalmitolein exhibited a coexisting gyroid (QGII) cubic phase with a lattice parameter of 104 and 106.6 Å, respectively. A summary of all phase changes observed from 0 to 15 mol% peptide is provided in Table 1.
Lipid | Bilayer thickness (Å) | Observed phase sequence |
---|---|---|
Monoolein | 32.1 | QDII → QDII + HII → QDII + HII + QGII |
Monopalmitolein | 30.3 | QDII → QDII + QGII |
Monovaccenin | 35.4 | QDII → QDII + HII |
Monoolein + 8% cholesterol | — | QDII → QDII + HII → QDII + HII + QGII |
Phytantriol | 29.6 | QDII → QDII + QGII |
Phytanoyl monoethanolamide | 25.8 | QDII → QDII + HII |
The binary system consisting of monoolein with 8 wt% cholesterol exhibited the QDII → HII transition at lower concentrations of gA′ than pure monoolein. An HII phase coexisting with the QDII phase was observed at 2 mol% gA′ and higher. By 15 mol% gA′, three coexisting phases were observed: QGII, QDII and HII phases, with the QGII lattice parameter measuring 101.0 Å. Whilst the presence of three coexisting phases is not common in pure lipid systems, the Gibbs phase rule allows for up to four coexisting phases for the quaternary system MO–cholesterol–gramicidin A′–H2O at a given temperature and pressure. The change in phase behaviour observed with increasing peptide concentration strongly indicates that gA′ is incorporated within the bilayer. This is consistent with past fluorescence studies and the very low solubility of the peptide in water.25
In phytantriol samples, a coexisting QGII phase with lattice parameter 104.8 and 104.4 Å at 3 and 5 mol% gA′, respectively, was evident. These peaks broadened with increasing peptide concentration such that determination of the lattice parameter was no longer possible, and resulted in a single very wide peak above 10 mol%. The breadth of these peaks is suggestive of a partially ordered system with a gradient of lattice parameters, consistent with a phase undergoing transition. With increasing peptide concentration phytanoyl monoethanolamide displayed the phase sequence QDII → QDII + HII (15 mol%).
The two main conformations adopted by gA′, the helical dimer and the double helix, are distinguishable by their characteristic CD spectra. It is believed that the active form in the bilayer is the helical dimer, although this issue is still debated and both conformations may be present in varying ratios.
For comparison with the CD spectra of gA′ in cubic phase lipids, we obtained spectra of the peptide in SDS micelles, where it is predominantly the helical dimer,31,32 and in THF, in which it is present as the double helix conformer.33 In summing these two spectra together, each with a 50% weighting, we obtain a third spectrum which closely resembles that of gramicidin A′ obtained in various lipidic mesophases. We therefore conclude that at the concentrations used in our CD study, gramicidin A′ is present as a mixture of the two conformations, with the magnitude of ellipticity at 220 nm serving as an indication of the ratio of one to the other.
The CD spectra of gA′ encapsulated in MAG lipid bilayers were all very similar (Fig. 8), and indicative of an approximately 50:
50 mixture of the helical dimer and the double helix conformations. This result is in agreement with recent observations from Patrick et al., 2016, which demonstrated that the peptide was present as a mixture of conformations under certain conditions.34 Note that for these samples the lipid absorbance caused saturation of the detector and precluded CD analysis below 215 nm. In phytantriol, however, a branched chain lipid, the CD spectra differed from those seen with the MAGS: the ellipticity at 220 nm is well below zero and overall the spectrum more closely resembled that obtained in THF. This suggests that the double helix is present in greater amounts in phytantriol. By fitting a weighted average of the spectra in SDS micelles (representing the helical dimer conformation) and the spectra in THF (representing the double helix) we approximated the contributions of each conformation to the CD spectra observed in the mesophase. For monoolein, the contributions are approximately 39% double helix, 61% helical dimer; and for phytantriol, 62% double helix, 38% helical dimer (see Fig. S2†). Whilst these numbers should not be considered to accurately reflect absolute percentages (the original spectra in THF and SDS micelles are assumed to contain 100% double helix and helical dimer respectively) they do give an indication of the relative amount of each conformation present. These results confirm the influence of lipid bilayer properties on the adopted conformation and will be discussed further in the following section.
The differences in mesophase behaviour of each of the different MAG lipids, may reflect differences in bilayer thickness and hence hydrophobic mismatch with the peptide. From the lattice parameters measured at 0 mol% peptide we can estimate the bilayer thickness via the infinitely periodic minimal surface model (IPMS) outlined by Turner et al.39Eqn (1) gives the lipid fraction (Φl) for a sample of known composition:
![]() | (1) |
![]() | (2) |
Given the lipid fraction, and the lattice parameter in Ångstroms (a) eqn (2) can be solved for lipid length (l). Here A0 is the dimensionless surface area, given by A0 = A/V2/3, where A is the surface area per unit cell of volume V = a3. The Euler characteristic χ provides a measure of the integrated Gaussian curvature per unit area of the surface. These values are specific to each cubic phase symmetry (A0 = 1.919 and χ = −2 for QDII). The lipid length is doubled to give the bilayer thickness values shown in Table 1, increasing in the order PE < PT < MP < MO < MV (this was not calculated for the cholesterol containing samples as the excess water point of this mixture is unknown).
A comparison of the channel length of the helical dimer gA′ conformation (believed to be the active form present within the bilayer, and the conformation responsible for HII formation40) with the bilayer thickness of each cubic phase MAG lipid suggests that hydrophobic mismatch (defined as the lipid length minus the peptide length) plays a significant role in determining lipid phase behaviour. The length of the helical dimer conformation is 25 Å, so that the hydrophobic mismatch for monovaccenin (−10.4 Å) is significantly larger than for monopalmitolein (−5.3 Å).
We find that the concentration at which the HII phase is formed is correlated with bilayer thickness and hence hydrophobic mismatch. The hydrophobic mismatch between the gA′ and the bilayer is greatest for MV, which has a calculated bilayer thickness of 35.4 Å in excess water. An HII phase is observed at relatively low peptide concentrations (3 mol% gA′) for this system. As the hydrophobic mismatch is reduced the peptide concentration at which the HII phase is observed increases. Thus, for MO, with a bilayer of 32.1 Å, an HII phase is not observed until 10 mol% gA′. For MP, which has the thinnest bilayer of the MAGs (30.3 Å) and hence a smaller hydrophobic mismatch, we note the absence of the HII phase even at the highest concentrations of peptide. Our results are in agreement with previous research indicating that greater hydrophobic mismatch leads to the formation of the HII phase at lower concentrations of gA′.25 The length of the gramicidin dimer appears to be a driver of phase behaviour for all lipids. This is supported by the observation that, despite a wide variation in the initial lattice parameter, the lattice parameters of all cubic phases tested approach a similar value at high peptide concentrations (Fig. 6 and 7). We suggest that peptide incorporation causes a thinning of the bilayer and concomitant reduction in lattice parameter, and that for all lipids this reduction continues until the cubic phase bilayer thickness matches the hydrophobic span of the peptide.
The ability of cholesterol to induce a coexisting HII phase at lower peptide concentrations in monoolein may be explained by its ‘rigidification’ or ordering effect on the bilayer.41 This effect has been observed in other studies42 at concentrations similar to those used here. We suggest that the increase in bilayer rigidity (or bending modulus) associated with the addition of cholesterol reduces the ability of the QDII phase to accommodate the gA′ by increasing its intrinsic curvature; triggering a phase transition to HII at lower concentrations of peptide.
The effects of variations in hydrophobic mismatch between peptides and membranes have been previously explored by Killian et al.43 on WALP peptides (peptides consisting of a repeating Leu–Ala chain capped on either end with Trp, designed to model transmembrane peptides) in diacylphosphatidylcholines of various lengths. The resultant phase of these systems is highly dependent on hydrophobic mismatch. Systems with a negative mismatch preferentially form HII phases, while small positive mismatches formed isotropic phases, and large positive mismatches formed bilayer structures. Similar behaviour was observed for gA′.44
Gramicidin A′ and similar peptides are proposed to promote HII formation by increasing the volume of the hydrocarbon region,31,45 and then occupying the positions in the HII phase in which the monolayer thickness is smallest.40 In the absence of peptide, the HII phase is generally not observed until much higher temperatures, e.g. around 90 °C for hydrated monoolein,46 due to the prohibitive steric strain imposed by the high curvature of the HII phase. These explanations conform well to our observations that at high peptide concentrations the resultant HII phases of MAG lipids all possess similar lattice parameters.
For monoolein, monoolein + 8 wt% cholesterol, monopalmitolein, and phytantriol samples, a QGII cubic was seen to coexist with a HII and/or a QDII cubic phase at 14–15 mol% peptide. This coexistence is unlikely to reflect a thermodynamically stable state, and likely results from the mixture of gramicidin A′ conformations present. The QDII → QGII transition was confirmed through the use of the Bonnet ratio: the predicted ratio between lattice parameters of coexisting cubic mesophases in excess water under equilibrium conditions.47 For example, in phytantriol with 3 mol% gA′ we observed both QDII (lattice parameter of 67.1 Å) and QGII (lattice parameter of 104.8 Å), the ratio of which (1.562) closely matches with that of the Bonnet ratio for coexisting QDII/QGII mesophases (1.576).
However, in comparing the disparate phase behaviour of monopalmitolein and phytantriol, two lipids with relatively similar bilayer thickness, it becomes clear that the bilayer thickness is not the only factor governing phase behaviour. Here we must also consider differences in the lateral pressure profile between the different lipid bilayers. For the MAG lipids, similarities in both the headgroup and the acyl chains will result in a very similar lateral pressure profile for these three lipids. However, previous work based on MD simulations has indicated that the local lateral pressure exerted within branched chain lipids is much higher than within unsaturated lipid chains. This has been demonstrated by Miklavcic et al. in their molecular dynamics study of branched chain archaeal lipid bilayers.27
To probe the effect of variations in lipid structure on the peptide conformation, CD spectra were obtained for gA′ in the three MAG lipids, MO + 8 wt% cholesterol, phytantriol, THF and SDS micelles. The spectra obtained for MO, MP, MV and MO + 8 wt% cholesterol are similar. For these unsaturated cubic phase lipids, the CD spectra resembles neither the helical dimer conformation in SDS micelles nor the double helix observed in THF, but appears to reflect a mixture of the two, with the proportion of the helical dimer conformation increasing in the order MO + 8 wt% cholesterol < MP < MO < MV. The proportion of the helical dimer conformation increases with bilayer thickness, while the reduction in helical dimer conformation caused by the addition of cholesterol is likely due to increased bilayer pressure, a known effect of cholesterol on the bilayer.48 The CD spectrum measured in phytantriol similarly represents a mixture of the two conformations. However, the drop in the ellipticity at 220 nm brings the spectrum closer to that found in THF solution, and hence is consistent with a larger proportion of the double helix conformation. An examination of the channel widths reveals that the double helix (β5.6) conformation is narrower (0.95 Å maximum pore radius) than the helical dimer (β6.3, 1.80 Å maximum pore radius).49
It should be noted that in our CD samples gramicidin A′ was added at a peptide to lipid ratio of 1:
3500, and at higher concentrations, such as those used in our SAXS study, peptide–peptide interactions may play a larger role in determining secondary structure. Even so, these results show that bilayer properties are capable of altering peptide conformation, and are consistent with the theory that peptides and proteins adopt narrower conformations in response to increased local lateral pressure within the bilayer.36
The finding that the branched chain lipid phytantriol stabilises the double helix conformation agrees well with our observations of cubic phase behaviour. The double helix is significantly longer (35 Å) than the helical dimer (25 Å),49 and as such should not promote an HII phase, but rather an isotropic or bilayer structure. This is consistent with our observations that branched chain lipids exhibit small changes in phase behaviour with gramicidin A′ addition compared to MAGs.
In examining the dissimilar phase behaviour of the MAG and isoprenoid type lipids subsequent to gA′ addition, we must also look at the intrinsic curvature of the bilayer systems involved. PT and PE, with lattice parameters much smaller than the MAG lipids tested, have an increased intrinsic curvature. It should be noted that although the total mean curvature of the surface is zero, the Gaussian curvature is everywhere non-positive and varies between zero at the flat points and most negative at the saddle points. With our CD experiments we observed that in PT, the double helix conformation is present in greater amounts than in the MAGs. Previous work by Kulkarni et al. has shown that the integral membrane protein bacteriorhodopsin may preferentially select a region of the bilayer based on Gaussian curvature, and will necessarily alter its tertiary structure to accommodate the principal curvature at these points.
The spontaneous curvature of lipid systems has also been shown to influence the uptake of hydrophobic antimicrobial peptides. Using EPR spectroscopy Cafiso et al. found that the binding of alamethicin to membranes consisting of DOPC/DOPE decreased with increasing spontaneous curvature, which was accomplished by increasing the DOPE content.50 The increase in binding free energy was found to be linearly dependent on the increase in DOPE mol fraction. This effect may be responsible for the exponential decay observed in our SAXS results; as curvature increases with increasing gA′ concentration, the ability of the membrane to bind further peptide is reduced. Furthermore, the decreased effect of gA′ on PT and PE lipid membranes may also stem from decreased uptake, an effect of the increased spontaneous curvature in these lipid systems.
Through the use of circular dichroism we have shown that peptides are able to adopt different conformations depending on the lipid mesophase environment. We suggest that this is due to differences in the lateral pressure profiles of these two lipid classes. Increasing the local lateral pressure within the bilayer may result in the peptide adopting a narrower conformation, consistent with previous observations.34
These results should be of use in the planning of in meso crystallization trials, where careful lipid selection is needed to retain desired protein conformations, as well as the design of peptide and protein based drug delivery and biosensor systems incorporating bicontinuous cubic phase materials.
Footnote |
† Electronic supplementary information (ESI) available: Plots of lattice parameters for top vs. bottom well (Fig. S1); the fitting parameters resulting from fitting a first order decay function to the plots of lattice parameter (Å) vs. gA′ concentration (mol%) (Table S1); comparison of the CD spectra of gA′ encapsulated in monoolein and phytantriol with a weighted, intensity adjusted combination of the spectra of gA′ in SDS micelles and THF (Fig. S2). See DOI: 10.1039/c6ra13658c |
This journal is © The Royal Society of Chemistry 2016 |