How much energy is stored in SARS‐CoV‐2 and its structural elements?

(R)evolution of Viruses: Introduction to biothermodynamics of viruses

As of 26 April 2024, the International Committee on Taxonomy of Viruses has registered 14690 virus species. Of these, only several dozen have been chemically and thermodynamically characterized. Every virus species is characterized by a specific empirical formula and thermodynamic properties - enthalpy, entropy and Gibbs energy. These physical properties are used in a mechanistic model of virus-host interactions at the cell membrane and in the cytoplasm. This review article presents empirical formulas and Gibbs energies for all major variants of SARS-CoV-2. This article also reports and suggests a mechanistic model of evolutionary changes, with the example of time evolution of SARS-CoV-2 from 2019 to 2024.

Iron‑sulfur clusters in viral proteins: Exploring their elusive nature, roles and new avenues for targeting infections

Viruses have evolved complex mechanisms to exploit host factors for replication and assembly. In response, host cells have developed strategies to block viruses, engaging in a continuous co-evolutionary battle. This dynamic interaction often revolves around the competition for essential resources necessary for both host cell and virus replication. Notably, iron, required for the biosynthesis of several cofactors, including iron‑sulfur (FeS) clusters, represents a critical element in the ongoing competition for resources between infectious agents and host. Although several recent studies have identified FeS cofactors at the core of virus replication machineries, our understanding of their specific roles and the cellular processes responsible for their incorporation into viral proteins remains limited. This review aims to consolidate our current knowledge of viral components that have been characterized as FeS proteins and elucidate how viruses harness these versatile cofactors to their benefit. Its objective is also to propose that viruses may depend on incorporation of FeS cofactors more extensively than is currently known. This has the potential to revolutionize our understanding of viral replication, thereby carrying significant implications for the development of strategies to target infections.

From genotype to phenotype with biothermodynamics: empirical formulas, biosynthesis reactions and thermodynamic properties of preproinsulin, proinsulin and insulin molecules

Insulin was discovered 100 years ago and has been well studied from the perspectives of life and biomedical sciences. This paper reports chemical and biothermodynamic properties of biosynthesis of insulin. This paper reports for the first time the molecular and empirical formulas, biosynthesis reactions, and thermodynamic properties of molecules and their biosynthesis for human preproinsulin, proinsulin, insulin chain A, insulin chain B, insulin, signal peptide and intermediate peptide (C-peptide). Based on these, metabolic reactions were formulated for conversion of preproinsulin to insulin and their thermodynamic feasibility was analyzed.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 strain wars continues: Chemical and thermodynamic characterization of live matter and biosynthesis of Omicron BN.1, CH.1.1 and XBC variants

SARS-CoV-2 has during the last 3 years mutated several dozen times. Most mutations in the newly formed variants have been chemically and thermodynamically characterized. New variants have been declared as variants under monitoring. The European Centre for Disease Prevention and Control has suggested the hypothesis that the new BN.1, CH.1.1 and XBC variants could have properties similar to those of VOC. Thermodynamic properties of new variants have been reported in this manuscript for the first time. Gibbs energy of biosynthesis, as the driving force for viral multiplication, is less negative for the new variants than for the earlier variants. This indicates that the virus has evolved towards decrease in pathogenicity, which leads to less severe forms of COVID-19.

Ghosts of the past: Elemental composition, biosynthesis reactions and thermodynamic properties of Zeta P.2, Eta B.1.525, Theta P.3, Kappa B.1.617.1, Iota B.1.526, Lambda C.37 and Mu B.1.621 variants of SARS-CoV-2

From the perspectives of molecular biology, genetics and biothermodynamics, SARS-CoV-2 is the among the best characterized viruses. Research on SARS-CoV-2 has shed a new light onto driving forces and molecular mechanisms of viral evolution. This paper reports results on empirical formulas, biosynthesis reactions and thermodynamic properties of biosynthesis (multiplication) for the Zeta P.2, Eta B.1.525, Theta P.3, Kappa B.1.617.1, Iota B.1.526, Lambda C.37 and Mu B.1.621 variants of SARS-CoV-2. Thermodynamic analysis has shown that the physical driving forces for evolution of SARS-CoV-2 are Gibbs energy of biosynthesis and Gibbs energy of binding. The driving forces have led SARS-CoV-2 through the evolution process from the original Hu-1 to the newest variants in accordance with the expectations of the evolution theory.

