Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A biomimetic approach to shielding from ionizing radiation: The case of melanized fungi

  • Thomas Vasileiou ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Software, Visualization, Writing – original draft

    Thomas.Vasileiou@esa.int

    Affiliation ESA - Advanced Concepts Team, European Space Research and Technology Centre (ESTEC), NL-2200AG Noordwijk, Netherlands

  • Leopold Summerer

    Roles Conceptualization, Project administration, Writing – original draft

    Affiliation ESA - Advanced Concepts Team, European Space Research and Technology Centre (ESTEC), NL-2200AG Noordwijk, Netherlands

Correction

31 Aug 2021: Vasileiou T, Summerer L (2021) Correction: A biomimetic approach to shielding from ionizing radiation: The case of melanized fungi. PLOS ONE 16(8): e0257068. https://doi.org/10.1371/journal.pone.0257068 View correction

Abstract

Melanized fungi have been shown to thrive in environments with high radionuclide concentrations, which led to the association of the pigment melanin with the protection against ionizing radiation. Several hypotheses regarding the function of melanin have been proposed. Yet, the exact mechanism behind the protective property of melanin is unclear and poorly explored. A better understanding of the mechanisms that are involved in increasing the tolerance of the organisms to ionizing radiation could lead to technology transfer to human-related applications. Effective protection from radiation is essential for human space flight in general and human missions beyond Low Earth Orbit specifically. In this paper, we follow a biomimetic approach: we test two of current hypotheses and discuss how they could be applied to radiation shield designs. First we focus on the interaction of melanin with high energy electrons, which has been suspected to reduce the kinetic energy of the electrons through a cascade of collisions, thus providing physical shielding. Second, we investigate if the spatial arrangement of melanin, organized as a thin film or a collection of hollow micro-spheres, affects its shielding properties. To this end, we measured experimentally and by numerical simulations the attenuation of β-radiation as pass through solutions and suspensions of melanin and contrasted the values to the ones of cellulose, a substance with similar elemental composition. Further, we investigate the spatial arrangement hypothesis using Monte Carlo simulations. In agreement with the simulations, our experiments indicated that melanin does not provide improved shielding in comparison to cellulose from β-radiation. However, our simulations suggest a substantial effect of the spatial arrangement on the shielding performance of melanin, a pathway that could be transferred to the design of composite radiation shields.

Introduction

Protection from ionizing radiation is one of the main challenges for human space flight, considering the biological stress radiation exerts on living organisms [14]. Risks associated with radiation exposure are present during the current manned low Earth orbit missions [5] and expected to be a major concern for future exploratory missions on the Moon and Mars [6]. Therefore, new approaches and materials able to address this problem and to provide effective shielding from ionizing radiation are urgently needed.

Environments with high background radiation levels can also be found on Earth. Yet, life is able to persist under these extreme conditions and living organisms have been found to adapt well and in some cases even thrive in such environments (e.g. the Chernobyl disaster site [7, 8]). Studying the mechanisms these organisms cope with ionizing radiation, may give inspiration for novel radiation shielding practises and materials. Interestingly, the percentage of fungi that synthesize the pigment melanin is substantially higher in the high background radiation areas than what is normally observed in places with typical ionizing radiation levels [8, 9]. This observation has lead to the hypothesis that melanin may provide a survival advantage in environments with ionizing radiation and has fueled many studies—from which we summarize the most relevant results in the following—confirming that indeed melanin plays a role in how fungi endure the ionizing radiation exposure.

Specifically, it have been demonstrated that melanized fungi exhibit increased proliferation when exposed to ionizing radiation in comparison to non-melanized strains or to low ionizing radiation environments [8, 10, 11]. Additionally, the protective effect of melanin can be transferred to organisms that do not produce the pigment; non-melanized fungal cells exhibited higher survival rates after irradiation, when melanin extracted from melanized cells was added to the culture medium [12]. The transfer of the protective effects of melanin has been also demonstrated on mice models; intravenous injection [1315] or ingestion [16] of melanin resulted in higher survival rates after exposure to γ-radiation in comparison to controls. Recently, the protective role of melanin was also demonstrated in fungal cells for deuteron radiation [17].

Still, how melanin is involved in the protection against the effects of ionizing radiation is poorly understood. The proposed mechanisms range from physical processes, like direct interaction with high-energy photons and electrons, to biochemical effects, like quenching of the cytotoxic free radical produced by radiation. Supporting the physical interaction hypothesis is the fact that fungal cells deposit melanin as a layer on the inner side of the cell membrane [18, 19], indicating that this layer might provide some kind of shielding. The word shielding is used here to indicate the physical interaction between radiation and melanin, even if the latter resides inside the cell. Equivalent shielding behaviour from cells has already been observed in the case of ultraviolet (UV) radiation: zebra-fish larvae employ an umbrella of melanocytes to protect the sensitive haematopoietic niche from UV-radiation [20] and human keratinocytes uptake synthetic melanin nano-particles and form a supranuclear cap to lessen UV damage [21]. While the UV absorption spectrum of melanin is well documented, the interaction of melanin with high-energy particles has so far only been hypothesized.

