COVID-19 pandemic

The pandemic caused by the new coronavirus SARS-CoV-2 appeared in Wuhan, the capital of the Chinese province of Hubei, in December 2019.

Immediately after China, SARS-CoV-2 infection hit Italy hard, causing thousands of victims to date. Globally, as of 27 August 2021, there have been 214,468,601 confirmed cases of COVID-19, including 4,470,969 deaths, reported to WHO (Source WHO).

At the beginning of last March 2020, a report by World Wildlife Fund (WWF) Italy “Pandemics and the boomerang effect of the destruction of ecosystems” Pratesi et al. (2020) aroused my curiosity because it proposed a connection between two apparently unrelated topics: “ecology and pandemics.” On the same issues, but from another point of view, an interesting video was published on March 20 2020 by the climate change scientist Antonello Pasini at the Italian National Research Council -CNR where it was speculated on a (impossible) cause-effect relationship between climate change and SARS-CoV-2 infection Pasini A and Italian National Research Council -CNR (2020).

So, both the WWF report and the CNR video discussed an impossible correlation between environment damage and impact (pollution, climate change) and the COVID-19 pandemic. In any case, the topic immediately aroused my interest as a theoretical physicist and since the end of February 2020, I was wondering how the SARS-CoV-2 new coronavirus could be directly affected by the environment and what effects this correlation could have on the COVID-19 pandemic.

The first step was therefore to study the thermodynamics of the new coronavirus virus in the environment.

The next step was to keep in mind that in the study of the thermodynamics of the new coronavirus SARS-CoV-2 in air, the presence of pollen, bacteria, other viruses, spores, dust, but more particularly, that of nanoparticles of health interest at the same scale threshold as the spike proteins of the pandemic virus, such as particulate matter, cannot be neglected.

This also applies to simulations and biologic computational models. Otherwise, we run into serious systematic errors.

The result of these observations was the theory of the energy landscape proposed here.

At the beginning of March 2020 I published my work as a preprint, wishing to immediately make my observations available to the community of colleagues, operators in the sector, doctors, journalists and social networkers and at the same time sent to the authorities of the Italian Government.

Now, after a year, I am truly happy and grateful to the many colleagues and research institutes who have found useful my methodological notes. See the position paper of the Italian Society of Environmental Medicine (SIMA) issued on March 16, 2020 and a survey conducted in the USA by the Harvard School of Public Health Wu, X., Nethery, RC, Sabath, MB, Braun, D. and Dominici, F. (2020). They really went beyond the scope of my work as they speculated on a correlation between air pollution and infection, according to which the areas most affected by the infection were also the most polluted. But our paper is not concerned with this. We just proposed some methodological notes which lead us to the formulation of a general mathematical apparatus, suitable to explain at the molecular level the energetic configurations of the quasispecies of the pandemic virus SARS-CoV-2 in its environment. We focus complexes between the viral particle and other objects in its environment, at the scale threshold of the spikes of the viral particle, which can be dangerous to human health, as particulate matter. Then, we wondered if such complexes can lead to the generation of more aggressive viral variants and how to predict their populations and energy configurations, in order to plan an adequate prophylaxis.

China and Italy

Let’s start on China and Italy the two nations tragically hit first by the epidemic. In the Chinese case, I will focus on Wuhan, the capital of Hubei Province, while in the Italian case I will focus on Northern Italy, with particular attention to the case of Lombardia. As an example, Lombardia has roughly the same population of Wuhan (about 10 million people) (The Government of Wuhan 2020; Setti et al. 2020; Regione Lombardia 2020). Both have been affected by the pandemic in the most devastating way. Both have suffered a severe environmental pollution situation in recent years, China with greater severity and impact on the health of population (Wang et al. 2019; CHINA National Environmental Monitoring Centre—CNEMC 2020; Rongshan Wu et al. 2018; Rohde and Muller 2015; He et al. 2017; Liang et al. 2019; Zhao et al. 2018; Liu et al. 2017; Anenberg et al. 2019; Zhu et al. 2018; Kuerban et al. 2020; Zhang and Cao 2015; Xu et al. 2020; Ren et al. 2017). Most papers analyze data from Chinese government agencies, as Wang et al. (2019), which refers to the Chinese Ministry of the Environment and Ecology, where we find the latest (2017) report of Ecology Ministry of Ecology of the People’s Republic of China (2017) translated in English (See Figs. 1 and 2).