COVID infection in 4 steps: Thermodynamic considerations reveal how viral mucosal diffusion, target receptor affinity and furin cleavage act in concert to drive the nature and degree of infection in human COVID-19 disease

We have developed a mechanistic model of SARS-CoV-2 and SARS-CoV infection, exploring the relationship between the viral diffusion in the mucosa and viral affinity for the angiotensin converting enzyme 2 (ACE2) target. Utilising the structural similarity of SARS-CoV and SARS-CoV-2 and a shared viral target receptor (ACE2), but a dramatic difference in upper or lower respiratory tract infectivity, we were able to generate insights into the linkage of mucosal diffusion and target receptor affinity in determining the pathophysiological pathways of these two viruses. Our analysis reveals that for SARS-CoV-2 the higher affinity of ACE2 binding, the faster and more complete the mucosal diffusion in its transport from the upper airway to the region of the ACE2 target on the epithelium. This diffusional process is essential for the presentation of this virus to the furin catalysed highly efficient entry and infection process in the upper respiratory tract epithelial cells. A failure of SARS-CoV to follow this path is associated with lower respiratory tract infection and decreased infectivity. Thus, our analysis supports the view that through tropism SARS-CoV-2 has evolved a highly efficient membrane entry process that can act in concert with a high binding affinity of this virus and its variants for its ACE2 which in turn promotes enhanced movement of the virus from airway to epithelium. In this way ongoing mutations yielding higher affinities of SARS-CoV-2 for the ACE2 target becomes the basis for higher upper respiratory tract infectivity and greater viral spread. It is concluded that SARS-CoV-2 is constrained in the extent of its activities by the fundamental laws of physics and thermodynamics. Laws that describe diffusion and molecular binding. Moreover it can be speculated that the very earliest contact of this virus with the human mucosa defines the pathogenesis of this infection.

XBB.1.5 Kraken cracked: Gibbs energies of binding and biosynthesis of the XBB.1.5 variant of SARS-CoV-2

The SARS-CoV-2 Hydra with many heads (variants) has been causing the COVID-19 pandemic for 3 years. The appearance of every new head (SARS-CoV-2 variant) causes a new pandemic wave. The last in the series is the XBB.1.5 "Kraken" variant. In the general public (social media) and in the scientific community (scientific journals), during the last several weeks since the variant has appeared, the question was raised of whether the infectivity of the new variant will be greater. This article attempts to provide the answer. Analysis of thermodynamic driving forces of binding and biosynthesis leads to the conclusion that infectivity of the XBB.1.5 variant could be increased to a certain extent. The pathogenicity of the XBB.1.5 variant seems to be unchanged compared to the other Omicron variants.

The SARS-CoV-2 Hydra, a tiny monster from the 21st century: Thermodynamics of the BA.5.2 and BF.7 variants

SARS-CoV-2 resembles the ancient mythical creature Hydra. Just like with the Hydra, when one head is cut, it is followed by appearance of two more heads, suppression of one SARS-CoV-2 variant causes appearance of newer variants. Unlike Hydra that grows identical heads, newer SARS-CoV-2 variants are usually more infective, which can be observed as time evolution of the virus at hand, which occurs through acquisition of mutations during time. The appearance of new variants is followed by appearance of new COVID-19 pandemic waves. With the appearance of new pandemic waves and determining of sequences, in the scientific community and general public the question is always raised of whether the new variant will be more virulent and more pathogenic. The two variants characterized in this paper, BA.5.2 and BF.7, have caused a pandemic wave during the late 2022. This paper gives full chemical and thermodynamic characterization of the BA.5.2 and BF.7 variants of SARS-CoV-2. Having in mind that Gibbs energy of binding and biosynthesis represent the driving forces for the viral life cycle, based on the calculated thermodynamic properties we can conclude that the newer variants are more infective than earlier ones, but that their pathogenicity has not changed.

Never ending story? Evolution of SARS-CoV-2 monitored through Gibbs energies of biosynthesis and antigen-receptor binding of Omicron BQ.1, BQ.1.1, XBB and XBB.1 variants