More specifically, it has been proposed that the π-electron rich oligomer units that compose melanin dissipate the energy of incident electrons in a controlled way [22]. In addition, it has been suggested that the interaction of melanin with photons through the Compton scattering mechanism—the inelastic scattering of a photon by charged particles—attenuates the photon energy and produces secondary electrons which melanin traps [10, 16]. The physical interactions also supported by the speculations that some fungi are able to utilize ionizing radiation as an energy source for metabolic processes, an ability termed radiotrophism. Radiotrophism was initially proposed by Zhdanova et al., after reporting hyphal growth of various fungi towards sources of radioactivity [7]. The implication of melanin in radiotrophism was suggested by Dadachova et al. [10], after demonstrating that ionizing radiation changes the electronic properties of melanin and that melanized fungi incorporate acetate faster under radiation exposure. We remark that acetate accumulation is indicative of the heterotrophic capability in photosynthetic bacteria. Two more studies provide support to the radiotrophism hypothesis: exposure to ionizing radiation reduces the adenosine triphosphate (ATP) levels only in melanized fungal cells, resembling utilization of ATP during the stage of simple sugar composition in photosynthesis [23]. In addition, significant up-regulation of ribosomal biogenesis genes has been reported in melanized yeast in carbon limited media after irradiation in comparison to melanin-deficient mutant [24], indicating that melanin may be able to harness the energy needed by the ribosomal biogenesis machinery.

Another hypothesis relates the spatial arrangement of melanin inside the fungal cells with its protective properties. In the case of the pathogenic fungus Cryptococcus neoformans, melanin is arranged in a spherical shape covering the inner surface of the cell membrane [18]. This arrangement has been suggested to increase the scattering of incident photons resulting in superior shielding [12]. Melanin hollow particles of roughly spherical shape, also referred to as “ghosts”, can be extracted from C. neoformans cells by digestion with acid, a process which leaves the original melanin structure intact [12, 18]. Irradiation of non-melanized fungal cells with a 137Cs source in the presence of intact C. neoformans ghosts resulted in higher survival rates than in the presence of crushed ghosts [12]. Moreover, higher attenuation of X-rays has been recorded for suspensions of C. neoformans ghosts in comparison to Sepia officinalis melanin, which forms smaller nano-particles [12, 25].

In the present study, we draw inspiration from the biological studies on melanized fungi and examine possible ways to transfer these principles into the design of shielding material for ionizing radiation. First, we test the potential of melanin as a shielding material against β-radiation and we examine the proposed role of melanin in “trapping” secondary Compton electrons. In contrast to previous studies, we separate the physical from any chemical or biological effects and we compare, experimentally and through numerical simulations, the transmitted electron energy through materials that contain or are free of melanin. In a second step, we consider the spatial arrangement hypothesis; we use numerical simulations to understand how the arrangement of different materials in a composite radiation shield affects its performance.

Materials and methods

Chemicals and sample preparation

We obtained the following chemical from Merck: melanin synthetic (M8631), melanin from S. officinalis (M2649) and ammonia solution (25 wt. %). We purchased cellulose nano-crystals (CNC, 12 wt. % aqueous gel) from the University of Maine, USA.

We prepared mixtures of synthetic and S. officinalis melanin in DI water and 2 mol L-1 ammonia solution. We note that the synthetic melanin is soluble to the ammonia solution, whereas all other combinations resulted in suspensions. We started by adding 20 mg of melanin in 1.14 ml of the solvent. In case of the solution, we vortex mixed for 1 min, we allowed for any undissolved particles to sink and we used the total volume of the supernatant for the experiments. For the suspensions, we added the ingredients in standard 15 ml falcon tubes and treated them in an ultrasonic bath (Branson 2510E-DTH, 100 W, 42 kHz) for 30 min.

We prepared CNC solution by diluting 167 mg of 12 wt. % CNC aqueous gel with 0.993 mL DI water, resulting in the same weight fraction as the melanin mixtures. Prior to the shielding experiments, we transferred all of the samples in the standard 12-well microtiter cell culture plates covered by a thin low-density poly-ethylene film to minimize evaporation.

Shielding experiments

We determined the shielding properties of samples using the following procedure; we placed the sample to be tested between the radioactive source and the detector (semiconductor spectrometer Cube 527, quasi-hemispherical 500 mm3 CdZnTe crystal, GBS Elektronik GmbH, Germany), in a vertical configuration as shown in Fig 1A. We irradiated the sample for 1000 sec with the 90Sr β-source. The specific source was selected because its energy spectrum overlaps with the Compton electron spectrum previous studies suggest an increased shielding effect is present (a detailed justification is given in S1 Text and a comparison of the spectra is plotted in S1 Fig). In all cases, the bottom of the well was in contact with the detector and we placed the source at the center of the well just above its edge. The detector was calibrated to the energy range of 10 keV to 2071 keV, which was separated into K = 1020 equally spaced bins. We remark that 10 keV is the minimum detectable energy. The bin counts were recorded by the WinSpec-I software (GBS Elektronik GmbH, Germany). Therefore, from each experiment we obtain a histogram of the deposited energy spectrum on the detector, which we describe by the energy mid-value of the bin, Eb,i, and the number of counts of the bin, ui, with i = 1, 2, …, K. We refer to the total number of counts, , as the histogram content. We estimated the absorbed dose at the detector by adding up all the energy contributions: (1) where md = 2.9 g is the mass of the detector crystal.

thumbnail
Fig 1. Setup schematic and simulation geometry.

(A) Illustration of the shielding experiments. (B) Simulated geometry for the shielding experiments in Geant4.

https://doi.org/10.1371/journal.pone.0229921.g001