Fig. 1
figure 1

Percentage of 338 Chinese cities with different concentrations of six major pollutants in 2017 [courtesy of the Chinese Ministry of the Environment and Ecology of Ecology and of the People’s Republic of China (2017)]. Note the high percentage values of PM10 and PM2.5

Fig. 2
figure 2

Air quality comprehensive index and primary pollutants of 74 Cities in 2017 [courtesy of the Chinese Ministry of the Environment and Ecology of Ecology and of the People’s Republic of China (2017)]. Note the preponderance of PM2.5 among the primary pollutants, also for the city of Wuhan

The Italian situation is summarized by the European Environment Agency [See Fig. 3 and 4 European Environment Agency (2019–2020)].

Fig. 3
figure 3

Italy—air pollution country fact sheet 2019 (courtesy of European Environment Agency 2019b)

Fig. 4
figure 4

The PM2.5 annual mean concentrations in 2017 in Italy compared to the EU annual limit value (25 µg/m3) (Left) The PM10 percentile 90.41 concentrations in 2017 in Italy compared to the EU daily limit value (50 µg/m3) (Right). Note that the highest values of Particulate Matter are reached in the regions of Northern Italy, particularly for PM 2.5. Courtesy of the European Environment Agency (2019b)

Note that the highest values of particulate matter are reached in the regions of Northern Italy (See Fig. 4), (See Fig. 5).

Fig. 5
figure 5

Contagion maps in Italy. Top left—situation on 25 April 2020. Top right—situation on 12 June 2020—center situation update March 11, 2021 [courtesy of Corriere della Sera della Sera (2021)]. If we compare these maps with the previous Fig. 4, we immediately notice that the areas most affected by the pandemic are the most polluted. Here, we find the same type of correlation studied by T.H. Chan School of Public Health at Harvard Wu et al. (2020)

Particulate matter

Among the particles of interest for human health at the scale threshold of the pandemic virus, we focus our study on particulate matter, whose damage to the human airways is well known and studied. The definition of particulate matter (hereinafter referred to as “PM”) and its health effects are well-known. We recall the definition given by the U.S. Environmental Protection Agency (EPA) United States Environmental Protection Agency-EPA (2020): PM is a mixture of solid particles and liquid droplets in air.

We then have:

PM10: particles, with diameters that are generally 10 μm and smaller.

PM2.5: fine particles, with diameters that are generally 2.5 μm and smaller.

Ultrafine particulate (UFPs) are defined as “respirable dust,” that is, capable of penetrating deeply into the lungs and undergo interstitialization, Friedlander et al. (2003)

We then have:

PM1, with a diameter of less than 1 μm;

PM0.1, with a diameter of less than 0.1 μm.

