Miriam
Simon‡
a,
Lauren
Matthews§
b and
Yeshayahu
Talmon
*a
aDept. of Chemical Engineering and The Russell Berrie Nanotechnology Institute (RBNI), Technion-Israel Institute of Technology, Haifa 3200003, Israel. E-mail: ishi@technion.ac.il
bESRF, The European Synchrotron, 71 avenue des Martyrs, CS 40220, 38043 Grenoble Cedex 9, France
First published on 25th July 2024
Self-assembly is a key process in forming biological materials. Especially the interaction between amphiphiles and polyelectrolytes has been widely investigated in recent years due to their potential application in industry and medicine, with a special focus on gene therapy. Accordingly, we investigated the formation of lipoplexes by mixing the cationic lipid DOTAP (1,2-dioleoyl-3-trimethylammonium-propane (chloride salt)) with different anionic polyelectrolytes (PE), such as NaPA (sodium polyacrylate), CMC (sodium carboxymethyl cellulose) with different degrees of substitution (DS, namely, different charge density), PSS (sodium polystyrenesulfonate) and DNA (deoxyribonucleic acid sodium salt). The goal of this project was to explore the influence of different system parameters, such as the charge ratio, CR = [+]/[−] = [DOTAP]/[PE], the charge density of the PE, or the type of PE on the morphology of the formed complexes. The investigation of these systems was performed by cryo-transmission electron microscopy (cryo-TEM), and with small-angle X-ray scattering (SAXS), to support our findings. In our experiments, we obtained a comprehensive picture of the formed lipoplexes, and how their structure depends on the different properties of the employed polyelectrolyte. Although the basic nanostructure of all complexes is lamellar, their detailed morphology depends strongly on parameters of the PE, e.g., the persistence length, charge density, or the polymer backbone. Understanding these specific interactions will allow the formation of more stable and optimized complexes as they are needed for drug or genetic material delivery.
Liposomes are widely used as delivery systems in cosmetics, such as sunscreen, creams, or hair products. For example, it was found, that the use of liposomal cream significantly contributed to repairing the surface lipid layer of the skin.5 In the food industry, liposomes may be used to encapsulate proteins, enzymes, vitamins, or flavors.6 Of course, liposomes are also a favored tool for drug delivery applications, as they can encapsulate hydrophilic drugs or bind hydrophobic drugs into or onto their membrane. Besides being carrier vehicles, liposomes also offer more advantages, such as protection of sensitive drug molecules and enhanced cellular uptake. New liposomal vaccines are also being developed. Here, the liposomes can act as an adjuvant, as well as a carrier for antigens and coadjuvants.7
For most applications, the lipids are mixed with other compounds to modify their properties with respect to rheology, liposome size, stability, zeta potential, or phase transition temperature.8,9 Here, especially the interactions between charged lipids and oppositely charged polyelectrolytes have been widely investigated in recent years, experimentally as well as theoretically,10 due to their many applications in industry and medicine, with a special focus on gene therapy.
One important potential application of lipid/polyelectrolyte complexes is in transfection.11 Transfection is the process of delivering genetic material into the cell nucleus to modify its genetic makeup. It is used as a therapy for cancer and a range of other genetic diseases. This approach is based on the complexation of DNA (deoxyribonucleic acid) with cationic lipids.12,13
Previous studies show that the addition of DNA to unilamellar, oppositely charged liposomes, induces disintegration of the liposomes and the formation of large multilamellar nanostructures instead, in which the polymer is sandwiched between two adjacent lipid bilayers. Alternatively, the DNA can be coated by lipid monolayers in an inverse hexagonal lattice.14,15 This was shown to occur for different lipids, such as DOTAP (1,2-dioleoyl-3-trimethylammonium-propane (chloride salt)), DOPC (dioleoylphosphatidylcholine), BFDMA (bis(n-ferrocenylundecyl) dimethylammonium bromide) and others, and was studied by small-angle X-ray scattering (SAXS), dynamic light scattering (DLS), cryo-TEM.16–18 Depending on the packing parameter of the lipid, also hexagonal structures in which thread-like lipid micelles are coated with DNA are possible as shown by Ewert et al.19
Lipids are available in many different forms and shapes. BFDMA, for example, contains a ferrocene group, whose oxidation state changes the net charge of the lipid. Using reduced BFDMA leads to multilamellar complexes, while aggregates with oxidized BFDMA are loose and disordered.20 Similarly, the mixing of two lipids, e.g., a charged and a neutral one, affects the structure of the formed complexes.21 Lamellar packing of lipid molecules was typical for the complexes formed from the cationic lipid-enriched mixtures, while inverted hexagonal arrays were found for the neutral lipid-enriched complexes.
Another very important parameter, influencing the transfection efficiency, is the cationic lipid/DNA ratio (charge ratio). However, the influence of the charge ratio on the complex formation seems to be very specific to the materials used.22,23
When charged lipids are mixed with other, oppositely charged polyelectrolytes, similar structures to the complexes with DNA are found, but their detailed morphology depends strongly on the respective system. Many factors, such as the type of lipid, the type of polyelectrolyte, and their mixing ratio are important. Also, external stimuli, such as the pH or temperature, can influence the aggregation behavior and with this the appearance of the resulting complexes. Some examples are summarized in the next few paragraphs:
Golan et al. employed sodium polyacrylate, (NaPA) as a polyelectrolyte with two different lipids (DOTAP and BFDMA) and found onion-like structures, very similar to the ones previously found when mixing the same lipids with DNA. They suggested, that the multilamellar aggregation is energetically favorable, and may be the preferred mode of aggregation for many of those systems.24
The effect of polyelectrolyte stiffness was also explored previously by adding the two very different polyelectrolytes, PAA and PSS (polystyrene sulfonate) to DOTAP and DDAB (didodecyldimethyl-ammonium bromide) systems. In all systems, multilamellar nanostructures were formed, but depending on the persistence length of the polyelectrolyte, very different spacings between the lamellae were found. The molecular weight (Mw) of the polyelectrolyte did not have a marked effect.25
Also, external stimuli were found to influence the lipid/polyelectrolyte complexes. Many polyelectrolytes are weak acids or bases, thus the pH of the solution affects their charge density, and, through it, their interaction potential. This is known for surfactant/polyelectrolyte systems,26 but has also been studied in lipid/polyelectrolyte complexes. If the pH is lowered below the pKa of the polyelectrolyte, no complexes are formed.25 Neumann, et al. showed how the lamellar structure of lipid/DNA films was sensitive to hydration and temperature.27
Obviously, the nanostructures of the complexes formed by lipids and polyelectrolytes also affect their macroscopic properties, be it rheological properties or the effectiveness of transfection. With respect to transfection, the scope of promising complexes for this application is not restricted to lipoplexes but can be extended to polyelectrolyte-DNA complexes, so-called polyplexes.28 Recently Rappoport et al. studied the complexation of a double hydrophilic block-copolymer with DNA and found a variety of very stable structures. At the right charge ratio, hexagonal structures can be formed (preferred for transfection), which are not affected by the interaction with blood proteins for several days.29
In our work we are interested in finding correlations of how the molecular constitution of a system affects its self-aggregation behavior in the liquid phase, and in how the formed aggregates lead to certain properties of the system and its functions. In the present study, we investigated the interactions of liposomes with selected polyelectrolytes. The cationic lipid DOTAP was mixed with different anionic polyelectrolytes (PEs), such as NaPA, CMC (sodium carboxymethyl cellulose) with different degrees of substitution (DS), namely, different charge density, PSS, and DNA. The goal of this project was to explore the influence of different system parameters, such as the charge ratio CR = [+]/[−] = [DOTAP]/[PE], the charge density of the PE or the type of PE on the morphology of the formed complexes. We studied these systems by cryo-transmission electron microscopy (cryo-TEM), complemented by small-angle X-ray scattering (SAXS) measurements.
From our experiments, we obtained a comprehensive picture of the formed complexes, and how their detailed structure depends on the different properties of the polyelectrolytes used. Although the basic nanostructures of the complexes are lamellar, their detailed morphology depends strongly on parameters like the persistence length, charge density, or polymer backbone diameter. These aggregates serve as a model system for lipoplexes as they might be used in gene therapy systems, where stable complexes with tunable properties are needed for optimized delivery.
To prepare mixed complexes, the DOTAP and polyelectrolyte stock solutions were mixed and diluted with water to obtain different charge ratios (CR = [+]/[−] = [DOTAP]/[PE]) at a constant DOTAP concentration of 1 mM (first DOTAP, then water was added, and last the PE). The complexed samples turned milky/turbid and sometimes phase-separated, showing a white precipitate settling at the bottom of the vial. The phase behavior was documented by photographs. For phase separated samples, only the supernatant was used for SAXS measurements and cryo-TEM imaging.
From previous studies, we know that the natural pH of the aqueous solutions we worked with in this study is around 7.25 Because the pKa of all the components in this study are well below 5, we did not control the pH of our solutions.
We imaged the specimens by an FEI (now Thermo Fisher Scientific) Talos F200C, FEG-equipped high-resolution TEM, operated at 200 kV. Specimens were equilibrated in the microscope below −178 °C in a Gatan 626 (Gatan, Pleasanton, CA) cryo-holder, and imaged using a low-dose imaging procedure to minimize electron-beam radiation-damage.33 Contrast was enhanced by a Volta phase plate. Images were recorded digitally by an FEI Falcon III, direct-imaging camera, and afterwards processed by the Photoshop software to enhance contrast.
We studied a series of DOTAP liposomes mixed with the biopolyelectrolyte CMC1.2 (MW = 250 kg mol−1, degree of substitution, DS = 1.2) in a wide range of charge ratios, both at polyelectrolyte and at lipid excess, using visual inspection, cryo-TEM, and SAXS. When DOTAP and CMC1.2 stock solutions are mixed, a precipitation in form of a white powder can be observed (see Fig. 2). The amount of precipitated material is highest close to the charge equilibrium, at CR = 1 and 0.5, and decreases in both directions of polyelectrolyte and DOTAP excess. The appearance of precipitation in such systems of oppositely charged materials at a charge ratio of 1 is not unusual, and can even be expected, as here the positive and negative charges completely neutralize each other, forming uncharged complexes, which lack the electrostatic repulsion needed for stabilization. It is however surprising to see precipitation occurring even at very low charge ratios (polyelectrolyte excess), which indicates a very strong interaction between the two materials, giving precipitate of charge ratio of 1, and leaving the supernatant with a charge ratio close to the specified global one. The supernatant of the samples remaining, once the precipitate has settled, ranges from clear (at high CR, hence DOTAP excess) over slightly blueish (CR = 1 and 0.5, charge equilibrium) to milky (low CR, PE excess). The color indicates differently sized complexes remaining in solution, the complex size increases with decreasing charge ratio.
Fig. 3 shows cryo-TEM images of three different DOTAP/CMC1.2 samples at charge ratios CR = 0.5, 0.2, and 0.1 all in the polyelectrolyte excess regime of the phase diagram. The specimens were prepared with the supernatant of the samples only, which had a bluish-to-milky color. As can be seen in Fig. 3, DOTAP and CMC1.2 form large complexes with a lamellar or multilayered inner structure at all charge ratios; no dependence on the charge ratio is visible in these images. Here it is important to note, that cryo-TEM images show only a projection of the sample, as we are looking at it in transmission mode. Lamellar stacks that are oriented perpendicular to the electron beam show the nice characteristic pattern seen in Fig. 3, while all other orientations show either blurred lamellae or no pattern at all. Hence, areas that do not show the lamellar pattern in the images, do not necessarily indicate lack of order, but most likely show different orientations of the lamellar stacks relative to the electron beam.
Further analysis of the multilayers (analyzed from 5 different images taken from different locations of the grid) yields an average spacing of 5.32 ± 0.07 nm (see ESI† for multilayer spacing analysis and FFT of the images). The dependence of the overall complex size on the charge ratio is not visible in the cryo-TEM images, due to the polydispersity of the complexes and their tendency to form large agglomerates.
In addition to the polyelectrolyte excess regime, we also investigated samples where CR ≥ 1, at lipid excess. Also here, the samples show multilayered complexes. With increasing DOTAP excess (higher CR), however, DOTAP vesicles can be found coexisting with the mixed complexes, as the available amount of polyelectrolyte is insufficient to sufficiently complex all available DOTAP (Fig. S4, ESI†).
While the cryo-TEM images give easily accessible qualitative information about the lamellar nature of the complexes, quantitative and especially statistically relevant information is harder to obtain. Thus, the same samples were also analyzed by small-angle X-ray scattering (SAXS), as shown in Fig. 4. All SAXS curves show a clear first-order correlation peak at q* ∼ 1.2 nm−1, this peak arises from the inner order of the complexes. In some curves a second-order lamellar peak is also visible, indicating a higher order of the multilayers. For lamellar structures the 1st order peak is expected at q*, the 2nd order peak at 2·q*, the 3rd order peak at 3·q*, and so on. The peaks found here match this pattern and therefore confirm the existence of the ordered lamellar structure we see in the cryo-TEM images. From the position of the 1st order peak, the repeating distance can be calculated via d = 2π/qmax ∼ 5.2 nm. This value is in good agreement with the layer thickness deduced from the cryo-TEM images. It can be seen that the position of the 1st order peak is unaffected by the CR, so the spacing between the multilayers does not change, as was also found in the cryo-TEM images. This finding seems logical since the spacing of the multilayers should depend mostly on the materials used and not on their mixing ratio.
Fig. 4 right is exemplarily shown in the ESI† as Fig. S6 with errors. It can be seen, that the error bars are smaller than the shown symbols in all cases, and all peaks, including higher order peaks, are significant compared to the error. Thus, errors are omitted for more clarity in other SAXS graphs.
However, some differences in the heights of the peaks can be observed (see normalized Fig. S5 (ESI†) for a better comparison of peak height). At CR = 10, a very high mixing ratio, the peak is very small and also somewhat broader than at the other mixing ratios. Here, there is only very little CMC1.2 present in the system, so not as many complexes can be formed, and their layers are less ordered. Whether excess DOTAP exists in the system cannot be easily seen from the scattering data. At a slightly lower CR of 2, the 1st order peak is already clearly visible, and the 2nd order peak is just emerging. At CR = 1, the 2nd order peak is clearly visible, but it is most pronounced at CR = 0.5, showing that at this charge ratio, the multilayers have reached their highest ordered state. The CR = 0.5 curve also shows a generally lower intensity than the other ones, which can be explained by the large amount of precipitation, which occurred in this sample. Less material is left in solution, which results in a lower scattering intensity.
When lowering the charge ratio even further to polyelectrolyte excess, it can be seen that the 2nd order peak is disappearing again (CR = 0.2). At CR = 0.1, there is a large excess of polyelectrolyte in the system, which is not participating in the complex formation, but adding to the background scattering.
One might wonder why the highest ordering is achieved at CR = 0.5, at PE excess, and not at the equimolar charge ratio of CR = 1. This can be explained by the degree of deprotonation of the PE, which could be controlled by the pH, but this parameter was already studied elsewhere,25 and would be beyond the scope of this manuscript. If the polyelectrolyte is not fully deprotonated at the given conditions, the actual charge ratio would be higher than the nominal one, meaning more PE would be needed to reach charge equilibrium. For that reason, experimentally, we find the highest ordering at CR = 0.5 and not CR = 1. This also explains the large amount of precipitate visually observed in the sample at CR = 0.5.
Samples were prepared analogously to the samples described above, mixing the different CMC solutions with the DOTAP stock solution, this time at a constant charge ratio of CR = 0.5, where the most ordered complexes had been found. Just as before, samples were examined by visual inspection, cryo-TEM, and SAXS.
From the visual appearance of the samples, it becomes obvious, that not the same amount precipitated from all samples (see Fig. S7, ESI†). The amount of precipitate decreases with increasing charge density (DS). In addition, the color of the supernatant changes from colorless to slightly blueish, which indicates an increase in particle size in the supernatant with increasing charge density of the polyelectrolyte.
Cryo-TEM images were taken of the three samples (Fig. 5). All images show large multilayered complexes. No systematic difference in the aggregate size could be determined by cryo-TEM because the aggregates are too large and too polydisperse in size to allow a reliable comparison between the three systems. From the images, it appears that the spacing of the multilayers within the structures does not change with DS. This is expected because the layer thickness depends on the polyelectrolyte backbone and on its flexibility, not on the charge density.
It should be added that all complexes shown in this work were extremely radiation sensitive, allowing only very short exposure times of ∼1 s in the low-dose imaging mode, where only ∼15 e− Å−2 were used. Higher electron doses resulted in destruction of their structure by electron-beam radiation-damage.37
A more quantitative comparison can be obtained by SAXS. In Fig. 6, we note the correlation peak, that arises from the ordered multilayered structure, increasing in intensity with increasing charge density (see normalized data in Fig. S8 for a better comparison of the peak height, ESI†). An increase in the peak area means, that more material exists in this conformation. The narrowing of the peak can be interpreted as a higher order in the system. In addition, the 2nd order peak of lamellar structures cannot be seen in the DS = 0.7 curve, emerging in the DS = 0.9 curve, and is clearly visible at DS = 1.2. It seems that the higher the charge density of the polyelectrolyte, the easier it is for the system to form ordered structures. Hence, the highest ordered systems can be found at the highest charge densities. In Fig. 6, it seems, that the 2nd order peak of the DS = 1.2 curve is shifted towards higher q, compared to the one of DS = 0.9. When looking at the numbers though, this difference is rather small (compare peak positions in Table S1, ESI†). Whether this difference actually reflects a slightly lower lamellar spacing of the DS = 1.2 complexes compared to the lower charge densities (5.33 vs. 5.39 nm), cannot be said for certain.
Spacing (from SAXS) (nm) | Spacing – DOTAP double layer (−3.72 nm) (nm) | Persistence length lP (nm) | Charge density (e− nm−1) | |
---|---|---|---|---|
PSS | 4.62 | 0.90 | 10 | 4 |
NaPA | 4.95 | 1.23 | 1.3 | 3.1 |
CMC | 5.33 | 1.61 | 6–16 | 0.7–1.2 |
DNA | 6.03 | 2.31 | 40–80 | 5.9 |
When analyzing the visual appearance (Fig. S9, ESI†) of the samples with different polyelectrolytes and DOTAP, at the same charge ratio of CR = 0.5, it is clearly seen, that the amount of precipitation decreases with increasing charge density, and also the color of the supernatant strengthens, indicating larger aggregates. This was already observed for the different CMCs (0.7–1.2 charges per nm), but this finding can now be generalized for different polyelectrolytes as well. The DNA sample (5.9 charges per nm) still shows a very small amount of precipitation, despite having the highest charge density, while the NaPA (3.1 charges per nm) and PSS (4 charges per nm) samples are completely homogeneous.38 Due to the double helical structure, DNA is also the stiffest polyelectrolyte, which indicates that this is also an important parameter for complex formation. Thus, the charge density directly influences the ability of the polyelectrolyte to form stable complexes with DOTAP, but a certain flexibility of the polyelectrolyte chain is also needed.
The cryo-TEM images of Fig. 7, show a range of different structures for different polyelectrolytes. Especially NaPA samples form highly ordered, spherical, multilayered complexes. In contrast to the complexes seen previously. These complexes are of very defined round shape and are made of very neatly organized layers. It seems that the high charge density of NaPA, together with its very flexible backbone (lP = 1.3–1.5 nm),39,40 allow it to organize into these well-formed structures, as it has also been observed before.25,41 It is interesting to note, that PSS with an even higher charge density but slightly lower flexibility (lP = 10 nm)38 does not form similarly ordered complexes. The backbones of the two other polyelectrolytes, CMC and DNA, are also much more rigid than that of NaPA, which results in less organized arrangements (lP(CMC) = 6–16 nm; lP(DNA) ∼ 40–80 nm).42,43 The comparison between DNA and CMC samples shows a slightly higher order for DNA, most likely because of the very high charge density. However, the rigid double-helical backbone of DNA prevents it from bending into uniform spherical multi-layered aggregates, like NaPA.
It also becomes apparent in these images, that the layer thickness of the formed complexes strongly depends on the polyelectrolyte. NaPA with the thinnest and most flexible backbone forms very narrow layers, while the double-helical DNA shows a much larger spacing.
This difference in layer thickness is confirmed by SAXS measurements (Fig. 8). The characteristic multilayer peak is shifting towards smaller q values with increasing polyelectrolyte diameter, which directly translates to larger repeat distances. The values calculated from the peak positions are: d = 4.6 nm for PSS, d = 4.95 nm for the narrow and very flexible NaPA molecules, d = 5.33 nm for CMC, which is composed of sugar units and therefore already is somewhat bulkier, and the largest spacing with d = 6.03 nm is found for DNA, due to its double-helical structure. 2nd order peaks are clearly visible in all of the samples (except for PSS, but this is most likely due to poor data quality) and appear at 2·q*, thereby confirming lamellar structures. The overall intensity of the CMC-curve is lower, due to the precipitation in that sample.
![]() | ||
Fig. 8 SAXS data of 1 mM DOTAP mixed with different polyelectrolytes, all at CR = 0.5. The peak position is shifting to smaller q, indicating larger spacings, from PSS to NaPA to CMC1.2 to DNA. |
When comparing the repeat distances of the multilayers as obtained from SAXS with the cryo-TEM images shown in Fig. 7, we are faced with a problem: the layers of the DOTAP/NaPA complex appear much thinner in the cryo-TEM image than the values obtained from the SAXS peak. This becomes even more clear when comparing the images and numbers with those of DOTAP/PSS complexes. As mentioned above, the scattering data is better suited for obtaining quantitative results and imaging is prone to artefacts. This is also the case here. Due to the high curvature of the NaPA layers in the complex and the transmission mode, the resulting image projection sees an overlay of multiple layers, looking much denser than they are in reality.
The highest ordered complexes are found at CR = 0.5, below the nominal point of charge equilibrium. This leads to the assumption, that the polyelectrolyte is not fully deprotonated at the studied conditions. The lipid/polyelectrolyte complexes will still achieve their highest ordering at the charge equilibrium, but experimentally we find this point to be at CR = 0.5 and not CR = 1. For transfection, however, usually higher charge ratios at lipid excess are employed.23 Presumably, the existence of free liposomes boosts the transfection efficiency,44 and the internal order of the complexes seems to be less important.
The complexes exhibit a higher degree of order the higher the charge density of the polyelectrolyte is. This was shown by employing CMC with different degrees of substitution: the height of the SAXS peak resulting from the lamellar ordering increased with increasing charge density. The very flexible and very highly charged NaPA even forms very symmetrical spherical multi-layered complexes, when interacting with DOTAP, showing that the interaction, not only depends on the mixing ratio and the charge density but also on the persistence length of the PE. If the polyelectrolyte is very flexible, like NaPA, it may bend around the previously existing vesicles and form multi-layered spheres. Stiffer polyelectrolytes break up the spherical shapes of the vesicles and form lamellar stacks that are randomly oriented within one complex. Fig. 9 shows a schematic of this situation.
Finally, the multilayer spacing is influenced by many parameters. By SAXS we found: d(PSS) = 4.6 nm < d(NaPA) = 4.95 nm < d(CMC) = 5.33 nm < d(DNA) = 6.03 nm. Naturally, a very bulky backbone like that of the double-helix DNA requires more space than the thin sodium polyacrylate. Hence, one could argue, that the backbone diameter of the polyelectrolyte defines the lamellar spacing. The persistence length could also be the dominating parameter since more flexible polymers could fold up more densely on top of the lipid bilayer, and therefore result in smaller lamellar spacings. But since we found the lamellar spacings with PSS to be smaller than with NaPA, both of these explanations cannot be the only factor here.
When the diameter of the pure lipid bilayer (3.72 nm) is subtracted from the lamellar spacing, it becomes clear, that the remaining space, which must be filled by the polyelectrolyte, is rather narrow (see Table 1). Geometrically it seems unlikely for any polyelectrolyte to be in a coiled or folded formation, independent of its persistence length. Hence, more factors, other than the persistence length must be influencing this situation. Possibly due to the slightly higher charge density of PSS, more electrostatic bonds are formed, freeing more Cl− ions (DOTAP counterions), which then enables a denser packing. Maybe other, hard-to-define parameters like hydration of the polyelectrolyte chain are the important parameters here. But on this, we can only speculate. Most likely it is a combination of all parameters mentioned above, that influences the multilayer formation, which means, that the lamellar spacing of lipid/polyelectrolyte complexes remains hard to predict.
Understanding the specific relations between lipids and oppositely charged polyelectrolytes, including the many very sensitive parameters that influence the final aggregate, is essential for the development of finetuned complexes as they are needed for optimized drug and genetic material delivery.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d4sm00489b |
‡ Present address: Bausch+Lomb GmbH, Brunsbüttler Damm 165/173, 13581 Berlin, Germany. |
§ Present address: ISIS Pulsed Muon and Neutron Source, Rutherford Appleton Laboratory, Chilton, OX11 7QY, UK. |
This journal is © The Royal Society of Chemistry 2024 |