To better assess the effect of melanin, we tested all samples against melanin-free controls, specifically wells with the same volume of solvent. In the following, we report the relative absorbed dose, (2) where Dc is the absorbed dose in the case of control. For no relative improvement to shielding RD → 1, whereas for perfect shielding RD → 0. The introduction of RD serves an alternative role; when determining D experimentally, there is an additional variability in the measurements of samples of same composition at different experimental campaigns, in between of which the radioactive source has to be unmounted from the setup for safe keeping. The variability can be attributed to minor changes on the geometry of the setup. The use of RD remedies the situation and allows to compare measurements from difference experimental campaigns, given that D and Dc were acquired during the same campaign.

Spatial arrangement

The original shielding experiments that support the spatial arrangement hypothesis were performed by Dadachova et al. [12] and compared the attenuation of X-rays between suspensions of C. neoformans ghosts and S. officinalis nanoparticles. We expand on the previous idea and we simulated a wider range of arrangements with the 90Sr and X-ray sources. Stated differently, we investigate if there is an optimal arrangement for a composite shield, made out of two materials with given atomic numbers, Z. To distinguish between the two materials, we loosely refer to them as low- and high-Z material. We characterize the composite shield the material volume ratio, RV, and the areal density, ρA; we define RV as the ratio of the volume of high-Z material over the total volume of the shield and ρA as the mass of the high-Z material per unit area of the shield.

For the simulations, we fix the values for ρA and RV, resulting in composite shields with constant mass and constant height per unit area for all possible arrangements. We simulated composite materials of two configurations: layered and mixture. In the layered configuration, which we referred to as “film”, the high-Z material forms a single sheet, which is either sandwiched between the low-Z material or placed at the outer facet of the composite, as shown in Fig 2A. The exact placement is described the relative position parameter hr ∈ [0, 1]. The mixture configuration is constructed by the repetition of a unit cell, a cube with side a of the low-Z material containing the high-Z material in the three arrangements shown in Fig 2B. In the “sphere” arrangement the high-Z forms a sphere with its center coinciding with the center of the unit cell. In the “packed sphere” and “ghost” configurations, the high-Z material is placed in a body-centered cubic arrangement (borrowing the term from atomic crystal characterization) of solid and hollow spheres, respectively. For the “ghost” configuration, the outer radius is equal to and the inner radius is defined by the RV parameter. To facilitate the comparison, we parameterized the lattice geometries by the equivalent radius, Req, the radius of a sphere with volume equal to the total volume of the high-Z material contained in the unit cell.

thumbnail
Fig 2. Schematic of the simulated spatial arrangements.

(A) Illustrations of the film spatial arrangement, at three relative positions (hr = 0, 0.5 and 1). The direction of the incoming radiation is indicated by the arrow. (B) Illustration of the lattice spatial arrangement for three configurations: sphere, packed sphere and ghost.

https://doi.org/10.1371/journal.pone.0229921.g002

We quantify the effectiveness of the difference arrangements by comparing the radiant fluence (3) after the shield, where M and Et,i is the total number and the energy of the transmitted particles and As is the shield area. For consistency, we report the relative radiant fluence, , which is normalized by the fluence Hc of the control sample. We use as control a homogeneous shield made of an ideal mixture of the two materials. Specifically, the density of the mixture is ρm = RV ρh + (1 − RV)ρl, where ρh and ρl are the densities of the high-Z and low-Z components respectively. The weight fraction is equal to for the high-Z and for the low-Z material.

Numerical simulations

We performed Monte Carlo simulations using the software Geant4 [26, 27], to gain insight on the theoretical values for H and D. All simulations were performed using the low-energy electromagnetic model (Livermore library) of Geant4 [28], with validity down to 250 eV.

For comparing with the shielding experiments, we simulated the exact same conditions and geometry, as shown Fig 1B. To avoid any systematic error, we model the 90Sr β-source using the energy spectrum taken from [29], which is plotted in panel (A) of S2 Fig. The source emits electrons at a right cone towards the sample, with apex semi-angle of 15.2°. For the simulations, we neglected the contribution from the infrequent γ-decay of 90Y. We modelled the solution as simple mixtures of elements, with the following simplifying assumption; the addition of melanin does not change the volume of the mixture (same volume as the solvent), but we added the contribution of the mass of the melanin to the density of the solution. For the suspensions we use the same modeling, since the melanin particles did not had the time to precipitate between the sonication and irradiation. In terms of the geometry, we model the sample and the detector (only the semiconductor crystal and the aluminum cap) and we determine the deposited energy, Ed, on the semiconductor crystal. We simulated 5 × 106 primary events and we summed up the contributions from the primary and all the secondaries particles in Ed. We discard events with Ed < 10 keV, which resides outside the range of our detector. The simulated dose was computed as (4) where P is the total number of events recorded by the detector. More information about the settings, geometries and the composition of materials used in the simulations are given in S1 Text.

The simulations to assess how the arrangement of different material affects the shielding properties of the composite were performed in a similar manner; in this case though we recorded the transmitted energy through the composite. We used melanin (synthetic, see S1 Table for elemental composition) and water as the high-Z and low-Z materials respectively. We fixed the geometric parameters of the composite shield at ρA = 3.37 mg cm−2 and RV = 0.234. We simulated various configuration for hr (ranging from 0 to 1) and Req (ranging from 16 nm to 8192 nm) with the 90Sr and the 40 kVp X-ray source. The energy spectrum of the X-ray photons was calculated using the method of interpolating polynomials for a tungsten anode source without filter, as described in [30]. The resulting spectrum is shown in panel (B) of S2 Fig.

