Photon attenuation parameters of non-essential amino acids using EPICS2017 library interpolations

Mass attenuation coefficients, effective atomic numbers, electron densities and energy absorption and exposure buildup factors for the non-essential and conditionally non-essential amino acids including alanine, arginine, asparagine, aspartic acid, cysteine, glutamic acid, glutamine, glycine, proline, serine, and tyrosine have been acquired using the latest evaluated photoatomic library of EPICS2017. The library was used by constructing an interpolation script that calculates for all photon attenuation parameters. Comparisons were made using alternative Monte Carlo simulation results for 15 energy points from 59.5 to 1333 keV, and in contrast with experimental works in literature. Good agreements for the mass attenuation coefficients were observed between EPICS2017-based values compared with Monte Carlo code and experimental results. Similar trends for the effective atomic numbers and electron densities were observed from EPICS2017 interpolation and from results found in literature. Conversely, buildup factors acquired by Geometric-Progression fitting parameters were reported in this work preliminarily for most of these biomolecules at different penetration depths. Overall, cysteine showed the most significant deviation among the other non-essential amino acids due to the presence of sulfur in its molecular structure.


Introduction
Protein is essential in various biological processes in the human body and plays an important role in the repair, growth, and maintenance of body cells. Proteins can be metabolized and broken down to their building blocks called amino acids (AAs). Each of these AAs can be subcategorized into essential (EAA) or non-essential (NEAA) based on the ability of the body to synthesize the particular AA [1]. In essence, NEAAs are those which can be produced by the human body, while EAAs are those obtained from diet [2]. With the vast roles of the NEAAs including gene expression, cell signaling, neurological processes, acid-base balance, food digestion and absorption, breakdown of nutrients, metabolism, fertilization and immunity, allosteric and epigenetic mechanisms, and tumorigenesis, researches have gained interest in targeting NEAAs for cancer therapy and field of oncology [3,4].
Radiotherapy in the treatment of malignancy uses high doses of radiation in order to target cancer cells and diminish the size of tumors and slow down the growth and proliferation of cells through DNA damage. To calculate the penetration and dose deposition of X-ray and gamma photons in a particular target, the photon mass attenuation coefficient ( ∕ ), effective atomic number, electron density, and buildup factor are important radiological parameters which can be used [5][6][7]. The former quantities are derived from the ∕ or the cross section. These parameters are significant in the field of radiology and cancer therapy for estimations and optimizations of doses received by the patient.
This work evaluates for NEAAs the photon attenuation parameters including mass attenuation coefficients, effective atomic numbers, electron densities, and exposure and energy absorption buildup factors, derived using EPICS2017 [23] library for interpolations. The new EPICS2017 library is the official evaluated photon and electron library of the latest ENDF/B-VIII [24]. To compare the mass attenuation coefficients obtained from the EPICS2017, the MC method was used as well as available experimental data found in literature. Exposure and energy absorption buildup factors for the NEAAs were reported in this work at multiple penetration depths using Geometric-Progression (G-P) fitting method.

