Abstract
The knowledge on the deposition and retention of the viral particle of SARS-CoV-2 in the respiratory tract during the very initial intake from the ambient air is of prime importance to understand the infectious process and COVID-19 initial symptoms. To give some light on that, we propose to use a modified version of a widely tested lung deposition model developed by the ICRP, in the context of the ICRP Publication 66, that provides deposition patterns of microparticles in different lung compartments. In the model, we mimicked the “environmental decay” of the virus, determined by controlled experiments related to normal speeches, by the radionuclide 11C that presents comparable decay rates. Our results confirm clinical observations on the high virus retentions observed in the extrathoracic region and the lesser fraction on the alveolar section (in the order of 5), which relevance is a subject to be investigated.
1 Introduction
Respiratory infections are relevant clinical conditions due to their diffusion and potentially severe consequences, such as presently observed for the SARS-CoV spread all around the world. Declared as a “global pandemic” by the World Health Organization (WHO) on March 11th, 2020, the new coronavirus disease 2019 (COVID-19) represents one of the greatest public health challenges in recent decades (1). Respiratory infections caused by viruses, mainly due to their high capacity of infection and spread, are a major cause of illness and death from a global perspective. The WHO has included coronaviruses into this group with the global alert caused by the SARS (Severe Acute Respiratory Syndrome) epidemic in 2003, caused by a newly described coronavirus (SARS-CoV). Since then, several other pathogens associated with acute respiratory system disorders have been described, such as more aggressive strains of Influenza, MERS-CoV in the Middle East, and also new types of coronavirus NL63 and KHU1 (2,3). Thus, understanding the modes of transmission of these emerging infectious diseases is a key factor both for the protection of health workers, who deal with infected people, as well as for the public, that sooner or later will be exposed to places where these agents circulate.
Airborne transmission is a key issue for the understanding of SARS-CoV-2 spread out, and it is particularly important for the community of healthcare professionals since they are more exposed to infected patients. However, it is not yet fully clarified the potential of contagion for the general public related to the high agglomeration of people, as observed in various modes of urban transport in large cities, especially in underdeveloped countries. In Italy (Bergamo city), a first study showed that RNA SARS-CoV-2 might be present in association with microparticles in the outdoor environment. Still, the detection of the virus itself was mostly inconclusive (4). Once in the free atmosphere, virus particles are submitted to conditions of natural denaturation or inactivation due to solar radiation, relative humidity, and air temperature (5,6). Such conditions result in dehydration of the viral particles derived from speech, sneezing, or coughing and, as a consequence, becoming an agglutination of viral particles with others suspended organic and inorganic molecules/microparticles in the air. Such interactions lead to changes in its size distribution and consequently effecting its diffusion/dispersion/deposition and residence time in the air. Unlike other stochastic models that determine the deposition fractions in the lung compartments from aerodynamic properties of aerosols, we present here a simulation of a model that considers the full biokinetics of the radionuclide 11C mimicking particles containing viruses. To perform that, we used a modified version of a widely tested lung deposition model developed by the International Commission on Radiological Protection (ICRP), in the context of the ICRP Publication 66 (Human Respiratory Tract Model for Radiological Protection-66) from the ICRP Task Group on Internal Dosimetry.
2 Methods
During the last decades, several studies have deeply investigated the dynamics of respiratory infections to gain information on effective treatment/prevention of these clinical events (7). Many different models have been developed till then. A respiratory infection model is a system that emulates the complexities observed on the relationship between the infectious agents and the host’s defenses. Several in vivo and in vitro respiratory models for humans and other vertebrates are available and could be easily reproduced. Mathematical models are also available and were created to numerically describe the principles and evolution of respiratory infections and their diffusion. All of them require specific inputs and have complexities related to the condition they attempt to emulate, therefore presenting both advantages and limitations.
2.1 The ICRP Model
The ICRP has developed models for aerosol pulmonary deposition, given the intensive use of natural and artificial radionuclides in the nuclear industry. Such activities range from uranium mining, where workers are exposed to dust containing naturally occurring radioactive aerosols (8), to sectors where there is the handling of wastes from nuclear power plants and related facilities (9,10). Such models depend on the aerosol size distribution of the radioactive aerosols suspended in the air, their chemical form, and its solubility and the corresponding biokinetic processes associated (11). In this context, extensive studies on radioactive aerosols biokinetics lead to successful practical application and have been particularly improved in calculating the internal dose of workers and individuals, especially those living in regions with high levels of natural radioactivity (12).
Herein we use a radionuclide-based model starting with the premise that the virus (SARS-Cov-2) has an “environmental half-life” due to its degradation process in the environment before the deposition in the respiratory tract. From this assumption, we draw a parallel with the half-life of a radionuclide, and thus being able to attribute to the virus an analog “constant of disintegration” related to the outdoor environment. The choice of the radionuclide was based on the study of Stadnytskyi et al. (13), that used a highly sensitive laser light scattering system to track the dynamics in the air of speech droplets generated by carriers of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Droplets were produced during the normal speaking condition. Speech droplets can be suspended in the air for a tenth of minutes and, therefore, can be a potential source of airborne virus transmission in low disturbed places. The settling down (vanishing) time scale of droplets before dehydration and attachment to other existing microparticles in the environment describes a behavior close to the disintegration curve of the artificial radioisotope 11C, which half-life is 20.33 minutes, Fig.1a. 11C is one of the most useful radionuclides employed for Positron Emission Tomography (PET) radiochemistry because its attachment to a biologically active molecule does not modify the biochemical properties of the inoculated compound (14). From the above, we elected the 11C to mimic the virus intake and deposition in the respiratory tract. Another reason relies on the fact that carbon molecules are one of the base constituents of cells and tissues; this allows very precise information on metabolism processes, receptor/enzyme function, and biochemical mechanisms. As a basic complementary input to the model, we used airborne SARS-CoV-2 size distribution data (15) in a stand of theoretical curves or data derived from other virus types as influenza, Fig. 1b. They obtained the size distribution of viral particles from a sequence of speeches in pre-sterilized gelatin filters installed inside a miniature cascade impactor during an air monitoring in the Renmin Hospital of Wuhan University - China, in February-March 2020. In their study, MMAD varied from 5.0 to 5.5 μm. Other studies of droplets produced during speeches with sustained vocalization found modes from 1.8 to 5.5 μm (16).
An imperative assumption to allow the use of the ICRP model, which is specially designed for radionuclides use, is to attribute a “decay property” to the virus. However, since the virus no longer is submitted to any “decay” after its intake by the human organism, the assumption of mimicking by the 11C stops at that stage. The number concentration of active virus particles inoculated in a body remains the same until the physiological responses start, and virus elimination/cell infections and replication takes place. This means that our estimates refer only to the initial deposition of the virus, and therefore the deposition fractions presented here were calculated for the first 1 minute after deposition.
The model used subdivides the respiratory system into two basic compartments: the extrathoracic region (composed of nose and trachea) and the thoracic region (composed of the bronchi to the alveolar sacs). The compartments AI, bb, BB, ET1, and ET2 represent the alveolar, bronchioles, bronchi, and extrathoracic regions, respectively, Fig. 1c. The lung is considered to have 16 generations where bifurcations occur in the bronchial and bronchioles tree. The BB region comprises the generations from 0 to 8 and the bb from 9 to 15. The rates of mechanical removal of a compartment are expressed in 1/time. Ciliary transport and retention time in the tissues migrating to the pulmonary lymph nodes are also considered in the model, as well as blood absorption from these compartments, causing then a dispute between mechanical removal and blood absorption. The compartments above act as source compartments in the infection of the lung tissues. All the calculations were performed by the AIDE (Activity and Internal Dose Estimates) software (17). The basic inputs of the model are presented in Table 1, where AMAD means “aerosol containing the viral load” median aerodynamic diameter, what implies that 50% of the viral load in the aerosol is associated with particles of aerodynamic diameter greater than the AMAD (used when deposition depends mainly on inertial impaction and sedimentation). At the same time, a compound refers to a material classified according to its rate of absorption from the respiratory tract to body fluids. In this case, deposited Type F materials are those that are readily absorbed into body fluids from the respiratory tract (Fast absorption), and f1 is the fractional absorption in the gastrointestinal tract.
2.2 SARS-Cov-2 Field Data
Because COVID-19 is highly associated through the respiratory airways, we also present data of nasopharyngeal swabs (1490) and bronco-tracheal aspirate (4) COVID-19 RT-PCR tests performed at Pedro Ernesto University Hospital (HUPE), one of the specialized units for COVID-19 treatment in the Rio de Janeiro City/Brazil. We used the data to observe the evolution of the pandemics moments just before and after the mandatory use of personal safety equipment (masks) has been implemented by the local government. The database for patients with COVID-19 started on March 29th (13th epidemiological week) up to July 12th (28th epidemiological week) when 939 patients were hospitalized and a total of 1,494 individuals were tested on this hospital. Among the realized tests, all the four bronco-tracheal aspirates resulted in negative for SARS-CoV-2, the remaining test results were presented in Table 2.
3 Results
Our results point to a far more relevant deposition fraction at the extrathoracic compartments of ET1 accounting to 47.4 and ET2 to 48.78 for the virus intake after 1 minute. In contrast, the sum for the bronchial and bronchioles compartments corresponded to 3.75 (Fig. 2). Data generated by the model is provided in the Supplementary Material.
4 Discussions
Our results point to a far more relevant deposition fraction at the extrathoracic compartments of ET1 accounting for 47.4% and ET2 to 48.78% for the virus intake after 1 minute of exposure. In contrast, the sum for the bronchial and bronchioles compartments corresponded to 3.75% (Fig. 2). For a comparison with a Stochastic Lung Deposition Model, as proposed by Madas et al. (18), they found 61.8% of the total inhaled mass fraction to be at the upper airways, ∼8.5% for the acinar airways and ∼5.5% for the bronchial compartment corresponding to a single inhalation. The differences between the two models can be attributed to the fact that Madas et al. (18) have used a mass size distribution of particles from influenza with modal values between 2-3 µm, derived from coughing of patients while we used a SARS-Cov-2 data in conditions of speeches with AMAD ∼ 5 µm. Our results stress the impact of the upper airways in the initial airborne virus retention in the respiratory tract since the total extrathoracic compartment may retain more than 96% of the virus load. Both models point to at least ∼4% of contribution to the most inner parts of the lung (bronchial and bronchioles). In the SARS-Cov-2 size distribution we have used, it is evident an existing viral load in the ultrafine particle size range (< 0.1 mm diameter). Their behavior differs from other groups of particles as the fine and coarse modes, since the virus attached to ultrafine particles may be deposit in the inner lung compartments by diffusion mechanisms. From controlled experiments, it is known that ultrafine particles peak deposition occurs in lung regions that encompass the transition zone between the conducting airways and the alveolar region (19). Therefore, since our model predicts a small fraction of particles containing virus reaching directly to the alveolar region its importance in the development of the disease should be considered since SARS-Cov-2 may bind to the cells at that compartment through ACE2 (Angiotensin-Converting Enzyme 2) receptors that are the host cell receptor responsible for mediating infection by SARS-CoV-2 (20). ACE2 plays an important role in the lungs protecting it from acute respiratory distress syndrome (ARDS) by breaking down Angiotensin II, which has inflammatory effects (21,22). Binding of the SARS-CoV-2 to ACE2 inhibits it and thus weakens its protective action on the organ (23).
The pathogenesis and virus transmission pathways are still being investigated and the already published results are, in fact, still under intense debate. That is the case for the airborne transmission of the virus. The analysis of the temporal dynamics of the SARS-Cov-2 infection suggests that after 2-3 days of the first symptoms the viral shedding starts, therefore favoring the pre-symptomatic/asymptomatic transmission of the virus (24). In this sense, understanding the first phases of the virus infection is paramount. Our results were derived from normal speech conditions as viral source-term and indoor conditions. In this specific case, we observed extrathoracic percentage infection levels in good agreement with clinical observations of patients that initially presented mild COVID-19 symptoms evolving to a more deteriorated health board (25). Although COVID-19 manifestation related to a minimum infectious dose, like that expected to the alveolar compartment, is so far not clear, and considering that: (a) none masks have 100% retention efficiency; (b) SARS-CoV-2 viral particles are found in aerodynamic diameters shorter than 0.1 μm in surveys; (c) model predicts around 5% direct virus insertion at the alveolar compartment, these facts when combined may explain why a fraction of the population even using masks get contaminated with “no apparent reason”. Respiratory tract models developed by the ICRP Task Groups have been largely used for safety and protection in nuclear activities in several countries reaching excellent performance and being validated by internal measurements. Their use for non-radioactive aerosols of biological and mineral origin and pollutants is an emerging topic and a potential to be explored.
Since ET1 and ET2 are key-compartments in the initial phase of SARS-Cov-2 deposition, we investigated the impact of the use of individual protection such as masks and face shields on the epidemiological data related to the COVID-19. Table 2 shows a survey developed by our working group from nasopharyngeal swabs and bronco-tracheal aspirate tests for COVID-19 by RT-PCR performed at Pedro Ernesto University Hospital (HUPE), in Rio de Janeiro, Brazil. The city of Rio de Janeiro in no time had adopted a full lockdown strategy, but just recommendations of social distance and personal care such as hands washing and alcohol use for personal and items/objects disinfection. However, COVID-19 data showed continuous prevalence forcing the local authorities to carry out a mandatory use of masks on April 23rd. Our results clearly suggest that the fraction of both hospitalized and total tested patients with SARS-CoV-2 detected from the RT-PCR tests exhibited a significant decrease when we observed an immediate drop in percentage from 85% to 60% (Figure 3). We do believe that the analyzed sample is representative of the city population since the Pedro Ernesto University Hospital (HUPE) receives patients from several health units spread all over the urban domain. We should also note that due to socio-economic reasons and also difficulties in the acquisition of high-quality safety equipment, most of the population made use of simple, cheap, and home-made masks and face shields. However, despite these circumstances, effective response against the pandemic was achieved.
In summary, respiratory tract models developed by the IAEA Group tasks have been largely used for safety and protection in nuclear activities in several countries, reaching excellent performance and being validated by in locus measurements. Their use for non-radioactive aerosols and microparticles of biological and mineral origin and pollutants is an emerging topic and a potential to be explored.
Data Availability
All data used to generate the results presented herein are provided in the supplementary material.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial frameworks that could be construed as a potential conflict of interest.
Author Contributions
All authors contributed equally to the development and text of the manuscript. Heitor Evangelista (Project coordinator, creator of the project), Luís Cristóvão Porto (medical survey), Sergio Gonçalves Junior (participation in the text content and discussions), Eduardo Delfino Sodré (participation in the text content and discussions), Juliana Nogueira (numeric model analysis), Cesar Amaral (participation in the text content and discussions), Marcio Cataldo Gomes da Silva (participation in the text content and discussions).
Funding
The authors received no funding for this research.
Acknowledgments
We greatly thank Dr. Luiz Bertelli for the availability of AIDE (Activity and Internal Dose Estimates) software and important comments on the numeric model and Roberta Priori for draw designs.