The health effects of PM are devastating, particularly in the lungs. PM2.5 is a thoracic powder capable of penetrating deeply into the lungs during breathing, while the ultrafine particulate matter (UFPs) is in all respects a “respirable dust,” which settles in the lungs, where the nanoparticles have the ability to penetrate the tissues and undergo interstitialization or to be absorbed directly into the bloodstream (Friedlander et al. 2003; Scientific Committee on Emerging and Newly Identified Health Risks-SCENIHR 2006; Scientific Committee on Health and Environmental Risks—SCHER 2005; European Environmental Agency (2019; Xing et al. 2016; MacNee et al. 2003). The studies on the damage caused by PM on the respiratory system generally take into account free radical peroxidation, imbalanced calcium homeostasis, inflammatory lesions, etc. In particular, the peroxidation of free radicals is one of the main factors causing oxidative DNA damage. Indeed, when damaged DNA is not repaired effectively over time, it can induce carcinogenesis, mutagenesis, and other irreversible damage (SCHER 2005; European Environment Agency 2019; Xing et al. 2016; MacNee et al. 2003). Moreover, oxidative stress within cells is an additional factor leading to transcription of pro-inflammatory genes (Xing et al. 2016; MacNee et al. 2003).

SARS-CoV-2

SARS-CoV-2 virus is a beta-coronavirus as SARS-CoV; its RNA genome is about 82% identical to the SARS-CoV (Zhang et al. 2020). SARS-CoV-2 has round or elliptic and often pleomorphic form, and a diameter of circa 60–140 nm (Cascella et al. 2020). Angiotensin converting enzyme 2 (ACE2) is the cell receptor for both the SARS-CoV and the new SARS-CoV-2 (Yan et al. 2020). Now, we know the X-ray structures of the unliganded SARS-CoV-2 Mpro and its complex with an α-ketoamide inhibitor (Zhang et al. 2020) and many other features. As virus SARS-CoV, the new virus SARS-CoV-2 also enters cells through an endocytic pathway (Hongliang Wang et al. 2008), as has recently been confirmed by many studies such as an interesting “snapshot” of COVID-19 by Blake Oberfield et al. (Harvard Medical School and MIT) (Oberfield et al. 2020), by Alexandra C. Walls et Al. (Washington University, Seattle, USA), Walls et al. (2020) and many others. With each passing day, new details emerge. Very recently (March, 2021), Wendy S. Barclay of Imperial College London hypothesized that new variants of the pandemic coronavirus such as the “Delta” variant, could specialize in an endosome-independent entry, exploiting a selective advantage in primary human airway epithelial cells, resulting from a higher differential expression of the cell surface protease TMPRSS2 (Peacock et al. 2021). But we will have to wait to know if the variants of the pandemic virus that specialize in endosome-independent entry are also expressed in the human population and, in this case, what will be their percentage ratio with respect to the viral species that follow the endocytic pathways. Furthermore we will have to study in detail the structure of SARS-CoV-2 spike proteins to understand what is the interaction of the spike with the nanoparticles of the environment and if this interaction can somehow influence and impact the cleavage sites of the protein and other crucial sites suitable to mediate host’s cell entry. Another interesting topic concerns the role of glycans, which camouflage the spikes, making the viral particles “stealth” for the host’s immune system. What are the interactions of glycans with particles and other objects in the environment? And beyond. Can nanoparticles of the environment at the same scale as spikes have the same “stealth” effect as glycans? Can they hide the pandemic virus from the immune system? We underline that the methodological observations carried out, also concern the studies that speculate on the endosome-independent cell entry, because they are affected by the same systematic errors already discussed. This confirms, in case there were still any doubts, the urgency to tackle the methodological problem once and for all. We will return to this fundamental topic soon. We acknowledge the excellent actions promoted by CERN community such as “CERN against Covid-19” CERN, Genève, CH (2020a), “Zenodo” (Open Science against COVID-19, which I also adhere to) CERN, Genève, CH (2020c), and the platform folding@home, which is a distributed computing project for simulating the dynamics of protein molecules CERN, Genève, CH (2020b).

Finally, I remember the very updated report of the Covid-19 Commission of the Accademia Nazionale dei Lincei, directed by Prof. Giorgio Parisi Accademia Nazionale dei Lincei (2021).

In the following figure, courtesy of Dimmock et al. (2016), we recall the main processes of entry of animal viruses genomes into cells (Fig. 6).

Fig. 6
figure 6

Entry of animal virus genomes into cells [courtesy of NJ Dimmock et al. Dimmock et al. (2016)]

In the present case, SARS-CoV-2 will exploit the process (b) of the previous Fig. 5 (receptor-mediated endocytosis), apart from the recent speculations of Wendy S. Barclay of Imperial College London concerning endosome-independent cell entry pathways of new variants of the pandemic virus (see above), which, however, has yet to be found and verified in the human population. Thus, SARS-CoV-2 is an enveloped virus that enters the cell in a vesicle (endocytosis). The vesicle subsequently merges with an endosome. The low pH environment inside the endosome promotes conformational changes in viral proteins that lead to the fusion of the lipid bilayer of the virus with the endosome allowing the viral genome and its associated proteins to enter the cytoplasm.

This is the big picture. However, in our case, things are a bit more bizarre and require a more advanced formalism.

SARS-CoV-2 and particulate matter

The generally accepted description of endocytosis, which we also find in the text by Dimmock et al. (2016) is wrong because nothing is specified about the environment, as if a viral particle could exist without its environment. Living matter cannot exist in a vacuum. Even more so the pandemic virus, given that it is an obligate parasite. Therefore it cannot exist in se ipsum. We are thus facing a complete violation, both logic and functional of the scientific method, that is involving a violations of the conditions of necessity and sufficiency. Remembering Bertrand Russell (1870–1970) and Karl Popper (1902, 1994), we can say that a theory or a model of living matter in a vacuum, leads to a nonsense, whose formal degradation is comparable with the brocardo “nihil sine causa datum est”; which  precipitates the theory (and the experiment that employs it) into a reductio ad absurdum. At this purpose, the present author is working on the so-called Russell categories or “Russellians”, defined to avoid these methodological nonsenses. Moreover, we noted that the same nonsense is achieved with mathematical simulations of both the background and the environment, offered by commercial software (sold with the equipment), and by taylor-made computational biology algorithms. Just on a computational level, this logical and methodological nonsense is a “blue error”, that is, a major school error, which becomes a serious systematic error that hits the statistics of an experiment to death.

So let’s begin our study of the pandemic virus in its air environment.

Consider a given concentration of viral particles per unit of air volume. As an example, among the many works, including Italian ones, we choose a recent paper by a Wuhan doctor, Liu et. al. (2020), about the concentration of different aerodynamic dimensions of SARS-CoV-2 aerosols in Wuhan hospitals during the lockdown, so as not to overestimate the concentrations.

Consider thus a spherical neighborhood (ball) \({\mathcal{B}}\left( r \right)\) of radius r > 0 around a viral particle such that r0 < r < pr0 where r0 is the average diameter of a viral particle and p = 1,…., \(L \in N\), where L is the side of a cubic box with a volume of 1 μm3. For each p > 0, there will be a ball \({\mathcal{B}}\left( r \right)\) in which a number N ≠ 0 of PM (particulate matter) particles will fall.

The following Fig. 7 shows a snapshot of a possible scenario of the ball (SARS-CoV-2 with the environment) for a low density of the UFP dust from 0.025 to 0.100 µm.

Fig. 7
figure 7

Snapshot of a simulation of a SARS-CoV-2 ball with its environment for a low density UFP particulate dust from 0.025 to 0.100 μm

As shown in Fig. 8 below, the endocytosis process captures the entire ball containing the virus with its micro-environment inside a vesicle, deforming the ball shape into an “ellipsoid” (see Fig. 8).

Fig. 8
figure 8

Snapshot of a simulation of the endocytosis of SARS-CoV-2 with its environment (ball) low density—UFP particulate dust 0.025 μm, 0.050 μm and − 0.100 μm

From the moment the endocytosis process takes place and the whole ball is incorporated into the vesicle, a thermodynamic isolated system is created inside the vesicle (See Fig. 9).

Fig. 9
figure 9

SARS-CoV-2 with its environment inside a vesicle, as an isolated thermodynamic system. Universe (SARS-CoV-2 with its environment) and the rest of the universe (environment). Low UFP dust 0.025–0.100 μm

A synergistic cooperation between SARS-CoV-2 and particulate matter

We see below how important it is to consider and study the environment where the pandemic virus combines with the host. Conversely, neglect of the environment leads to nonsenses and serious systematic errors throughout the pandemic COVID-19 research chain, as we said earlier.

Consider that the new coronavirus SARS-CoV-2 is an obligate parasite and our species is its host, that is, the virus needs us to be able to replicate itself.

Let’s focus on the parasite’s attack strategies in its hunting environment, where the host lives.

The adaptation of the parasite to the host’s environment conditions its attack strategy also at the molecular level. The global spread of the pandemic virus SARS-CoV-2 tells us that the parasite specialized in exploiting to its own advantage the environment. At this point, recall that the atmospheric environment is not an empty box as is currently modeled, but it contains a huge amount of items, pollen, bacteria, other viruses, spores, dust, but more particularly, that of nanoparticles of health interest at the same scale threshold of the spike proteins of the pandemic virus, as particulate matter. An important point to emphasize is that none of these objects should and can be overlooked, because each component of the environment is not additive, but has a non-linear dependence between the other components and the pandemic virus, being the environment in the presence of the virus (i.e. the thermodynamic universe) a complex system that cannot be reduced.

To get an idea of how the parasite specialized in exploiting the host's environment to its advantage, let's consider a case belonging to human technology. This example can help us to have a qualitative and quantitative idea of the dynamics we are trying to understand.

Much human research on nanotechnologies has sought to exploit the ability of certain compounds such as mesoporous silica nanomaterials (MSNs) to penetrate human cells through endocytosis, to perform intracellular detection and monitoring functions or to release genetic material or drugs inside of the cell Trewyn et al. (2008).

As we will see, this is precisely the behavioral strategy that we expect from the interaction of the pandemic virus with the objects in the environment. Among the set of competitive processes that the pandemic virus can carry out by exploiting the environment, the one that has the greatest selective advantage, as we will demonstrate shortly, is the process of endocytosis through nanoparticles. Likewise, in the host’s environment, at the threshold scale of the spikes, particulate matter is, as we have seen, the one that has the most devastating impact on health. Consider therefore the effect of nanoparticles of particulate matter, which we indicate as “PM” (See Fig. 8). Endocytosis through PM is performed by the parasite through a complex \(\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right]\), where we have chosen to indicate the equilibrium complex between the pandemic and PM with the operator symbol “\( \oplus\)”. With this “direct sum” symbol we mean an equilibrium complex, where all the interactions between the components of the thermodynamic universe have stabilized. We will see shortly that this occurs inside a vesicle. Among these we consider an annihilation operator Aleph ℵ=ℵ(Q) on a certain cell set |a>, such that, for a certain element K(Q) (belonging to the thermodynamic universe), the average destructive action against the cell set is maximum, per volume, sample, per unit of time. We will now see that Q is the order parameter of our model. That is <ℵ {[SARS-CoV-2]⊕[K]}> is maximum if the K(Q) component coincides with [PM], that is for Q ≈ 1. We can take a 1mm3 cubic sample of human lung or upper airway epithelial cells. Now we will see how this complex will work out and what its selective strength or advantage will be.