Statistical analysis

The dose accumulation and the radiant fluence can be modelled as a compound Poisson process: D and H are the sum of independent and identically distributed (i.i.d.) random variables, where the numbers of terms to be added follows the Poisson distribution. In general, the energy per particle, Ep,i, which is either deposited on the detector or exiting the shield, is the i.i.d. random variable and the number of incident particles in a given time interval follows Pois(λ), with λ the event rate. The probability distribution of Ep,i derives from the distorted spectrum of the radiation source, as it passes through the shielding material. The parameter λ is related to the decay rate of the radioactive source, but it is also affected by the shielding material through the absorption and generation of particles.

In the case of the shielding experiments, the theoretical mean and the variance of D, denoted by μD and respectively, are given by and , where E denotes the expected value. Therefore, Eqs (1) and (4) constitute appropriate estimators of μD, the first one in case the histogram of the energy spectrum is available, whereas the second if each of the individual realizations of Ep,i are registered. The respectful estimators for in case of histograms and individual recording are (5) (6) These estimators also indicate that doses obtained from two separate experiments can be combined by averaging the means and the variances of the two. As P → ∞, as a consequence to the Berry–Esseen theorem [31]. For our calculation of D, the individual terms in the summations are in the order of 106, thus we treat D as normally distributed. We apply similar reasoning and treatment for the calculation of H.

The shielding material can change both the energy spectrum and the rate of the incident particles on the detector. Comparison between different shielding approaches is performed in three steps: first we employ the Anderson–Darling test [32, 33] to detect differences in the energy spectra, then we apply the binomial test to compare the rate of the two event or the content of the two histograms [33]. Finally we apply the Z-test to compare the dose. In all of the previous tests, the null hypothesis states that the random variables were drawn from the same distribution. We accompany the p-values from the hypothesis tests with estimated effect sizes. To quantify differences between two energy spectra we calculate the Kolmogorov–Smirnov score, DKS, namely the maximum absolute difference between the empirical cumulative distribution functions computed from the histogram data [33]. For the rate comparisons, we estimate Cohen’s w [34] and we report the Z-test score (7) for the dose comparisons, where the subscripts 1 and 2 indicate the samples to be compared.

In the case of RD and RH, we follow a similar procedure. As random variables, RD and RH have ratio distribution as the quotient of two normally distributed random variables. A normal approximation of the ratio distribution is valid under some assumptions [35, theorem1]. Specifically, the normal approximation can be justified if both variables in the ratio are strictly positive and the root square sum of their inverse coefficient of variation is smaller than 0.1, which translates to in the case of D. The variables and denote the mean and standard deviation of Dc. In the case of RD, the mean of the approximation is given by Eq (2) and the standard deviation by [35]. These criteria were always valid in our case, therefore we use the Z-test to detect differences in the values of RD and RH. Moreover, we assume that histograms of the same sample that recorded during different experimental campaigns differ mainly on the content and not on the shape. This implies that the shape of histograms between different samples can be compared by the Anderson–Darling test. For the content though, we apply Fisher’s exact test under the null hypothesis that the ratio of the histogram content between sample and control remains constant. For Fisher’s exact test, the associated Cohen’s w effect size is reported.

For the comparison of multiple experiments, first we assess the hypothesis that all incident particles have the same energy spectrum using the k-sample Anderson–Darling test. To compare the rates, in the case of D and H we use Pearson’s χ2-test for goodness of fit under the null hypothesis that the categorical data has the same frequency; each separate experiment constitutes a category and the histogram content the corresponding frequency. In the case of RD and RH, we use χ2-test for independence, where the content of the control histogram is also used. In both cases, we report Cohen’s w. If any of these two tests reject the null hypothesis, we perform pairwise comparisons, similar to the case of two samples. The hypothesis rejection rates for the pairwise comparisons are adjusted using the Bonferroni correction.

In the following, all quantities are reported as the expected value with the 99% confidence interval (CI). The p-values correspond to two-sided tests and are reported as continuous variables. We indicate p-values smaller than 0.001 as <0.001. We signify rejection of the null hypothesis with one, two or three asterisks at confidence level of α = 0.05, 0.01 and 0.001 respectively. We note that for the k-sample Anderson–Darling, the provided interpolating formula for the p-values range from 0.001 to 0.25 and as a result the reported values are restricted to this range.

Results

Melanin shielding from β-radiation

The experimentally measured values for RD for the different mixtures are shown in Fig 3A, along with the numerically simulated values. All samples exhibit measurable RD, namely the deposited energy on the detector was substantially lower in comparison to control. An example of a spectrum recorded during the experiments is shown in Fig 3B for a S. officinalis melanin and water suspension, where it is contrasted to the simulated one. A comparison of the previous spectrum to its control is shown in S3 Fig. All spectra are normalized to represent the estimated probability density function (PDF) for better inspection (area under the curve is equal to 1). In comparison to the simulated one, the recorded spectrum from the experiments has more content on the low energy region; nonetheless, the simulations captured the overall trend and the discrepancy is attributed to neglected detector dynamics, as discussed in detail in S1 Text and S4 Fig.

thumbnail
Fig 3. Shielding effectiveness of melanin mixtures.

