Dimitris G.
Mintis
ab,
Dimitra-Danai
Varsou
bc,
Panagiotis D.
Kolokathis
bc,
Andreas
Tsoumanis
ab,
Georgia
Melagraki
d,
Johannes P.
Seif
e,
Alejandro J.
del Real
ef,
Iseult
Lynch
*bg and
Antreas
Afantitis
*abc
aNovaMechanics Ltd., Nicosia 1070, Cyprus. E-mail: afantitis@novamechanics.com; Tel: +357 99048039
bEntelos Institute, Larnaca 6059, Cyprus
cNovaMechanics MIKE, Piraeus 18545, Greece
dDivision of Physical Sciences and Applications, Hellenic Military Academy, Vari 16672, Greece
eIDENER, 41300 Sevilla, Spain
fDepartment of Systems and Automation, University of Seville, 41092 Seville, Spain
gSchool of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham B15 2TT, UK
First published on 7th July 2025
This work presents the development of LungDepo, a freely accessible online web application hosted on the Enalos Cloud platform (https://enaloscloud.novamechanics.com/proplanet/lungdeposition or https://enaloscloud.novamechanics.com/insight/lungdeposition/), which enables users to simulate particle deposition in the human lung across the head airway, the tracheobronchial and the alveolar (pulmonary) regions based on the aerodynamic diameter of inhaled particles. The LungDepo web-based tool offers two modelling approaches: the International Commission on Radiological Protection (ICRP) model and the multiple path particle dosimetry (MPPD) model, to simulate regional particle deposition and assess inhalation toxicity risk. In addition to modeling the fractional regional particle deposition based on the two implemented models, LungDepo offers the computation of the regional deposited mass (mg) and the computation of the regional lung surface area covered by deposited particles (m2) providing statistical analyses of the relative contributions of coarse (2.5–10.0 micron), fine (<2.5 micron) and ultrafine (<0.1 micron) particles. The web-based tool includes predefined modelling scenarios for various particle-bound substances or co-pollutants such as polycylic aromatic hydrocarbons (PAHs) or per- and polyfluoroalkyl substances (PFAS) and micro-sized engineered materials, along with provision for the inclusion of user-defined particle size distributions. The integration of an application programming interface (API) and the development of an intuitive graphical user interface (GUI) facilitate data exchange and integration with external environmental and regulatory toxicity models by simplifying complex tasks and broadening the tool's applicability for diverse users in environmental research, monitoring, and regulatory activities.
Environmental significanceEvery human breath may expose individuals to a toxic cocktail of ultrafine particles and engineered nanomaterials that are silently deposited in the human lungs. In this study, LungDepo is introduced as a free-to-use web application that uses dosimetry models to quantify the deposition of airborne particles within the human respiratory tract. The tool provides statistical analyses of the contributions of coarse (2.5–10.0 μm), fine (<2.5 μm), and ultrafine (<0.1 μm) particles. Designed for non-programmers, LungDepo is accessible to a wide audience and can potentially support broader regulatory frameworks for assessing human inhalation toxicity and safeguarding public health. |
Understanding the deposition of inhaled particulate matter in the different regions of the respiratory tract is crucial for evaluating, quantifying and predicting inhalation toxicity. The chemical composition of PM is primarily carbonaceous with contributions from inorganic ions (including sulphates, nitrates and metal oxides) and volatile and semi-volatile organic compounds such as polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs), and per- and polyfluoroalkyl substances (PFAS),27–29 which can alter biological activity and contribute to carcinogenic, neurotoxic, and endocrine-disrupting effects.30–34 Fine and ultrafine particles remain airborne longer, penetrate deeper into the lungs, and often carry highly reactive transition metals and oxidized organics, enhancing their toxicity.35 Furthermore, microparticles or micro-sized engineered materials and their agglomerates with aerodynamic diameters between 0.1 and 100 μm such as zinc oxide (ZnO),36,37 iron oxide (Fe2O3),38 titanium dioxide (TiO2),39 aluminum oxide (Al2O3),40 copper oxide (Cu2O),41 cobalt oxide (Co3O4)41 and nickel hydroxide (Ni(OH)2)42 may pose additional inhalation hazards due to their ability to penetrate deep into the respiratory tract. Engineered micro-sized materials are applied in various industrial and consumer products.39 Consequently, it is crucial to evaluate the effects of inhalation exposure, particularly in terms of pulmonary toxicity, that may arise as a result of their very small dimensions.39
From a morphological perspective, the human respiratory system consists of two distinct structural regions. The first consists of the extrathoracic airways, termed as the head airways (HA) consisting of the nose, mouth, pharynx and larynx. The second consists of the tracheobronchial (TB) airway trees and the alveolar (AL) regions43 (Fig. 1D). Deposition in the HA and TB regions act as barriers to safeguard the AL region, which constitutes the air–blood barrier, from potential irritation or harmful particles.44 In the AL region, where the airway walls lacks a mucus layer, deposited insoluble particles are effectively cleared by the alveolar macrophages.43 However, soluble particles deposited in the AL region may dissolve and pass through the alveolar membranes, entering the systemic circulation.43,45 Numerous factors affect the regional deposition of particles within the human lungs, primarily dependent on an individuals' pulmonary physiology such as breathing patterns and lung geometry.46,47 Additionally, the physicochemical characteristics of inhaled particles, such as their size, size distribution, density, shape, charge, surface properties, hygroscopicity (ability to absorb moisture from the surrounding air) and the interactions between particles and pulmonary surfactant play a critical role.35,48–50 Of these factors, the size of inhaled particles is considered as the key predominant factor affecting the fate and regional deposition of particles in the human lungs.51 Coarser size particles (PM10) tend to lodge in the throat and/or the bronchi, whereas fine and ultrafine particles can travel deeper into the airways of the respiratory tract, reaching the alveoli, binding with proteins that support (opsonising) or hamper (disopsonising) recognition by macrophages and influence crossing of anatomical barriers, and penetrating into the bloodstream.2,8,52
The experimental determination of particle deposition within the human respiratory tract poses significant challenges in both in vivo and in vitro studies. In vivo studies are ethically problematic, while in vitro models have significant limitations in accurately replicating realistic physiological conditions.53 To overcome these limitations, particle dosimetry models have been extensively developed over the past decades to calculate the regional deposition patterns of inhaled aerosols.54–56 It is still challenging however, to accurately measure and quantify analytically the deposition of particles in the respiratory tract due to the intricate interplay between the structural morphology of lungs (that differs between individuals) and the hydrodynamic flow field within the airways, which constantly undergoes dynamic changes.44 According to Hofmann,57 existing particle dosimetry models are categorized based on lung morphometry and mathematical modelling techniques into semi-empirical models,58–60 one-dimensional cross-section or “trumpet” models,61–64 mechanistic models,65–70 and computational fluid and particle dynamics (CFPD)-based models.71–75 Among the most widely used particle dosimetry models are two semi-empirical models that combine (first-principle) mechanistic frameworks with experimental (empirical) data:57,58,76 the ICRP66 model developed by the International Committee of Radiological Protection (ICRP) in 1994 (ref. 58) (which has since been updated to ICRP130),59 and the NCRP model developed by the National Council on Radiation Protection and Measurements (NCRP) in the USA, in 1997.60 Additionally, mechanistic models such as the deterministic asymmetric branching model, also known as the multiple pathway particle dosimetry (MPPD) model, introduced in 2001,77,78 and the stochastic asymmetric branching model, referred to as IDEAL (inhalation, deposition and exhalation of aerosols in the lung), developed between 1990 and 1992,69,70,79 which derive deposition fractions directly from classical flow equations. The clearance mechanisms of particles are not discussed here, as they are beyond the scope of the present work.
The objective of this work is to introduce a publicly accessible web application, LungDepo, designed to facilitate users including policy-makers, public health professionals and risk assessors in employing the ICRP and MPPD models for the quantification of particle deposition within the three anatomical compartments of the human lungs: the head airways (HA), the tracheobronchial region (TB), and the alveolar region (AL). Designing and developing a web application in the context of inhalation toxicity is essential for providing accessible and efficient tools that can be utilized by a wide range of users without the need for advanced mathematical or programming skills. Well-designed graphical user interfaces (GUI) are crucial for enabling users with limited technical expertise to navigate and apply complex models for assessing inhalation risks. This accessibility promotes broader use of the models, enabling researchers, healthcare professionals, and policymakers to efficiently conduct evaluations and make informed decisions. Additionally, incorporating application programming interfaces (API) capabilities enhances the flexibility of the web application, allowing it to integrate with external datasets, automate processes, and support the implementation of customized workflows. This facilitates collaboration across different research domains and accelerates innovation in the design of safer materials and substances and evaluation of the potential impacts of mitigation measures to reduce exposures to particles. LungDepo is hosted on the Enalos Cloud platform facilitating its integration with other freely accessible web applications within the framework of nano-informatics, including several human health related models such as the integrated approach to testing and assessment (IATA) for lung exposure and toxicity.80–82 The Enalos Cloud platform also offers access to several other models including NanoBioAccumulate,83 which simulates nanomaterial uptake and elimination in aquatic and soil invertebrates; MicroPlasticFate,84 which provides the dynamic and steady-state modeling of the environmental fate of micro- and nanoplastics across key environmental scales and compartments; and UANanoDock,85 which predicts protein adsorption onto nanoparticles. Beyond nano-informatics, the Enalos Cloud platform also hosts web applications in cheminformatics,86,87 and machine learning,88,89 providing a versatile infrastructure for developing, deploying, and accessing predictive models in materials and life sciences.
The LungDepo web application presented in this study is specifically designed to offer a user-friendly GUI that supports the application of complex models for computing the mass deposition (of particles themselves or based on concentrations of adsorbed co-pollutants), the lung surface area covered by deposited particles and the relative contributions of inhaled particulates classified by their size (coarse, fine and ultrafine) in the various regions of the human lungs. As regulators and industry increasingly shift towards the implementation of safe and sustainable by design (SSbD) principles for the design of chemicals and materials, the LungDepo web application provides valuable insights into how new materials can be engineered to reduce harmful effects on human health as well as visualizing the impacts of altering the particle size distributions through specific emissions mitigation measures. Additionally, LungDepo features predefined scenario models to simulate the deposition of several airborne substances and micro-sized engineered materials in the human lungs depending on the breathing frequency and concentration of airborne particles.
This model provides calculations for the deposition efficiency of inhaled particles across three distinct regions of the respiratory tract: the head airways (HA), the tracheobronchial region (TB), and the alveolar region (AR). The deposition efficiency of particles across these three distinct regions is calculated as follows:
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
Fig. 2, presents the analytical calculation of the deposition efficiency (fraction) in the three regions of the human lung as a function of the aerodynamic particle size, ranging from 0.001 μm to 100 μm, using the simplified equations proposed by Hinds44 for the ICRP model. The x-axis, which represents the aerodynamic diameter of particles, is displayed on a logarithmic scale. As shown in Fig. 2, the ICRP model suggests that particles with aerodynamic diameter of less than 10 μm are inhaled at nearly 100%. In contrast, particles larger than 10 μm in diameter exhibit a significantly reduced inhalable fraction, becoming lodged within the head airways of the human lung. Ultrafine particles with aerodynamic diameter of less than 0.1 μm (PM0.1) exhibit significant deposition efficiencies, reaching up to 50% in the alveolar region and approximately 30% in the tracheobronchial regions.
![]() | ||
Fig. 2 Deposition efficiency (fraction) computed based on the simplified equations proposed by Hinds44 for the ICRP model, plotted against the aerodynamic particle diameter in the range 0.001–100 μm. |
Fig. 3, presents the analytical calculation of the deposition efficiency (fraction) for particles deposited in the three regions of the human lungs using the MPPD model version 3.04, developed by Applied Research Associates Inc.78,93 based on the conditions described in the previous paragraph. The analysis reveals three notable discrepancies between the outcomes of the ICRP model and the MPPD model. Firstly, the ICRP model estimates the AL deposition efficiency for ultrafine particles to be approximately 50%, whereas the MPPD model yields a lower estimate of approximately 30%. Secondly, regarding the deposition efficiency for ultrafine particles in the TB region, the MPPD model estimates the efficiency to be approximately 50%, contrasting with the ICRP model's lower estimate of around 30%. According to the ICRP model, the inhalation efficiency for particles with an aerodynamic diameter up to 100 μm declines to approximately 60%. In contrast, the MPPD model suggests that all particles within this diameter range (100%) continue to be inhaled.
![]() | ||
Fig. 3 Deposition efficiency (fraction) computed using the MPPD model version 3.04, developed by Applied Research Associates Inc.,78,93 plotted against the aerodynamic particle diameter in the range 0.001–100 μm. |
![]() | (5) |
The total deposited surface area of particles, Adeposition, is calculated as:
![]() | (6) |
The number Ni of particles deposited for a given particle aerodynamic diameter, Dp,i is given by:
![]() | (7) |
![]() | (8) |
The surface area, Ai, assuming the particles are spherical, is given by:
Ai = π·D2p,i | (9) |
Fig. 4 presents the graphical user interface (GUI) of the LungDepo web application. In this example, the chosen model is the ICRP, accompanied by a predefined scenario for the particle size distribution of PFAS in particulate matter collected from a coastal area, as reported in the study by Lin et al.91 Once the model and scenario related to particle size distribution are selected, two intuitive graphical plots are generated. The first plot represents the normalized probability density of the aerodynamic diameter size distribution of the particulate matter. The second plot shows the deposition efficiency, expressed as a fraction, across the relevant range of particle sizes. In the example shown, the concentration of exposure to particulate matter is set at 6.644 × 10−7 mg m−3, with a respiratory volume rate reflective of a female in a sitting position, set at 0.39 m3 h−1. The exposure duration is defined as 3600 seconds. By selecting the ‘enable advanced options’ check box, the user can specify the density of particles, which is set to 1 g cm−3 in this example. Also, the user can define the total surface area of the HA, TB and AL regions, measured in m2. Upon selecting the ‘calculate dose’ button, users can obtain the lung regional deposited doses of particles, measured in milligrams (mg) and as a percentage (%). Here, the percentage represents the fraction of the total deposited mass that occupies the surface area of a specific lung region relative to the total mass deposited across the entire lung. Additionally, users can retrieve the deposited dose of particles per regional lung surface area (mg m−2), the surface area (of the lung region) covered by deposited particles (m2) and the ratio of the surface area covered by deposited particles to the total regional lung area, expressed as percentage (%). The output also includes the relative contributions of particles classified by size in the three different regions of the human lung, measured in percentage. This classification includes coarse particles larger than 2.5 μm, fine particles less than 2.5 μm, and ultrafine particles smaller than 0.1 μm. Additionally, there is an option to download the results via the ‘Download Dose Data’ button.
The user may choose from a dropdown menu that includes the ICRP and MPPD models. Additionally, for the scenario pertaining to particle size distribution, users can select from a range of predefined examples representing various substances that may be constituents of particulate matter. These substances include91 perfluorobutanoic acid (PFBA),91 perfluorobutane sulfonate (PFBS),91 corn starch,97 chitosan,98 glycerol,99 zinc oxide (ZnO),100 acetic acid,101 silane-based102 and siloxane-based103 substances, 2-octenylsuccinic anhydride,104 sodium alginate,105 alkyl ketene dimers (AKDs)106 and polycyclic aromatic hydrocarbons (PAHs)107 substances as reported in previous studies. Users have the option to upload a custom scenario as well pertaining to the size distribution of the substance of interest within particulate matter. Additionally, for the respiratory volume rate, users may choose from predefined scenarios established by the ICRP model, which takes into account gender as well as varying activity levels, such as sitting, light exercise, or heavy exercise.43 To support broader applicability, LungDepo includes a step-by-step user guide, which is directly available within the web application interface, explaining how to upload a custom particle size distribution in .txt format and how to manually define all relevant input parameters, including particle concentration, respiratory volume rate, exposure duration, particle density, and lung region surface areas to ensure that users can apply easily LungDepo to a wide variety of compounds and exposure scenarios.
The LungDepo web application incorporates a representational state transfer (REST) application programming interface (API) that enhances both its functionality and usability. This is shown in Fig. 5. The functionality of the ‘POST’ and ‘GET’ requests for the endpoints https://enaloscloud.novamechanics.com/proplanet/apis/lungdeposition, https://enaloscloud.novamechanics.com/proplanet/apis/lungdeposition/scenario and https://enaloscloud.novamechanics.com/proplanet/apis/lungdeposition/respiratoryVolumeRate, respectively, were tested using Postman (see Fig. S1, S3 and S4†). The same tests were performed using Swagger. The API documentation is publicly accessible via the Swagger user interface (see Fig. 5) at https://enaloscloud.novamechanics.com/proplanet/swagger-ui/index.html#/ and the underlying specification at https://enaloscloud.novamechanics.com/proplanet/apis/swagger.json. Fig. S2† demonstrates the use of Swagger to perform a ‘/lungdeposition’ POST request, using the same input data as shown in Fig. S1† (which used Postman), confirming consistency of the output across tools. This testing ensured that the expected responses were returned for each request type. In section S2 of the ESI,† robustness tests for the ‘/lungdeposition’ POST query are presented by submitting extreme input scenarios. Specifically, two high-value and two low-value test cases were constructed using exaggerated concentrations, nanoparticle densities and tidal surface areas across the alveolar, tracheobronchial, and head airway regions, using both the ICRP and MPPD models implemented in LungDepo. All tests were performed using Postman and, in every case, the API successfully returned results, confirming the robustness and reliability of its implementation under a wide range of input conditions. A REST API allows the LungDepo web application to interact with external systems and services in a standardized way, facilitating the exchange of data between the application and other tools, databases, or platforms. This means that users can access the application's features programmatically, enabling integration with broader research ecosystems and allowing for automation of tasks, such as data input, analysis, and reporting. Additionally, the REST API provides scalability and flexibility, allowing different users – whether researchers, developers, or public health professionals – to customize their workflows by integrating the web application into their existing systems. For instance, researchers can automate the analysis of large-scale datasets to simulate particle deposition under various conditions or integrate LungDepo with other modeling tools to enhance the breadth of their inhalation toxicity studies.
![]() | ||
Fig. 5 LungDepo is available through a REST API to enable programmatic access and integration with other models, tools and into an IATA for hazard and risk assessment. |
Furthermore, a REST API facilitates collaboration by enabling multiple users to access the same application from different locations or systems, allowing for shared research efforts and cross-platform compatibility. This makes the application more adaptable to different research needs, supporting a wide range of use cases – from routine regulatory assessments to advanced scientific research in inhalation toxicology. By providing a flexible and interoperable interface, a REST API significantly expands the reach and impact of the LungDepo web application, helping users make more informed decisions in evaluating inhalation toxicity, evaluating the impact of exposure mitigation measures such as air filtration, and designing materials by application of the SSbD principles. Integrating the LungDepo web application into the Enalos Cloud Platform further enhances its capabilities by allowing users to access other web applications hosted on the same platform, including the integrated approach to testing and assessment (IATA) for lung exposure and toxicity.80–83 This integration fosters a collaborative environment where different tools can be used together to provide comprehensive insights and support multifaceted research efforts.
Lin et al.91 conducted extensive work quantifying the (normalized) particle size distribution of particle-bound PFAS and the sum concentrations of gaseous and particulate PFAS emitted from various reference sites, including wastewater treatment plants (WWTPs), landfills, coastal areas and natural reserve sites. They investigated collectively 49 PFAS and reported the particulate concentrations of each, measured in picograms of PFAS per cubic meter in the air (pg m−3). The total concentrations were designated as Σ PFAS. The normalized distribution of total PFAS (Σ PFAS) on the different size fractions in the particulate matter, along with the total particulate concentrations of PFAS from different sampling locations as reported by Lin et al.,91 have been included as predefined scenarios in the LungDepo web application. Among the predefined scenarios the normalized distributions of PFBA and PFBS are included, as obtained from the study of Lin et al.91
Fig. 6A, presents the predicted mass deposition rate (pg h−1) of inhaled PFAS at or near a landfill site using the ICRP and MPPD models in LungDepo as compared to those calculated by Lin et al.91 also using the ICRP model, while Fig. 6B presents the corresponding predicted mass deposition percentage (%) within the three distinct regions of the human lung calculated from the LungDepo model (this was not calculated by Lin et al.,91 and is one of the new features of LungDepo). The corresponding figures for predicted exposures based on PFAS emissions at coastal and natural reserve sites and at WWTPs are provided in the supplemental information file (Fig. S9–S18†). These calculations were performed using the ICRP and MPPD dosimetry models integrated within the LungDepo web application, alongside the findings reported in the study by Lin et al.91 to validate and verify the accuracy of the calculation outputs generated by the LungDepo web application. Lin et al.91 used the simplified equations proposed by Hinds,44 similar to our approach in this study, to calculate the mass deposition rate of inhaled PFAS based on the ICRP model. As can be observed in Fig. 6A, the mass deposition rate calculated by the LungDepo web application is consistent with that reported by Lin et al.91 thereby verifying the calculations performed in relation to the ICRP dosimetry model. The total mass of inhaled PFAS deposited, as estimated by both the ICRP and MPPD models, is approximately 110 pg h−1. Specifically, the ICRP model predicts that the contributions to the HA, TB and AL region are roughly 77%, 5% and 18%, respectively. The MPPD model predicts that the contributions to the HA, TB and AL region are around 66%, 13% and 21%, respectively. These discrepancies are attributed to the different deposition efficiencies predicted by the two models, as shown in Fig. 2 and 3, respectively. The particle size distribution of the PFAS, as reported in the study by Lin et al.,91 spanned from 0.08 to 25 μm. The distribution exhibited a peak at an aerodynamic diameter of 0.78 μm that corresponded to an inhalation probability of around 17%. It is evident from Fig. 2 and 3, that for particles with an aerodynamic diameter of approximately 0.78 μm, the ICRP model predicts slightly higher deposition efficiency in the HA, estimated at around 8%, compared to the 5% predicted by the MPPD model. Conversely, the deposition efficiency in the TB region, according to the ICRP model, is lower at approximately 1.7%, whereas the MPPD model predicts a value of around 5%. The differences predicted between the MPPD and ICRP models are minor and reflect the improved lung geometry in the MPPD model.
![]() | ||
Fig. 6 Predictions of PFAS mass deposition at a landfill site in Hong Kong (based on input data from Lin et al.), computed using the ICRP and MPPD models integrated within LungDepo. (A) Deposition in pg h−1 compared with Lin et al.; (B) deposition in % of Σ PFAS, comparing ICRP and MPPD results.91 |
The total mass deposition and the contributions of different particle sizes of the particles to which the PFAS are bound, as calculated here using input data (on PFAS concentration and particle size distribution) retrieved from the study of Lin et al.,91 are found to be consistent and comparable to the findings reported in the study of Guo et al.90 Guo et al.,90 reported that the mass deposition of eight PFAS within the human lungs, in the urban atmosphere of Shanghai, China is approximately 361 pg h−1, indicating that the coarse PFAS particles contribute the most to the particle and PFAS load in the HA region. They also reported that PFAS particles with an aerodynamic diameter of less than 2.1 μm account for approximately 72.5% of the particle deposition in the AL region. This observation is comparable with the results presented here that demonstrated that fine and ultrafine particles contribute approximately 75–80% to the total deposition in the AL region.
The LungDepo web application offers further advanced analytical features for calculating the regional percentage contributions of the three distinct particle size categories (coarse, fine and ultrafine particles) as well as the computation of the deposited dose per regional area (mg m−2) and evaluation of the percentage of the regional surface area on which particles have deposited relative to the total regional area (as a percentage) which can also be thought of as lung surface area in which effective lung function (e.g., oxygen exchange) will likely be reduced. As can be seen in Fig. 7, both ICRP and MPPD models predict that coarse particles predominate in the HA region. In more detail, the MPPD model predicts this percentage to be slightly higher at around 89% of the total particle load, whereas ICRP predicts that coarse particles constitute around 74%. In the TB region, the MPPD model indicates that fine particles predominate corresponding to around 55% of the total particle load, whereas ICRP model indicates equal contribution between fine and coarse particles. In the AL region, both models predict that fine particles predominate this region with percentages around 60–65% and ultrafine particles contributing to around 15% of the total particle load here.
![]() | ||
Fig. 7 Contributions (expressed as percentages of the Σ PFAS) of inhaled Σ PFAS associated with particles of different sizes (collected from a landfill site in Hong Kong by Lin et al.91) as calculated using the ICRP and MPPD models integrated within the LungDepo web application. |
The predicted deposition dose per lung regional area of people in close proximity to the landfill, using the ICRP model, is found to be 5.03 × 10−5 mg m−2 in the HA region, 2.19 × 10−7 mg m−2 in the TB region and 2.34 × 10−9 mg m−2 in the AL region. The corresponding computations using the MPPD model yielded values of 4.24 × 10−5 mg m−2 in the HA region, 5.84 × 10−7 mg m−2 in the TB region and 2.80 × 10−9 mg m−2 in the AL region. Although the predicted particles' coverage in the AL region is relatively low, in the order of 10−9 mg m−2, as observed by both models, the computed mass of PFAS deposited in the AL region, which is found to be 20–25 pg h−1, should be of concern regarding potential toxicological effects and health implications. While there are as yet limited health-based guidance values or limit values for individual or total PFAS concentrations in air, some suggestions are that annual exposure to PFOA (for example) should be limited to 0.0053 μg m−3.108 The European Food Safety Authority has advised a tolerable weekly intake (TWI) for total PFAS via food of 0.63 ng kg−1 from food and water for example.109 In parallel, the U.S. Environmental Protection Agency (EPA) has set maximum contaminant levels (MCLs) in drinking water at or near zero for several PFAS compounds, including PFOA and PFOS, reflecting increasing concern over their health impacts and a shift toward more protective regulatory thresholds.110
Using the probability density size distribution as reported by Lv et al.107 and the concentrations from the sum of 12 PAHs measured by Voliotis et al.29 during cold (January–March 2013) and warm (May–July 2013) seasons at an urban traffic site in Thessaloniki (northern Greece), the mass deposition rate (ng h−1) of particle-bound PAHs is calculated. The respiratory volume rate is set to 1.7 m3 s−1 to be consistent with the study of Voliotis et al.29 This volumetric rate corresponds to light activity levels in adults. The mass deposition rate is computed using the ICRP and MPPD models integrated within the LungDepo web application. The results were compared to those reported by Voliotis et al.,29 who utilized the MPPD model with the stochastic lung option, applying a functional residual capacity (FRC) of 3389 mL. In contrast, the MPPD model integrated within the LungDepo web application is fitted based on a Yeh/Schum symmetric lung geometry model with an FRC of 3300 mL. Nevertheless, as presented in Fig. 8A and B respectively, the predicted outcomes from the MPPD model using LungDepo and the outcomes predicted by Voliotis are in close agreement during both cold and warm seasons. This comparison further assures the validity of the implementation of the MPPD model integrated within the LungDepo web application. The predicted outcome based on the ICRP model overestimates the mass deposition of particle-bound PAHs. This discrepancy as discussed previously is attributed to the different computations of deposition efficiency as evidenced in Fig. 2 and 3.
![]() | ||
Fig. 8 Predictions of the mass deposition in human male adult lung measured in ng h−1 of PAHs during the cold season (A) and warm season (B) at an urban traffic site as computed using the ICRP and MPPD models integrated within LungDepo and compared against the previous study29 (whose data was used as the input to the LungDepo model).29 |
As can be seen in Fig. 8A and B, the mass deposition of particle-bound PAHs during the cold season is significantly higher, predicted to be approximately 18–20 ng h−1, compared to the warm season, where deposition levels are predicted to be around 6–8 ng h−1. This observation aligns with previous studies that reported higher concentrations of PAHs during cold seasons due to the increase from seasonal emission sources such as domestic heating etc.118–122 and lower concentrations of PAHs during warm seasons due to the rise of temperature and concentrations of atmospheric oxidants such as ozone and radicals.123–126
A case study presented by Tsiros et al.131 is used to examine the robustness and accuracy of the LungDepo web application to predict the mass of TiO2 deposited in the three different human lung regions, following the exposure of female and male workers to 22 nm ultrafine TiO2 particles under light exercise during an 8 hour work shift. Tsiros et al.,131 employed the ICRP model using the Hinds approach,44 and the MPPD model for assessing the risk exposure of female and male workers exposed to a concentration of 5.85 mg m−3 of TiO2 ultrafine particles, considering a particle density (assuming spherical particles) of approximately 4.26 g cm−3 and a count median diameter of the particles (CMD) of 0.02 μm. Fig. 9A and B, present the mass deposition measured in mg calculated using the ICRP and MPPD models integrated within the LungDepo web application, following the exposure of a male or female, respectively to 22 nm ultrafine TiO2 particles during an 8 hour work shift under light exercise. The ICRP model predicts that the largest mass of deposited TiO2 microparticles is in the AL region, estimated at approximately 33 mg for males and 28 mg for females. This finding underscores the significant risk of pulmonary inflammation that may result from the 8 hour (workday) exposure to TiO2 ultrafine particles, to both males and females. This observation is in perfect agreement with the calculations conducted by Tsiros et al.,131 using the ICRP model. Using the MPPD model integrated within the LungDepo web application, it is observed that the predicted mass of TiO2 deposited in the alveolar (AL) region is also considerable, approximately 20 mg for males and 16 mg for females however, the largest accumulation of TiO2 is predicted to be in the TB region. This discrepancy between the two models can be described through the calculation of deposition efficiency for particles measuring 0.022 μm, as illustrated in Fig. 2 and 3. It is noteworthy to state that although Tsiros et al.,131 used the same lung geometry for MPPD model as used in this study, they considered variable exposure (30 minutes of active emissions per hour) in contrast to the MPPD model used in this study which was fitted based on a constant exposure scenario. This distinction largely explains the differences observed in the predictions yielded by the MPPD model integrated within the LungDepo web application and the predictions yielded by the MPPD model utilized by Tsiros et al.131
Using the advanced functionalities of the LungDepo web application, it has been observed that both the ICRP and MPPD models predict the surface area occupied by deposited particles per the regional area to be around 2–3% in the AL region, around 250% in the TB region and around 2500% in the HA region. These findings clearly indicate that inhalation of 22 nm ultrafine TiO2 particles for 8 hours work shift, could potentially lead to significant pulmonary inflammation and other adverse health implications for both males and females.
Future extensions of LungDepo could include parameterisation of different lung sizes/breathing rates and functional residual capacities to cover also children and vulnerable populations (e.g., pregnant women, the elderly, people with asthma, emphysema or other lung-related conditions, and obesity) which will influence both exposure and susceptibility to additional health effects arising from air pollution.132,133
The enhanced user interaction and flexibility of the LungDepo web application is affirmed by its user-friendly GUI, which facilitates users with no programming knowledge to perform calculations using the ICRP and MPPD models for prediction of the mass deposition of inhaled particles and aerosols in the human lungs. The tool also features REST API integration that facilitates interaction with other web applications and software such as for integrated exposure, hazard and risk assessment, thereby promoting efficient data exchange and enhancing interoperability. This integration opens new opportunities for collaborative research and a thorough evaluation of inhalation toxicity, ultimately contributing to the design and development of inherently safer materials through application of the SSbD principles.
In this study, clearance mechanisms of deposited particles were not considered. To more accurately assess the risk from inhaled particles to human lungs, future work will integrate comprehensive physiologically-based pharmacokinetic (PBPK) models within LungDepo, such as the one presented in the work of Tsiros et al.131 The inclusion of additional variants of MPPD models fitted to alternative lung geometric models, such as the Yeh/Schum 5-Lobe, the stochastic lung, the Weibel etc. will also be considered. Finally, the variable exposure conditions featured in the MPPD model version 3.04, developed by Applied Research Associates Inc.78,93 will also be integrated within LungDepo in future upgrades.
Footnote |
† Electronic supplementary information (ESI) available: ESI file is included, providing a comprehensive analysis of the outcomes generated using the LungDepo web-tool. See DOI: https://doi.org/10.1039/d5en00299k |
This journal is © The Royal Society of Chemistry 2025 |