RNA viruses exhibit a great tendency to mutate. Mutations occur in the parts of the genome that encode the spike glycoprotein and less often in the rest of the genome. This is why Gibbs energy of binding changes more than that of biosynthesis. Starting from 2019, the wild type that was labeled Hu-1 has during the last 3 years evolved to produce several dozen new variants, as a consequence of mutations. Mutations cause changes in empirical formulas of new virus strains, which lead to change in thermodynamic properties of biosynthesis and binding. These changes cause changes in the rate of reactions of binding of virus antigen to the host cell receptor and the rate of virus multiplication in the host cell. Changes in thermodynamic and kinetic parameters lead to changes in biological parameters of infectivity and pathogenicity. Since the beginning of the COVID-19 pandemic, SARS-CoV-2 has been evolving towards increase in infectivity and maintaining constant pathogenicity, or for some variants a slight decrease in pathogenicity. In the case of Omicron BQ.1, BQ.1.1, XBB and XBB.1 variants pathogenicity is identical as in the Omicron BA.2.75 variant. On the other hand, infectivity of the Omicron BQ.1, BQ.1.1, XBB and XBB.1 variants is greater than those of previous variants. This will most likely result in the phenomenon of asymmetric coinfection, that is circulation of several variants in the population, some being dominant.

Why doesn't Ebola virus cause pandemics like SARS-CoV-2?

Ebola virus is among the most dangerous, contagious and deadly etiological causes of viral diseases. However, Ebola virus has never extensively spread in human population and never have led to a pandemic. Why? The mechanistic biophysical model revealing the biothermodynamic background of virus-host interaction) could help us to understand pathogenesis of Ebola virus disease (earlier known as the Ebola hemorrhagic fever). In this paper for the first time the empirical formula, thermodynamic properties of biosynthesis (including the driving force of virus multiplication in the susceptible host), binding constant and thermodynamic properties of binding are reported. Thermodynamic data for Ebola virus were compared with data for SARS-CoV-2 to explain why SARS-CoV-2 has caused a pandemic, while Ebola remains on local epidemic level. The empirical formula of the Ebola virus was found to be CH_1.569O_0.3281N_0.2786P_0.00173S_0.00258. Standard Gibbs energy of biosynthesis of the Ebola virus nucleocapsid is -151.59 kJ/C-mol.

Strain wars 3: Differences in infectivity and pathogenicity between Delta and Omicron strains of SARS-CoV-2 can be explained by thermodynamic and kinetic parameters of binding and growth

In this paper, for the first time, empirical formulas have been reported of the Delta and Omicron strains of SARS-CoV-2. The empirical formula of the Delta strain entire virion was found to be CH1.6383O0.2844N0.2294P0.0064S0.0042, while its nucleocapsid has the formula CH1.5692O0.3431N0.3106P0.0060S0.0043. The empirical formula of the Omicron strain entire virion was found to be CH1.6404O0.2842N0.2299P0.0064S0.0038, while its nucleocapsid has the formula CH1.5734O0.3442N0.3122P0.0060S0.0033. Based on the empirical formulas, standard thermodynamic properties of formation and growth have been calculated and reported for the Delta and Omicron strains. Moreover, standard thermodynamic properties of binding have been reported for Wild type (Hu-1), Alpha, Beta, Gamma, Delta and Omicron strains. For all the strains, binding phenomenological coefficients and antigen-receptor (SGP-ACE2) binding rates have been determined and compared, which are proportional to infectivity. The results show that the binding rate of the Omicron strain is between 1.5 and 2.5 times greater than that of the Delta strain. The Omicron strain is characterized by a greater infectivity, based on the epidemiological data available in the literature. The increased infectivity was explained in this paper using Gibbs energy of binding. However, no indications exist for decreased pathogenicity of the Omicron strain. Pathogenicity is proportional to the virus multiplication rate, while Gibbs energies of multiplication are very similar for the Delta and Omicron strains. Thus, multiplication rate and pathogenicity are similar for the Delta and Omicron strains. The lower number of severe cases caused by the Omicron strain can be explained by increased number of immunized people. Immunization does not influence the possibility of occurrence of infection, but influences the rate of immune response, which is much more efficient in immunized people. This leads to prevention of more severe Omicron infection cases.

Strain wars 5: Gibbs energies of binding of BA.1 through BA.4 variants of SARS-CoV-2

This paper reports, for the first time, standard Gibbs energies of binding of the BA.1, BA.2, BA.3, BA.2.13, BA.2.12.1 and BA.4 Omicron variants of SARS-CoV-2, to the Human ACE2 receptor. Variants BA.1 through BA.3 exhibit a trend of decreasing standard Gibbs energy of binding and hence increased infectivity. The BA.4 variant exhibits a less negative standard Gibbs energy of binding, but also more efficient evasion of the immune response. Therefore, it was concluded that all the analyzed strains evolve in accordance with expectations of the theory of evolution, albeit using different strategies.

