Inward Bound - The T'ai Chi Corner

Inward Bound - The T'ai Chi Corner
Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Mon Feb 22, 2021 1:40 p

Spring is very close, and this following is a very good way, especially for those confined indoors who wish a relaxing full body exercise in a small space, for your mind, body and soul to greet spring, as wood is the element of spring.

It is performed and narrated by a very accomplished and knowlegeable young Chinese woman named Mimi Kuo Deemer who has a very pleasant voice to listen to as she guides you through the routine, which is 17.09 minutes long if you do the full routine, or you can extract portions like the deer frolic if you have less time, all of which can be done in a very small space as you are not stepping any distance:

Qigong for the Wood Element: Nourishing Practices for Body, Mind and Spirit

Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sat Sep 11, 2021 1:40 p


ACS Nano. 2020 Nov 25 : acsnano.0c08484.

Published online 2020 Nov 25. doi: 10.1021/acsnano.0c08484

PMCID: PMC7724984

PMID: 33236896

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author


Coronavirus disease 2019 (COVID-19), due to infection by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is now causing a global pandemic.

Aerosol transmission of COVID-19, although plausible, has not been confirmed by the World Health Organization (WHO) as a general transmission route.

Considering the rapid spread of SARS-CoV-2, especially nosocomial outbreaks and other superspreading events, there is an urgent need to study the possibility of airborne transmission and its impact on the lung, the primary body organ attacked by the virus.

Here, we review the complete pathway of airborne transmission of SARS-CoV-2 from aerosol dispersion in air to subsequent biological uptake after inhalation.

In particular, we first review the aerodynamic and colloidal mechanisms by which aerosols disperse and transmit in air and deposit onto surfaces.

We then review the fundamental mechanisms that govern regional deposition of micro- and nanoparticles in the lung.

Focus is given to biophysical interactions between particles and the pulmonary surfactant film, the initial alveolar-capillary barrier and first-line host defense system against inhaled particles and pathogens.

Finally, we summarize the current understanding about the structural dynamics of the SARS-CoV-2 spike protein and its interactions with receptors at the atomistic and molecular scales, primarily as revealed by molecular dynamics simulations.

This review provides urgent and multidisciplinary knowledge toward understanding the airborne transmission of SARS-CoV-2 and its health impact on the respiratory system.


Coronavirus disease 2019 (COVID-19), due to infection of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently causing a global pandemic, with more than 53 million confirmed cases and 1 million deaths, as of November 15, 2020, in more than 200 countries, areas, and territories in the world. 1

Given the current world population of 7.8 billion, approximately 1 out of every 150 humans on earth has or had been infected with COVID-19.

The incidence of COVID-19 shows strong ethnic disparities.

Emerging epidemiological data suggest that black communities are affected disproportionately hard by COVID-19.

Although they account for only 13% of the United States population, African Americans constitute 24% of COVID-19 deaths, nearly twice of what would be expected based on their share of the national population. 4

In some regions of the United States, such as Chicago, the death rate of African Americans has been reported to be six times higher than that of the white population. 5

Similarly, it was found that indigenous peoples of the Pacific, including Native Hawaiians and Pacific Islanders, have a significantly higher COVID-19 infection rate than other ethnic groups. 6

These health disparities among racial and ethnic groups are often related to different living conditions, work circumstances, and underlying health conditions including diabetes, heart disease, and lung disease, all of which are closely linked to more severe cases of COVID-19. 7,8

All available evidence suggests that SARS-CoV-2 is highly contagious and quickly spreads in our communities. 9

To date, the confirmed modes of SARS-CoV-2 transmission include respiratory droplets, direct (person-to-person) and indirect (fomite) contacts, as well as scarce reports of fecal–oral transmission. 10

However, the latest research suggests that fomite transmission is unlikely to be a major route of transmission as attempts to culture SARS-CoV-2 from surfaces were largely unsuccessful. 11,12

On the other hand, airborne transmission of SARS-CoV-2, although plausible, has not been confirmed by the World Health Organization (WHO) as a general transmission route. 10

However, on October 5, 2020, the United States Centers for Disease Control and Prevention (CDC) updated their guideline and acknowledged that airborne transmission of SARS-CoV-2 could occur under special circumstances that include enclosed spaces, inadequate ventilation, and prolonged exposure to events involving heavy breathing, such as singing and exercising. 13

Given the rapid spread of the coronavirus, especially nosocomial outbreaks and other superspreading events, there is an urgent need to carefully access the possibility of airborne transmission, especially the colloidal and aerodynamic mechanisms of aerosol dispersion and deposition, and the molecular interactions between SARS-CoV-2 and receptors. 14−22

An open letter by 239 researchers from 32 countries, published on July 6, 2020, challenged the WHO’s view on aerosol transmission. 23

Virus transmission via aerosols can be fundamentally different from transmission via respiratory droplets in terms of their colloidal and aerodynamic mechanisms. 24

Respiratory droplets are emitted when an infected individual coughs or sneezes.

Due to their large size (usually >5–10 μm), these droplets in general have a very short lifetime in air and a limited transmission distance of <2 m, although some studies suggest that these droplets could travel as far as 8 m. 25

Aerosols, often called droplet nuclei, are airborne particles or droplets much smaller than 5 μm. 24

It should be noted that the aerodynamic size of 5 μm is not a strictly defined threshold, but a historical value that provides a loose definition of aerosols.

In contrast to respiratory droplets, it appears that coughs or sneezes are unnecessary for producing infectious aerosols, and simple exhalation is sufficient. 26,27

In general, most aerosols in exhaled breath are found to be are smaller than 4 μm, with a median between 0.7 and 1 μm. 24,28

Also different from respiratory droplets, aerosols are too small to rapidly settle out of the atmosphere under the influence of gravity. 22

Respiratory droplets in the size range between 60 and 100 μm, although could be carried more than 6 m by sneezing, generally fall to the ground within two meters. 24,29

Depending on the environmental conditions, virus particles (i.e., virions) may bind to aerosols and thus travel over a significantly longer distance (>2 m) and remain floating in air for a much longer time period (up to hours) than respiratory droplets.

It is now known that air pollution is an associated risk factor of COVID-19.

An increase of 1 μg/m3 in PM2.5 is associated with an 8% increase in the COVID-19 death rate. 32

In addition, it was found that outside the optimal relative humidity (RH) range of 40–60%, the viability of influenza virus in droplets increased both at higher (>60%) and lower (<40%) RH. 33

A similar humidity dependence has been found for the transmission of SARS-CoV-2. 34

E-cigarette aerosols have also been identified as suspicious anthropogenic aerosols potentially capable of spreading SARS-CoV-2. 35,36

Aerosol transmission was reported to be responsible for a superspreading event of SARS-CoV in a housing block in Hong Kong in 2003. 37,38

The 50-story building had 342 confirmed cases of SARS and 42 deaths.

The report identified defects in the wastewater plumbing system as a transmission mode within the building, which facilitated the transport of virus-laden aerosols through empty U-bends in bathrooms.

In a superspreading event at a church in Skagit County, WA, 45 out of 60 singers of a choir were diagnosed with COVID-19 with two deaths, after only a 2.5 h rehearsal in which nobody present was coughing or sneezing or appeared ill. 39

Another superspreading event reported in China involved 10 individuals from three families in the same restaurant. 40

This study found that the virus was able to travel far enough to infect individual members of the three families.

It appears that virus transmission in this outbreak cannot be explained by droplet transmission alone and consequently might involve aerosol transmission. 40

In fact, although available evidence remains scarce, SARS-CoV-2 has been detected in hospital air. 41,42

A recent study suggested that in every hour COVID-19 patients can exhale millions of SARS-CoV-2 RNA copies into air. 27

Given a study demonstrating that SARS-CoV-2 remained viable in aerosols for at least 3 h with only limited reduction in infectious titer, 43 there is an urgent need to understand the likelihood of SARS-CoV-2 transmission via aerosols.

In this paper, we review the complete pathway of airborne transmission of SARS-CoV-2 from aerosol dispersion in air to subsequent biological uptake after inhalation.

The review is organized in three main sections.