First, let's look at some characteristics that will be useful to us when we will deal with entropy of configurations. From the parasite's point of view, this strategy reveals intelligence, adaptability. Can we talk about intelligence for a virus? Maybe this is our best chance. While from the host's point of view, the [SARS-CoV-2] \( \oplus\)[PM] complex is a devastating weapon, not only because it introduces lethal poisons into cells, weakening the organism and irreparably damaging cellular structures and architectures (i.e., the pulmonary epithelial cells), but also because it exploits the ability of the particulate dust to promptly alter molecular structures and functions and induce genetic modifications. In the end, it is as if the host is being bombed by a “smart bomb.”. From our point of view as researchers, this strategy expresses a strong selective advantage, therefore a necessity, vital for the obligate parasite, because it maximizes the reproductive efficacy of the viral species. Later we will develop a formalism (the energy landscape) that will help us to express these concepts in mathematical language.

Our “smart bombs” are thus \(\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right]\) complexes entropy, close packed, highly cooperative and synergistic structures. As we will see, the dynamics leading to the formation of the complexes will be governed by a control parameter. In this way, the spikes of the new coronavirus will strictly interact with the particulate dust. Since the system is isolated inside a vesicle, that is the thermodynamic universe (See Fig. 9), all the interaction processes (mechanical-thermal-chemical) between the viral particle and the particulate dust will go to thermodynamic equilibrium, thus stabilizing the compound in a synergistic cooperative single entity, that is the complex.