(A) Shielding effectiveness of melanin mixtures in ammonia (Amm) and water, for synthetic melanin (Mel#1) and melanin of S. officinalis (Mel#2). For comparison, slurry of cellulose nano-crystals (CNC) was also tested. Experimental measurements are shown in circles and numerical simulation results in squares. (B) Comparison between experiment and simulation of the estimated probability density function of the deposited energy spectrum. Results for the S. officinalis and water suspension.

https://doi.org/10.1371/journal.pone.0229921.g003

Statistical comparison among the experiment recordings revealed difference in both the energy spectra and the rate of incident particles (k-sample Anderson–Darling p < 0.001, χ2-test p < 0.001, Cohen’s w = 0.0094). Similar results were obtained from comparing the simulated dose (k-sample Anderson–Darling p < 0.001, χ2-test p < 0.001, Cohen’s w = 0.0193). Differences were also detected among the experiment controls (k-sample Anderson–Darling p < 0.001, χ2-test p < 0.001, Cohen’s w = 0.1073). The pairwise comparisons for the experimental and simulation data are presented in Table 1. The correspondence between samples and controls, along with the experimental campaign they were recorded in, is summarized in S2 Table. We point out that the reported comparisons are with respect to RD in the case of the experiments and to D in the case of simulations and experiment controls.

For the melanin mixtures, we detected no noticeable influence of the melanin type (synthetic or S. officinalis) on RD, but the observed differences seem to be connected to the solvent. The measured RD for CNC solution is lower than the ones for the water-melanin suspensions, although the simulated value show no detectable discrepancy. For the experimental data, the variations in RD seem to be associated to the deposited energy spectrum, whereas for the simulations the variations seem to stem from the histogram content. In particular, the experimental values of DKS are higher than the simulated ones for comparisons between the same samples, whereas Cohen’s w follows the opposite trend. The highest values for θ were recorded among the experiment controls, confirming that direct comparison of D between different experimental campaigns is not straightforward. In accordance to our assumption, histograms recorded during different experimental campaigns (Ammonia#1 vs Water#1 and Water#1 vs Water#2) have an order of magnitude higher values for Cohen’s w than histograms recorded during the same campaigns (Ammonia#1 vs Water#2), while the values of DKS are comparable for all sample comparisons.

Spatial arrangement of melanin affects shielding effectiveness

Fig 4A presents the simulated RH for the 90Sr source. The results for the X-ray source are shown in S5 Fig, because RH was identical for all arrangements within the estimated uncertainties. To amplify the effect of the spatial arrangement, we repeated the simulations with tungsten (W, high-Z) and poly-ethylene (PE, low-Z) instead, two materials that have been proposed as components for radiation shields [36]. For the W-PE composite, RV was kept as previously and ρA = 38.71 mg cm−2, resulting in a shield with same height per unit area as in the melanin-water case. The results of the simulations are presented in Fig 4B and 4C.

thumbnail
Fig 4. Comparison between different arrangements.

(A) Relative radiant fluence for different spatial arrangements for the melanin-water composite and the 90Sr source. The relative film position is marked on the bottom x-axis and the equivalent radius for the lattice configurations at the top logarithmic x-axis. (B) Relative radiant fluence for the W-PE composite and the 90Sr source. The CI intervals are smaller than the line width. (C) Relative radiant fluence for the W-PE composite and the 40 kVp X-ray source. The CI intervals are smaller than the line width.

https://doi.org/10.1371/journal.pone.0229921.g004

For the melanin-water composite and the 90Sr source, the film configuration showed noticeable trend between RH and hr, whereas among the lattice configurations the effect was negligible. The most favorable configuration with respect to shielding is the film with hr = 0. The results for the W-PE composite with the 90Sr source reveal the tendency of RH to decrease with decreasing hr and Req. Surprisingly, for the X-ray source the different configurations resulted in well separated regions. In this case, the sphere configuration demonstrated RH > 1, practically performing worst that the ideal mixture of the two material. The diminished performance of the sphere configuration can be attributed to the fact that for perpendicular incidence of the incoming photons, there are direct trajectories through the material that do not intersect with the high-Z material. Overall, the film configuration with hr = 0 followed by the ghost with Req ≈ 1 μm appear to provide the highest shielding in all cases.

The absence of detectable differences in RH for melanin-water composite with the X-ray source is not in agreement with the shielding results from the Dadachova et al. experiments. Consequently, we examined if the discrepancy can be attributed to differences in ρA. In the original experiments, three different amounts of melanin were added to 96-well plates, in particular 30 mg, 50 mg and 100 mg. For a typical well diameter of 6.94 mm for the 96-well plate, ρA is roughly equal to 82 g cm-2, 137 g cm-2 and 274 g cm-2. We repeated the simulations for the aforementioned values of ρA for two configurations: melanin ghosts and film with hr = 0. We model ghosts as perfect hollow spheres with an outer diameter of 3 mm and wall thickness of 200 nm, values that closely resemble C. neoformans ghosts [18]. We assume that the ghosts are closely packed with the in-between space filled with water, as in Fig 2B, resulting in RV = 0.234. Table 2 holds the simulated values for RH for the two arrangements, marked with the subscripts f for the film and g for the ghosts, along with p-values from the statistical testing and the estimated effect sizes. In accordance to the previous results, we did not observe any detectable differences in the RH values, the histogram shapes or the histogram content for the X-ray source. In contrary, the ghost arrangement has a substantially lower RH for the 90Sr source up to ρA of 137 mg cm-2, but not for ρA = 274 mg cm−2. The difference in RH values seems to arise from the number of incoming particles on the detector and not from the shape of the energy distribution; θ shows stronger correlation to Cohen’s w than to DKS.

Discussion

Our study investigated possible bioinspired mechanisms to improve current radiation shielding techniques. In addition, through this study we test some of the main speculation regarding the high radiation resistance of the melanized fungus C. neoformans. With respect to interaction of melanin with β-radiation, our experiments indicate that melanin does not exhibit improved shielding; the recorded RD for melanin was comparable to cellulose, a substance with similar chemical composition to melanin. This fact is further supported by the good agreement between experiments and numerical simulations. The physics models in the numerical simulations treat the materials as a mixture of elements or isotopes [26], without taking into account the molecular structure. Thus, the molecular structure of melanin seems to be irrelevant for its shielding capabilities. Moreover, the chemical composition and the pH of the solvent do not affect the shielding capabilities of melanin, as comparable RD was measured for synthetic and S. officinalis melanin in water and ammonia solvents. The interaction of melanin with Compton recoil electrons has previously been speculated to explain the reduced generation of free radical species from the radiolysis of water [22]. Given that melanin acts also as free radical scavenger [37], identifying the exact source of the observed higher survival rates in biological systems poses a challenging task.

Still melanin may be of interest for shielding applications in space, in an indirect way. Fungal-based biocomposites have been proposed as building and shielding materials for habitat structures on the Moon and Mars [38]. The in-situ production of these materials relies on the resistance of the fungi to the extreme radiation environment of space. Hence, melanized fungi are potential candidates for the production of such biocomposites. Besides, we emphasise that this paper does not address the use of melanin as a possible radioprotector, namely a molecule able to reduce the radiation toxicity and mitigate the health risks from human space flight [4].

The simulation results confirm the spatial arrangement hypothesis; the arrangement of melanin in ghosts may provide increased shielding. Although the arrangement of melanin on the cell membrane will marginally reduce the absorbed dose for highly energetic particles at the nucleus of the cell, the effect will be cumulative for the fungal colony consisting of multiple cells. A lattice of melanin ghosts will feature consistent shielding for the colony cells in an anisotropic radiation environment. Most likely the protective effect of melanin can be attributed to various factors; nonetheless, accumulation of melanin in separate structures inside the cytoplasm instead of a diffused state may have a beneficial effect for shielding. Finally, the discrepancy between our simulations and the Dadachova et al. experiments may have arisen from geometrical variations or differences in the modeling and quantity of the solvents.

The superior performance of composite shields has already been identified by previous studies [36, 3941] and can be understood from the perspective of dose enhancement effects taking place at the interface of dissimilar materials [4244]. Most of these studies though approach the problem from another perspective: they either compare the composite shield to aluminum, a commonly used and well-studied material, or they compare different composites materials with one another. Here, we are interested in investigating if there is an optimal geometric arrangement from a radiation shield made of two materials. Our results suggest that spatial arrangement alone is able to reduce the absorbed dose, even up to 50% in comparison to the ideal mixture, just by layering the high-Z material at the incident plane of the radiation. On the other hand, the source of radiation plays an important role; for the 90Sr source, hr and Req play significant role in the dose reduction, whereas for the X-ray source the arrangement itself is the main factor. Moreover, as the shield thickens, which is reflected by an increase in the value of ρA, the effect of the spatial arrangement seem to diminish, as shown in Table 2. Therefore, for composite shield design, the radiation source and the ρA should be taken into account. Although the current study has only examined the reduction in the radiant fluence through the shield, neglecting other design considerations such as manufacturability and cost, it provides an overview on the main design consideration for composite radiation shields.

Supporting information

S1 Text. Detailed description of simulation parameters and β-source selection.

https://doi.org/10.1371/journal.pone.0229921.s001

(PDF)

S1 Fig. Comparison of Compton electron and 90Sr spectra.

Comparison of Compton electron energy spectrum for the 137Cs and the 90Co to the β-spectrum of 90Sr. Compton electron spectra were calculated as as described in [45].

https://doi.org/10.1371/journal.pone.0229921.s002

(TIF)

S2 Fig. Spectra used in Geant4 simulations.

(A) Energy spectrum for the 90Sr source. (B) Energy spectrum for the 40 kVp X-ray source.

https://doi.org/10.1371/journal.pone.0229921.s003

(TIF)

S3 Fig. Example of recorded sample and control spectra.

Comparison of the recorded spectrum for the S. officinalis and water suspension to its control (water only).

https://doi.org/10.1371/journal.pone.0229921.s004

(TIF)

S4 Fig. Simulated 90Sr spectra for two detector models.

(A) Comparison between the simulated detector spectra, using the source spectrum from [29] (Devaney) or the Geant4 radioactive decay module (Geant4), for a detector that adds the contribution of the secondary particles to the primary and the experimentally recorded spectrum (exp). (B) Simulated spectra for a detector that registers each particle separately. The experimentally recorded spectrum is the same as in panel (A).

https://doi.org/10.1371/journal.pone.0229921.s005

(TIF)

S5 Fig. Comparison of different arrangements.

Relative radiant fluence for different spatial arrangements for the melanin-water composite and the X-ray source. The relative film position is marked on the bottom x-axis and the equivalent radius for the lattice configurations at the top logarithmic x-axis.

https://doi.org/10.1371/journal.pone.0229921.s006

(TIF)

S1 Table. Composition of the simulated materials.

Elemental composition used in the numerical simulations for the synthetic and the S. officinalis melanins, and the cellulose.

https://doi.org/10.1371/journal.pone.0229921.s007

(PDF)

S2 Table. Correspondence between samples and controls.

Summary of the samples with the corresponding controls and experimental campaigns.

https://doi.org/10.1371/journal.pone.0229921.s008

(PDF)

Acknowledgments

We thank Alessandra Costantino and Michele Muschitiello from the Co-60 Facility (ESTEC, the Netherlands) for the valuable input in the shielding experiments and the handling of the radioactive sources. We thank Alan Dowson from the Life Science Laboratory (ESTEC, the Netherlands) for his support with the sample preparation.

References

  1. 1. White RJ, Averner M. Humans in space. Nature. 2001;409(6823):1115. pmid:11234026
  2. 2. Durante M, Cucinotta FA. Physical basis of radiation protection in space travel. Reviews of Modern Physics. 2011;83(4):1245.
  3. 3. Cucinotta FA, To K, Cacao E. Predictions of space radiation fatality risk for exploration missions. Life sciences in space research. 2017;13:1–11. pmid:28554504
  4. 4. Cortese F, Klokov D, Osipov A, Stefaniak J, Moskalev A, Schastnaya J, et al. Vive la radioresistance!: converging research in radiobiology and biogerontology to enhance human radioresistance for deep space exploration and colonization. Oncotarget. 2018;9(18):14692. pmid:29581875
  5. 5. Cucinotta FA. Space radiation risks for astronauts on multiple International Space Station missions. PLoS one. 2014;9(4):1–14.
  6. 6. Cucinotta FA, Hu S, Schwadron NA, Kozarev K, Townsend LW, Kim MHY. Space radiation risk limits and Earth-Moon-Mars environmental models. Space Weather. 2010;8:n/a–n/a.
  7. 7. Zhdanova NN, Tugay T, Dighton J, Zheltonozhsky V, Mcdermott P. Ionizing radiation attracts soil fungi. Mycological research. 2004;108(9):1089–1096. pmid:15506020
  8. 8. Dadachova E, Casadevall A. Ionizing radiation: how fungi cope, adapt, and exploit with the help of melanin. Current opinion in microbiology. 2008;11(6):525–531. pmid:18848901
  9. 9. Dighton J, Tugay T, Zhdanova N. Fungi and ionizing radiation from radionuclides. FEMS microbiology letters. 2008;281(2):109–120. pmid:18279333
  10. 10. Dadachova E, Bryan RA, Huang X, Moadel T, Schweitzer AD, Aisen P, et al. Ionizing radiation changes the electronic properties of melanin and enhances the growth of melanized fungi. PloS one. 2007;2:e457. pmid:17520016
  11. 11. Shuryak I, Bryan RA, Nosanchuk JD, Dadachova E. Mathematical modeling predicts enhanced growth of X-ray irradiated pigmented fungi. PloS one. 2014;9(1):e85561. pmid:24454887
  12. 12. Dadachova E, Bryan RA, Howell RC, Schweitzer AD, Aisen P, Nosanchuk JD, et al. The radioprotective properties of fungal melanin are a function of its chemical composition, stable radical presence and spatial arrangement. Pigment cell & melanoma research. 2007;21(2):192–199.
  13. 13. Kunwar A, Adhikary B, Jayakumar S, Barik A, Chattopadhyay S, Raghukumar S, et al. Melanin, a promising radioprotector: Mechanisms of actions in a mice model. Toxicology and applied pharmacology. 2012;264(2):202–211. pmid:22968190
  14. 14. Rageh MM, El-Gebaly RH, Abou-Shady H, Amin DG. Melanin nanoparticles (MNPs) provide protection against whole-body γ-irradiation in mice via restoration of hematopoietic tissues. Molecular and cellular biochemistry. 2015;399(1-2):59–69. pmid:25300618
  15. 15. Schweitzer AD, Revskaya E, Chu P, Pazo V, Friedman M, Nosanchuk JD, et al. Melanin-covered nanoparticles for protection of bone marrow during radiation therapy of cancer. International Journal of Radiation Oncology, Biology, Physics. 2010;78(5):1494–1502. pmid:20421152
  16. 16. Revskaya E, Chu P, Howell RC, Schweitzer AD, Bryan RA, Harris M, et al. Compton scattering by internal shields based on melanin-containing mushrooms provides protection of gastrointestinal tract from ionizing radiation. Cancer Biotherapy and Radiopharmaceuticals. 2012;27(9):570–576. pmid:23113595
  17. 17. Pacelli C, Bryan RA, Onofri S, Selbmann L, Shuryak I, Dadachova E. Melanin is effective in protecting fast and slow growing fungi from various types of ionizing radiation. Environmental microbiology. 2017;19(4):1612–1624. pmid:28127878
  18. 18. Eisenman HC, Nosanchuk JD, Webber JBW, Emerson RJ, Camesano TA, Casadevall A. Microstructure of cell wall-associated melanin in the human pathogenic fungus Cryptococcus neoformans. Biochemistry. 2005;44(10):3683–3693. pmid:15751945
  19. 19. Prados-Rosales R, Toriola S, Nakouzi A, Chatterjee S, Stark R, Gerfen G, et al. Structural characterization of melanin pigments from commercial preparations of the edible mushroom Auricularia auricula. Journal of agricultural and food chemistry. 2015;63(33):7326–7332. pmid:26244793
  20. 20. Kapp FG, Perlin JR, Hagedorn EJ, Gansner JM, Schwarz DE, O’Connell LA, et al. Protection from UV light is an evolutionarily conserved feature of the haematopoietic niche. Nature. 2018;558(7710):445–448. pmid:29899448
  21. 21. Huang Y, Li Y, Hu Z, Yue X, Proetto MT, Jones Y, et al. Mimicking melanosomes: Polydopamine nanoparticles as artificial microparasols. ACS central science. 2017;3(6):564–569. pmid:28691067
  22. 22. Schweitzer AD, Howell RC, Jiang Z, Bryan RA, Gerfen G, Chen CC, et al. Physico-chemical evaluation of rationally designed melanins as novel nature-inspired radioprotectors. PloS one. 2009;4(9):e7229. pmid:19789711
  23. 23. Bryan R, Jiang Z, Friedman M, Dadachova E. The effects of gamma radiation, UV and visible light on ATP levels in yeast cells depend on cellular melanization. Fungal biology. 2011;115(10):945–949. pmid:21944206
  24. 24. Robertson KL, Mostaghim A, Cuomo CA, Soto CM, Lebedev N, Bailey RF, et al. Adaptation of the black yeast Wangiella dermatitidis to ionizing radiation: molecular and cellular mechanisms. PloS one. 2012;7(11):e48674. pmid:23139812
  25. 25. Dadachova E, Casadevall A, inventors; Melanin nanoshells for protection against radiation and electronic pulses. US8586090B; 2014.
  26. 26. Agostinelli S, Allison J, Amako Ka, Apostolakis J, Araujo H, Arce P, et al. GEANT4—a simulation toolkit. Nuclear instruments and methods in physics research section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2003;506(3):250–303.
  27. 27. Allison J, Amako K, Apostolakis J, Araujo H, Arce Dubois P, Asai M, et al. Geant4 developments and applications. IEEE Transactions on Nuclear Science. 2006;53(1):270–278.
  28. 28. Allison J, Amako K, Apostolakis J, Arce P, Asai M, Aso T, et al. Recent developments in Geant4. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2016;835:186–225.
  29. 29. Devaney JJ. Beta spectra of 90Sr and 90Y. Los Alamos National Lab; 1985.
  30. 30. Boone JM, Seibert JA. An accurate method for computer-generating tungsten anode x-ray spectra from 30 to 140 kV. Medical Physics. 1997;24(11):1661–1670. pmid:9394272
  31. 31. Korolev V, Shevtsova I. An improvement of the Berry–Esseen inequality with applications to Poisson and mixed Poisson random sums. Scandinavian Actuarial Journal. 2012;2012(2):81–105.
  32. 32. Scholz FW, Stephens MA. K-sample Anderson–Darling tests. Journal of the American Statistical Association. 1987;82(399):918–924.
  33. 33. Porter FC. Testing consistency of two histograms. arXiv preprint arXiv:08040380. 2008;.
  34. 34. Cohen J. Statistical power analysis for the behavioral sciences. 2nd ed. Lawrence Erlbaum Associates; 1988.
  35. 35. Díaz-Francés E, Rubio FJ. On the existence of a normal approximation to the distribution of the ratio of two independent normal random variables. Statistical Papers. 2013;54(2):309.
  36. 36. Steffens M, Hepp F, Höffgen SK, Krzikalla P, Metzger S, Pellowski F, et al. Characterization of novel lightweight radiation shielding materials for space applications. IEEE Transactions on Nuclear Science. 2017;64(8):2325–2332.
  37. 37. Różanowska M, Sarna T, Land EJ, Truscott TG. Free radical scavenging properties of melanin: interaction of eu-and pheo-melanin models with reducing and oxidising radicals. Free Radical Biology and Medicine. 1999;26(5-6):518–525. pmid:10218640
  38. 38. Lakk H, Krijgsheld P, Montalti M, Woesten H. Fungal based biocomposite for habitat structures on the Moon and Mars. In: 69th International Astronautical Congress; 2018. Available from: http://www.esa.int/gsp/ACT/doc/HAB/ACT-RPR-HAB-2018-IAC-Fungal_Biocomposite.pdf.
  39. 39. Fan WC, Drumm CR, Roeske SB, Scrivner GJ. Shielding considerations for satellite microelectronics. IEEE Transactions on Nuclear Science. 1996;43(6):2790–2796.
  40. 40. Spieth BD, Qassim KS, Pittman RN, Russell DA. Shielding electronics behind composite structures. IEEE Transactions on Nuclear Science. 1998;45(6):2752–2757.
  41. 41. Nambiar S, Yeow JTW. Polymer-composite materials for radiation protection. ACS Applied Materials & Interfaces. 2012;4(11):5717–5726.
  42. 42. Long DM, Millward DG, Wallace J. Dose Enhancement Effects in Semiconductor Devices. IEEE Transactions on Nuclear Science. 1982;29(6):1980–1984.
  43. 43. Murthy M, Lakshmanan A. Dose enhancement due to backscattered secondary electrons at the interface of two media. Radiation research. 1976;67(2):215–223. pmid:781720
  44. 44. Ibarmia S, Eck J, Ivanchenko V, Lavielle D, Rivera A, Cueto J, et al. Experimental dose enhancement in multi-layer shielding structures exposed to high-energy electron environments. IEEE Transactions on Nuclear Science. 2013;60(4):2486–2493.
  45. 45. Blumenthal GR, Gould RJ. Bremsstrahlung, synchrotron radiation, and Compton scattering of high-energy electrons traversing dilute gases. Reviews of Modern Physics. 1970;42(2):237.