First, we review the aerodynamic and colloidal mechanisms by which aerosols disperse and transmit in air and deposit onto surfaces.

Second, we focus on lung deposition of the aerosols and their interaction with the first-line host defense pulmonary surfactant (PS) film at the alveolar–capillary barrier.

Third, we summarize the current understanding of molecular interactions between SARS-CoV-2 and receptors.

This review provides urgent and multidisciplinary knowledge in understanding the airborne transmission of SARS-CoV-2.


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sat Sep 11, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Aerosol Dispersion and Deposition

Basic Considerations

The dispersion of infectious aerosols in air and their subsequent deposition onto surfaces are physical processes that can be modeled by treating air as a continuous medium that exchanges mass, energy, and momentum with the aerosol particles and associated virions.

As physical processes, they can potentially be controlled through an engineering approach, for example, by designing interior spaces with controlled ventilation and humidity or by engineering surface coatings that disrupt virus transmission.

In this section, we review the literature on the transport of infectious aerosols in the atmosphere and their deposition onto surfaces.

In particular, our review is intended to facilitate development of realistic models that capture all relevant physical effects and are parametrized by experimental measurements.

First, we develop some basic considerations concerning SARS-CoV-2 virions and aerosols.

Aerosols are commonly defined as solid particles or liquid droplets with diameters between a few nanometers and micrometers.

This range includes the smallest e-cigarette aerosols, with characteristic size D = 150 nm, to the largest PM2.5 particles, with D = 5 μm.

Given the typical diameter of a SARS-CoV-2 virion at approximately Dv = 70 nm, 44 the size ratio Dv/D can vary from Dv/D ≈ 0.5 to Dv/D ≈ 0.01.

The typical number of virions per aerosol is also an important quantity.

We are not aware of any study that has sought to directly measure this quantity for aerosols laden with SARS-CoV-2, which must depend on the biophysical process generating infectious aerosols.

Still, some idea of this number can be obtained from other quantitative data on SARS-CoV-2.

According to Bar-On et al., 45 COVID-19 patient sputum contains 106–1011 viral RNAs per milliliter, although this number “can overestimate infectious virions”.

If sputum matter is aerosolized as 5 μm particles, then the average number of virions per aerosol, using 1011 virions/mL, is no greater than 50.

The actual number is probably much less than this upper bound.

Furthermore, previous studies have sought to characterize the number of virions per aerosol for other viruses.

For instance, in a study of artificially generated, infectious aerosols carrying the MS2 bacteriophage, the authors estimated ∼10–4 virions per aerosol particle. 46

Since this is a sample-averaged quantity, it gives little insight into the probability distribution of the number of adsorbed virions for an individual aerosol.

Under the reasonable a priori assumption that adsorption of an individual virion by an aerosol is an independent random event, the conditional probability of an aerosol with at least one adsorbed aerosol carrying two or more virions is negligible.

However, this assumption may be incorrect if there are physical mechanisms that endow a small number of aerosols with a large number of virions, while leaving the rest of the aerosols without virions.

In an older study, artificially generated aerosols with adsorbed adenovirus virions were individually resolved with electron microscopy. 47

It was found that the vast majority (67%) of aerosol particles had only one virion, 11% had two virions, etc.

Overall, it appears to be most likely that each infectious aerosol carries one to a few SARS-CoV-2 virions.

A recent review of the soft matter science of COVID-19 broadly concurs with the analysis here. 48

Generation and Morphology of Infectious Aerosols

Conceptually, we can distinguish several possible modes of association between virions and aerosols.

For liquid aerosols, virions could be encapsulated within the aerosol or adsorbed at the aerosol surface.

For solid aerosols, virions could be adsorbed to the aerosol surface or, if the aerosol is porous, contained within the aerosol structure.

Moreover, the prevailing mode of association is clearly influenced by the biophysical process that generates the infectious aerosol.

Exhalation from infected individuals is obviously the most important process for generation of aerosols laden with SARS-CoV-2 virions.

Coughing, sneezing, breathing, and human speech can produce liquid droplets with a broad size range, from <0.8 to 10 μm.

The droplets thus generated can contain virions, along with salts, biopolymers, and other biomolecules.

Upon exposure to air, the fate of a droplet depends on its size and ambient conditions, especially RH.

It is commonly assumed that large droplets (>5 μm) will quickly sediment to the ground, while smaller droplets (i.e., aerosols) will remain aloft for a long period.

​Figure11a shows a high-speed image of a sneeze; large droplets with ballistic trajectories are clearly visible (green), while a moist and buoyant “puff” (red) carries small droplets (discussed in detail below). 51

Aerosols in the size range between 0.5 and 5 μm are known to be the most breathable fraction of all aerosols (see Dispersion of Bioaerosols in Air for details).

It should be noted that droplet size is a dynamic quantity, as expelled droplets will start to evaporate if the ambient RH is <100%. 52

This will reduce the size of droplets, prolonging their time in air.

The droplet settling time therefore largely depends on initial size and RH. 53

Since the pioneering work of Maxwell in 1877, many studies have sought to model, using continuum microscopic theory and numerical simulations, the coupled heat and mass fluxes in the vicinity of an evaporating droplet as well as the changing size of the droplet. 54−60

The coupling of the two fluxes through the latent heat and the moving boundary condition makes this a challenging problem in transport theory. 61

To simplify the modeling, various assumptions are commonly made, for example, that the droplet radius can be regarded as quasi-static or stationary for determination of the instantaneous heat and mass fluxes and that the surface of the droplet is in thermodynamic equilibrium.

Recently, the problem of droplet evaporation has received fresh impetus from the COVID-19 pandemic; for instance, it was recently estimated that at 50% RH, droplets with the initial size smaller than 40 μm, released from a typical human height, will dry out before hitting the ground. 62

Droplets that dry out in air produce solid residues, or “droplet nuclei,” as shown schematically in ​Figure11b. 63

Owing to their small size, these dry particles can potentially remain aloft for a long time.

Using an aerodynamic particle sizer for droplets produced by coughing, Yang et al. determined that the dry nuclei size distribution peaked at 1 μm, consistent with a much earlier work from Duguid. 50,64

Drying of the droplet can potentially lead to gelation of mucin and other proteins contained in the droplet, and virions may be contained within the gel network. 65

We are not aware of any attempt to directly image droplet nuclei obtained from air in order to determine their morphology or microstructure.

For instance, it would be useful to know whether virions are typically contained within droplet nuclei or are carried as surface adsorbates.

However, there have been some experimental studies of the drying of model virus-laden solutions that were aerosolized and deposited on solid surfaces.

In one such study, Vejerano and Marr determined via fluorescence microscopy that phi6 virions, used as a surrogate for influenza, were homogeneously dispersed within desiccated droplets (​Figure11c). 65

In another such study, three phi6-laden solutions — saliva, SM buffer, and water — were deposited on a surface and left to dry.

Interestingly, the virions in saliva exhibited the greatest viability after 14 h, even though the saliva droplets completely evaporated in that time for all tested RH conditions, leaving a thin dry deposit.

It was hypothesized that biopolymers in saliva may enhance virion survival.

Necessarily, these studies give little information on the interaction of virions and aerosols of comparable, that is, submicron, size.

However, based on experimental and simulation studies involving nanoparticles and virions, it is plausible that SARS-CoV-2 virions and nanoaerosols bind to form a loose cluster. 47,67,68

For instance, the Cyro-TEM image in ​Figure11d provides evidence of binding between vesicular stomatitis virions and small nanoparticles.

Recently, this conceptual model of SARS-CoV-2 virions adsorbed to “carrier” nanoaerosols has motivated the development of nanofiber filters for capture and filtration of ambient nanoaerosols. 69

Infectious aerosols can also be generated by medical procedures, such as intubation, 70 and even by flushing a toilet containing infectious material. 71

The relevance of these other mechanisms to COVID-19 transmission is poorly understood, but could easily be related to the superspreading events listed in the Introduction.

Dispersion of Bioaerosols in Air

The dispersion of aerosols has a strong dependence on ambient airflow, even when the aerosols are released only meters away from the ground, as in human exhalation or in medical procedures.

