Takuya
Suguri
ab and
Bradley D.
Olsen
*a
aDepartment of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA 02139, USA. E-mail: bdolsen@mit.edu; Tel: +1 (617) 715-4548
bYokkaichi Research Center, JSR Corporation, 100 Kawajiricho, Yokkaichi, Mie, Japan
First published on 5th March 2019
Bioconjugates made of the model red fluorescent protein mCherry and synthetic polymer blocks show that topology, i.e. the BA, BA2, ABA and ABC chain structure of the block copolymers, where B represents the protein and A and C represent polymers, has a significant effect on ordering transitions and the type and size of nanostructures formed during microphase separation. ABA and ABC type block copolymers were synthesized by using two site-specific bioconjugation reactions: the thiol–ene reaction with a cysteine on mCherry and maleimide functionalized polymers, and the sortase A ligation reaction with an LPETG sequence at the C-terminus on mCherry and a triglycine functionalized polymer. The phase behaviors of mCherry–poly(N-isopropylacrylamide) (PNIPAM) and mCherry–(PNIPAM)2 show that the shapes of the phase diagrams are similar overall, but mCherry–(PNIPAM)2, i.e. BA2 type, yields a narrower domain spacing than mCherry–PNIPAM, i.e. BA type. PNIPAM–mCherry–PNIPAM (ABA type) shows only lamellar phases in the range of conditions under which ordered structures appear. PDMAPS–mCherry–PNIPAM (ABC type) shows an ordered structure across the widest range of conditions in the four bioconjugates and also the widest range of different nanodomain structures. The phase behavior of the ABC type implies that the repulsive interaction between two water-soluble coil polymers can be a key factor in enhancing the self-assembly of globular protein–polymer block copolymers.
Because the self-assembly of protein–polymer block copolymers is very different from traditional block copolymer systems, understanding the physics of this self-assembly process remains a challenge.25–27 One important parameter that has an impact on block copolymer self-assembly is the topology of the block, for example changes between BA, ABA, and BA2 copolymers. Although there are many reports about topology effects on polymer–polymer block copolymers including branch and triblock types of coil–coil28–32 and rod–coil block copolymers,33–38 the effect on protein–polymer block copolymers has not been clarified. Recently, we have reported the topological effects on the self-assembly of globular protein–coil protein fusion proteins22,39 and protein polymer block copolymers.40 In this work, changing the elastin-like polypeptide (ELP) tails from the N to C termini of mCherry has no impact on the ability to form well-ordered nanostructures.39 On the other hand, a bola fusion protein, which has two ELP tails on both N and C termini of mCherry, showed a lower order–disorder transition concentration than the linear fusion protein, which has one EPL tail on the N terminus of mCherry.22 This work provided important new insight into the effect of chain topology in the self-assembly of block copolymers including globular proteins. However, fusion proteins have some limitations of topological and functional designs because globular and coil proteins can be connected to each other only at the N and C termini, and only linear primary chain structure is possible. Furthermore, new genes are necessary for each design to modify or introduce functions derived from the coil protein, such as stimuli responsiveness. In addition to these fusions, changing the bioconjugation site in mCherry–PNIPAM bioconjugates has been explored. Changing the bioconjugation site showed qualitatively similar phase diagrams, indicating that self-assembly is robust with respect to changes in the conjugation site. However, differences in domain spacing were observed, suggesting changes in protein orientation within the lamellar phase.40
The protein–polymer block copolymer approach is a promising approach to expand the capability to obtain topological designs and functional designs. To obtain well-defined bioconjugates, site-specific modification is one of the most important tools. To control topologies of polymers, living polymerizations such as RAFT polymerization are also necessary.41 Two or more orthogonal site-specific reactions provide the ability to access a wider variety of bioconjugates; however, this remains challenging.42
In this study, the effect of topology on the self-assembly of protein–polymer bioconjugates is investigated by comparing four types of protein–polymer bioconjugates, namely BA, BA2, ABA and ABC types, where B represents the protein and A and C represent polymers. In order to examine the effect of branching architecture, 1 and 2-arm PNIPAMs are respectively introduced on mCherry, i.e. mCherry–PNIPAM and mCherry–(PNIPAM)2. In addition, in order to examine the effect of ABA and ABC triblock types, PNIPAM–mCherry–PNIPAM and PDMAPS–mCherry–PNIPAM are also studied. The use of multiple conjugations necessitated the use of two orthogonal site-specific bioconjugations, motivating the development of RAFT polymers functionalized for sortase ligation. Finally, the microphase separation behaviour of each bioconjugate is investigated using small-angle X-ray scattering (SAXS) and depolarized light scattering (DPLS).
To obtain BA and BA2 type copolymers, a site-specific thiol–maleimide coupling between a cysteine on the mutated mCherry (mCherryS131C) and maleimide–PNIPAM was used (Scheme 1). To obtain one and two arm PNIPAMs, two types of reversible addition fragmentation chain transfer (RAFT) agents containing a protected maleimide group (pMI-EMP and pMI-(EMP)2) were used. After polymerization, the maleimide groups were thermally deprotected to yield maleimide end-functionalized 1-arm and 2-arm PNIPAMs. Each conjugate was designed to have an approximately equal volume fraction of mCherry and PNIPAM, because previous results suggested that an approximately equal volume fraction of each polymer would enhance self-assembly and show a minimum order–disorder transition concentration.46 Bioconjugation reactions were performed under mild conditions similar to literature procedures.21
For ABA and ABC copolymers, two distinct bioconjugation sites are necessary, and ABC type copolymers need two orthogonal site-specific chemistries to attach the two different polymers. Cysteine is selected as the first reacting group for thiol–ene reactions, which is on the one side of the β-barrel structure; the mutation site is the same as mCherryS131C. The second group is the LPETG sequence, installed at the C-terminus on the other side of the β-barrel structure for Sortase A mediated ligation (SrtA ligation). Sortase A from Staphylococcus aureus catalyzes the cleavage reaction between threonine and glycine at an LPXTG recognition sequence on the first molecule, where the variable residue X is most commonly E, followed by the formation of an amide bond between the carboxyl group of threonine and a second molecule with an N-terminal nucleophilic group, typically an oligoglycine motif. SrtA ligation was selected for the second bioconjugation reaction because of its high specificity, mild reaction conditions, ease of introducing the LPETG sequence into recombinant proteins, and compatibility with synthetic polymers. In addition, a histidine tag (His-tag) could be introduced in mCherry downstream from the LPETG sequence to facilitate purification. To introduce N-terminal oligoglycine in polymers, a new RAFT agent with the tert-butyloxycarbonyl (Boc) protected triglycine (Boc-GGG-EMP) was synthesized (Scheme 2). HO-EMP was synthesized according to previously reported procedure (details in ESI†). HO-EMP was reacted with Boc-GGG-OH in the presence of N-(3-dimethylaminopropyl)-N′-ethylcarbodiimide hydrochloride and 4-dimethylaminopyridine to give Boc-GGG-EMP. The Boc moiety was removed after RAFT polymerization.
ABA and ABC type triblock copolymers were synthesized in two steps (Scheme 3). First, PNIPAM and PDMAPS with maleimide end groups were reacted with a cysteine on mCherryS131C–LPETGG–His6 to give PNIPAM–mCherry–LPETGG–His6 and PDMAPS–mCherry–LPETGG–His6. After both bioconjugates were purified, they were reacted with triglycine–PNIPAM in the presence of sortase A. Both polymer chains for each triblock were designed to have approximately equal volume fraction in the final block copolymer. Unreacted bioconjugates (PNIPAM–mCherry–LPETGG–His6 and PDMAPS–mCherry–LPETGG–His6) and sortase A which has a His-tag were removed using Ni-NTA. Bioconjugates (PNIPAM–mCherry–PNIPAM and PDMAPS–mCherry–PNIPAM) and excess polymer (GGG-PNIPAM) were collected as flow-through. Then, FPLC purification was used to isolate PNIPAM–mCherry–PNIPAM (ABA) and PDMAPS–mCherry–PNIPAM (ABC), respectively. The mCherry function after purification of the triblock type was relatively lower than that of the diblock type: BA = 94%, BA2 = 90%, ABA = 74%, and ABC = 85% based on UV-vis absorbance measurements (Fig. S17, S20 and S22†). These differences suggest that sortase ligation or the longer two step conjugation and purification process results in some small loss of protein function. The SDS-PAGE shows that the bioconjugates are pure, indicating that the self-assembly behavior is due to the bioconjugates, not due to a mix of native or partially modified and fully modified proteins. The molecular properties of all four protein–polymer block copolymers are summarized in Table 1.
![]() | ||
Scheme 3 Synthesis of (a) PNIPAM–mCherry–PNIPAM and (b) PDMAPS–mCherry–PNIPAM by thiol–ene bioconjugation reaction and sortase A ligation reaction. |
Bioconjugates | Type | mCherry | Maleimide–PNIPAM | Maleimide–PDMAPS | Triglycine–PNIPAM | Polymer volume fractiona | |||||
---|---|---|---|---|---|---|---|---|---|---|---|
Reaction site | M n | PDI | M n | PDI | M n | PDI | ϕ maleimide | ϕ triglycine | ϕ total | ||
a Density of PNIPAM 1.05, PDMAPS 1.37 and mCherry 1.35 g cm−3 and a molecular weight of mCherry 28.2 kDa were used to calculate the volume fraction of polymers in bioconjugates. | |||||||||||
mCherry–PNIPAM | BA | Cys | 29.1 | 1.16 | — | — | — | — | 0.606 | — | 0.606 |
mCherry–(PNIPAM)2 | BA2 | Cys | 28.6 | 1.15 | — | — | — | — | 0.599 | — | 0.599 |
PNIPAM–mCherry–PNIPAM | ABA | Cys/LPETG | 17.8 | 1.05 | — | — | 17.4 | 1.07 | 0.313 | 0.318 | 0.631 |
PDMAPS–mCherry–PNIPAM | ABC | Cys/LPETG | — | — | 19.2 | 1.20 | 17.4 | 1.07 | 0.303 | 0.322 | 0.625 |
Comparison of mCherry–PNIPAM (BA type) and mCherry–(PNIPAM)2 (BA2 type) block copolymer self-assemblies in concentrated solution shows that both materials exhibit a similar phase behavior. Phase identifications were performed using small-angle X-ray scattering (SAXS) and depolarized light scattering (DPLS), and each identified phase is summarized in phase diagrams (Fig. 1a and b). The phase diagrams of both mCherry–PNIPAM and mCherry–(PNIPAM)2 showed a quite similar general shape (Fig. 1c and d). In both bioconjugates, disordered phases are observed at concentrations below 40 wt% at low temperature. Increasing temperature above 35 °C yields a disordered micellar phase due to the LCST behavior of PNIPAM46 which causes desolvation of this block in the copolymer upon heating. Increasing concentration to 40 wt% yields a lamellar phase with peaks at q* and 2q*. At 50 wt%, both bioconjugates show well-ordered lamellar phases, and the peaks broaden with further increasing concentration to 60 wt%. This suggests that the degree of order decreases. Only in mCherry–PNIPAM, lamellar–lamellar coexistence is observed at 45 and 50 wt% (see also Fig. S25†). Increasing the temperature at concentrations above 40 wt% yielded a hexagonal phase with peaks at q*, √3q* and √4q* for a narrow range of temperature and concentration. In addition, in both bioconjugates, birefringent lamellae are observed at temperatures below 25 °C (mCherry–PNIPAM) or 30 °C (mCherry–(PNIPAM)2). In the case of mCherry–(PNIPAM)2, a birefringent lamellar phase is observed from 50 wt% to 60 wt%, while mCherry–PNIPAM shows birefringent lamellae from 40 wt% to 60 wt%. While these results might imply that the BA2 type shows relatively weaker ordered structure than the BA type, it has also previously been shown that liquid crystalline ordering can emerge in protein–polymer conjugates at high concentrations,43 and it is possible that the propensity for such ordering is stronger in the linear BA conjugate.
In a previous comparison of bola and linear block copolymers using mCherry and ELP fusion proteins,22 the bola type which has two ELP tails showed a much wider concentration range for ordering than the linear fusion copolymer. While the linear mCherry–ELP fusion protein self-assembled over a similar concentration range to the linear mCherry–PNIPAM bioconjugate (above 40%), ordered structures were observed for the bola fusion for concentrations above 30 wt%. Comparing the BA type copolymer in this work and linear fusion in the previous work, both block copolymers have similar structures, which are globular protein–coil block copolymers, despite differences in the chemistry of the coil block. Furthermore, both phase diagrams suggest a similar phase behavior. In other words, these results imply that the difference between PNIPAM and ELP has only a minor impact on the phase behavior of these two polymers. In contrast, BA2 type and bola type show different behaviors for the ELP. Because ELP and PNIPAM behave similarly in the linear system, this difference in the effect of topology between PNIPAM and ELP may arise from the disparities in the way the individual chains are connected to mCherry: there is a single branch point for PNIPAM, while the ELP is connected in two different locations on the same side of the protein.
A theoretical study of the BAn Gaussian coil block copolymer systems50,51 and rod–coil block copolymer systems including coil branching52,53 predicts that phases could be shifted by the molecular architectural change; however, this is not observed with these protein copolymers. Recent experimental results show a topology effect on self-assembly of polyhedral oligomeric silsesquioxane (POSS)–polystyrene (PS) conjugates, which suggests that introduction of two PS tails with half-length of one tail affects the boundaries between different ordered phases, order–disorder transition temperatures and packing configurations of the functional POSS cages.54 Although both BA and BA2 in this study show similar boundaries and order–disorder transition temperatures, the difference between topologies affected the packing configurations, namely domain spacings and cross-sectional areas.
The domain spacings of both mCherry–PNIPAM and mCherry–(PNIAPM)2 vary with temperature and concentration over the range of 24.8–30.8 nm in the case of mCherry–PNIPAM and 20.2–24.2 nm in the case of mCherry–(PNIAPM)2 (Fig. 2a and b). In mCherry–PNIPAM at 50 wt% and 10 °C, the domain spacing is 27.4 nm, which is nearly two times as large as the calculated domain spacing (ESI†).
![]() | ||
Fig. 2 Domain spacing calculated by primary peaks of SAXS for (a) mCherry–PNIPAM and (b) mCherry–(PNIPAM)2 and (c) cartoon schematics of lamellar phase. |
This value suggests bilayer lamellar morphologies, as previously reported.45 In mCherry–(PNIPAM)2, the smaller domain spacing (e.g. 21.7 nm at 50 wt% 10 °C) suggests that branched architectures result in narrower polymer domain morphologies (Fig. 2c) because of the lower hydrodynamic volumes, as observed in other studies.37,38 For this to occur while conserving the volume of each component, there must be a larger average lateral spacing between proteins. To estimate lateral spacing, cross-section area54 can be calculated using domain spacings based on SAXS curves (ESI†). The cross-section area corresponding one bioconjugate of mCherry–PNIPAM and mCherry–(PNIAPM)2 at 50 wt%, 10 °C is calculated to be 6.43 nm2 and 7.97 nm2 based upon SAXS curves, which suggest that diameters of them are 2.90 nm and 3.22 nm, respectively. These estimated diameters are reasonable considering the previous study54 and the calculated diameter at 50 wt% based upon mCherry β-barrel structure of 3.15 nm (ESI, Fig. S29†).
Both linear and branched copolymers show a similar dependence of the domain spacing on temperature and concentration. For temperatures below 30 °C, the nanostructures at 40 wt% show the smallest domain spacing; at higher temperatures, there is little dependence of the domain spacing on concentration. As temperature increases, domain spacing initially decreases as the solvent quality for the PNIPAM block grows worse and the PNIPAM chain contracts, and then increases when the solvent changes from good to poor for PNIPAM, presumably due to an increased strength of segregation. Over the range of conditions, mCherry–(PNIAPM)2 shows smaller domain spacings than mCherry–PNIPAM.
Compared to BA and BA2 diblock copolymers, ABA triblocks show a higher critical concentration for ordering. PNIPAM–mCherry–PNIPAM undergoes order–disorder transitions at 50 wt%, while diblock type bioconjugates undergo transitions at 40 wt% (Fig. 3a and b). The SAXS peaks of PNIPAM–mCherry–PNIPAM at concentrations above 50 wt% are broader than for diblock types. At concentrations below 45 wt%, increasing temperature yields a disordered micellar phase, which is similar to the behavior observed in BA and BA2 type copolymers. On the other hand, in the range of ordered concentration and temperature, only lamellae are observed, unlike mCherry–PNIPAM and mCherry–(PNIAPM)2 which also show hexagonal phases. In addition, unlike BA and BA2 type, PNIPAM–mCherry–PNIPAM never shows birefringent phases. Theoretical studies of ABA type Gaussian coil block copolymers55 and comparison between rod–coil and coil–rod–coil56 provide insights into the self-assembly behaviors of our bioconjugates, which suggested that the lamellar phase of a coil–rod–coil becomes more unstable than a coil–rod. Since mCherry has a large diameter compared to the rod, it is not straightforward to explain this self-assembly behavior using previous theory, but our results suggest similar behaviors.
![]() | ||
Fig. 3 Representative pattern (a), phase diagram (b), domain spacing (c), and cartoon schematics (d) of PNIPAM–mCherry–PNIPAM. |
The relatively small domain spacings for PNIPAM–mCherry–PNIPAM change only slightly with varying concentration and temperature in the range 20.2–21.7 nm (Fig. 3c). The domain spacing of the lamellar morphology at 55 wt%, 10 °C is 20.2 nm, which is nearly 1.4 times larger than the calculated domain spacing for a rigid block copolymer. This value suggests a monolayer lamellar phase, unlike a bilayer lamellar phase in the case of the diblock type (Fig. 3d). This is consistent with the fact that the mCherry block has PNIPAM conjugated to both ends, preventing bilayer association.
In contrast to PNIPAM–mCherry–PNIPAM, PDMAPS–mCherry–PNIPAM shows ordered structures at the lowest concentration of any of the four copolymers studied here. This ABC copolymer shows ordered structure at all concentrations studied above 30% (Fig. 4a and b). In the concentration range of 30–35 wt%, a lamellar morphology was observed at 10–25 °C. Increasing temperature in the concentration range of 30–35 wt% yields phase transitions to disordered and disordered micellar phases at temperatures above 35 °C due to the LCST behavior of PNIPAM. Increasing concentration to 40 wt%, hexagonal and lamellar ordering within a perforated lamellar structure are observed. The perforated lamellar phase is characterized by SAXS peaks at (100), (001), (110), (002), and (003) (e.g. 45 wt%, 10 °C) or (100), (001), (002), (210), (300), (310) and (003) (e.g. 50 wt%, 20 °C) (Fig. S26(b)†). Increasing concentration from 40 wt% to 60 wt% yields phase transitions from perforated lamellar to lamellar to hexagonal phases. Over the concentration range of 40–50 wt%, increasing temperature yields phase transitions from lamella to perforated lamellar to hexagonal to cubic phases including co-existence regions. Changes in the experimental conditions change the form factor of the cylinder nanodomains, resulting in changes in the relative peak intensities for the hexagonal phase: 1: 2q*: √7q* (e.g. 40 wt%, 30 °C), 1: √3q*: √7q* (e.g. 60 wt%, 10–40 °C), 1: √3q*: 2q*: √7q* (e.g. 45 wt%, 30 °C) with higher-order peaks (3q*, √12q*, v13q*) at some conditions (Fig. S26(a)†). The cubic phase is characterized by SAXS peaks at q*, √2q*, √3q* √5q* and √6q*, matching the pattern commonly observed for BCC sphere phases formed by coil–coil block copolymers. Phase transitions from lamellar to hexagonal to cubic are possible in theory,57,58 and experimental studies have been reported as a function of changing concentration, salt concentration, or polymer volume fraction.59–63 Therefore, there is precedence in the literature for such transitions, which can be driven in this system by changes in volume fraction and the LCST/UCST behaviors of PNIPAM and PDMAPS which changes polymer hydration as a function of temperature. PDMAPS–mCherry–PNIPAM shows ordered structures across the widest range of concentration in four bioconjugates in this study. In addition, PDMAPS–mCherry–PNIPAM shows a larger number of higher order peaks, q*: 2q*: 3 q*: 4 q* (lamellar) and q*: √3q*: 2q*: √7q*: 3q*: √12q*: v13q* (hexagonal), than the other three types of bioconjugates. Higher order peaks and ordered structures across the widest range of concentration suggest that ABC type block copolymers (e.g. PDMAPS–mCherry–PNIPAM) can achieve well-controlled nanostructured protein materials driven in part by repulsive interactions between two different water-soluble polymers.
The domain spacings of PDMAPS–mCherry–PNIPAM change from 27.3 nm to 32.9 nm depending on the concentrations and temperatures (Fig. 4c). However, there is no obvious trend for the domain spacings depending on the concentrations or temperatures, potentially due to the fact that many phase boundaries are observed. However, the domain spacings are typically twice the calculated domain spacing, suggesting bilayer type arrangements (Fig. 4d), which is expected due to the asymmetry of the block copolymer.
In previous studies, it was revealed that electrostatic segregation of the protein block from a polymer block contributes significantly to the driving force for microphase separation, and the ionic polymer block results in much weaker segregation: mCherry–PNIPAM showed the lowest CODT but mCherry–PDMAPS showed the highest CODT.44 In this study, PDMAPS–mCherry–PNIPAM shows ordered structure across the widest range of concentration of the four bioconjugates (Fig. 5). Because the protein–polymer repulsion for PDMAPS and mCherry is relatively weak, this comparison implies that the repulsive interaction between two coil polymers is one of the key factors enhancing the self-assembly of ABC triblock copolymers (e.g. PDMAPS–mCherry–PNIPAM). The use of triblocks therefore provides a powerful method to induce ordering in protein–polymer bioconjugate systems that may be weakly ordered or disordered.
![]() | ||
Fig. 5 The range of concentration and temperature yielding ordered structures for the four types of bioconjugates. |
Despite their different topologies, the phase diagrams of mCherry–PNIPAM (BA) and mCherry–(PNIPAM)2 (BA2) showed a quite similar general shape. However, the BA2 type showed smaller domain spacing than the BA type, suggesting that branching architectures result in narrower polymer domain morphologies and larger cross-section area. PNIPAM–mCherry–PNIPAM (ABA) showed weaker ordered structures than other bioconjugates: only lamellar ordered structures were observed at concentrations above 50 wt%. In contrast, PDMAPS–mCherry–PNIPAM (ABC) shows ordered structure across the widest range of conditions in the four bioconjugates. Our previous studies showed that (1) volume fraction of the polymer46 and (2) electrostatic segregation of the protein block from the polymer block44 contribute significantly to the driving force for microphase separation. This study further reveals that (3) topology is also a key factor in governing microphase separation and (4) the repulsive interaction between the two coil polymers of ABC type bioconjugates is the main factor enhancing the self-assembly of globular protein materials. Therefore, the ABC topology provides a practical approach to enhance microphase separation of a variety of proteins, allowing nanostructural control in protein materials.
The molar masses and dispersities (Đ) of synthesized PNIPAMs were analyzed by gel permeation chromatography using a Waters GPC system equipped with a Waters 1515 Isocratic HPLC Pump, Waters 2487 Dual λ Absorbance Detector, and two ResiPore columns, 7.5–300 mm (Agilent Technologies) in N,N-dimethylformamide (DMF) with 0.02 M LiBr as the mobile phase. The signals were collected from a miniDAWN TREOS MALS detector and a Wyatt Optilab T-rEX refractometer. A dn/dc value of PNIPAM = 0.0761 mL g−1 was used to obtain absolute molar masses. The molar masses and dispersity of synthesized PDMAPS were analyzed by gel permeation chromatography using an Agilent Technologies 1260 Infinity system with two Aquagel columns in 0.5 M NaCl (aq) with 0.02% sodium azide as the mobile phase. The signals were collected from a Wyatt DAWN HELEOS II multi-angle light scattering detector and a Wyatt Optilab T-rEX refractometer. A dn/dc value of PDMAPS in 0.5 M NaCl solution of 0.1423 mL g−1 was used to obtain absolute molar masses.
After bioconjugation with PNIPAM, the bioconjugates were purified by 1 M ammonium phosphate precipitation to remove unreacted mCherry (three times) and then by NiNTA chromatography to remove the excess PNIPAM.
After bioconjugation with PDMAPS, the bioconjugates were purified by Ni–NTA chromatography to remove the excess PDMAPS. The eluent containing free protein and bioconjugate was dialyzed into 10 mM Tris–Cl pH 7.0 buffer for FPLC purification. Bioconjugates were purified using a HiTrap Q HP 5 mL anionic exchange column, 10 mM Tris–Cl pH 7.0 as buffer A and 10 mM Tris–Cl 0.5 M NaCl pH 7.0 as buffer B in an AKTA pure FPLC system. The results of SDS-PAGE are shown in Fig. S16 and S19.†
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8py01228h |
This journal is © The Royal Society of Chemistry 2019 |