Omicron BA.2.75 Subvariant of SARS-CoV-2 Is Expected to Have the Greatest Infectivity Compared with the Competing BA.2 and BA.5, Due to Most Negative Gibbs Energy of Binding

Omicron BA.2.75 may become the next globally dominant strain of COVID-19 in 2022. The BA.2.75 sub-variant has acquired more mutations (9) in spike protein and other genes of SARS-CoV-2 than any other variant. Thus, its chemical composition and thermodynamic properties have changed compared with earlier variants. In this paper, the Gibbs energy of the binding and antigen-receptor binding rate was reported for the BA.2.75 variant. Gibbs energy of the binding of the Omicron BA.2.75 variant is more negative than that of the competing variants BA.2 and BA.5.

Strain Wars: Competitive interactions between SARS-CoV-2 strains are explained by Gibbs energy of antigen-receptor binding

Since the beginning of the COVID-19 pandemic, SARS-CoV-2 has mutated several times into new strains, with an increased infectivity. Infectivity of SARS-CoV-2 strains depends on binding affinity of the virus to its host cell receptor. In this paper, we quantified the binding affinity using Gibbs energy of binding and analyzed the competition between SARS-CoV-2 strains as an interference phenomenon. Gibbs energies of binding were calculated for several SARS-SoV-2 strains, including Hu-1 (wild type), B.1.1.7 (alpha), B.1.351 (beta), P.1 (Gamma), B.1.36 and B.1.617 (Delta). The least negative Gibbs energy of binding is that of Hu-1 strain, -37.97 kJ/mol. On the other hand, the most negative Gibbs energy of binding is that of the Delta strain, -49.50 kJ/mol. We used the more negative Gibbs energy of binding to explain the increased infectivity of newer SARS-CoV-2 strains compared to the wild type. Gibbs energies of binding was found to decrease chronologically, with appearance of new strains. The ratio of Gibbs energies of binding of mutated strains and wild type was used to define a susceptibility coefficient, which is an indicator of viral interference, where a virus can prevent or partially inhibit infection with another virus.

Strain wars 2: Binding constants, enthalpies, entropies, Gibbs energies and rates of binding of SARS-CoV-2 variants

SARS-CoV-2 virus is the cause of COVID-19 pandemic and belongs to RNA viruses, showing great tendency to mutate. Several dozens of mutations have been observed on the SARS-CoV-2 virus, during the last two years. Some of the mutated strains show a greater infectivity and are capable of suppressing the earlier strains, through interference. In this work, kinetic and thermodynamic properties were calculated for strains characterized by various numbers and locations of mutations. It was shown that mutations lead to changes in chemical composition, thermodynamic properties and infectivity. Through competition, the phenomenon of interference of various SARS-CoV-2 strains was explained, which results in suppression of the wild type by mutant strains. Standard Gibbs energy of binding and binding constant for the Omicron (B.1.1.529) strain were found to be Δ_BG⁰ = −45.96 kJ/mol and K_B = 1.13 ∙ 10⁺⁸ M⁻¹, respectively.

Atom counting method for determining elemental composition of viruses and its applications in biothermodynamics and environmental science

Quantitative physicochemical perspective on life processes has been a great asset, in bioengineering and biotechnology. The quantitative physicochemical approach can be applied to practically all organisms, including viruses, if their chemical composition and thermodynamic properties are known. In this paper, a new method is suggested for determining elemental composition of viruses, based on atom counting. The atom counting method requires knowledge of genetic sequence, protein sequences and protein copy numbers. An algorithm was suggested for a program that finds elemental composition of various viruses (DNA or RNA, enveloped or non-enveloped). Except for the nucleic acid, capsid proteins, lipid bilayer and carbohydrates, this method includes membrane proteins, as well as spike proteins. The atom counting method has been compared with the existing molecular composition and geometric methods on 5 viruses of different morphology, as well as experimentally determined composition of the poliovirus. The atom counting method was found to be more accurate in most cases. The three methods were found to be complementary, since they require different kind of input information. Moreover, since the 3 methods rest on different assumptions, results of one model can be compared to those of the other two.