As a reasonable first approximation, aerosols can be regarded as “point-like” or passive tracer particles that simply follow fluid pathlines.

(Recall that a streamline is a curve tangent to the flow field at a particular moment in time, while a pathline is the trajectory followed by an individual parcel of fluid.)

Thus, significant effort has been invested in understanding air flow in complex geometries characteristic of interior and exterior built environments.

Meter-scale air flows in these environments are typically turbulent.

Moreover, temperature variations (e.g., induced by operation of medical devices or the heat output of human beings) can drive the natural convection of air.

Due to the complexity of the governing physics, computational fluid dynamics (CFD) packages like Fluent or Comsol are typically used to solve the energy equation (for the air temperature) and the Navier–Stokes equations (for the air velocity) in a meshed spatial domain. 72

Direct numerical simulation (DNS) of the Navier–Stokes equations is usually computational prohibitive, requiring an approximate approach.

One such approach is to consider the Reynolds-averaged Navier–Stokes (RANS) equations.

These equations arise from decomposition of the Navier–Stokes equations into time-averaged and fluctuating components and subsequent time averaging of the equations.

As a result, one obtains equations without an explicit time derivative, but with a new contribution to the fluid stress, the Reynolds stress.

This stress is not known a priori, as it represents fluctuating components; one must choose a model (or “closure”) for the Reynolds stress, such as the k-epsilon model.

Another approximate approach, distinct from RANS, is large eddy simulation (LES).

A recent paper by Blocken discusses the relative advantages of RANS and LES in modeling interior air flows. 73

As an example of CFD used to model bioaerosol dispersion, Buchanan and Dunn-Rankin considered an operating room with heat generated by personnel and medical devices, which were meshed as solid objects. 74

They used Fluent to solve the energy balance equation and the RANS equations with k-epsilon closure.

Aerosol particles were assumed to have velocities equal to the local air flow velocity, plus a stochastic component to model the effect of air flow velocity fluctuations.

In another study, Li et al. considered a SARS superspreading event in a Hong Kong hospital ward, in which one patient infected 138 other patients, visitors, and medical personnel. 75

Airborne transmission was suspected in this event.

Accordingly, the study combined inspections and measurements of the on-site ventilation system, some limited on-site experimental measurements of aerosol dispersion, and CFD modeling of the ventilation and air flow.

In their CFD modeling, they included convection due to heat from patients, lighting, and heating systems and meshed the detailed geometry of the ward.

Patients were included as rectangular prisms, and bioaerosols were assumed to disperse as passive tracers.

The simulated aerosol distribution displayed reasonable agreement with measured concentrations.

In these single-particle models, one can obtain some additional realism without much additional model complexity by incorporating particle inertia.

Specifically, one can add a differential equation for the particle velocity, allowing a particle to exhibit relative motion with respect to a background flow field.

Particle density is also easily incorporated.

Although the governing equations for particle motion are simple, tracers with inertia can have rich physics, even from a fundamental point of view.

For instance, Vilela and Motter considered heavy inertial tracers in model open background flows. 76

Despite the expectation that aerosols can escape open flows in finite time, they found that the presence of vortices in the background flow can lead to permanent trapping of the aerosols.

Concerning CFD of droplet dispersion in ventilated rooms, Sun and Ji modeled droplets produced by coughing as heavy tracers with inertia and simultaneously solved differential equations for the mass and energy fluxes from the evaporating droplets. 77

The background flow in a meshed model ventilated room was solved with RANS.

There may be other physical effects, microscopic in origin, that cause particles to move relative to background flows.

For instance, phoretic forces in concentration and/or temperature gradients could contribute to particle motion.

Although some recent papers have considered phoretic effects in turbulent flows, 78,79 this topic is largely unexplored, and further work is needed to assess the relevance of phoresis in the specific context of air flows.

The effect of nonspherical aerosol shape in turbulent flows is also a vigorous area of research. 80

Concerning Brownian motion of aerosols, one study has concluded that it is generally dominated by turbulent diffusion, that is, random motion from turbulent velocity fluctuations. 81

So far, our discussion has considered the motion of single droplets in ambient air flows determined by room ventilation and temperature conditions.

However, recent research has shown that, during sneezing and coughing, droplets are emitted in turbulent gas clouds, that is, discrete puffs or continuously emitted jets. 82,83

(More precisely, jets are continuously provided momentum by the emitter, while puffs have a conserved quantity of momentum.)

The momentum, temperature, humidity, and chemical composition of the cloud can have a significant effect on droplet trajectories and evaporation dynamics. 25

The RH in a cloud is typically much higher (∼85%) than the RH in interior spaces (∼40%, although this can vary according to season and climate), prolonging droplet lifetime.

Moreover, the infectivity of respiratory viruses is apparently lower under conditions of intermediate ambient RH (i.e., 40–60%) than high ambient RH (>60%); the humid microenvironment provided by the cloud could potentially promote infectivity before it disperses. 33,84

Concerning temperature, the jet or puff emitted from the body will generally be warmer than ambient conditions, providing some initial positive buoyancy to the cloud.

Bourouiba et al. used high-speed photography to directly image puffs, including suspended droplets, emitted by adult human subjects during coughing and sneezing events (Figure​Figure22a). 82

Large droplets were observed to move ballistically, with little effect from the surrounding cloud, while smaller droplets (visible as mist) tended to move with the cloud over long distances.

Puff clouds initially evolved as cones with some upward curvature due to buoyancy, consistent with an established model of self-similar growth of turbulent jets by entrainment of ambient fluid.

An analogue experiment involving the injection of fresh water into salty water allowed the authors to probe the crossover between jet-like and puff-like cloud dynamics.

Moreover, a theoretical model was developed to describe fallout of droplets as the puff entrains fluid and decelerates (​Figure22b).

While Bourouiba et al. used analytical modeling, CFD modeling can be extended to include turbulent jets and puffs ejected during coughing and sneezing.

Redrow et al. used CFD to model the dispersion of particles in an initially buoyant jet emitted from the mouth of a meshed “mannequin”, using initial conditions determined from cough flow rate data from the National Institute of Occupational Safety and Health (NIOSH). 85

They found that the maximal particle concentration appeared 0.65 m from the mouth of the mannequin (​Figure22c).

In the same work, Redrow et al. sought to extend droplet evaporation models to account for the local temperature and chemical composition of ejected gas clouds.

A similar study was recently carried out by Pendar and Pascoa, motivated by the need to determine social distancing guidelines in the context of COVID-19 (​Figure22d,e). 86

Indeed, CDC guidelines, including the familiar six-foot rule, are still based on the 1930s-era concept of an isolated respiratory droplet settling under gravity, despite the wealth of contemporary literature reviewed above.

In a recent work, several investigators in the airborne transmission community have sought to develop an integrated framework that consolidates recent insights concerning ejection, dispersion, and inhalation of respiratory droplets, and which can provide concrete guidance for social distancing (in both space and time) via simple model equations. 51

Modeling of Aerosol Deposition onto Surfaces

Infectious aerosols in the atmosphere can deposit onto surfaces as fomite contaminants.

Physical contact of a potential host with fomites can lead to infection hours or even days after deposition. 43

Moreover, infectious aerosols that are inhaled must travel through lung airways and ultimately deposit on an interior surface in order to transmit virions (for details see Future Directions in Modeling Aerosol Dispersion and Deposition).

In both scenarios, an aerosol and associated virions must move in near-surface fluid flow, subject to sedimentation, diffusion, phoretic forces, and other effects, and adsorb to the surface, as shown schematically in ​Figure33a.

An understanding of this physical process could facilitate development of mitigation strategies, for example, surface coatings that hinder deposition or reduce the viability of adsorbed virions.

Empirically, it is well-established that fluid flows tend to “stick” to solid surfaces, that is, the fluid velocity obeys a no-slip boundary condition on surfaces.

Accordingly, air flow will slow down in the vicinity of surfaces, especially within surface boundary layers.

Within boundary layers, which include the laminar sublayer of turbulent flow, the flow can be regarded as viscous.

Accordingly, the physics of deposition can be significantly different from the physics of dispersion in the atmosphere.

