I.
Orue
a,
L.
Marcano
b,
P.
Bender
c,
A.
García-Prieto
de,
S.
Valencia
f,
M. A.
Mawass
f,
D.
Gil-Cartón
g,
D.
Alba Venero
h,
D.
Honecker
i,
A.
García-Arribas
be,
L.
Fernández Barquín
c,
A.
Muela
ej and
M. L.
Fdez-Gubieda
*be
aSGIker, Universidad del País Vasco – UPV/EHU, 48940 Leioa, Spain
bDpto. Electricidad y Electrónica, Universidad del País Vasco – UPV/EHU, 48940 Leioa, Spain. E-mail: malu.gubieda@ehu.eus
cCITIMAC, Universidad de Cantabria, 39005 Santander, Spain
dDpto. Física Aplicada I, Universidad del País Vasco – UPV/EHU, 48013 Bilbao, Spain
eBCMaterials, Basque Center for Materials, Applications and Nanostructures, UPV/EHU Science Park, 48940 Leioa, Spain
fHelmholtz-Zentrum Berlin für Materialien und Energie, Albert-Einstein-Str. 15, 12489 Berlin, Germany
gStructural Biology Unit, CIC bioGUNE, CIBERehd, 48160 Derio, Spain
hISIS, STFC Rutherford Appleton Laboratory, Chilton, Didcot OX11 0QX, UK
iInstitut Laue-Langevin, 38042 Grenoble, France
jDpto. Inmunología, Microbiología y Parasitología, Universidad del País Vasco – UPV/EHU, 48940 Leioa, Spain
First published on 26th February 2018
Magnetospirillum gryphiswaldense is a microorganism with the ability to biomineralize magnetite nanoparticles, called magnetosomes, and arrange them into a chain that behaves like a magnetic compass. Rather than straight lines, magnetosome chains are slightly bent, as evidenced by electron cryotomography. Our experimental and theoretical results suggest that due to the competition between the magnetocrystalline and shape anisotropies, the effective magnetic moment of individual magnetosomes is tilted out of the [111] crystallographic easy axis of magnetite. This tilt does not affect the direction of the chain net magnetic moment, which remains along the [111] axis, but explains the arrangement of magnetosomes in helical-like shaped chains. Indeed, we demonstrate that the chain shape can be reproduced by considering an interplay between the magnetic dipolar interactions between magnetosomes, ruled by the orientation of the magnetosome magnetic moment, and a lipid/protein-based mechanism, modeled as an elastic recovery force exerted on the magnetosomes.
Among all the MTB species, Magnetospirillum gryphiswaldense strain MSR-1 is frequently used as a model strain of MTB due to its ease of cultivation in the laboratory. Magnetospirillum cells are helical and contain a variable number of ∼40–45 nm sized cuboctahedral magnetite (Fe3O4) magnetosomes arranged in a chain. The biomineralization of the magnetite nanoparticles and the assembly in a chain are complex and comprise different steps.4,13,14 Firstly, the magnetosome membrane is formed from the invagination of the cytoplasmic membrane comprised with proteins.15 Secondly, iron is taken up from the environment and transported into the magnetosome vesicle where nucleation in magnetite nanocrystals occurs,16–18 and finally, magnetosomes are assembled forming the chain.4,19–21 The present work is devoted to shed light on this last step. In particular, here we address the underlying mechanisms that determine the arrangement of the magnetosomes and consequently the geometry of the chain.
The arrangement of magnetosomes in a chain results from the interplay between the active assembly mechanism mediated by proteins and the magnetic dipolar interactions between the nearest magnetite nanoparticles.22–24 As shown by cryotomography imaging on Magnetospirillum, bundles of cytoskeletal filaments intervene in the chain assembly.15,25,26 These filaments, formed by the actin-like protein MamK, traverse the cell and position the chain in the middle of the cell. Another important protein involved in the chain formation is MamJ. One possible function of MamJ is to connect the magnetosome membrane to the cytoskeletal filament. As magnetosomes get closer together, magnetic dipolar interactions arise.
Here we show by electron cryotomography imaging that magnetosome chains of M. gryphiswaldense are not straight lines but appear slightly bent, a fact that has been also observed in previous studies.14,26 The implications of the chain shape on the magnetic response of magnetosome chains have then been addressed with complementary techniques performed on a set of bacterial arrangements: (i) small angle neutron/X-ray scattering (SANS/SAXS) on a bacterial colloid, (ii) macroscopic magnetometry on 3D and 2D fixed arrangements of randomly distributed and aligned bacteria, and (iii) X-ray photoemission electron microscopy (XPEEM) on an individual chain of magnetosomes extracted from bacteria.
Two main findings are achieved from this work. Firstly, the equilibrium magnetic moment of the magnetosomes is tilted 20° out of the [111] crystallographic axis, as concluded from the magnetic analysis. Secondly, the tilting of the magnetic moment is the key to understand the helical-like shape of the magnetosome chains. Indeed, the experimental chain patterns imaged by electron cryotomography are accurately reproduced by counterbalancing the magnetic dipolar interactions between magnetosomes, strongly affected by the orientation of the individual magnetic moments, and a lipid/protein-based mechanism, modeled as an elastic recovery force exerted on magnetosomes.
Besides the basic interest of this study, a precise knowledge of the mechanisms determining the chain shape is decisive in applications of MTB such as biological micro-robots. Precisely, recent studies propose MTB to be exploited as motors with embedded source propulsion secured by their flagella and embedded steering system to control the directional motion provided by their magnetosome chain.7
Fig. 1 Electron cryotomography (ECT) of magnetosome chains. (a) and (b) are ECT images of two different bacteria. Top: 3D tomograms reconstructed from the tilt-series of images obtained by ECT. Videos of the tomograms of the two chains viewed from different perspectives (as consecutive slices in XY, XZ and YZ plane orientations, and as an electron density map) can be found in the ESI (movies S1 and S2†). The insets show the Z-projection of the tomograms. Bottom: XY, YZ and XZ slices of the tomograms shown on top. The three slices are taken along the marked lines and cross on one of the magnetosomes. (c) ECT images of extracted magnetosomes. Top: Reconstructed tomograms. Bottom: Central XY, YZ and XZ slices of the magnetosome tomogram shown on top. (d) Schematic representation of a cuboctahedral magnetosome. |
ECT images of a magnetosome extracted from M. gryphiswaldense reveal the morphology of the magnetosome (Fig. 1c). The figure on the top is a reconstructed 3D tomogram, and below, the XY, YZ and XZ central sections of the tomogram are shown. These images evidence that magnetosomes have a strongly faceted morphology similar to the truncated octahedral one, which is adopted as a reference (Fig. 1d), with a mean size of 40 nm. In this truncated octahedron the 〈001〉 crystallographic axes define the growth directions of the square faces and the 〈111〉 axes those of the hexagonal faces. The magnetosome in the image has two neighbours (only partially observed in the image). The three magnetosomes have self-assembled so that their hexagonal faces are towards each other. In the same way, when magnetosomes are forming the chain inside the cell, numerous studies in Magnetospirilla show that they align with their hexagonal faces towards each other,27,28i.e., along one of their 〈111〉 crystallographic directions, as sketched in Fig. 5a and d. The [111] direction along which magnetosomes align defines the so-called chain axis.
The magnetic properties of the bacterial colloid have been further investigated via small angle X-ray scattering (SAXS) and polarized small angle neutron scattering (SANS). In a SAXS/SANS experiment the X-ray/neutron beam hits the sample and the scattered beam impacts on a 2D detector perpendicular to the incoming beam (Fig. 2b). The intensity I as a function of the scattering vector gives pure structural (Inuc(q)) information from SAXS/SANS experiments, or, additionally, magnetic (Icross(q)) information for SANS depending on the magnetic field on the sample and the polarization of the incoming and scattered neutron beam (Fig. 2c). Inuc(q) and Icross(q) are extracted from the 2D scattering patterns (Fig. 2c) as explained in the Methods section.
Fig. 2d shows the radially averaged 1D SAXS intensity I(q) of the colloidal dispersion measured in zero field. An indirect Fourier transform of I(q) results in the pair distance distribution function P(r) as displayed in Fig. 2e. The extracted distribution function exhibits three distinct peaks. The first peak has its maximum at about 25 nm and is nearly bell-shaped. The comparison of this peak with the pair distance distribution function of a homogeneous sphere with diameter DSAS = 48 nm shows qualitatively very good agreement. Hence, we surmise that the first peak corresponds to the nuclear scattering of the individual magnetosome of an average size of about 48 nm. This particle size is considerably larger than the one obtained by ECT (≈40 nm). The reason for this difference is that ECT is only sensitive to the core of the magnetosome (the magnetic nanoparticle), while the small angle scattering (SAS) techniques are sensitive to both the core and the surrounding lipid bilayer membrane. In particular for neutrons, the scattering cross section of H is very large, which means that organic materials with many H-atoms (such as the lipid membrane in our case) can generate large scattering signals.29 Considering a membrane thickness of ≈4 nm, the core diameter obtained by SANS is 48 − 8 = 40 nm, in agreement with the value obtained from ECT. The maxima of the second and third peaks of the distribution function are at about 75 and 125 nm, and both peaks are also nearly bell-shaped. The positions of the peaks agree well with the expected positions of the center of mass of the next neighbors in a chain-like structure of DSAS = 48 nm sized particles with a center-to-center distance of d ≈ 50 nm, close to the d ≈ 60 nm distance measured by ECT.
Fig. 2d additionally shows the structural 1D SANS intensities Inuc(q) detected for field strengths of μ0H = 2 mT and 1 T applied perpendicularly to the neutron beam. Both Inuc(q)'s are virtually identical, which means that already at 2 mT the bacteria are fully aligned in the field direction. Consequently an indirect Fourier transform of both nuclear scattering intensities results in superimposed pair distance distribution functions, whose second peak (r ∼ 75 nm) is about twice the height of the first peak (Fig. 2e). This verifies that the bacteria align with the chain axis parallel to the field, considering that within the chain each particle is surrounded on average by two neighbors and thus the probability to find a scatterer at this position is two times the probability to find a scatterer within the primary particle (first peak). To investigate if structural alignment equals magnetic saturation the cross-terms Icross(q) were analyzed.
The two cross-terms Icross(q) detected at 2 mT and 1 T display the same functional form (Fig. 2d) and accordingly the extracted distribution functions are qualitatively similar (Fig. 2e). The observation of only one peak can be explained by the fact that the cross terms depict the correlation between nuclear scattering length density and magnetization Mz(r) along the axis perpendicular to the applied field. In all cases the shape of the distribution function is comparable to the distribution function of a single sphere with a homogeneous scattering length density. This verifies that the particles are homogeneously magnetized (i.e. Mz(r) = Mz) and thus can be regarded as single-domain particles. However, the absolute values of P(r) detected at 2 mT are over the whole r-range systematically reduced by a factor of 0.83 compared to 1 T. Assuming that at μ0H = 1 T the system is saturated in field direction (i.e. Mz = Ms), this means that at 2 mT the magnetization in field direction amounts to Mz = 0.83Ms and that the magnetic saturation is achieved by a coherent rotation of the spins within the individual nanoparticles towards the field direction. These two processes (chain alignment followed by a coherent rotation of the magnetosome spins) explain the isothermal magnetization measurement of the colloid shown in Fig. 2a. In the latter, in agreement with the SANS result, at 2 mT Mz = 0.83Ms, and the coherent rotation of magnetosome spins would start at ≈15 mT, where Mz = 0.9Ms corresponds to a misalignment of the magnetic moment with the chain axis of θ ≈ 25°, where Mz = Mscosθ.
The hysteresis loops of randomly arranged and oriented bacteria are shown in Fig. 3a and b, respectively. The hysteresis loops of oriented bacteria have been measured at different angles with respect to between 0° and 90° in steps of 20 degrees for 3D arrangement and at discrete angles, 0°, 45° and 90°, for 2D arrangement. These measurements evidence that bacteria are highly anisotropic magnetic objects, since their hysteresis loops depend strongly on the direction of the applied field. This is not surprising, as one would expect that magnetosome chains behaved as a good compass, with a single magnetic easy axis oriented along the chain axis, which, as noted previously, is coincident with one of the magnetosome 〈111〉 crystallographic axes. However, as previously observed in Magnetospirillum magneticum AMB-1,30 the hysteresis loops of 3D and 2D bacterial arrangements do not correspond to those expected for a uniaxial Stoner–Wohlfarth model31 along the chain, since the hysteresis loop perpendicular to the chain axis is not anhysteretic as expected. This could be attributed to the misalignment of the magnetic moment already detected by SANS measurements.
More information on the magnetism of magnetosome chains has been gathered from the theoretical modelling of the hysteresis loops. Previous studies on M. gryphiswaldense show that the magnetic anisotropy of magnetosome chains can be described as a superposition of the cubic magnetocrystalline anisotropy of magnetite, with four equivalent easy axes directed along the 〈111〉 crystallographic directions, and a uniaxial shape anisotropy directed along the chain axis (thus parallel to one of the 〈111〉 axes) that originates from intra-chain dipolar interactions and dominates the magnetic response of the chain.32,33
Following this argumentation, as a first approximation, our approach to simulate theoretically the experimental hysteresis loops consists of considering the magnetosome chain as a collection of independent magnetic dipoles, where the equilibrium orientation of each magnetic moment is obtained by minimizing the single dipole energy density E, calculated as the sum of three contributions: (i) cubic magnetocrystalline energy of magnetite with anisotropy constant Kc, (ii) uniaxial anisotropy energy along the chain axis ([111] direction defined as ûuni = û111 in Fig. 1d) with anisotropy constant Kuni, and (iii) Zeeman energy in an external magnetic field :
E = Ecubic + Kuni[1 − (ûuni·ûm)2] − μ0MH(ûH·ûm) | (1) |
Misalignments of the chains with respect to the aligning field occurring during sample preparation and detected by SANS experiments have been considered by including a Gaussian angular distribution of the chain axes. We have tested three angular distributions around the chain axis: 15°, 25°, and 35°, in agreement with misalignment values reported previously.28,30 The data used for the simulations are the magnetocrystalline anisotropy constant (Kc = −11 kJ m−3) and magnetization (Ms = 48 × 104 A m−1) of magnetite and an effective uniaxial anisotropy constant along the chain axis Kuni = 12 kJ m−3 to account for both shape and magnetic interaction anisotropies. Even though the simulations are good enough at the perpendicular orientation, i.e., with the applied magnetic field perpendicular to the bacterial chain axis, they clearly deviate at smaller angles as shown in the polar representation in Fig. 3c, where the experimental and calculated reduced remanent magnetization and coercivity of bacteria in 2D and 3D arrangements are plotted for different orientation angles between 0° and 90°.
Since the misalignment of the bacteria due to sample preparation does not reproduce the hysteresis loops, we have tried a different approach in which we have considered that there is a tilting of the magnetization with respect to the chain axis inherent to the magnetosome. This approach allows the accurate reproduction of the experimental hysteresis loops for both 2D and 3D arrangements by setting the uniaxial easy axis ûuni at 25° with the chain axis in eqn (1), in the plane containing the chain axis ([111]) and the [100] directions. By setting this angle, the effective easy axis, and as a consequence the equilibrium magnetic moment, is found to lie 20° out of the chain axis, in close agreement with the 25° tilting observed in the previous SANS analysis of the colloid. Note that there is a good resemblance with experiment even for an angular dispersion as small as 15° (Fig. 3a and b). We used the same data for the simulations of the hysteresis loops at all orientation angles and for both 3D and 2D arrangements, Kc = −11 kJ m−3, Ms = 48 × 104 A m−1 and Kuni = 12 kJ m−3. The resemblance with experiment is more evident in the polar plots of the reduced remanent magnetization and coercivity shown in Fig. 3c.
The tilting of the magnetosome magnetization with respect to the chain axis is attributed to the competition between the shape anisotropy of the magnetosome, which presents a well faceted morphology as shown in Fig. 1c, and the magnetic interactions trying to align the magnetic dipoles along the [111] direction parallel to the chain axis. In fact, electron holography experiments on individual magnetosomes from Magnetovibrio blakemorei MV-1 clearly show that the magnetization direction of the particle is tilted with respect to the [111] crystallographic direction towards a long dimension of the particle, consistent with shape anisotropy dominating the magnetic state of the crystal.36
XPEEM measurements have been performed at room temperature on extracted magnetosomes due to the high absorbing power of the whole bacteria. Magnetosomes were deposited onto a conductive Si substrate either randomly (Fig. 4c) or oriented along an aligning field (Fig. 4d). As expected for magnetite, the Fe L3 XMCD signal shows three peaks (Fig. 4b): two minima at 707.4 eV and 709.6 eV and a maximum at 708.6 eV. The two minima correspond to the Fe2+ and Fe3+ ions occupying tetrahedral sites, and the maximum corresponds to the Fe3+ ions occupying octahedral sites. The magnetic contrast images were recorded at the Fe L3 resonance energy (709.6 eV) of magnetite. The set of magnetic images displayed in Fig. 4c and d show the space-resolved dichroic images obtained at different values of an in-plane magnetic field. The dichroic signal yields the magnetic moment, so that the hysteresis loops of selected areas of the magnetic images (enclosed in rectangles in Fig. 4c and d) have been obtained by plotting the corresponding dichroic signal as a function of the applied magnetic field (Fig. 4a).
The hysteresis loops of the randomly deposited magnetosomes enclosed in the region marked with the black rectangle in Fig. 4c and of a chain parallel to the applied field, marked with the green rectangle in Fig. 4c, are shown in Fig. 4a.
Magnetosomes deposited under an applied aligning field form longer chains. The XPEEM image together with the SEM image of one of these chains is shown in Fig. 4d. This chain is clearly not a straight line but is rather a zigzag, formed by segments oriented at different angles with the applied field. Two of these segments, enclosed in blue rectangles in Fig. 4d, form 60° with the applied field, and another two segments, in red rectangles, form 90°. The corresponding hysteresis loops are shown in Fig. 4a.
As shown in Fig. 4a, the experimental XPEEM loops can be accurately simulated to the theoretical model developed previously for the 3D and 2D chain arrangements (eqn (1)) with a smaller tilting angle of the effective easy axis (15°) and a larger anisotropy constant Kuni = 16 kJ m−3 as compared to the SQUID loops (Kuni = 12 kJ m−3). This is attributed to the larger distance between magnetosomes in chains inside the bacteria (d = 60 nm) than in chains of extracted magnetosomes (d = 50 nm). Indeed, for a 40 nm sized particle, the dipole pair potential energy is given by ∼μ0m2/4πd3 = 0.75 eV (m = MsV being the particle magnetic moment and d = 60 nm), while for the chains of extracted magnetosomes the dipolar energy from the neighbors at d = 50 nm would increase up to ∼1.35 eV.
Fig. 5 Equilibrium configuration of the magnetosome chain. (a) Schematic representation of two competing mechanisms: magnetic force pushing to align magnetosome magnetic moments along the stray field lines from neighboring particles, and lipid/protein-based mechanism modeled as an elastic recovery force acting perpendicularly to the chain axis, where the chain axis is the [111] crystallographic direction along which magnetosomes align, as highlighted in the figure. (b) Experimental reconstructions obtained from ECT imaging of the magnetosome chains shown in Fig. 1a and b. A zoom-in of the first magnetosome chain reconstruction highlights the deviation from a straight line. (c) Two stable solutions for the chain patterns obtained as explained in the text. A potential filament has been drawn as a guide for the eye. A video of the experimental reconstruction of the chain on the left (chain (a) in Fig. 1) together with the corresponding simulated chain viewed from different perspectives can be found in the ESI (movie S3†). (d) Schematic representation of the magnetic dipoles and the three independent variables used in the simulation: radial (ρ) and azimuthal (ϕ) coordinates for the magnetosome positions, and azimuthal orientation (φ) of the magnetic dipoles. θ is the polar angle of the magnetic dipoles, fixed to 20°. |
With this aim, we have developed an approach to explain the shape of the magnetosome chains that consists of quantifying the total energy of the chain by including the magnetostatic interactions between nanoparticles, and the contribution of the lipid/protein-based architecture embedding the magnetosome chain, modeled as a spring-like elastic energy. We focus on the stable geometry of chain arrangements and assume that close loop configurations such as rings or 3D clusters are avoided by the cytoskeleton inside bacteria. The same approach has been used by other authors,22,38 but in the present case, as proposed previously, the magnetization of each magnetosome is tilted 20° out of the chain axis.
Here we will summarize the main steps followed to calculate the total energy of the chain, but a detailed description can be found in the ESI.†Fig. 5a shows a section of a magnetosome chain composed of three magnetosomes. The particle in the centre is subjected to the stray magnetic field produced by the two neighbors (the dotted red lines in Fig. 5a represent the stray field produced by the particle at the bottom). We implicitly assume that the local torque that tends to align neighboring dipoles is counter-balanced by the cytoskeleton, since cryotomography and TEM images reflect that the hexagonal faces are face to face.27,28 The particle in the middle will tend to align its magnetization along the stray field lines and will undergo a magnetic force that compels it to shift towards the direction of its own magnetic moment. The magnetostatic energy associated with the magnetic force on magnetosomes is then implemented as the sum of the dipolar pair potential energies between nearest neighboring particles.
On the other hand, magnetite crystals are embedded in a lipid/protein-based architecture composed, among others, of the magnetosome membranes, the cytoskeletal filaments, and the proteins connecting the magnetosome membranes with the cytoskeleton. The contribution to the total chain energy of this lipid/protein-based architecture has been modeled based on two simplifying assumptions: firstly, that the distance between particles, d, is fixed, so that chains can bend or twist but cannot stretch; and secondly, that the net force on the magnetosomes can be ascribed to an elastic recovery force with elastic constant k acting perpendicularly to the chain axis, as shown in Fig. 5a.
The most stable chain configurations have then been predicted by minimizing the total energy of the chain considering three independent variables per particle, namely the radial (ρ) and azimuthal (ϕ) coordinates for the magnetosome positions, and azimuthal orientation (φ) of the magnetic dipoles (see sketch in Fig. 5d). The results have then been compared to the 3D reconstructions obtained from the chains imaged by ECT (Fig. 5b).
Considering a particle diameter D = 40 nm, distance between them d = 60 nm, and k = 0.1 × 10–3 N m−1, two stable solutions give the patterns shown in Fig. 5c, in excellent match with the two experimental chain reconstructions shown in Fig. 5b. They are helical-like shaped chains with slowly bending axes and a different pitch, given by ≈2L/n, L = Nd being the chain length, N the number of magnetosomes per chain, and n = 1 and 2, respectively. The energy difference between the two chain shapes is less than 1% of the total energy, hence they are approximately equally stable. The agreement between the simulated and the experimental chains is even more explicit in a video of images of both chains taken from different perspectives shown in the ESI (movie S3†).
A zoom-in of a section of one of the simulated chains shows the magnetosome dipole moments. The azimuthal increment of the individual dipole orientations (Δφ) between consecutive positions along the chain increases with n approximately as:
(2) |
The role of k in chain geometry is to scale up or down the radial positions of magnetosomes along the chain, so that larger k favours configurations approaching straight lines. In our case, k = 0.1 × 10–3 N m−1 results in assembling forces of the order of 0.6 pN for radial displacements of ≈10 nm as obtained from the simulations, similar to the force generated by the actin filament,39 and far below recent estimations of the fracture limit of the actin-based scaffolding filaments (≈30 pN) as reported by Körnig et al.40
On the basis of our results, a mechanism of chain formation is proposed. Our findings reveal that the chain shape is mostly driven by passive spontaneous magnetostatic effects triggered by the intrinsic anisotropy of magnetosomes, ultimately defined by their morphology. The magnetosome morphology is regulated during the biomineralization process. This is a genetically controlled process which involves a specific set of about 30 mam (magnetosome membrane) and mms (magnetic particle-membrane specific) genes.4 The product of the mamJ gene is an acidic protein that connects the magnetosome membrane to the filament. The deletion of mamJ results in bacteria that produce magnetosomes arranged in compact three-dimensional clusters instead of arranged in chains.26 The fact that they form clusters and not closed rings is consistent with our conclusions, namely the tilting of the magnetosome magnetic moment with respect to the chain [111] axis, since according to previous studies, more than four magnetosomes tend to form closed rings when their magnetic moment is parallel to the [111] axis.38 In the same way, when depositing the magnetosomes onto a 2D surface, magnetosomes tend to self-assemble in a close-packing configuration,41 and the tilting could also be behind the zigzag configuration of the chains observed by SEM in oriented magnetosomes. But, why would this helical-like shape benefit the bacteria? M. gryphiswaldense are long cells, easily reaching several microns long. They need a high magnetic torque to overcome the drag forces and orientate along the Earth's magnetic field. As a consequence, their magnetosome chain needs to be long to maximize the chain net magnetic moment. Indeed, their chain is frequently composed of more than 20 magnetosomes, which brings the total length of the chain to 1.5 μm or more, about 50–60% of the bacterial length. Such an object should necessarily be bent in order to accommodate this spiral-shaped microorganism. Thus a helical-like shape fits better, but it only changes slightly the total magnetic moment, hence hardly affecting the magnetic orientation. A genetic control of the magnetosome shape towards an energetically optimum chain arrangement is also observed in Magnetobacterium bavaricum MYR-1, a species which synthesize bullet-shaped magnetosomes arranged into bundles of magnetosome chains. The magnetosome magnetization within MYR-1 magnetosomes is parallel to the chain axis, which coincides with the long axis of the magnetosomes. Unexpectedly, the latter is parallel to the [100] crystallographic axis, a magnetically hard axis in magnetite, rather than along the [111] magnetocrystalline easy axis, due to compromise effects of shape anisotropy and intra-chain and intra-bundle interactions.42
In summary, our finding sheds light on the understanding of the magnetosome chain assembly during the biomineralization process of MTB, which may influence their potential future applications as biological micro-robots. Indeed, one of the major technical issues in the development of interventional platforms for guiding drug-loaded MTB is the directional magnetic field strength that needs to be produced at the human scale to induce sufficient directional torque on the chain of magnetosomes. A good knowledge of the magnetic configuration of the magnetosome chain will lead to more efficient and less costly drug-delivery platforms that may benefit a larger population.
For the alignment and 3D reconstruction of the tilt-series, we used IMOD software.46 We employed the Au fiducial markers during the alignment process, and 3D reconstruction was carried out by weight back-projection followed by a reconstruction algorithm named the Simultaneous Iterative Reconstruction Technique (SIRT). The resulting tomograms were visualized with ImageJ47 as a sequence of cross sectional slices in different plane orientations. Tomograms were processed using a median filter and visualized as 3D electron density maps using UCSF Chimera48 software.
The 2D bacterial configurations were prepared by depositing five 5 μL drops of the bacterial suspension (109 cells per mL) by the drop coating method41,49 onto a Si substrate. For the oriented bacteria, the deposition was done under an external applied magnetic field of 0.5 T. To obtain homogeneous samples, infrared radiation was used during the deposition aimed at accelerating the drying and minimising the surface tension. The resulting samples were finally oriented at different angles with respect to the applied magnetic field.
The 3D bacterial configurations were prepared by resuspending 250 μl of a bacterial colloid (1011 cells per mL) in 750 μl of an agar solution (2% agar and 98% water) at 60 °C to maintain the solution in a liquid state. To align the bacteria, a uniform magnetic field of 1 T was applied. After 3 minutes, the field was turned off, and the sample was cooled using liquid nitrogen until the temperature reached around 0 °C. This caused the agar to solidify, trapping the bacteria, and keeping this solid state at room temperature.
The application of the POLARIS mode enabled us to detect the non-spin flip (nsf) intensities , , where the superscripts indicate the polarization of the incoming neutron beam and the scattered neutrons with regard to the applied field direction, respectively (“+”: parallel, “−”: antiparallel). Defining x as the direction of the neutron beam and z as the direction of the applied magnetic field at the sample position the nsf intensities can be written as:56
(3) |
Here, Θ is the angle between the scattering vector and the magnetic field and the terms represent the Fourier transforms of the magnetization in the y, z direction. The superscript * indicates the complex-conjugated quantities and the constant bH = 2.7 × 10−15 m/μB−1, with μB being the Bohr magneton. To investigate the alignment of the bacteria in field direction we analyzed the 1D nuclear cross sections Inuc(q)∝|Ñ|2 and the 1D nuclear magnetic cross-terms . The purely nuclear scattering intensities were determined by integration of the nsf intensities in 10° sectors around Θ = 0° (, eqn (3)) and the cross terms Icross(q) by integration of in 10° sectors around Θ = 90° (, eqn (3)). The difference between the two nsf cross sections yields information on the polarization-dependent nuclear-magnetic terms. This difference allows one to highlight weak magnetic contributions relative to strong nuclear scattering (or vice versa). From the determined cross sections we extracted the underlying pair distance distribution functions in the same manner as for SAXS. In H2O, there is a contrast match of bacteria in the dispersion, such that only the magnetosomes and their arrangement are probed by SANS.
Magnetic imaging of the isolated magnetosomes was performed at room temperature by means of photoelectron emission microscopy (PEEM) by using X-ray magnetic circular dichroism (XMCD) as a magnetic contrast mechanism. Measurements were carried out at the UE49_PGM SPEEM beamline at Helmholtz-Zentrum Berlin. A magnetic solenoid attached to the sample holder allowed the application of an in-plane pulsed magnetic field ranging ±0.1 T to saturate the sample. During imaging the magnetic field range was restricted to ±27.5 mT. The incoming photon energy was tuned to the Fe L3 resonance (709.6 eV) to obtain element specific XMCD images of the magnetosomes as a function of the applied magnetic field. At each field a sequence of images was acquired with incoming circular polarized radiation (90% of circular photon polarization) with left (σ−) and right helicity (σ+), respectively. To improve statistics we acquired up to 400 images per helicity and field (3 s of exposure time). After normalization to a bright-field image, the sequence was drift-corrected, and frames recorded with the same helicity were averaged.
Footnote |
† Electronic supplementary information (ESI) available: TEM image of oriented bacteria, details on the calculation of the hysteresis loops, calculation of the chain energy, and 2D projection of the magnetic dipoles. Movies S1 and S2: Videos of the tomograms of the three chains in Fig. 1a and b viewed as consecutive slices in XY, YZ, and YZ plane orientations, and as an electron density map. Movie S3: Video of the experimental reconstruction of the chain in Fig. 1a (corresponding to n = 1 in Fig. 5b) together with the corresponding simulated chain (Fig. 5c) viewed from different perspectives. See DOI: 10.1039/C7NR08493E |
This journal is © The Royal Society of Chemistry 2018 |