For the perspective (and health) of the host, complexes \(\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right]\) thus achieve a synergistic and highly cooperative smart optimization of the cell demolition strategies, which exceed the individual potentials of both viral and particulate particles, as soon as they are considered as independent forcing agents on the host cell. With the previous formalism we can say that ℵ{[SARS-CoV-2]⊕[PM]} ≥ ℵ[SARS-CoV-2]+ ℵ[PM] per volume and sample, per unit of time

The presence of multiple non-covalent interactions between SARS-CoV-2 and PM is a source of cooperativity between them. When these interactions cooperate in the isolated system of the vesicle, a single stable conformation is created: a minimum entropy complex given by SARS-CoV-2 and PM: \(\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right]\) . The complex created is therefore much stronger than might be expected from the sum of their individual strengths because we are dealing with complex systems.

Let’s see from the theoretical point of view how the complex formation mechanism can be represented and described by a mathematical model.

The energy landscape of [SARS-CoV-2] \( \oplus\) [PMs] complexes

The theoretical structure that is the most appropriate to describe the complexes \(\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right]\) is that of the “rugged funnel” energy landscape by Frauenfelder (2005), Frauenfelder et al. (1991), Frauenfelder (1997), Wolynes and Luthey-Schulten (1997).

Molecular biophysical interactions between viral particles and PM dust are mediated by SARS-CoV-2 proteins (e.g., spike proteins); therefore, the energy rugged funnel is certainly the best possible model.

Postulate 1. Rugged energy landscapes

The energy landscape of SARS-CoV-2 and PM dust is  rugged due to the possibility of continuous contacts and interactions, among which we consider the non-covalent interactions that combine cooperatively to create complexes. We call them “complex interactions.” When non-complex interactions are made, the energy contributions are random.