For instance, Brownian diffusion may play a significant role near surfaces even when it is irrelevant in the bulk fluid.

A point-like description of the particle may no longer be valid, as the near-surface flow may vary on length scales comparable to the particle size.

Moreover, the inertia of the fluid and/or of the particle can potentially be disregarded in these regions, owing to the reduction of the flow speed and, consequently, of the fluid and particle Reynolds numbers.

In the remainder of this section, we review some of the literature on modeling deposition of aerosols.

Much of the groundwork on deposition in the viscous sublayer of turbulent flow, similar to what would be expected in fomite deposition, was laid by McLaughlin, Ahmadi, and co-workers. 87,88

An influential early model considered the deposition of a spherical particle via Brownian dynamics. 87

The mean background flow was modeled as a linear shear profile.

The dynamics of the particle was modeled with a Langevin equation involving an inertial term, viscous drag terms, and a Brownian fluctuating term.

The effect of turbulent fluctuations was crudely modeled with a wall-normal contribution to the flow velocity that varied quadratically with distance from the wall.

The authors found that Brownian motion dominated the deposition process close to the wall.

We highlight this particular model as capturing the most basic physical effects affecting deposition and as a point of departure for development of more realistic models.

Limitations of the model include the neglect of aerodynamic (i.e., fluid-mediated) interactions of the particle with the wall; these interactions lead to variation of the drag coefficient and noise term with distance from the wall. 89

The particle is considered to be point-like, that is, gradients in the background flow velocity in the vicinity of the particle do not affect the particle velocity.

Phoretic, electrostatic, and other forces are not included.

Inertia of the fluid, which can give rise to wall-normal forces, is neglected.

The particle shape is assumed to be spherical, and the background flow is severely approximated.

Finally, if the virion size is comparable to the aerosol size (recalling the earlier discussion of relative size), the multibody dynamics of aerosol and associated virion or virions must be considered.

Many of these limitations were addressed in subsequent papers that also worked within the Brownian dynamics framework.

For instance, Fan and Ahmadi considered deposition of ellipsoidal particles. 90

The model includes the anisotropic drag on the particle, that is, the fact that drag is smaller for motion along the major axis than for motion along a minor axis.

The particle orientation evolves according to coupled Langevin equations, with the orientation represented with quaternions.

A term representing the shear-induced lift force (which arises from fluid inertia) is included in the Langevin equation.

This addition is somewhat ad hoc, in the sense that the lift force is represented by a phenomenological term in force balance equations and does not emerge from a solution of the dynamics of the suspending fluid.

With the growth of computing power, a recent trend is to simulate the background flow in the viscous sublayer with CFD by DNS.

For instance, Zhang and Ahmadi used a pseudospectral DNS method to generate time-dependent instantaneous velocity fields in vertical and horizontal ducts. 91

The motion of a spherical aerosol was again modeled with a Langevin equation, in this case with a Reynolds number-dependent drag term.

Another paper from the same group combined CFD simulation of velocity and temperature fields with a model for thermophoretic force on a sphere. 92

Marchioli, Soldati, and co-workers have also been very active in applying DNS, especially pseudospectral methods, to aerosol deposition. 93

In particular, their work highlights the key role of coherent flow structures like quasi-streamwise vortices in segregating and trapping aerosols, depositing them onto walls, and re-entraining them from walls (​Figure33b).

In the specific context of bioaerosol inhalation and infection, various studies have addressed aerosol motion and deposition in the confined spaces of lung alveoli. 94−100

Small aerosols (<1–2 μm) preferentially deposit in alveoli (see Future Directions in Modeling Aerosol Dispersion and Deposition for details). 101

Flows in these spaces generally possess low Reynolds numbers, with the Reynolds number varying characteristically between 0.01 and 10. 102

Alveoli expand and contract during breathing, making the flow and confining geometry time dependent.

Balashazy et al. modeled this flow analytically using a superposition of model uniform and radial flow components in an idealized hemispherical alveolus geometry.

Aerosol trajectories were calculated via Brownian dynamics and included the effects of particle inertia, Brownian motion, and sedimentation.

They found that the orientation of the alveolus with respect to gravity can affect the deposition efficiency for larger aerosols. 97

Similarly, Haber et al. applied a mixed numerical and analytical approach to solve for the flow in a model alveolus in the shape of a truncated spherical cavity. 94

Additionally, they solved for the deterministic trajectories of massless particles (i.e., passive tracers) and found that they followed complex, chaotic, and effectively stochastic pathlines.

Similar results for flow streamlines and particle trajectories were obtained by Sznitman, as shown in ​Figure33d. 102

A comprehensive review of aerosol delivery in the respiratory tract is provided by Longest and Holbrook. 103

Future Directions in Modeling Aerosol Dispersion and Deposition

We suggest two promising directions for simulation of aerosol deposition.

The Brownian dynamics method can be combined with the boundary element method (BEM) for the solution of Stokes flow (the viscous limit of the Navier–Stokes equations) in complex confined geometries.

Specifically, the BEM can capture the aerodynamic effect of complex particle shape and aerodynamic interactions with confining surfaces. 104

A recent work by Gubbiotti et al. demonstrates how to consistently model fluctuating Brownian motion in confinement, suggesting combination of this Brownian dynamics approach with the BEM. 105

Since the particle surface is meshed in the BEM, anisotropic and “patchy” or surface-heterogeneous forces can be included in a straightforward manner.

The multibody dynamics of several particles, representing an aerosol and associated virions, can be modeled in this approach.

However, the BEM is most suitable to linear partial differential equations; for instance, it is not straightforward to include the effect of fluid inertia.

Another promising approach is the lattice Boltzmann method (LBM) with fluid–structure interaction between the fluid (modeled as an LBM fluid) and moving particles. 106,107

Here, the effect of fluid inertia naturally arises from the underlying fluid lattice dynamics.

The LBM has been successfully applied to airflow analysis in the alveolar region of the lung, 108 electrophoretic deposition of aerosols onto a sensor surface, 109 and filtration of aerosols by fiber arrays (such as in N95 masks), 110 as shown in ​Figure33c.


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sat Sep 11, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Aerosol Deposition in the Lung

Regional Deposition in the Lung

The respiratory system is the first and primary body region affected by SARS-CoV-2.

Inhaled virus-laden aerosols deposit in different regions of the respiratory system, depending on the physicochemical properties of the aerosols, including their aerodynamic size, density, shape, charge, hydrophobicity, and hygroscopicity.

Among these influencing factors, the aerodynamic size plays a predominant role in determining the regional deposition of inhaled aerosols (​Figure44b). 112,113

As shown in ​Figure44a, the lung, a key infection site, features a branched architecture that consists of 23 generations of dichotomous branching, from the relatively thick trachea (generation 0), with an average diameter of 15–25 mm, to small airways with diameters of <2 mm, ending up with 300–600 million of terminal alveoli (generation 23), with an average diameter of only 200 μm per alveolus. 114,115

The primary mechanisms of aerosol deposition in the lung include inertial impaction, gravitational sedimentation, and Brownian diffusion as well as, to a lesser extent, interception (applied to elongated aerosols such as fibers) and electrostatic precipitation (applied to charged aerosols).112,114,116,117

Inertial impaction mainly affects large aerosols >5–10 μm.

These large aerosols primarily deposit in the nasopharynx and the tracheobronchial region of the lung (first 10 generations) where the air flow rate is high and the flow condition is predominantly turbulent.

The degree of particle deposition by inertial impaction is quantified by the dimensionless Stokes number (Stk), defined by an equation where, ρp and dp are the density and aerodynamic diameter of the airborne particles, respectively, u and μ are the mean velocity and dynamic viscosity of the air flow, respectively, and D is the characteristic diameter of the airway.

Gravitational sedimentation mainly affects medium-sized aerosols between 1 and 5 μm.

These medium-sized aerosols mostly deposit in small airways (11th–16th generations) where the flow rate is low and the particles have a relatively long retention time in this region.

The terminal settling velocity (Vts) of the particles is expressed by an equation where, g is the gravitational acceleration.

