Imola Cs.
Szigyártó
a,
Judith
Mihály
a,
András
Wacha
a,
Dóra
Bogdán
ab,
Tünde
Juhász
a,
Gergely
Kohut
ac,
Gitta
Schlosser
c,
Ferenc
Zsila
a,
Vlada
Urlacher
d,
Zoltán
Varga
a,
Ferenc
Fülöp
e,
Attila
Bóta
a,
István
Mándity
*ab and
Tamás
Beke-Somfai
*af
aInstitute of Materials and Environmental Chemistry, Research Centre for Natural Sciences, H-1117 Budapest, Hungary. E-mail: beke-somfai.tamas@ttk.mta.hu; mandity.istvan@ttk.mta.hu
bDepartment of Organic Chemistry, Faculty of Pharmacy, Semmelweis University, H-1092 Budapest, Hungary
cInstitute of Chemistry, Eötvös Loránd University, H-1117 Budapest, Hungary
dInstitute of Biochemistry, Heinrich-Heine University, 40225 Düsseldorf, Germany
eMTA-SZTE Stereochemistry Research Group, Institute of Pharmaceutical Chemistry, University of Szeged, H-6720 Szeged, Hungary
fDepartment of Chemistry and Chemical Engineering, Physical Chemistry, Chalmers University of Technology, SE-41296 Göteborg, Sweden
First published on 15th June 2020
Self-assembling peptides offer a versatile set of tools for bottom-up construction of supramolecular biomaterials. Among these compounds, non-natural peptidic foldamers experience increased focus due to their structural variability and lower sensitivity to enzymatic degradation. However, very little is known about their membrane properties and complex oligomeric assemblies – key areas for biomedical and technological applications. Here we designed short, acyclic β3-peptide sequences with alternating amino acid stereoisomers to obtain non-helical molecules having hydrophilic charged residues on one side, and hydrophobic residues on the other side, with the N-terminus preventing formation of infinite fibrils. Our results indicate that these β-peptides form small oligomers both in water and in lipid bilayers and are stabilized by intermolecular hydrogen bonds. In the presence of model membranes, they either prefer the headgroup regions or they insert between the lipid chains. Molecular dynamics (MD) simulations suggest the formation of two-layered bundles with their side chains facing opposite directions when compared in water and in model membranes. Analysis of the MD calculations showed hydrogen bonds inside each layer, however, not between the layers, indicating a dynamic assembly. Moreover, the aqueous form of these oligomers can host fluorescent probes as well as a hydrophobic molecule similarly to e.g. lipid transfer proteins. For the tested, peptides the mixed chirality pattern resulted in similar assemblies despite sequential differences. Based on this, it is hoped that the presented molecular framework will inspire similar oligomers with diverse functionality.
Accordingly, for progress towards the controlled formation of oligomers with membrane activity, here we focused on short β-peptide sequences. The obtained β3-peptides readily form oligomers which are soluble in water and in hydrophobic organic solvents such as octanol and can insert into a lipid bilayer.
Accordingly, five β3-hexapeptides (1–5, named here as “RS” set) were designed and synthesized (Fig. 1). Initially, peptide 1 was synthesized containing aromatic residues and leucines on the apolar side, whereas glutamines were chosen for the polar side. To enhance membrane activity and result discrete small oligomers, instead of infinite nanofibers,69 the N-terminal was kept more hydrophobic with the first two residues being homochiral for all sequences. The phenyl and 2-naphthyl-homoalanine residues were selected to test anchoring properties upon potential membrane interactions. Sequence 2 was designed to be less hydrophobic by introducing charged residues on the hydrophilic side but maintaining a neutral total charge with only [R]-β3-homoleucines on the apolar side. To increase water solubility, three additional sequences were synthesized (3–5). The peptide synthesis were performed using solid phase peptide synthesis (SPPS) and Fmoc strategy, in a continuous-flow apparatus, whereas peptide purity was determined and analyzed using HPLC-MS as described in the Materials and methods section of the ESI (Fig. S1–S5).†
Solubility test of the designed sequences has shown that 1, the most hydrophobic sequence, was soluble only in hexafluoroisopropanol (HFIP), whereas 2 could be dissolved in 10% acetic acid (Fig. 1). The more hydrophilic sequences 3–5 were highly water-soluble.
FTIR spectroscopy can be used for characterization of β-peptide sequences.70,71 Similarly, to IR spectra of natural proteins and peptides, the amide I and II absorption bands, originating mainly from CO stretching and N–H bending vibrations (∼1651 and ∼1540 cm−1) of the peptide backbone were considered.72–74 The amide I vibration depends on the backbone structure and is commonly used for secondary structure analysis by fitting with components bands or by decomposing into basis spectra.75
The most intense contribution of deconvoluted amide I band of 1 in HFIP is at 1648 cm−1 which could be assigned to a random coil conformation (Table S1†). For peptide 2 in 10% acetic acid solution the unordered structure was dominant too, and only a very small amount of sheet-like conformation was observed (Table S1†). The band maximum at 1646 cm−1 in the spectra of 3 in water (Fig. 2A) referred to a random coil structure. The band component at ∼1630 cm−1 could be assigned to the vibration of the guanidinium group (νs(CN3H6+)) of the arginine side chain.75 Similar structure was observed for foldamer 4 in water with small shoulders at 1628 and 1688 cm−1 which was attributed to sheet-like conformation (Fig. 2A). In the aqueous solution of 5 the random coil structure was also dominant, while an additional band at 1690 cm−1 could be assigned to turn-like conformations.76,77 Among these, the small bands ∼1567 and 1534 cm−1 originating from deprotonated glutamic acid carboxylate and protonated lysine side chains70,75,78,79 together gave the possibility of salt bridge formation. In the deconvoluted spectra of all three water-soluble foldamers (3–5) (Fig. S6A–C†) a band at 1658 cm−1 with medium intensity was observed, which was previously identified as originating from H-bond formation between the β-peptides.20 These results suggest that random coil structure predominates in “RS” peptides while there is a small amount of sheet-like structures with the presence of inter- and intramolecular H-bonds.
To further study secondary structure, CD spectroscopy was also applied for all water-soluble foldamers (3–5). However, in contrast to α-peptides, for β-peptides, helical and sheet secondary structures often produce similar spectra.80 Consequently, absolute structure determination of these new heterochiral molecules is challenging. Thus, CD spectra were used to compare foldamer conformations to each other (Fig. 2B), rather than to identify individual secondary structures. In each case, positive signals were observed below 206 nm, whereas slightly weaker negative peaks were detected at longer wavelengths. Nevertheless, some differences, between these spectral signatures should be noted, such as the intensity, shape and spectral position of the bands. The CD curve of peptide 4 exhibits similar curve pattern to that of alternating heterochiral trans-2-aminocyclopentanecarboxylic acid (trans-ACPC) hexamer which forms an elongated strand-like structure.25 Peptide 3 shows similar properties, however with lower signal intensities, suggesting a less ordered structure. The CD curve of 5 is similar to those of elongated strand-like β-peptides composed of cis-ACPC with monotonic chirality.69
To obtain further structural information, NMR measurements were also carried out for peptides 3–5. In line with expectations, the dynamic nature of these assemblies resulted in low signal resolution of the 1H NMR spectra (Fig. S8†), which is also a typical property for non-helically ordered structures.25,69 The complete signal assignation of the backbone protons using 2D NMR (Fig. S9 and S10†) measurements cannot be made unequivocally due to the similar chemical shift of the second and fourth amino acids. Thus, we performed DOSY experiments for these peptides (Fig. S11–S13†) which proved more informative. The results revealed the formation of small assemblies in aqueous solution. The average number of molecules in the assemblies were determined as detailed previously20 using 1,4-dioxane as an internal standard.81,82 Based on molar mass calculation of peptides 3–5 it results in oligomers with low number of peptide molecules, with estimated tetrameric to decameric composition (Fig. S11–S13†).
The morphology and the size of the assemblies were investigated by using transmission electron microscopy (TEM) and dynamic light scattering (DLS) techniques. The freshly dissolved peptide 5 showed the formation of rod-or needle-like assemblies in the range of 60–90 nm size (Fig. 3). Detailed investigation of these particles shows that this morphology is formed from smaller individual bundles of ∼8–10 nm in size. Assemblies similar to those of peptide 5 were reported earlier for longer polyglutamine sequences with β-hairpin and for peptides of amyloid β.83–85 This not only suggests capacity of the method to identify these oligomeric species, but also that combination of residues with hydrophobic and hydrogen bonding character can result similar discrete assemblies for otherwise distinct peptidic compounds. DLS measurements also support the size range of formations observed on the TEM images (Fig. S14†).
To support the experimental observations and provide molecular level insight, molecular dynamics (MD) simulations were performed on self-assembly of peptides 3–5. For this we employed our recently developed parameters for β-peptide backbone dihedrals, which enables compounds to explore non-standard conformations as well.86 Besides simulations on the monomer state, spontaneous assembly of eight individual peptides was also tested in explicit water where the computations started with eight individual molecules placed in random orientation at the corners of a cube of 4.5 nm edge length (for more details see the Materials and methods section of the ESI†). Interestingly, after very short simulation time (∼100 ns), the peptide molecules quickly assembled, forming a sandwich-like octameric bundle (Fig. 4). Similar assemblies for peptides 3 and 4 were also obtained (Fig. S15 and S16†). These oligomers consist of two nearly parallel layers (I and II), each of them built from peptides aligned parallel, where the layers are stabilized separately with intermolecular hydrogen bonds. The alignment of the peptide backbones in the two layers is nearly orthogonal (Fig. 4A). The two layers are formed by keeping the Leu side chains in the inner hydrophobic core, while the external ones are hydrophilic with Lys and Glu side chains. The two layers face each other with the hydrophobic side inwards, “protecting” the Leu side chains from the polar solvent (Fig. 4B and C).
The secondary structure of the peptides in the bundle can be assessed by conformational clustering. For analysis, the part of the trajectory corresponding to the aggregated state has been split into eight parts after 400 ns, corresponding to the individual peptides. These sub-trajectories were concatenated, giving a set of frames sampling the secondary structures of all peptides. Cluster analysis was performed87 which identified the most common fold as an elongated conformation termed as zig-zag conformation,59,88,89 accounting for nearly 53% population of the conformations explored throughout the simulation (Fig. S17 and S18†).
The charged Lys and Glu side chains can form salt bridges (Fig. S19 and Table S2†) which are also supported by FTIR results. There are no salt bridges found between chains on opposite (I and II) layers of the aggregate, nonetheless, besides intramolecular salt bridges, almost all neighboring chains in the same layer form intermolecular Lys–Glu side chain hydrogen bonds as well. This suggests that the separate layers are rather stable, with ∼10 intermolecular peptide H-bonds and several intermolecular salt bridges per layer formed between the Lys–Glu side chains.
Fig. 5 DSC thermograms of pure DPPC and DPPC/peptide (3–5) systems collected at 68 mM total lipid concentration and the peptide:lipid molar ratio was 1/10. |
To better understand β-peptide–lipid interactions, we investigated the two major regions of the lipid bilayer using IR spectroscopy. The band position of the symmetric (νsCH2) at 2853 cm−1 and antisymmetric (νasCH2) at 2922 cm−1 methylene stretching vibrations are a sensitive parameter to lipid conformational order and acyl chain packing.75,95 No significant changes were observed in the CH2 region of DOPC and DOPC/DOPG liposomes when peptides 1 and 2 were added (data not shown), suggesting that these peptides do not perturb the hydrophobic lipid chains in the membrane. The phosphate group vibrations with antisymmetric (νasPO2−) and symmetric (νsPO2−) stretching modes are at around 1240–1230 cm−1 and 1087 cm−1, respectively. The R–O–P–O–R′ stretching vibration generally appears as a shoulder near the νsPO2− band at ∼1066 cm−1. For 1 and 2, νasPO2−, the most sensitive band for H-bonding interactions, showed a shift towards lower wavenumbers compared to the pure lipid systems (Fig. S7A†). This indicates that peptides 1 and 2 are situated on the surface of the lipid bilayer, perturbing the primary hydration shell of the lipid molecules, and inducing a conformational change of the lipid headgroup moiety. The amide I band, tracked regularly for natural peptides,96–98 was also analyzed (Fig. S7B†). In the presence of DOPC liposomes, the main band around 1648 cm−1 assigned to the random coil structure remained dominant for 1 and 2 and a new band component appeared around 1688 cm−1. More pronounced changes were observed with the negatively charged DOPC/DOPG liposomes. The spectrum showed a more ordered structure with a pronounced band at 1688 cm−1 and a small one at 1629 cm−1, which could be assigned to a sheet-like conformation. The most intense contribution of amide I band of 2 in the presence of DOPC/DOPG liposomes is shifted to 1653 cm−1 and two shoulders could be observed at 1671 and 1638 cm−1, which may again suggest sheet-like conformations.
Based on IR spectra of peptides 1 and 2, the more water-soluble RS compounds (3–5) were investigated only in the presence of DOPC/DOPG liposomes. Interestingly, for peptide 5 a shift was observed for the acyl chain vibrations in both the antisymmetric (νasCH2) (from 2922.1 to 2924.1 cm−1) and the symmetric (νsCH2) (from 2852.7 to 2853.7 cm−1) methylene stretching vibrations compared to that of pure lipids (Fig. 6A). The extent of these shifts agrees well in magnitude with those observed for other membrane inserted peptides.99 These results suggest that peptide 5 interacts with the lipid acyl chain region by inserting into the lipid bilayer, which cannot be observed for peptides 3 and 4. For the latter two, the antisymmetric and symmetric PO2− stretching vibrations showed a shift towards lower wavenumbers (from 1243 to 1233 cm−1 and from 1091 to 1087 cm−1), which refers to the “on surface” localization of these peptides (Fig. 6B). The intensity of the phosphate-diester stretching band (R–O–P–O–R′) also decreased as observed for 1 and 2 in liposome systems (Fig. S7A†). Regarding the amide I band of peptide 3 (Fig. 6C) in the presence of DOPC/DOPG liposomes a broadening of the band was observed with a small shoulder at 1620 cm−1 which may refer to some sheet-like secondary structure. A similar broad band was also seen for 4 in the presence of liposomes with a shoulder around 1665 cm−1 which may be allocated to intramolecular hydrogen bonding. When comparing the water-soluble RS compounds, the highest change in peptide secondary structure was presented by peptide 5 upon interaction with the lipid bilayer (Fig. S6D–F†). The relative amount of random coil structure in the peptide 5–liposome system decreased, while the appearance of the band component at ∼1675 cm−1 strongly indicates an increase in intermolecular H-bonding.
In overall, these observations from FTIR spectra confirmed that the peptide 5 with increased amount of intermolecular H-bonds is situated below the head group region. This effect was observed neither for peptides 1 and 2, nor for 3 and 4, as these were located mainly near to or embedded into the lipid head group regions.
Further on, the presence of DOPC/DOPG liposomes caused a substantial blue shift in the CD spectrum of peptide 3 (Fig. 6D). Both the zero-crossing point and the λmin were displaced to lower wavelengths by 6 and 4–5 nm, respectively, suggesting a stronger effect of the membrane on peptide conformation. The negative band of 4 in the presence of liposomes totally overlaps with that obtained in water and its minimum shows only a small, 2 nm red shift (Fig. 6D). Such a slight alteration of the CD profile does not indicate gross conformational change in the peptide. When compared to CD spectra of peptide 5 in water, the zero-crossing point remains at 203 nm also in the lipid system, whereas the positive band exhibits a moderate red shift. The negative band shows two minima at 210 and 217 nm, which cannot be seen in aqueous solution. These changes may correlate with the increase of strand-like conformation content induced by the membrane interaction. However, the lower signal intensity due to the presence of liposomes makes direct comparison of the aqueous and membrane phases difficult. Thus, conformational variations of foldamer 5 were further investigated in solvents with different polarity (Fig. S23†). The minimum of the negative band detected at 215 nm in water is blue shifted to 207 nm in methanol, and the same shift can be observed for the zero-crossing point as well. On the other hand, this blue shift in octanol is smaller, resulting in a shallow negative minimum at ∼212 nm. Octanol is widely used to mimic membrane environment and based on the slight difference observed in the shape of CD curves for water and octanol, it is proposed that the peptide conformations in these solvents lie close to each other. The most significant shift in methanol suggests larger structural changes, which may indicate appearance of monomers. To support this, MD calculations were also performed using methanol as solvent. Eight individual peptide molecules were placed in random orientation using the same setting parameters as in water simulations (for more details see the Materials and methods section of the ESI†). In contrast to the aqueous phase, oligomer formation did not occur in methanol over the first μs of the simulation time. Note that formation of dynamic assemblies in water occurred already at the first 100 ns of the simulation time. To further confirm monomer state, the simulation was continued until a full 2 μs of simulation time, under which no assembly formation was observed as demonstrated by the lower number of intermolecular hydrogen bonds (Fig. S20†).
To gain molecular level insight to membrane interactions, 5 was subjected to further MD simulations in the presence of a lipid bilayer. The aim of these simulations was to investigate how peptide 5 could be embedded in the lipid chain region as suggested by experimental results detailed above. Accordingly, simulations were started with peptide 5 centered into the middle of the bilayer. The applied atomic restraints during the minimization and equilibration process are summarized in Table S3.† To estimate the optimal extent of self-assembly, three parallel simulations were performed, with a monomeric, tetrameric and an octameric form as starting structure. Considering the hydrophobic surrounding, the oligomeric form was assumed to be similar to the water-soluble two-layered assemblies but having the leucine on the exterior side and the charged residues on the interior side. For the monomeric and tetrameric state a rapid, spontaneous reorientation occurred in the first ∼100 ns (Fig. 7A), after which both forms are situated mainly in the vicinity of the bilayer surface, where they are stabilized over the remaining simulation time (900 ns). In contrast, the octameric form remained in the hydrophobic interior of the lipid acyl chains (Fig. 7A), with very little variation throughout the entire simulation time. This foldamer oligomer also retains connection with both ends of the lipid bilayer, by having several neighboring lipids bent from the original membrane surface. An overall thinning of the membrane from both sides is caused in the vicinity of the oligomer (Fig. 7B). The experimental and simulation limitations prevent the quantification of the number of peptides present in aqueous or lipid bilayer inserted oligomers. Nevertheless, based on the above it is concluded that an octameric association of peptide 5 is already large enough to properly shield the charged residues in the depth of the lipid bilayers. The overall structure resembles to the water-soluble one, namely it is also a two-layered assembly with several H-bonds between the individual peptides in each layer. Conformational analysis of the amino acid residues indicates similar torsional angle distributions for the aqueous and bilayer forms, though with a somewhat altered ratio (Fig. S18†). However, in contrast to the water soluble form, the angle between the plane of layer I and layer II shows much less variation in the lipid bilayer, its value stays ∼60°–80° throughout the entire simulation time (Fig. S21†). The lipid bilayer structure also contains numerous salt bridges formed by the lysine and glutamate side chains in the interior of the assembly. Note, that when comparing relative average positions of the monomeric, tetrameric and octameric forms, a qualitative mechanism on membrane insertion can be outlined (Fig. 7A). Accordingly, as the number of peptides in the oligomers increase, the position of peptide 5 gets closer to the center of the lipid bilayer. In the case of the octameric form, the observed stability in the depth of the membrane will also be supported by the apolar butyryl groups on the N-terminal of peptide 5. Furthermore, the large, ∼60°–80° angle between the two layers allows for such an orientation where all apolar N-terminals are placed towards the hydrophobic core of the bilayer (Fig. 7B). The presence of oligomers inside the lipid bilayer is also in line with the sequential build-up of the molecule. Since the alternating sequence prevents helix formation, the charged residues would fall on the same side in a more elongated conformation of a monomer, where Glu and Lys are likely to form a salt bridge in the hydrophobic membrane environment. However, a single Glu–Lys salt bridge will be still a rather polar region and the monomer state cannot shield that from the lipid bilayer. Thus, based on the IR and DSC results, the MD simulations and the sequence of peptide 5, it can be concluded that monomeric form inside the bilayer is unlikely.
In aqueous solutions, ANS, a well-known hydrophobic probe, shows a weak fluorescence signal with an emission maximum at ∼530 nm, which is blue shifted in hydrophobic environment to below 500 nm.105 Model membranes or micelles with net positive or neutral head-groups bind ANS buried at the hydrophobic–hydrophilic edge, as reflected by the emission maximum of 480 nm for the PC liposomes (Fig. 8A). Similarly, proteins with exposed hydrophobic patches with a nearby positively charged moiety where ANS fits well, exhibit remarkably blue-shifted emission maxima, too. When tested with the assemblies of 3–5, a dose-dependent shift of the maximum in the 465–470 nm range were detected upon consecutive addition of the peptides. This suggests the presence of an optimal ANS binding site where the probe is exposed to the solvent/buried from the solvent to a comparable extent, as in the case of lipid bilayers or suitable protein patches.
In contrast to the minor changes in emission maximum, larger differences were detected for the signal intensities as observed for peptide 5 in comparison to peptides 3 and 4 (Fig. 8B). Considering the sequential differences of peptides, this can presumably be attributed to the replacement of Arg in 3 and 4 to Lys in 5. This change might manifest in the effect observed as the positively charged peptide side chains could contribute to the binding via electrostatic attraction to the ANS sulfonyl group, or via cation–π interaction formed between cationic side chains and ANS aromatic rings.106 Alternatively, the higher capacity of Arg over Lys to form multiple thus stronger H-bond networks might result in higher oligomerization tendency, leading to the formation of bigger, less structured aggregates with enhanced ANS binding properties. In addition, the oligomeric state might also affect ANS binding, resulting in intensity variations, supported by the enhanced signal background arising from light scattering of larger particles observed in TEM images. Moreover, peptide assemblies and proteins in partially unfolded states, decorated with exposed hydrophobic patches, are characterized by remarkably enhanced ANS fluorescence, while well-folded structure with a more buried hydrophobic interior exhibit less pronounced ANS signals.107 Based on these considerations, the high signal intensities detected for 3 and 4 are likely indicative of a less folded, loosely packed overall structure, whereas the relatively low ANS binding intensity observed for 5 suggests more organized, folded associates, in agreement with the IR and TEM results. Considering the signal enhancement, titration data showed a nearly linear dose-dependence on foldamer concentration (Fig. 8B). This data may also suggest a critical association concentration (cac) below 20 μM as ANS binding was detected even at concentrations as low as 20–50 μM, indicating presence of foldamer assemblies capable of effectively hosting guest molecules at low peptide levels.
Host–guest interaction was further tested using the probe ThT, which is sensitive for the rigidity of the microenvironment.100 ThT alone also exhibited self-fluorescence,108 the intensity of which increased upon addition of peptides 3–5 (Fig. 8C). In consideration with the structure and morphology of the assemblies derived from our MD simulations and TEM images, ThT could bind to the side cavity formed between the two layers. This arrangement could provide a favourable surface for the dye to align in an extended conformation, which is the supposed prerequisite for the induced fluorescence. Further on, this is in line with the preferred ThT binding site, which was concluded earlier to be enriched in hydrophobic and aromatic side chains, but lacking many charged groups.100 Analysing the ThT signal, the induced intensity was the lowest for 4 and the highest for 3 which can be explained in terms of slightly different morphology of the assemblies, or sequential environment. For amyloid structures, the induced signal was proven to be proportional to the concentration, allowing fibril quantification based on the intensity, at least up to 10 μM amyloid.108 In comparison, the dose dependence curves show saturation type behaviour at higher foldamer levels, indicating a dissociation constant of the interaction in the micromolar range (Fig. 8C). Altogether, fluorescence-based findings suggest similar host capacity for 3–5 indicating the formation of assemblies. Fluorescent probes differing in hydrophobicity and charge reported variations in binding capacities, which are consistent with sequential and spatial preferences of the particular dye. The rather linear dependence of the probe interactions indicates no critical association concentration for the foldamers in the 20–250 μM range suggesting that the observed assemblies are formed readily at low concentrations.
Finally, to test these assemblies toward potential applicability as hosts or transporters of hydrophobic biomolecules or drugs, we employed peptide 5 and the planar small molecule, pyrene, which readily inserts into lipid membranes due to its hydrophobic nature.109–111 UV-Vis measurements indicated that the pyrene molecule is inserted into the foldamer assembly similarly as it inserts into lipid bilayers (Fig. 9A). The broad band with a maximum at 360 nm assigned to pyrene aggregates present in water,112 almost vanished upon titration with peptide 5 and the shape of the spectra in the lower (300–350 nm) wavelength range resembled that of monomeric pyrene inserted in a lipid bilayer.113,114 The reduction of the absorption band intensity at 360 nm has already initiated with the very first addition of 5 with as low as 5 μM concentration (Fig. 9A inset). This observation suggests that oligomerization of compound 5 takes place already in the low micromolar range. This is also supported by the signal enhancement results during ANS titration for 5 (Fig. 8B). Thus, qualitatively both measurements point towards a critical association concentration (cac) to be in this range, estimated to be ∼5–20 μM. However, taking into consideration that our peptides show some similarity to sheet-like amyloid peptide systems, which are known to have complex oligomerization processes,115–117 we cannot exclude that multiple aggregation states, starting potentially well below micromolar concentrations, occur also for the currently presented compounds.
Insertion of pyrene into the bundles formed by 5 has also been investigated by MD simulations. A single pyrene molecule was simulated with the solvated octameric oligomer of 5. According to the simulation results, the pyrene molecule interacts closely with peptide 5 in most of the total 1 μs simulation time. For a more in-depth analysis, the oligomer was sectioned into four regions (as detailed in the ESI†) and the position of the guest molecule was followed by measuring pair distances with each section. Accordingly, the fully hydrophobic pyrene is located either in the internal leucine core, providing an optimal hydrophobic environment, or near the core but closer to the butyryl groups of the N-terminus, partially exposed to the solution (Fig. 9B and S22†). In this way a small bundle is capable of hosting water-insoluble compound and no significant change in the oligomeric state and peptide conformation was found during the simulation. Only some re-ordering of the chains was observed whereas the number of intermolecular hydrogen bonds remained nearly the same. These suggest that the two-layered assemblies can be sufficiently dynamic to allow rearrangements upon incoming guest molecules and decrease or negate potentially arising steric strains within the assembly.
An extended Method section is reported in ESI.†
Footnote |
† Electronic supplementary information (ESI) available: Experimental procedures and data: peptide purification and characterization, FTIR, NMR, MS, DLS, MD, DSC, TEM and CD measurements. See DOI: 10.1039/d0sc01344g |
This journal is © The Royal Society of Chemistry 2020 |