Theory
The mass attenuation coefficient ∕ ( cm 2 g −1 ) for a medium is described by the Beer-Lambert law or Equation (1), where I 0 and I are beam incident and transmitted intensities, respectively, is the density (g cm −3 ) of the medium, and x is the thickness or length (cm) [25].
The law of mixtures shown in Equation (2) gives the total mass attenuation coefficients for a multi-element material, which describes it as a linear weighted combination of the (μ∕ρ) i of its constituent elements, where w i is weight fraction of the ith element [26].
The mass attenuation coefficient is subsequently proportional to the total atomic cross section σ T ( cm 2 atom −1 ) as shown in Equation (3), where N A is the Avogadro's number, n i is the number of atoms in the molecular formula, and A i is the atomic mass of the ith element.
The total atomic cross section for each element (σ T ) i is the sum of the partial cross sections for photoelectric, coherent, incoherent, and pair production, described in Equation (4).
The electronic cross section is a quantity described by Equation (5), where f i is the atom fraction of each element and Z i is the atomic number of each element [27]. It is related to effective atomic number Z eff that is a unitless quantity given by the ratio of total atomic cross section and electronic cross section shown in Equation (6). On the other hand, the electron density N eff (electrons g −1 ) is related to Z eff by Equation (7). Both the Z eff and N eff are widely used for the calculations of doses in multi-element materials by radiation therapy [28,29].
The Beer-Lambert law in Equation (1) is only valid for narrow-beam geometry, thin absorbing materials, and monoenergetic photons. In the case of divergences from such conditions, a correction factor called buildup factor B(E, x) accounts for the effects of scattering described in Equation (8). It is a function of both photon energy and penetration depth in terms of the mean free path (mfp) [30]. The mfp is the inverse of the linear attenuation coefficient μ as shown by Equation (9).
Two classifications of the quantity B(E, x) are referred to as the energy absorption buildup factor (EABF) and the exposure buildup factor (EBF), where the quantity of interests are absorbed energy and exposure, respectively [31]. For both classifications of B(E, x), the ANSI/ANS-6.4.3 [32] contains a compilation of G-P fitting parameters (a, b, c, d, X k ) for 23 elements in the photon energy range between 0.015 and 15 MeV, for depth of up to 40 mfp. The use of G-P fitting parameters is an accurate method of estimating buildup factors for different materials at a broad range of energies [33][34][35].
There are three steps for obtaining the EABF and EBF using the ANSI/ANS-6.4.3 G-P fitting parameters: (1) determination of equivalent atomic number Z eq for the material, (2) evaluation of the G-P fitting parameters for the material, and (3) estimation of EABF and EBF using the parameters [30]. The Z eq represents an equivalent atomic number for multi-element materials. It is obtained by considering the ratio of scattering and total attenuation, which monotonically decreases as a function of atomic number [35]. This is described by Equation (10), where R is the ratio of the ∕ of Compton scattering and the total ∕ , for the material. The R 1 and R 2 are the ratios of two successive elements with atomic numbers Z 1 and Z 2 . The Z eq is used to obtain the G-P fitting parameters of the material. This is done by considering the tabulated G-P fitting parameters of the elements in ANSI/ANS-6.4.3 as illustrated in Equation (11), where F represents each of the five parameters (a, b, c, d, X k ).
Afterward, the term K(E, x) is calculated using the solved G-P fitting parameters through Equation (12). This term represents the photon dose multiplication and change in the shape of the spectrum, with penetration from that at 1 mfp [32]. The buildup factors B(E, x) representing EABF or EBF are then estimated using Equations (13) or (14).

Non-essential amino acids
All NEAAs investigated include alanine, arginine, asparagine, aspartic acid, cysteine, glutamic acid, glutamine, glycine, proline, serine, and tyrosine. The chemical constituents of the biomolecules are shown in Table 1. For alanine, arginine, glycine, and proline, experimentally determined mass attenuation coefficients were previously evaluated by Pawar and Bichile [21]. For asparagine, aspartic acid, glutamine, values were recently evaluated by Kore and Pawar [19], and for serine by Gaikwad et al. [16].