Postulate 2. The “principle of minimal frustration”

The rugged landscape is not globally flat but has a preferred direction of flow. There is an overall regular slope of the energy landscape, which is a symptom of cooperativity and synergy. The interactions and conformational energies between SARS-CoV-2 and PM dust are highly stabilizing, until reaching the minimum of entropy.

Consider therefore the rugged funnel landscape diagram in the following Fig. 10.

Fig. 10
figure 10

Energy rugged funnel landscape for [SARS-CoV-2] \( \oplus\) [PMs] complex formation

We introduce the order parameter “Q” as the percentage of complex-like interactions between SARS-CoV-2 and PM dust, according to the different size and density of PM (as an example, in the scenario of the previous section, we considered a med-hi density for three dust sizes, from 0.025 to 0.100 μm) (See Figs. 8, 9).

The functional dependence of the thermal average energy will be \(\overline{E}\left( Q \right)\), the ruggedness will be \(\sqrt {{\Delta }E^{2} \left( Q \right)}\); the density of states will be \({\Omega }\left( {E,Q} \right)\); the entropy will be \(S\left( {E,Q} \right)\) and the local glass transition temperature will be \(T_{g} \left( Q \right)\).

In the “ball” stages, as the interaction between the viral particles and dust are random (Q ≈ 0), there will be many states and the ensemble average is most appropriate. The hopping rate R between micro-states at each stage is

$$R \approx e^{{ - {\Delta }E^{2} /k_{B} T^{2} }}$$

where kB is the Boltzmann constant and T is the absolute temperature.

The probability distribution of energies in the interval δE at any position in the funnel is a Gaussian centered about the mean energy

$$P\left( E \right)\delta E = \frac{1}{{\sqrt {2\pi \Delta E^{2} \left( Q \right)} }}\exp \left\{ { - \frac{{\left( {E - \overline{E}\left( Q \right)} \right)^{2} }}{{2\Delta E^{2} \left( Q \right)}}} \right\}\delta E$$

Now, we will define the configuration space.

Consider a spike protein of SARS-CoV-2. If we consider γ dust configurations per spike in the ball state (Q ≈ 0), then the total number of configurations for a viral particle with N spikes will be \({\Omega }_{0} \left( Q \right) = \gamma^{N} .\)

When the system approaches the state of the complexes, the number of configurations will decrease, until the only configurations will be those of the complexes (Q ≈ 1).

If \({\Omega }_{0} \left( Q \right)\) is the number of structures that have a fraction Q of the contacts similar to those in the complexes, then an approximation will be

$${\Omega }_{0} \left( Q \right) = \gamma^{{*N\left( {1 - Q} \right)}}$$

where γ* is an experimental value for compact configurations. In the case of protein folding, we have γ* ≃ 1.5 (Frauenfelder et al. 1991; Frauenfelder 1997; Wolynes and Luthey-Schulten 1997).

The entropy will decrease as the complex is approached.

$$S_{0} \left( Q \right) = k_{B} \log \Omega_{0} \simeq k_{B} N\left( {1 - Q} \right)\log \gamma *$$

So that in the minimum of the rugged funnel (See Fig. 10 above) we would have compact, closed packed configurations, which precisely are our complexes

$$\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right].$$

The density of conformational states with energy E and the order parameter Q will be

$${\Omega }\left( {E,Q} \right) = {\Omega }_{0} \left( Q \right)P\left( E \right)$$

Such that the total entropy will be

$$S\left( {E,Q} \right) = S_{0} \left( Q \right) - k_{B} \frac{{\left( {E - \overline{E}\left( Q \right)} \right)^{2} }}{{2\Delta E^{2} \left( Q \right)}} - k_{B} \log \left( {\frac{{\sqrt {2\pi \Delta E^{2} \left( Q \right)} }}{\delta E}} \right)$$

where δE is large relative to the spacings between energy levels, but small relative to \(\sqrt {{\Delta }E^{2} \left( Q \right)}\). At thermodynamic equilibrium (inside the vesicle), the most probable energy will be

$$E_{{{\text{best}}}} \left( Q \right) = \overline{E}\left( Q \right) - \frac{{\Delta E^{2} \left( Q \right)}}{{k_{B} T}}$$

The entropy associated with the most probable energy will be

$$S\left( {E_{{{\text{best}}}} ,Q} \right) = S_{0} \left( Q \right) - \frac{{\Delta E^{2} \left( Q \right)}}{{2k_{B} T^{2} }}$$