Brownian diffusion is the mechanism that controls the motion of small aerosols mainly <1 μm in diameter.

These small aerosols can penetrate into the alveolar region where the air flow rate can be approximated as negligible.

Deposition by Brownian diffusion is proportional to the Brownian diffusion coefficient (DB) defined by an equation where k is Boltzmann’s constant, T is the absolute temperature, and c is the Cunningham correction factor.

Depending on their specific deposition regions, inhaled aerosols experience two different clearance mechanisms in the lung.

Aerosols larger than 5 μm, including most respiratory droplets, are deposited in the upper regions of the respiratory tract, from which they are removed by the mucociliary escalator and eventually spit out or swallowed (​Figure44c).

Mucociliary clearance is regulated by a dual-layered airway surface liquid composed of a ∼10 μm-thick aqueous periciliary liquid (PCL) plus a 10–70 μm-thick viscoelastic mucus gel layer atop the PCL. 119,120

This clearance mechanism is chemically nonspecific and highly efficient, with 80–90% of inhaled aerosols being excreted from the lung within 24 h. 121,122

In contrast, aerosols smaller than 5 μm, especially those smaller than 3 μm, 112 are capable of penetrating to the distal alveolar regions of the lung, at which the primary clearance mechanism is macrophage clearance with the assistance of the immunological surfactant proteins (​Figure44d). 122−124

Macrophage clearance is found to be most effective for aerosols in the size range 0.5 and 5 μm,125 which happens to be the most breathable fraction of all aerosols. 112,122


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sun Sep 12, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Aerosol Deposition in the Lung, continued ...

Interaction with the Pulmonary Surfactant Film

In contrast to direct and fomite transmissions, by which SARS-CoV-2 first infects the nasal cavity, replicates itself in the upper airways, and eventually propagates into the lung periphery, 126 SARS-CoV-2 transmission mediated by aerosols may directly infect the lung interior, 20,127 thus leading to a rapid onset of the most severe third-phase alveolar infection of COVID-19. 128,129

As illustrated in ​Figure 55, once reaching the alveolar regions, SARS-CoV-2 attacks both type I and type II alveolar epithelial cells, promotes lung flooding, and causes the life-threatening acute respiratory distress syndrome (ARDS). 130−132

SARS-CoV-2 attacks cells via the angiotensin-converting enzyme 2 (ACE2) receptor (see Pulmonary Surfactant Biomolecular Corona for detailed mechanisms of virus interaction with ACE2 receptors), 133 which is highly expressed on the surface of type II cells. 134

Among the most important physiological functions of these cells is the production of PS, a lipid–protein complex that lines the entire air–water surface of alveoli as a thin film. 135

It is composed of mostly lipids (∼90% by weight) and four surfactant associated proteins, that is, SP-A, B, C, and D (∼10% by weight). 135

A closer evaluation of the lipid portion of PS reveals that the majority of the lipids are phospholipids that contribute to its high surface activity (i.e., the ability to reduce surface tension). 136,137

SP-B and SP-C are hydrophobic proteins embedded in the phospholipid matrix that work with phospholipids to achieve low surface tensions during normal tidal breathing. 138

SP-A and SP-D are hydrophilic proteins that assist in innate immune responses, including macrophage clearance, 139 although SP-D is only partially associated with the PS membranes. 140

The PS film serves two major physiological functions.

The first function of the PS is innate immunity, that is, serving as the initial alveolocapillary barrier and first-line host defense against aerosols and pathogens that reach the deep lung.

The immunological function of the PS is mostly regulated by the two hydrophilic surfactant proteins, SP-A and SP-D.

It has been reported that macrophage phagocytosis of nanoparticles was triggered by SP-A. 124

In addition, SP-D bound to the spike protein of SARS-CoV triggered macrophage phagocytosis of the virus. 143

The second function of the PS is its biophysical function that allows for normal lung mechanics and prevents lung collapse by reducing the alveolar surface tension.

With surface tension reduction to near-zero levels, inhalation–exhalation cycles are maintained at minimum mechanical energy expense. 135

Depending on their physicochemical properties, inhaled particles could interfere with the physiological function of PS via two major perspectives.

First, particles may directly inhibit the biophysical function of the PS film at the air–water surface of alveoli (​Figure66). 144−146

By binding and denaturing hydrophobic surfactant proteins (SP-B and/or SP-C), particles cause surfactant inhibition, quantitatively measured by increasing surface tension and PS film compressibility during respiratory cycles, 146,147 thereby leading to alveolar collapse and respiratory failure, which is a major symptom of ARDS. 135

A wealth of in vitro and in vivo studies have demonstrated that the biophysical function of PS can be inhibited by various particles, including nanoparticles, 144−148 nanoenabled sprays and drugs,149,150 PM2.5,151−153 and pathogens such as the H1N1 influenza virus. 154

Second, particles may interfere with PS metabolism and homeostasis in the alveolar hypophase.

It was found that nanoparticles caused surfactant degradation from the more surface active large aggregates to less surface active small aggregates. 144

The poor surface activity of small aggregates is attributed, at least in part, to depletion and/or denaturation of the surfactant proteins, especially SP-B and/or SP-C, which are membrane associated proteins that promote vesicle fusion to form more surface active large aggregates. 155,156

In COVID-19-induced ARDS, damage to type II cells inevitably compromises PS synthesis.

The combination of activated macrophages and large numbers of neutrophils recruited to the lung could lead to seven times more phospholipid clearance than the normal lung. 157

This increased phospholipid clearance could, over time, lead to smaller pools of functional surfactant and hence further impair its biophysical function of reducing the alveolar surface tension. 157

Consequently, there is a strong rationale for using exogenous surfactant replacement as a supportive therapy to treat COVID-19 patients, especially those critically ill patients requiring mechanical ventilation.

In a recent report of a small-scale clinical trial, intratracheal administration of Curosurf, a porcine surfactant preparation, to mechanically ventilated patients appeared to promote positive clinical outcomes in four out of five patients. 158


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sun Sep 12, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Aerosol Deposition in the Lung, concluded ...

Pulmonary Surfactant Biomolecular Corona

Upon contacting the alveolar lining fluid, nanoparticles immediately adsorb biomolecules to form the so-called PS biomolecular corona. 144,145,159

This PS corona formed on inhaled aerosols can be considered as a counterpart of the well-recognized protein corona formed on nanoparticles entering the bloodstream. 160−163

It is generally accepted that formation of the biomolecular corona is an entropy-driven process, in which the decrease in entropy of the adsorbed biomolecules is more than compensated for by the increase in entropy of the released water molecules. 160−163

Numerous studies have shown that it is the biophysicochemical properties of the biomolecular corona, rather than that of the pristine particles, that determines the biological identity of the particles entering the human body. 162

The PS corona determines the biodistribution, bioreactivity, biopersistence, retention, translocation, and clearance of inhaled particles. 164,165

Hence, it is important to better understand the effect of this PS biomolecular corona on virus-laden aerosols and how the PS corona affects cellular uptake and cytotoxicity of the virus-laden aerosols.

A number of in vitro studies have demonstrated the formation of PS corona on nanoparticles. 144−146,166

For example, pre-incubating various nanoparticles, for example, single-walled carbon nanotubes, 167 polystyrene, 168 and silica nanoparticles, 169 in PS-containing bronchoalveolar fluid significantly altered the particle uptake by alveolar epithelial cells and macrophages.

The specific chemical composition and molecular conformation of the PS corona depend on the physicochemical characteristics of the particles, such as their size, shape, charge, hydrophobicity, and agglomeration state. 144−147,159,166

Among these properties, hydrophobicity of the nanoparticles appears to play a predominant role in affecting the PS corona (​Figure 66).

However, because of technical difficulties in determining the hydrophobicity of nanoparticles, there are only scarce studies that quantitatively evaluate the effect of particle hydrophobicity on PS biophysical properties 147 and corona formation. 159

A recent invention of an easy-to-use optical method that is capable of measuring the surface free energy of micro- and nanoparticles, as a quantitative measure of hydrophobicity, may be used for quantitatively evaluating the effects of particle hydrophobicity on the biomolecular corona. 170

