Mutagenic Distinction between the Receptor-Binding and Fusion Subunits of the SARS-CoV-2 Spike Glycoprotein and Its Upshot.

We observe that a residue R of the spike glycoprotein of SARS-CoV-2 that has mutated in one or more of the current variants of concern or interest, or under monitoring, rarely participates in a backbone hydrogen bond if R lies in the S1 subunit and usually participates in one if R lies in the S2 subunit. A partial explanation for this based upon free energy is explored as a potentially general principle in the mutagenesis of viral glycoproteins. This observation could help target future vaccine cargos for the evolving coronavirus as well as more generally.

Conserved High Free Energy Sites in Human Coronavirus Spike Glycoprotein Backbones.

Methods previously developed by the author are applied to uncover several sites of interest in the spike glycoproteins of all known human coronaviruses (hCoVs), including SARS-CoV-2 that causes COVID-19. The sites comprise three-dimensional neighborhoods of peptides characterized by four key properties: (1) they pinpoint regions of high free energy in the backbone whose obstruction might interrupt function; (2) by their very definition, they occur rarely in the universe of all gene-encoded proteins that could obviate host response to compounds designed for their interference; (3) they are common to all known hCoV spikes, possibly retaining activity in light of inevitable viral mutation; and (4) they are exposed in the molecular surface of the glycoprotein. These peptides in SARS-CoV-2 are given by the triples of residues (131, 117, 134), (203, 227, 228), and (1058, 730, 731) in its spike.

Backbone Free Energy Estimator Applied to Viral Glycoproteins.

Earlier analysis of the Protein Data Bank derived the distribution of rotations from the plane of a protein hydrogen bond donor peptide group to the plane of its acceptor peptide group. The quasi Boltzmann formalism of Pohl-Finkelstein is employed to estimate free energies of protein elements with these hydrogen bonds, pinpointing residues with a high propensity for conformational change. This is applied to viral glycoproteins as well as capsids, where the 90th+ percentiles of free energies determine residues that correlate well with viral fusion peptides and other functional domains in known cases and thus provide a novel method for predicting these sites of importance as antiviral drug or vaccine targets in general.

Model of Morphogenesis.

We build a theoretical model of morphogenesis. This model describes cell fate in the developing organism using the notion of epigenetic code of each cell. Namely, given the epigenetic spectra of a cell and its neighboring cells, we can determine the corresponding cell event it will perform.

Theory of Morphogenesis.

A model of morphogenesis is proposed based on seven explicit postulates. The mathematical import and biological significance of the postulates are explored and discussed.