The free energy will be

$$F\left( {Q,T} \right) = E_{{{\text{best}}}} \left( Q \right) - TS\left( {E_{{{\text{best}}}} ,Q} \right) = \overline{E}\left( Q \right) - \frac{{\Delta E^{2} \left( Q \right)}}{{2k_{B} T}} - TS_{0} \left( Q \right)$$

Definitely, free energy has relative minima either near the complex states (Q ≈ 1), or near the Ball states (Q ≈ 0). If we neglect the entropy of the complex states, its free energy will be equal to its internal energy. We thus have states of minimal entropy where the complexes \(\left[ {{\text{SARS}} - {\text{CoV}} - 2} \right] \oplus \left[ {{\text{PM}}} \right]\) are stabilized in compact, closed packed structures.

Finally, the local glass transition temperature will be

$$T_{g} \left( Q \right) = \sqrt {\frac{{\Delta E^{2} \left( Q \right)}}{{2k_{B} S_{0} \left( Q \right)}}}$$

Recall that a transition to a “glassy state” occurs at the temperature where so few states are available such that the system remains frozen in one of few distinct states (Wolynes and Luthey-Schulten 1997; Mezard et al. 2004; Fischer and Hertz 1991; Parisi 1988); ; ; (See Fig. 10 above).

In conclusion, we have complex [SARS-CoV-2] \(\oplus\) [PM], stabilized by strongly cooperative and synergistic conformation interactions and energies.

Evolution of SARS-CoV-2: mutations and quasi-species

The evolution of the pandemic virus depends on mutations that give the parasite a selective advantage. A large viral population of SARS-CoV-2 with different genomic sequences is generally called a “quasi-species” or “quasispecies” (Dimmock et al. 2016).

In our landscape model, we can draw a lower risk margin to measure the ability of a quasispecies of SARS-CoV-2 to create a devastating impact on the host.

We can say that a quasispecies of SARS-CoV-2 have sufficient selective advantage to create a devastating impact on the host if the following conditions on entropy, energy and conformations are met, respectively, as follows:

  1. 1.

    \(T\left( Q \right) \le T_{g} \left( Q \right)\) (temperature of viral quasispecies)

  2. 2.

    the most probably energy of quasispecies is \(E_{{{\text{best}}}} \left( Q \right) = \overline{E}\left( Q \right) - \frac{{\Delta E^{2} \left( Q \right)}}{{k_{B} T}}\)

  3. 3.

    the energy variability of quasispecies is given by \(\sqrt {\Delta E^{2} \left( Q \right)}\) (ruggedness)

All quasispecies that do not meet these conditions do not pose a danger to the host. In particular, the temperature condition (1) means that a dangerous quasispecies must be composed of complexes [SARS-CoV-2] \(\oplus\) [PM]. It is clear that the most dangerous quasi-species are found in the neighborhood of Q ≈ 1 because they are the most stable and resistant. They therefore have the maximum mean destructive action on human health <ℵ(Q ≈ 1)>. Conversely, we must not fear the quasi-species around Q ≈ 0.

In any case, (3) defines the variability of each quasi-species such as

$${\text{Ruggedness}} = \sqrt {\Delta E^{2} \left( Q \right)}$$

where

$$\sqrt {\Delta E^{2} \left( Q \right)} \gg \delta E$$

where δE is large with respect to the distances between the energy levels, but small with respect to the ruggedness. Finally, our concept of selective advantage will be W = W(S,ΔE,Q) ≈ Tg(Q)2.

Landscape prophylaxis

For effective prophylaxis, we can avoid studying stochastic states (Q ≈ 0), focusing only on the viral quasi-species with

  1. 1.

    \(T\left( Q \right) \le T_{g} \left( Q \right)\)

  2. 2.

    Q ≈ 1

Conclusions and discussion

In this work, some methodological notes have been addressed concerning the research on the COVID-19 pandemic caused by the new coronavirus SARS-CoV-2.