Other than direct experimentation, coarse-grained molecular dynamics (CGMD) simulations have been used to study the detailed chemical composition and molecular conformation of the PS biomolecular corona.

​Figure 77 shows the PS corona formed on two representative nanoparticles, that is, the hydrophilic silver nanoparticle and the hydrophobic polystyrene nanoparticle. 159

It is clear that the PS corona formed on both nanoparticles consists of a lipoprotein complex that contains both lipids and proteins.

The specific structure of the corona and its molecular conformation are affected by the hydrophobicity of the nanoparticle, with a lipid bilayer formed on the hydrophilic silver nanoparticle and a lipid monolayer formed on the hydrophobic polystyrene nanoparticle.

The two hydrophobic surfactant proteins, SP-B and SP-C, appear to mediate the formation of the corona by directly interacting with lipids in the corona, with SP-B orienting in parallel to the lipid bilayer and the helical fragment of SP-C adopting a transmembrane orientation.

The hydrophilic SP-A was also found associated with the PS corona, mainly via its carbohydrate recognition domain, albeit to a lesser extent of association than SP-B and SP-C.

Although the detailed structure and composition of the PS biomolecular corona are yet to be discovered, both MD simulations and mass spectrometry studies suggest that this corona consists of primarily phospholipids, 167,171 hydrophobic surfactant proteins, 144,145 and hydrophilic surfactant proteins. 123,172,173

A fine balance between electrostatic and hydrophobic interactions may be involved in the formation of the PS biomolecular corona. 145

The PS corona provides the particles with a physicochemical barrier against the environment, reduces the hydrophobicity of the pristine particles, and enhances biorecognition of the particles. 159


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sun Sep 12, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Future Directions in Understanding Aerosol Deposition in the Lung

Although plausible, direct clinical evidence for lung infection by virus-laden aerosols is still lacking.

A common objection against the likelihood of aerosol transmission is the seemingly insufficient viral concentration to cause an infection, after individual airborne viral emissions in the form of a sneeze or a cough, or after continuous airborne viral emissions in the form of speaking or breathing.

The aerosol volumes generated by an individual sneeze and cough are estimated at 1000 and 100 nL, respectively. 174

The aerosol volumes generated by continuous talking and breathing are estimated at 10 and 1 nL/min, respectively. 174

Given the estimated viral load of SARS-CoV-2 in saliva at 1000 virions per nL, 174 a COVID-19 patient could generate at least 1000 airborne SARS-CoV-2 virions per minute.

This rough estimation is in general agreement with experimental measurements of airborne SARS-CoV-2 virions exhaled by COVID-19 patients. 27

After aerosol emission from the reservoir (i.e., an infectious individual) and before aerosol deposition in the deep alveolar regions of a susceptible host, the viral concentration would be continuously diluted during aerosol transmission, decay, and nonspecific deposition. 175

Only a very small fraction of the virus-laden aerosols, likely those within the most breathable size range between 0.5 and 5 μm, may carry SARS-CoV-2 deep to the terminal alveoli.

It is unclear if these amounts of virus are enough to cause lung infection.

The infectious dose of SARS-CoV-2, that is, the average number of virions necessary to cause an infection for COVID-19, is still unknown, but it is estimated at ∼1000 virions by analogy to influenza and other coronaviruses.

If so, there appears to be a high likelihood of airborne transmission of SARS-CoV-2, especially in scenarios involving indoor environments, poor ventilation, prolonged exposure to events involving heavy breathing, or aerosol-generating medical procedures.

It is unknown how interactions with the PS system may affect the behavior of virus-laden aerosols in the lung.

On the one hand, the PS biomolecular corona formed on aerosols reaching the alveolar regions would mediate macrophage clearance of these virus-laden aerosols and their interactions with the ACE2 receptor expressed on the surface of alveolar type II cells.

On the other hand, injury of the surfactant-generating type II cells would worsen lung homeostasis and lead to ARDS symptoms such as alveolar flooding and leakage of the alveolar-capillary barrier, which would further induce surfactant inhibition.

Hence, in spite of its failure in past ARDS trials, exogenous surfactant therapy may become beneficial for patients with COVID-19-induced ARDS, in which damage of the type II cells is an established pathogenesis. 129

Molecular Interactions of SARS-CoV-2 with the ACE2 Receptor

Modeling SARS-CoV-2 Spike Protein and Virion

Fundamental studies about the molecular structure of the spike (S) protein of SARS-CoV-2 using experimental and simulation methods are crucial to the understanding of the virus transmission mechanism in the lung and could contribute to corresponding therapeutic development.

Cryogenic electron microscopy (cryo-EM) 176−179 has provided precise measurements of the structure of the full-length S protein in its trimeric form and the S protein-ACE2 complex with angstrom resolution.

Computer simulations, such as atomistic molecular dynamics (MD), 180−184 have become a crucial approach to complement experimental studies and provide insight into the structural dynamics of the S protein and its interactions with ligands.

Atomistic MD simulations compute trajectories and positions of N interacting atoms in a system by numerically solving Newton’s equations of motion. 185,186

Simulations can probe molecular structures with atomistic resolutions and sample MD from subnanoseconds to microseconds and elucidate molecular interactions. 180−182,187−194

It is worth noting that the computational load using atomistic simulations to study a large-sized system in a long process is prohibitively expensive, particularly with an explicit water model. 181,187,189,195

Based on the current computational capabilities, it is highly challenging to simulate the entire SARS-CoV-2 virion using atomistic MD simulations. 181,196

An atomistic model of a virion such as SARS-CoV-2 and influenza with the diameter of 50–200 nm in water involves more than 200 million atoms, which presents a daunting computational load, requiring about 250,000 processing cores on a high-performance computational cluster. 181,182,187

A recent simulation of four glycosylated S proteins anchored into a patch of viral membrane and embedded in an aqueous solvent with 0.15 M NaCl involved 4.1 million atoms. 182

A viral membrane system of a fully glycosylated full-length S protein with palmitoylation can be built using Membrane Builder of the website: CHARMM-GUI.180 Amaro and co-workers performed atomistic MD simulations, 181 which consisted of ∼1.7 million atoms, to study the structural dynamics of a full-length glycosylated S protein in explicit water.

The protein structure with atomistic resolution can be retrieved from the Protein Data Bank, and the glycomic data about glycan sites of the S protein can be found in literature. 197,198

Molecular interactions between biomolecules and water can be described using common force field parameters, such as Chemistry at Harvard Macromolecular Mechanics (CHARMM),199−201 Assisted Model Building with Energy Refinement (AMBER),202−204 and Optimized Potentials for Liquid Simulations-All Atom (OPLS-AA). 205,206

To achieve higher efficiency, an alternative method is to use coarse-grained (CG) models instead of full-atom models in MD. 207

In CGMD simulations, the resolution representation is reduced, where a number of atoms are simplified as a CG bead and the intramolecular interaction functions are simplified. 207,208

The general MARTINI force field 209 can be used for CG simulations of proteins and lipids.

More precise atomistic or CG models of SARS-CoV-2 and their associated force field parameters, particularly regarding glycans, are still under development, and related information is available in literature 196,210,211 and the open database GitHub.


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sun Sep 12, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

SARS-CoV-2 and Its S Protein

The SARS-CoV-2 genome encodes approximately 25 proteins, including the S protein, which interacts with human ACE2 receptor in the initial stage of infection, two proteases, which cleave viral and human proteins, the RNA polymerase for synthesizing viral RNA, and the RNA-cleaving endoribonuclease (​Figure 88a). 212

The 23 nm long S protein grows on the virus surface as an outward-facing molecular “spike”, which is responsible for human infection (​Figure 88b). 176,178,179,181,212

The S protein, as a class I fusion trimeric protein, has a single polypeptide chain, which can be divided into three topological domains: the head, stalk, and cytoplasmic tail (CT). 182

Cryo-EM experiments and atomistic MD simulation showed that each stalk has three flexible hinge structures, named the hip, knee, and ankle, that link the S protein to the viral membrane (​Figure88c). 182

The hinges provide the heavily glycosylated head with pronounced orientational flexibility, which is critical for the SARS-CoV-2 S protein to scan and interact with the host cell membranes. 182

