A neutral evolution test derived from a theoretical amino acid substitution model

The Spike Protein of SARS-CoV-2 Is Adapting Because of Selective Pressures

The global scale of the COVID-19 pandemic has demonstrated the evolution of SARS-CoV-2 and the clues of adaptation. After two years and two months since the declaration of the pandemic, several variants have emerged and become fixed in the human population thanks to extrinsic selective pressures but also to the inherent mutational capacity of the virus. Here, we applied a neutral substitution evolution test to the spike (S) protein of Omicron’s protein and compared it to the others’ variant of concern (VOC) neutral evolution. We carried out comparisons among the interactions between the S proteins from the VOCs (Alpha, Beta, Gamma, Delta and Omicron) and the receptor ACE2. The shared amino acids among all the ACE2 binding S proteins remain constant, indicating that these amino acids are essential for the accurate binding to the receptor. The complexes of the RBD for every variant with the receptor were used to identify the amino acids involved in the protein—protein interaction (PPI). The RBD of Omicron establishes 82 contacts, compared to the 74 of the Wuhan original viral protein. Hence, the mean number of contacts per residue is higher, making the contact thermodynamically more stable. The RBDs of the VOCs are similar in sequence and structure; however, Omicron’s RBD presents the largest deviation from the structure by 1.11 Å RMSD, caused by a set of mutations near the glycosylation N343. The chemical properties and structure near the glycosylation N343 of the Omicron S protein are different from the original protein, which provoke reduced recognition by the neutralizing antibodies. Our results hint that selective pressures are induced by mass vaccination throughout the world and by the persistence of recurrent infections in immunosuppressed individuals, who did not eliminate the infection and ended up facilitating the selection of viruses whose characteristics are different from the previous VOCs, less pathogenic but with higher transmissibility.

Neutral evolution test of the spike protein of SARS-CoV-2 and its implications in the binding to ACE2

As the SARS-CoV-2 has spread and the pandemic has dragged on, the virus continued to evolve rapidly resulting in the emergence of new highly transmissible variants that can be of public health concern. The evolutionary mechanisms that drove this rapid diversity are not well understood but neutral evolution should open the first insight. The neutral theory of evolution states that most mutations in the nucleic acid sequences are random and they can be fixed or disappear by purifying selection. Herein, we performed a neutrality test to better understand the selective pressures exerted over SARS-CoV-2 spike protein from homologue proteins of Betacoronavirus, as well as to the spikes from human clinical isolates of the virus. Specifically, Tyr and Asn have higher occurrence rates on the Receptor Binding Domain (RBD) and in the overall sequence of spike proteins of Betacoronavirus, whereas His and Arg have lower occurrence rates. The in vivo evolutionary phenomenon of SARS-CoV-2 shows that Glu, Lys, Phe, and Val have the highest probability of occurrence in the emergent viral particles. Amino acids that have higher occurrence than the expected by the neutral control, are favorable and are fixed in the sequence while the ones that have lower occurrence than expected, influence the stability and/or functionality of the protein. Our results show that most unique mutations either for SARS-CoV-2 or its variants of health concern are under selective pressures, which could be related either to the evasion of the immune system, increasing the virus' fitness or altering protein - protein interactions with host proteins. We explored the consequences of those selected mutations in the structure and protein - protein interaction with the receptor. Altogether all these forces have shaped the spike protein and the continually evolving variants.

On the Importance of Asymmetry in the Phenotypic Expression of the Genetic Code upon the Molecular Evolution of Proteins

The standard genetic code (SGC) is a mapping between the 64 possible arrangements of the four RNA nucleotides (C, A, U, G) into triplets or codons, where 61 codons are assigned to a specific amino acid and the other three are stop codons for terminating protein synthesis. Aminoacyl-tRNA synthetases (aaRSs) are responsible for implementing the SGC by specifically amino-acylating only its cognate transfer RNA (tRNA), thereby linking an amino acid with its corresponding anticodon triplets. tRNAs molecules bind each codon with its anticodon. To understand the meaning of symmetrical/asymmetrical properties of the SGC, we designed synthetic genetic codes with known symmetries and with the same degeneracy of the SGC. We determined their impact on the substitution rates for each amino acid under a neutral model of protein evolution. We prove that the phenotypic graphs of the SGC for codons and anticodons for all the possible arrangements of nucleotides are asymmetric and the amino acids do not form orbits. In the symmetrical synthetic codes, the amino acids are grouped according to their codonicity, this is the number of triplets encoding a given amino acid. Both the SGC and symmetrical synthetic codes exhibit a probability of occurrence of the amino acids proportional to their degeneracy. Unlike the SGC, the synthetic codes display a constant probability of occurrence of the amino acid according to their codonicity. The asymmetry of the phenotypic graphs of codons and anticodons of the SGC, has important implications on the evolutionary processes of proteins.