From a careful examination of the scientific literature in subiecta materia, we have noticed that the current research and characterization of SARS-CoV-2 pandemic virus lack a correct thermodynamic description. In some cases, as in the works referred to here, thermodynamics is completely lacking. In other cases, an attempt has been made to calculate and measure some thermodynamic quantities, such as free energy and temperature, but in the laboratory, i.e. under conditions other than the experimental ones, or (this is the worst error) by algorithmic simulations based on manually adjusted boundary conditions or copied from “similar” conditions, but different from the experimental ones and taken in different spaces and times. In other cases, the measurement conditions were very far from equilibrium, or not reproducible, or vitiated by inadequate statistics. A review of these cases will be the subject of a future paper. But the most worrying thing from a methodological point of view is the study of the viral particle in vacuum, which causes a chain of logical absurdities and, therefore, systematic errors that can cast doubt on the results of an entire research.

The viral particle is described and represented in se ipsum, i.e., independently of the environment, as if the viral particle could exist without its environment, in a vacuum (nihil sine causa datum est). Not only in theoretical descriptions and experimental results, but also in biological and biophysical-computational simulations, unfortunately, the origin of this abomination is well known and must be sought in x-ray crystallography of proteins. From here, the crystallographic technique and the software, commercial or tailor-made, of molecular and macro-molecular simulations are unfortunately taken and exported tout-court in the physical world of living matter. From here it is easy for many to build toy models of the dynamics of a viral particle in an empty box. In some cases, the box is also missing, that is, the space-time references, so that the dynamics of the pandemic virus is also found free from a particular region of space-time. But we know, and the present research bears witness to it, that the distribution of the pandemic also depends on geography and on the particular climatic conditions. So, the methodological problems are therefore many. In the present work, we have focused on the thermodynamic environment of the pandemic virus. This seems like a good start. We therefore studied the thermodynamics of SARS-CoV-2 in its “hunting” environment, from air transport to endosome and endosome-independent cell entry pathways. In the study of the thermodynamics of the air environment of the pandemic virus SARS-CoV-2 in air, the presence of pollen, bacteria, other viruses, spores, dust, but more particularly, that of nanoparticles of health interest at the same scale threshold as the spike proteins of the pandemic virus, such as particulate matter, cannot be neglected. This work therefore starts from a comparative study of the air environments in China and Italy, the first countries affected by the infection. Currently, a correlation between the spread of infection and pollution is still very controversial. But our paper is not concerned with this. We propose some methodological notes which lead us to the formulation of a general mathematical apparatus (an energy landscape theory), suitable to explain at the molecular level the energetic configurations of the quasi-species of the pandemic virus SARS-CoV-2 in its environment. We focus on complexes between the viral particle and particulate matter in its environment at the scale threshold of the spikes of the viral particle. Then, we wondered if such complexes can lead to the generation of more aggressive viral variants and how to predict their populations and energy configurations, in order to plan an adequate prophylaxis. To do this we trace on the energy landscape some quantities of interest as a function of entropy, glass temperature, energy and the number of configurations, as a function of a control parameter Q that measures the percentage quantity of non-covalent interactions between spikes and nanoparticles which generate stable complexes. In this way we have a quantitative and qualitative measure of the most dangerous quasi-species. A useful quantity in this sense is that of the selective advantage W, which is proportional to the square of the glass temperature, and the average destructive action ℵ(Q) on human lung and upper airway epithelial cells. We can therefore trace a useful methodological protocol, as follows.

After a correct sampling of the organic materials and the environmental background, in the same space–time intervals and in the same physical conditions, the thermodynamics of the system can be correctly defined according to the observations outlined here. Therefore, viral quasi-species can be analyzed as a function of energies and configurations, from which it is possible to trace an experimental landscape for each micro-area of study, which will guide us in the analysis and simulations, as explained in this work.

The first thing to do is to create a diagram of the viral quasi-species identified so far as a function of the average (free) energies measured in the same environmental conditions as the background. Each geographical micro-region must have its own diagram. For example, we will have the Lombardy diagram for a certain time interval, the Wuhan diagram and so on. For each quasi-species, we will have to identify the danger on the basis of the selective advantage W and the average destructive action ℵ(Q) on human lung and upper airway epithelial cells, calculated point-to-point by the tables generated and updated from time to time by the health institutes and scientific committees of each country. The intersection of the diagram with the most dangerous species W and the average destructive action ℵ(Q) on human lung and upper airway epithelial cells, will provide us with some points in the neighborhood of the control parameter Q ≈ 1. These will be our experimental minima of the free energy, which will help us to build our experimental landscape, which will serve as a condition for outline to proceed with the simulations and models of the evolution of the disease. From which it is easy to calculate the various coefficients such as contagion indices and so on.

We reserve further information and updates.