Each monomer of the S protein is composed of two subunits (S1 and S2), including a total of about 22 N-glycosylation sites, in three different domains (head, stalk, and CT) (​Figure 99a).181,198

The S1 subunit contains the N-terminal domain (NTD) and the receptor binding domains (RBD) that interact with ACE2, and the S2 subunit includes the fusion peptide (FP), central helix (CH), connecting domain (CD), heptad repeat 2 (HR2), and transmembrane (TM) domains. 181

The S2 subunit is responsible for membrane fusion, which is key to viral infection. 176,178,179,181

The viral infection process is triggered when the RBD in the S1 subunit of the S protein binds to the peptidase domain (PD) of ACE2, resulting in subsequent exposure and cleavage of the S2 subunit by host proteases. 176,178,179,213

Compared to SARS-CoV, the SARS-CoV-2 S protein has a furin cleavage site between S1 and S2 subunits (S1/S2) (​Figure 99a,b).181

The S1/S2 site can be cleaved by serine proteases, including trypsin, cathepsins, elastase, the host type 2 transmembrane serine protease, and plasmin, facilitating virus entry into epithelial cells. 133,181

Another proteolytic cleavage is found at site S2 that releases the FP, which penetrates the host cell membrane for fusion. 181

Experimental characterizations 197 showed that the head and stalk domains are extensively covered with glycans (Figure​ 99c–f).

Viral glycosylation plays a wide range of roles in viral pathobiology, such as mediating protein folding and stability and shaping viral tropism and immune recognition. 197

Experimental studies using site-specific mass spectrometric approaches revealed site-specific glycan signatures characteristic of a natively folded trimeric spike. 197

For S proteins, no mannose clusters are detected compared with LASV GPC and HIV-1 Env.

It is worth noting that a recent study using atomistic MD simulations and experiments 181 has shown that besides its functionality of shielding the S protein from antibodies, N-glycans at positions N165 and N234 can modulate the conformational dynamics of the spike’s RBD and further enhance the RBD-ACE2 binding.

SARS-CoV-2 is more contagious than other coronaviruses such as SARS-CoV, most likely due to the stronger binding affinity of its S protein, compared with that of SARS-CoV, to ACE2 receptors, even though SARS-CoV-2 shares ∼80% sequence identity of the SARS-CoV genome.

Composed of three RBDs, which can take different conformations (closed or open) (​Figure10), the S protein head plays a key role in binding receptors. 178,182,214,215

Cryo-EM measurements showed that in the closed conformation, the RBDs are shielded by the N-terminal domains (Figure​ 10 a–c), whereas in the open conformation, one RBD is exposed upward in a receptor-accessible conformation (Figure 10 d–f). 176,178,182

Atomistic MD simulations indicated that an increased number of interdomain salt bridges in the closed state more likely reduces RBD’s mobility compared to in the open state. 216

Free energy landscape analysis also suggested that the salt bridge network is critical to the transition between closed and open states. 216

Moreover, recent studies showed that conformational changes from the open to the closed state can be manipulated by the removal of N-glycans at sites N165 and N234 of the S protein. 181


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sun Sep 12, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Interactions between SARS-CoV-2 S Protein and the ACE2 Receptor

SARS-CoV-2 enters host cells, such as type II alveolar epithelial cells, by binding its S protein with ACE2 receptors expressed on the surface of the host cells. 176,217

Experimental studies reported that the ectodomain of the S protein binds to the ACE2 PD with a dissociation constant (Kd) of ∼15 nM, which demonstrates a higher binding affinity compared to SARS-CoV. 176

The SARS-CoV-2 RBD consists of twisted five-stranded antiparallel β sheets (β1–7), short connecting helices, and loops that form the core (​Figure 11). 218

There are nine cysteine residues in the RBD and eight of them form four pairs of disulfide bonds.

The β sheet structure is stabilized by three cysteine pairs (Cys336-Cys361, Cys379-Cys432, and Cys391-Cys525) in the core.

The ACE2 is a dimer accommodating the RBD in its peptidase domain.

Two lobes of the N-terminal peptidase domain of the ACE2 form the peptide substrate binding site between them.

A concave outer surface in the receptor-binding motif (RBM) of the SARS-CoV-2 RBD binds to the N-terminal helix of the ACE2 (Figure 11). 218

The ACE2 PD uses the α1-helix as its primary contact to bind RBD, though the ACE2 α2-helix and the linker of the β3- and β4-sheets also help. 179,218

Studies of interactions among the amino acids at the interface between the SARS-CoV-2 RBD and the ACE2 PD are crucial to the understanding of virus transmission in the lung and thus for inhibitor design.

Atomistic MD simulations suggested that the S protein-ACE2 binding is governed by hydrogen bonding, electrostatic interactions, and hydrophobic interactions. 181,183

There are three unique stable hydrogen bonds between SARS-CoV-2 S protein’s residues (Tyr449, Gln493, and Gln498) and ACE2 residues (Asp38, Glu35, and Lys353) at the binding interface (​Figure 12a,b), compared to SARS-CoV (Figure 12 c,d). 183

Residues (Tyr473, Ala475, and Phe486) of SARS-CoV-2 form hydrophobic interactions with ACE2 residues (Figure​ 12 b). 183

In the middle region of the binding interface, besides hydrophobic interactions between Phe456 and Tyr489 of the S protein and Thr27 of ACE2, there is a strong salt bridge between Lys417 of the S protein and Asp30 of ACE2 (​Figure 12 b), which increases the S protein-ACE2 electrostatic interactions and therefore further stabilizes the complex structure. 183

Both the S protein and ACE2 contain disulfide bonds between cysteine residues.

Atomistic MD simulation and free energy computation showed that the binding affinity between the S protein and ACE2 becomes weaker when the disulfides of ACE2 are reduced to sulfydryl groups. 219

However, the S protein-ACE2 binding is insensitive to the reduction of disulfide in the S protein’s RBD. 219

Atomistic MD simulations showed that the absolute value of binding free energy between SARS-CoV-2 S protein and human ACE2 is 9.5 kcal/mol greater than that with pangolin ACE2.

Human ACE2 is among the most polymorphous gene and has 317 missense single-nucleotide variations (SNVs). 220

Different SNVs exhibit different binding affinity with the S proteins. 220

Han and Kral performed atomistic MD simulations to study the binding mechanism of peptide inhibitors of SARS-CoV-2. 221

Their study demonstrated that the pronounced structural stability of α-helical peptide bundle extracted from the ACE2 PD is critical to their binding specificity to the RBD of SARS-CoV-2. 221

Their work also suggested that nanoparticle carriers with many such peptides can provide multivalent binding and full coverage to SARS-CoV-2 receptors. 221

Recent CGMD simulations elucidated that the unique polybasic cleavage site of R682RAR685 (i.e., S1/S2 on Figure​ 9 a) of SARS-CoV-2 S protein can enhance the binding affinity of the S protein’s RBD with ACE2 via electrostatic interactions and hydration bonding, although they are about 10 nm away from the RBD (Figure 13). 222

The simulations further demonstrated that the binding of a negatively charged tetrapeptide (GluGluLeuGlu) to one of the three polybasic cleavage sites, which are positively charged, reduces the RBD-ACE2 binding strength by 34%. 222

Consistent with these findings, experiments showed that the mutation of the T678NSPRRAR685 residues in S protein to a variant of T678-IL-R685 lowered the SARS-CoV-2 transduction efficiency in human ACE2-expressing baby hamster kidney cells. 178

To predict immunogenic regions, Colombo and co-workers performed an energetic coupling analysis on a solvated fully glycosylated S protein using atomistic MD simulation trajectories. 184

Their structure-based computational approach showed antibody recognition sites (epitopes) can be correlated to localized regions, which exhibit only low-intensity energetic coupling with the rest of the structure and display large functional conformational transitions. 184

The three most antigenic regions identified were: the upward part of RBDs, including the N-terminal domain, and the central/C-terminal part of the S1A domain (Figure 14 a,b), which corresponds to the binding region of the monoclonal antibody CR3022, the binding region of the new antibody 484, and the antigen-binding fragments of COV57, respectively. 184