EPICS2017 interpolation
EPICS2017 [23] holds the most recent photoatomic library and is an official library in the latest ENDF/B-VIII release [24]. EPICS2017 is downloadable from the IAEA NDS website in both the ENDL and ENDF data formats. For this work, the version in the ENDF-6 data format was used. The EPICS2017 library was extracted previously for the development of a simple spreadsheet-based program for photon shielding and dosimetry applications [36]. Our previous program was built using MS Excel spline interpolation extensions and additional features via the VBA language. In this present work, we implement the interpolations using an executable Windowsbased program.
It is worth mentioning that the predecessors of EPICS2017 are interpolable using the log-log interpolation law. However, due to several users mis-interpolating by using lin-lin, it was decided to convert the new EPICS2017 into a linearized library. Hence, the interpolation in this work was based on the recommended (linlin) law [37] and was applied to all partial cross sections, including coherent scattering, incoherent scattering, pair production, triplet production, and photoelectric separately for each subshell. The total atomic cross section was obtained by summing all relevant partial cross sections. Notably, due to the linearization, it is now possible to interpolate the total atomic cross sections also available from the library [38].
The modeled geometry is shown in Fig. 1 comprising of a detector covered with a lead shield. This shield was 5-cm thick and had a beam entry diameter of 4 mm for achieving narrow-beam geometry. Each amino acid sample was modeled as a cylinder of 2-cm diameter. Length of each sample was in the range conforming to Nordfors criteria (2 ≤ ln(I 0 ∕I ) ≤ 4) [41]. The surface-to-surface distance between the lead shield and the amino acid sample was fixed at 10 cm. Photon source was modeled as a monoenergetic pencil beam with 4-mm diameter. The starting energies simulated were 59. 5 The [T-Track] tally with energy binning was used to obtain uncollided fluences inside the detector cell. Each run was achieved at a total of 10 7 simulated source photons per photon energy for each NEAA sample.

Results and discussion
An evaluation of the photon attenuation parameters for the NEAAs was made using the EPICS2017 photoatomic library. The results for the mass attenuation coefficients are shown in Tables 2, 3, 4 and 5, which are from EPICS2017 interpolation, from comparative PHITS simulations, and from available experimental values found in literature [16,19,21]. It is seen that the PHITS simulation results agree well with EPICS2017 with only marginally lower values. Nonetheless, the difference is owing to the slight dissimilarity in photoatomic libraries, deviation from narrow   [21] (alanine, arginine, glycine, and proline) and Gaikwad et al. [16] (serine) showed agreement with the theoretical results. Those obtained by Kore and Pawar [19] for asparagine, aspartic acid, and glutamine were slightly lower than the theoretical results but were nevertheless found to be in similar trends.
Effective atomic numbers and electron densities obtained using EPICS2017 interpolations are shown in Figs. 2 and 3. In the lower energy region in which photoelectric effect is dominant, the values are highest with a negative slope toward the Compton dominated region. For Compton dominated region, both plots tend  decrease buildup factors within the medium. EBF and EABF were previously reported for EAAs and some conditional NEAAs [6,42]. For some EAAs, buildup factors were described with an inverse relation to effective atomic numbers Z eff [6]. In this study, the three highest Z eff and three lowest showed this inverse relationship with EBF and EABF. This was assessed from the NEAAs at 1 to 40 mfps. It is highlighted that all NEAAs are chemically composed of H, C, N, and O atoms with the exception of cysteine that also contains S atoms. This contributes to cysteine's larger effective atomic numbers and electron densities at the lower and higher portions of the energy range considered. Consequently, cysteine also demonstrated the lowest values for both types of buildup factors at a large range of photon energies.

Conclusion
Photon attenuation parameters of eleven NEAAs have been determined using interpolation method based on the new EPICS2017 library. The EPICS2017 is the current latest library for electron-photon transport and is officially included in ENDF/B-VIII. The EPICS2017 has been interpolated in this work in the recommended lin-lin law using a Windows-based program. Obtained mass attenuation coefficients of NEAAs were compared with MC method results and experimental values in literatures. An agreement with simulation and experimental values was found. Effective atomic numbers and electron densities were also in agreement with those in literature based on MCNP5 calculations. Buildup factors for energy absorption and exposure quantities have been reported via G-P fitting parameters with the use of the EPICS2017 library. Cysteine was found to have a significantly larger effective atomic number which translates to significantly lower buildup factor values, owing to the S atom comprising its chemical structure. These results can serve in the fields of radiology and medical dosimetry. Furthermore, the interpolation method using EPICS2017 was found to be a useful and recent alternative to the generally used XCOM-NIST.

Conflict of interest
The authors declare no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.