From the tertiary structure, the boundary of the N-terminal region locates close to the furin-targeted motif RRAR (Figure 14 c), which is critical to preactivation of the S protein through proteolysis.184

The β-sheet at the initial boundary of the C-terminal domain in the S2 subunit is near the fusion peptide (Figure​ 14 c). 184

It was also found that glycans with stronger energetic coupling protect underlying peptidic epitopes, such as N165-linked oligosaccharide. 184

In contrast, glycans with weaker energetic coupling can enhance antibody recognition, for example, an oligosaccharide fragment bound to N234,184 an observation which was also reported in the work of Amaro and co-workers.181


Site Admin
Posts: 38771
Joined: Thu Aug 30, 2018 1:40 p

Re: Inward Bound - The T'ai Chi Corner

Post by thelivyjr » Sun Sep 12, 2021 1:40 p

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions, continued ...

Yi Y. Zuo, corresponding author, William E. Uspal, corresponding author and Tao Wei corresponding author

Future Directions in Molecular Dynamics Simulations of Virus–Receptor Interactions

A comprehensive understanding of the structural dynamics of SARS-CoV-2 S protein and virion, along with their interactions with the ACE2 receptor at the atomistic and molecular levels, is critical to the development of vaccines and therapeutics.

Atomistic MD simulations, which are based on the well-calibrated force field parameters and experimental measurements of protein structure, have been extensively employed to offer insights into the SARS-CoV-2 biomolecules with atomistic resolution.

However, it presents a major challenge to perform atomistic MD simulations to model the entire virion in an explicit water environment involving hundreds of millions of atoms with a prolonged time scale, due to the heavy computational load.

Most recent simulations only focused on the atomistic modeling of a fully glycosylated S protein without and with the viral membrane instead of an entire virion.

More efficient CGMD simulations, which can present molecular structural changes, are still under development, and reliable force field parameters are not yet available.

Complementary to experimental measurements, atomistic MD simulations have revealed structures of the S protein, including its subunits, RBDs, plausible active sites, and location of surface glycans, structural stability, and its interaction with receptors.

Simulations showed that a SARS-CoV-2 S protein displays an extremely high rotational structural flexibility at the connection regions between the headgroup and cell membranes, which helps the virus locate and interact with the ACE2 acceptor of human cells.

Analysis of the binding free energy showed that the binding affinity between SARS-CoV-2 and the ACE2 receptor is higher than that with SARS-CoV due to its domain structural stability, strong hydrophobic and electrostatic interactions, and unique hydrogen bonds.

Extraordinary efforts have been exerted to predict active sites of the S protein for the purpose of controlling the viral binding.

Glycans covering the cell surface were found not only to protect the virus from antibody recognition but also to contribute to the viral binding with ACE2.

It is noteworthy that several critical issues remain unknown for an in-depth understanding of the virus transmission via aerosols.

It is still unknown how the S protein interacts with aerosol surfaces and biomembranes at the atomistic and molecular levels, which is a critical issue that determines the affinity and transmission of the virus via aerosols.

The effect of water content on the S protein structure on the aerosol surfaces and their correlation with virus survival time need further investigation.

MD simulations hold promise for all of these investigations.

Concluding Remarks

Engineering, physical, and chemical sciences have much to contribute toward the global fight against the COVID-19 pandemic.

In particular, aerosols are strongly suspected to play a significant role in the rapid spread of the disease.

Key steps in the hypothesized airborne transmission pathway, from the initial aerosol ejection by a cough or exhalation to deposition of the virus-laden aerosols in the lung and to the binding of the S protein of SARS-CoV-2 to the ACE2 receptor on the surface of a host cell, are largely governed by physicochemical mechanisms, that is, the exchange of mass, energy, and momentum between virions, aerosol particles, and the surrounding environment.

Viewing COVID-19 transmission through the lens of bioaerosol mechanics may lead to the development of mitigation approaches, such as rational ventilation protocols, antiviral materials and coatings, and even new therapeutic interventions, such as the surfactant therapy to alleviate COVID-19-induced ARDS symptoms.

Although we have no expectation of a “magic bullet” against COVID-19, widespread adoption of physically and chemically motivated and justified mitigation strategies would effectively reduce the virus reproduction rate R0, giving relief to stressed medical systems and healthcare providers and saving hundreds of thousands of lives.

Indeed, after the first few months of the pandemic, most governments and individuals have come to understand and accept the effectiveness of simple face coverings, which may greatly reduce the chance of airborne transmission.

The pandemic has also revealed gaps in our understanding about particulate matter and bioaerosol mechanics.

For instance, dispersal of infectious virus-laden aerosols in air and their deposition onto general surfaces and various regions of the lung are multiscale and multiphysical processes.

Accordingly, it is challenging to identify and select all relevant contributing physicochemical processes and resolve them under the framework of one unified computational or analytical model.

On the one hand, the CFD approach has demonstrated success in modeling aerosol transmission by air flow in built environments.

On the other hand, MD simulations have demonstrated particular impacts on revealing specific virus–receptor interactions at the atomistic and molecular scales.

Effective multiscale interfacing and exchanging parameters between the micro-to-macroscale CFD modeling and the atomistic-to-molecular scale MD simulations, however, have proven to be technically challenging.

Several key fundamental questions remain open.

For example, correlations between the virus viability and the S protein denaturation on the aerosol surface are still unknown.

The effects of water content and physicochemical properties of the aerosols on the S protein structural dynamics and virion packing remain largely undefined.

A comprehensive understanding of the complete pathway for the airborne transmission of COVID-19 is yet to be established.

Although plausible, direct clinical evidence for lung infection by virus-laden aerosols is still lacking.

It is also unknown whether the virus-laden aerosols, especially those within the most breathable size range between 0.5 and 5 μm, can carry SARS-CoV-2 deep to the terminal alveoli.

However, if this transmission pathway does exist, it would bypass the mucociliary clearance and incubation period of the virus in the upper airways and thus cause direct detrimental effects on the alveolar regions of the lung, which would greatly change the progression of the disease.

This could also in part explain the widely differential symptoms and clinical outcomes for COVID-19 patients even in a seemingly homogeneous population.

In the alveolar region, interactions with the endogenous PS would determine the subsequent macrophage clearance of the virus-laden aerosols and interactions between the virus and ACE2 receptors expressed on the surface of alveolar type II cells.

Synergetic effects between direct surfactant inhibition and reduced surfactant metabolism would worsen lung homeostasis and cause alveolar collapse and instability, thus, in the most severe cases, leading to the ARDS symptoms of pulmonary edema and alveolar flooding.

Hence, there is a strong rationale for surfactant replacement as a plausible supportive therapy to treat COVID-19 patients.

We hope this review provides useful multidisciplinary knowledge that helps understanding and ultimately contributes to curbing the airborne transmission of COVID-19.


We thank Drs. Fred Possmayer and Neil Uspal for valuable discussion about our paper.

This research was supported by the National Science Foundation (NSF) grant numbers CBET-1604119 and CBET-2011317 to Y.Y.Z., the American Chemical Society (ACS) Petroleum Research Fund grant number 60809-DNI9 to W.E.U., and NSF grant numbers CBET-1831559 and CBET-1943999 to T.W.



COVID-19 coronavirus disease 2019, is the disease that is causing a global pandemic

SARS-CoV-2 severe acute respiratory syndrome coronavirus 2, is the virus that causes COVID-19

ACE2 angiotensin converting enzyme 2, is the human receptor responsible for SARS-CoV-2 entry into cells

RBD receptor binding domain, is a key part of the viral spike protein that allows it to dock to the receptor for cell entry

PS pulmonary surfactant, is a lipid–protein complex that is synthesized by the alveolar type II epithelial cells and is responsible for innate immunity and surface tension reduction in the lung

MD molecular dynamics, is a computer simulation method for studying the physical movements and interactions between atoms and molecules

CFD computational fluid dynamics, is a branch of fluid mechanics concerned with numerical methods for the solution of physical problems involving fluid flow


The authors declare no competing financial interest.


Post Reply