Analyzing host-viral interactome of SARS-CoV-2 for identifying vulnerable host proteins during COVID-19 pathogenesis

CAPZB mRNA is a novel biomarker for cervical high-grade squamous lesions

This study aimed to evaluate the potential of capping protein (actin filament) muscle Z-line subunit β (CAPZB) messenger ribonucleic acid (mRNA) levels as a biomarker for distinguishing low-grade squamous intraepithelial lesions of the cervix (LSIL) from high-grade squamous intraepithelial lesions of the cervix (HSIL). We collected a total of 166 cervical exfoliated cells and divided them into five groups based on histopathological results. Each sample was divided into two portions, one for fluorescence in situ hybridization (FISH) detection and the other for bisulfite sequencing polymerase chain reaction (BSP) detection. We found that FISH detection of CAPZB mRNA mean fluorescence intensity (MFI) and BSP detection of CAPZB deoxyribonucleic acid (DNA) percentage of methylation rate (PMR) performed as biomarkers for distinguishing HSIL from LSIL, with an area under the receiver operating characteristic curve (AUC), sensitivity, specificity and cut-off value of 0.893, 81.25%, 80.39% and 0.616, 0.794, 64.06%, 81.37% and 0.454, respectively. Furthermore, FISH detection of CAPZB mRNA exhibited a greater AUC (0.893) for the detection of HSIL than the CAPZB DNA methylation method (0.794), indicating the CAPZB mRNA levels can be used as a biomarker for assessing cervical lesions.

PyPCN: protein contact networks in PyMOL

MOTIVATION

Protein contact networks (PCNs) represent the 3D structure of a protein using network formalism. Inter-residue contacts are described as binary adjacency matrices, which are derived from the graph representation of residues (as α-carbons, β-carbons or centroids) and Euclidean distances according to defined thresholds. Functional characterization algorithms are computed on binary adjacency matrices to unveil allosteric, dynamic, and interaction mechanisms in proteins. Such strategies are usually applied in a combinatorial manner, although rarely in seamless and user-friendly implementations.

RESULTS

PyPCN is a plugin for PyMOL wrapping more than twenty PCN algorithms and metrics in an easy-to-use graphical user interface, to support PCN analysis. The plugin accepts 3D structures from the Protein Data Bank, user-provided PDBs, or precomputed adjacency matrices. The results are directly mapped to 3D protein structures and organized into interactive diagrams for their visualization. A dedicated graphical user interface combined with PyMOL visual support makes analysis more intuitive and easier, extending the applicability of PCNs.

AVAILABILITY AND IMPLEMENTATION

https://github.com/pcnproject/PyPCN.

Strategies and Trends in COVID-19 Vaccination Delivery: What We Learn and What We May Use for the Future

Vaccination has been the most effective way to control the outbreak of the COVID-19 pandemic. The numbers and types of vaccines have reached considerable proportions, even if the question of vaccine procedures and frequency still needs to be resolved. We have come to learn the necessity of defining vaccination distribution strategies with regard to COVID-19 that could be used for any future pandemics of similar gravity. In fact, vaccine monitoring implies the existence of a strategy that should be measurable in terms of input and output, based on a mathematical model, including death rates, the spread of infections, symptoms, hospitalization, and so on. This paper addresses the issue of vaccine diffusion and strategies for monitoring the pandemic. It provides a description of the importance and take up of vaccines and the links between procedures and the containment of COVID-19 variants, as well as the long-term effects. Finally, the paper focuses on the global scenario in a world undergoing profound social and political change, with particular attention on current and future health provision. This contribution would represent an example of vaccination experiences, which can be useful in other pandemic or epidemiological contexts.

SARS-COV-2 Coronavirus Papain-like Protease PLpro as an Antiviral Target for Inhibitors of Active Site and Protein-Protein Interactions

The papain-like protease PLpro of the SARS-CoV-2 coronavirus is a multifunctional enzyme that catalyzes the proteolytic processing of two viral polyproteins, pp1a and pp1ab. PLpro also cleaves peptide bonds between host cell proteins and ubiquitin (or ubiquitin-like proteins), which is associated with a violation of immune processes. Nine structures of the most effective inhibitors of the PLpro active center were prioritized according to the parameters of biochemical (IC 50) and cellular tests to assess the suppression of viral replication (EC 50) and cytotoxicity (CC 50). A literature search has shown that PLpro can interact with at least 60 potential protein partners in cells, 23 of which are targets for other viral proteins (human papillomavirus and Epstein-Barr virus). The analysis of protein-protein interactions showed that the proteins USP3, UBE2J1, RCHY1, and FAF2 involved in deubiquitinylation and ubiquitinylation processes contain the largest number of bonds with other proteins; the interaction of viral proteins with them can affect the architecture of the entire network of protein-protein interactions. Using the example of a spatial model of the PLpro/ubiquitin complex and a set of 154 naturally occurring compounds with known antiviral activity, 13 compounds (molecular masses in the range of 454-954 Da) were predicted as potential PLpro inhibitors. These compounds bind to the "hot" amino acid residues of the protease at the positions Gly163, Asp164, Arg166, Glu167, and Tyr264 involved in the interaction with ubiquitin. Thus, pharmacological effects on peripheral PLpro sites, which play important roles in binding protein substrates, may be an additional target-oriented antiviral strategy.

Computational analysis of the sequence-structure relation in SARS-CoV-2 spike protein using protein contact networks

The structure of proteins impacts directly on the function they perform. Mutations in the primary sequence can provoke structural changes with consequent modification of functional properties. SARS-CoV-2 proteins have been extensively studied during the pandemic. This wide dataset, related to sequence and structure, has enabled joint sequence-structure analysis. In this work, we focus on the SARS-CoV-2 S (Spike) protein and the relations between sequence mutations and structure variations, in order to shed light on the structural changes stemming from the position of mutated amino acid residues in three different SARS-CoV-2 strains. We propose the use of protein contact network (PCN) formalism to: (i) obtain a global metric space and compare various molecular entities, (ii) give a structural explanation of the observed phenotype, and (iii) provide context dependent descriptors of single mutations. PCNs have been used to compare sequence and structure of the Alpha, Delta, and Omicron SARS-CoV-2 variants, and we found that omicron has a unique mutational pattern leading to different structural consequences from mutations of other strains. The non-random distribution of changes in network centrality along the chain has allowed to shed light on the structural (and functional) consequences of mutations.

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.

Structural analysis of SARS-CoV-2 Spike protein variants through graph embedding

Since December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected almost all countries. The unprecedented spreading of this virus has led to the insurgence of many variants that impact protein sequence and structure that need continuous monitoring and analysis of the sequences to understand the genetic evolution and to prevent possible dangerous outcomes. Some variants causing the modification of the structure of the proteins, such as the Spike protein S, need to be monitored. Protein contact networks (PCNs) have been recently proposed as a modelling framework for protein structures. In such a framework, the protein structure is represented as an unweighted graph whose nodes are the central atoms of the backbones (C- α ), and edges connect two atoms falling in the spatial distance between 4 and 7 Å. PCN may also be a data-rich representation since we may add to each node/atom biological and topological information. Such formalism enables the possibility of using algorithms from graph theory to analyze the graph. In particular, we refer to graph embedding methods enabling the analysis of such graphs with deep learning methods. In this work, we explore the possibility of embedding PCN using Graph Neural Networks and then analyze in the embedded space each residue to distinguish mutated residues from non-mutated ones. In particular, we analyzed the structure of the Spike protein of the coronavirus. First, we obtained the PCNs of the Spike protein for the wild-type, α , β , and δ variants. Then we used the GraphSage embedding algorithm to obtain an unsupervised embedding. Then we analyzed the point of mutation in the embedded space. Results show the characteristics of the mutation point in the embedding space.

Analysis and identification of potential type II helper T cell (Th2)-Related key genes and therapeutic agents for COVID-19

COVID-19 pandemic poses a severe threat to public health. However, so far, there are no effective drugs for COVID-19. Transcriptomic changes and key genes related to Th2 cells in COVID-19 have not been reported. These genes play an important role in host interactions with SARS-COV-2 and may be used as promising target. We analyzed five COVID-19-associated GEO datasets (GSE157103, GSE152641, GSE171110, GSE152418, and GSE179627) using the xCell algorithm and weighted gene co-expression network analysis (WGCNA). Results showed that 5 closely correlated modular genes to COVID-19 and Th2 cell enrichment levels, including purple, blue, pink, tan and turquoise, were intersected with differentially expressed genes (DEGs) and 648 shared genes were obtained. GO and KEGG pathway enrichment analyses revealed that they were enriched in cell proliferation, differentiation, and immune responses after virus infection. The most significantly enriched pathway involved the regulation of viral life cycle. Three key genes, namely CCNB1, BUB1, and UBE2C, may clarify the pathogenesis of COVID-19 associated with Th2 cells. 11 drug candidates were identified that could down-regulate three key genes using the cMAP database and demonstrated strong drugs binding energies aganist the three keygenes using molecular docking methods. BUB1, CCNB1 and UBE2C were identified key genes for COVID-19 and could be promising therapeutic targets.

A systems biology approach for investigating significantly expressed genes among COVID-19, hepatocellular carcinoma, and chronic hepatitis B

Background

Worldwide, COVID-19's death rate is about 2%, considering the incidence and mortality. However, the information on its complications in other organs, specifically the liver and its disorders, is limited in mild or severe cases. In this study, we aimed to computationally investigate the typical relationships between liver-related diseases [i.e., hepatocellular carcinoma (HCC), and chronic hepatitis B (CHB)] and COVID-19, considering the involved significant genes and their molecular mechanisms.

Methods

We investigated two GEO microarray datasets (GSE164805 and GSE58208) to identify differentially expressed genes (DEGs) among the generated four datasets for mild/severe COVID-19, HCC, and CHB. Then, the overlapping genes among them were identified for GO and KEGG enrichment analyses, protein-protein interaction network construction, hub genes determination, and their associations with immune cell infiltration.

Results

A total of 22 significant genes (i.e., ACTB, ATM, CDC42, DHX15, EPRS, GAPDH, HIF1A, HNRNPA1, HRAS, HSP90AB1, HSPA8, IL1B, JUN, POLR2B, PTPRC, RPS27A, SFRS1, SMARCA4, SRC, TNF, UBE2I, and VEGFA) were found to play essential roles among mild/severe COVID-19 associated with HCC and CHB. Moreover, the analysis of immune cell infiltration revealed that these genes are mostly positively correlated with tumor immune and inflammatory responses.

Conclusions

In summary, the current study demonstrated that 22 identified DEGs might play an essential role in understanding the associations between the mild/severe COVID-19 patients with HCC and CHB. So, the HCC and CHB patients involved in different types of COVID-19 can benefit from immune-based targets for therapeutic interventions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43042-022-00360-3.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Network for network concept offers new insights into host- SARS-CoV-2 protein interactions and potential novel targets for developing antiviral drugs

SARS-CoV-2, the causal agent of COVID-19, is primarily a pulmonary virus that can directly or indirectly infect several organs. Despite many studies carried out during the current COVID-19 pandemic, some pathological features of SARS-CoV-2 have remained unclear. It has been recently attempted to address the current knowledge gaps on the viral pathogenicity and pathological mechanisms via cellular-level tropism of SARS-CoV-2 using human proteomics, visualization of virus-host protein-protein interactions (PPIs), and enrichment analysis of experimental results. The synergistic use of models and methods that rely on graph theory has enabled the visualization and analysis of the molecular context of virus/host PPIs. We review current knowledge on the SARS-COV-2/host interactome cascade involved in the viral pathogenicity through the graph theory concept and highlight the hub proteins in the intra-viral network that create a subnet with a small number of host central proteins, leading to cell disintegration and infectivity. Then we discuss the putative principle of the "gene-for-gene and "network for network" concepts as platforms for future directions toward designing efficient anti-viral therapies.

Disease spreading modeling and analysis: a survey

MOTIVATION

The control of the diffusion of diseases is a critical subject of a broad research area, which involves both clinical and political aspects. It makes wide use of computational tools, such as ordinary differential equations, stochastic simulation frameworks and graph theory, and interaction data, from molecular to social granularity levels, to model the ways diseases arise and spread. The coronavirus disease 2019 (COVID-19) is a perfect testbench example to show how these models may help avoid severe lockdown by suggesting, for instance, the best strategies of vaccine prioritization.

RESULTS

Here, we focus on and discuss some graph-based epidemiological models and show how their use may significantly improve the disease spreading control. We offer some examples related to the recent COVID-19 pandemic and discuss how to generalize them to other diseases.

Beyond COVID-19 Pandemic: Topology-aware optimization of vaccination strategy for minimizing virus spreading

•

Paper discusses the relevance of the adoption of ad-hoc vaccination strategies.

•

Paper shows how to evaluate the impact of different vaccination strategy by considering network-based models.

•

Tailored interventions, e.g., vaccination, applied on central nodes of these networks may efficiently stop the propagation of an infection.

•

The way node "centrality" is defined is the key to curb infection spreading.

The mitigation of an infectious disease spreading has recently gained considerable attention from the research community. It may be obtained by adopting sanitary measurements (e.g., vaccination, wearing masks), social rules (e.g., social distancing), together with an extensive vaccination campaign. Vaccination is currently the primary way for mitigating the Coronavirus Disease (COVID-19) outbreak without severe lockdown. Its effectiveness also depends on the number and timeliness of administrations and thus demands strict prioritization criteria. Almost all countries have prioritized similar classes of exposed workers: healthcare professionals and the elderly, obtaining to maximize the survival of patients and years of life saved. Nevertheless, the virus is currently spreading at high rates, and any prioritization criterion so far adopted did not account for the structural organization of the contact networks.

We reckon that a network where nodes are people while the edges represent their social contacts may efficiently model the virus’s spreading. It is known that tailored interventions (e.g., vaccination) on central nodes may efficiently stop the propagation, thereby eliminating the “bridge edges.” We then introduce such a model and consider both synthetic and real datasets. We present the benefits of a topology-aware versus an age-based vaccination strategy to mitigate the spreading of the virus. The code is available at https://github.com/mazzalab/playgrounds.

Identifying Potential Mitochondrial Proteome Signatures Associated with the Pathogenesis of Pulmonary Arterial Hypertension in the Rat Model

Pulmonary arterial hypertension (PAH) is a severe and progressive disease that affects the heart and lungs and a global health concern that impacts individuals and society. Studies have reported that some proteins related to mitochondrial metabolic functions could play an essential role in the pathogenesis of PAH, and their specific expression and biological function are still unclear. We successfully constructed a monocrotaline- (MCT-) induced PAH rat model in the present research. Then, the label-free quantification proteomic technique was used to determine mitochondrial proteins between the PAH group (n = 6) and the normal group (n = 6). Besides, we identified 1346 mitochondrial differentially expressed proteins (DEPs) between these two groups. Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to analyze the mainly mitochondrial DEPs' biological functions and the signal pathways. Based on the protein-protein interaction (PPI) network construction and functional enrichment, we screened 19 upregulated mitochondrial genes (Psmd1, Psmc4, Psmd13, Psmc2, etc.) and 123 downregulated mitochondrial genes (Uqcrfs1, Uqcrc1, Atp5c1, Atp5a1, Uqcrc2, etc.) in rats with PAH. Furthermore, in an independent cohort dataset and experiments with rat lung tissue using qPCR, validation results consistently showed that 6 upregulated mitochondrial genes (Psmd2, Psmc4, Psmc3, Psmc5, Psmd13, and Psmc2) and 3 downregulated mitochondrial genes (Lipe, Cat, and Prkce) were significantly differentially expressed in the lung tissue of PAH rats. Using the RNAInter database, we predict potential miRNA target hub mitochondrial genes at the transcriptome level. We also identified bortezomib and carfilzomib as the potential drugs for treatment in PAH. Finally, this study provides us with a new perspective on critical biomarkers and treatment strategies in PAH.

Differential Co-Expression Network Analysis Reveals Key Hub-High Traffic Genes as Potential Therapeutic Targets for COVID-19 Pandemic

Background

The recent emergence of COVID-19, rapid worldwide spread, and incomplete knowledge of molecular mechanisms underlying SARS-CoV-2 infection have limited development of therapeutic strategies. Our objective was to systematically investigate molecular regulatory mechanisms of COVID-19, using a combination of high throughput RNA-sequencing-based transcriptomics and systems biology approaches.

Methods

RNA-Seq data from peripheral blood mononuclear cells (PBMCs) of healthy persons, mild and severe 17 COVID-19 patients were analyzed to generate a gene expression matrix. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules in healthy samples as a reference set. For differential co-expression network analysis, module preservation and module-trait relationships approaches were used to identify key modules. Then, protein-protein interaction (PPI) networks, based on co-expressed hub genes, were constructed to identify hub genes/TFs with the highest information transfer (hub-high traffic genes) within candidate modules.

Results

Based on differential co-expression network analysis, connectivity patterns and network density, 72% (15 of 21) of modules identified in healthy samples were altered by SARS-CoV-2 infection. Therefore, SARS-CoV-2 caused systemic perturbations in host biological gene networks. In functional enrichment analysis, among 15 non-preserved modules and two significant highly-correlated modules (identified by MTRs), 9 modules were directly related to the host immune response and COVID-19 immunopathogenesis. Intriguingly, systemic investigation of SARS-CoV-2 infection identified signaling pathways and key genes/proteins associated with COVID-19's main hallmarks, e.g., cytokine storm, respiratory distress syndrome (ARDS), acute lung injury (ALI), lymphopenia, coagulation disorders, thrombosis, and pregnancy complications, as well as comorbidities associated with COVID-19, e.g., asthma, diabetic complications, cardiovascular diseases (CVDs), liver disorders and acute kidney injury (AKI). Topological analysis with betweenness centrality (BC) identified 290 hub-high traffic genes, central in both co-expression and PPI networks. We also identified several transcriptional regulatory factors, including NFKB1, HIF1A, AHR, and TP53, with important immunoregulatory roles in SARS-CoV-2 infection. Moreover, several hub-high traffic genes, including IL6, IL1B, IL10, TNF, SOCS1, SOCS3, ICAM1, PTEN, RHOA, GDI2, SUMO1, CASP1, IRAK3, HSPA5, ADRB2, PRF1, GZMB, OASL, CCL5, HSP90AA1, HSPD1, IFNG, MAPK1, RAB5A, and TNFRSF1A had the highest rates of information transfer in 9 candidate modules and central roles in COVID-19 immunopathogenesis.

Conclusion

This study provides comprehensive information on molecular mechanisms of SARS-CoV-2-host interactions and identifies several hub-high traffic genes as promising therapeutic targets for the COVID-19 pandemic.

Mutation-Induced Long-Range Allosteric Interactions in the Spike Protein Determine the Infectivity of SARS-CoV-2 Emerging Variants

The emergence of a variety of highly transmissible SARS-CoV-2 variants, the causative agent of COVID-19, with multiple spike mutations poses serious challenges in overcoming the ongoing deadly pandemic. It is, therefore, essential to understand how these variants gain enhanced ability to evade immune responses with a higher rate of spreading infection. To address this question, here we have individually assessed the effects of SARS-CoV-2 variant-specific spike (S) protein receptor-binding domain (RBD) mutations E484K, K417N, L452Q, L452R, N501Y, and T478K that characterize and differentiate several emerging variants. Despite the hundreds of apparently neutral mutations observed in the domains other than the RBD, we have shown that each RBD mutation site is differentially engaged in an interdomain allosteric network involving mutation sites from a distant domain, affecting interactions with the human receptor angiotensin-converting enzyme-2 (ACE2). This allosteric network couples the residues of the N-terminal domain (NTD) and the RBD, which are modulated by the RBD-specific mutations and are capable of propagating mutation-induced perturbations between these domains through a combination of structural changes and effector-dependent modulations of dynamics. One key feature of this network is the inclusion of compensatory mutations segregated into three characteristically different clusters, where each cluster residue site is allosterically coupled with specific RBD mutation sites. Notably, each RBD mutation acted like a positive allosteric modulator; nevertheless, K417N was shown to have the largest effects among all of the mutations on the allostery and thereby holds the highest binding affinity with ACE2. This result will be useful for designing the targeted control measure and therapeutic efforts aiming at allosteric modulators.

Exploiting the molecular basis of age and gender differences in outcomes of SARS-CoV-2 infections

Motivation: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (coronavirus disease, 2019; COVID-19) is associated with adverse outcomes in patients. It has been observed that lethality seems to be related to the age of patients. While ageing has been extensively demonstrated to be accompanied by some modifications at the gene expression level, a possible link with COVID-19 manifestation still need to be investigated at the molecular level. Objectives: This study aims to shed out light on a possible link between the increased COVID-19 lethality and the molecular changes that occur in elderly people. Methods: We considered public datasets of ageing-related genes and their expression at the tissue level. We selected human proteins interacting with viral ones that are known to be related to the ageing process. Finally, we investigated changes in the expression level of coding genes at the tissue, gender and age level. Results: We observed a significant intersection between some SARS-CoV-2 interactors and ageing-related genes, suggesting that those genes are particularly affected by COVID-19 infection. Our analysis evidenced that virus infection particularly involves ageing molecular mechanisms centred around proteins EEF2, NPM1, HMGA1, HMGA2, APEX1, CHEK1, PRKDC, and GPX4. We found that HMGA1 and NPM1 have different expressions in the lung of males, while HMGA1, APEX1, CHEK1, EEF2, and NPM1 present changes in expression in males due to ageing effects. Conclusion: Our study generated a mechanistic framework to clarify the correlation between COVID-19 incidence in elderly patients and molecular mechanisms of ageing. We also provide testable hypotheses for future investigation and pharmacological solutions tailored to specific age ranges.

Comparative analysis of human coronaviruses focusing on nucleotide variability and synonymous codon usage patterns

The prevailing COVID-19 pandemic has drawn the attention of the scientific community to study the evolutionary origin of Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2). This study is a comprehensive quantitative analysis of the protein-coding sequences of seven human coronaviruses (HCoVs) to decipher the nucleotide sequence variability and codon usage patterns. It is essential to understand the survival ability of the viruses, their adaptation to hosts, and their evolution. The current analysis revealed a high abundance of the relative dinucleotide (odds ratio), GC and CT pairs in the first and last two codon positions, respectively, as well as a low abundance of the CG pair in the last two positions of the codon, which might be related to the evolution of the viruses. A remarkable level of variability of GC content in the third position of the codon among the seven coronaviruses was observed. Codons with high RSCU values are primarily from the aliphatic and hydroxyl amino acid groups, and codons with low RSCU values belong to the aliphatic, cyclic, positively charged, and sulfur-containing amino acid groups. In order to elucidate the evolutionary processes of the seven coronaviruses, a phylogenetic tree (dendrogram) was constructed based on the RSCU scores of the codons. The severe and mild categories CoVs were positioned in different clades. A comparative phylogenetic study with other coronaviruses depicted that SARS-CoV-2 is close to the CoV isolated from pangolins (Manis javanica, Pangolin-CoV) and cats (Felis catus, SARS(r)-CoV). Further analysis of the effective number of codon (ENC) usage bias showed a relatively higher bias for SARS-CoV and MERS-CoV compared to SARS-CoV-2. The ENC plot against GC3 suggested that the mutational bias might have a role in determining the codon usage variation among candidate viruses. A codon adaptability study on a few human host parasites (from different kingdoms), including CoVs, showed a diverse adaptability pattern. SARS-CoV-2 and SARS-CoV exhibit relatively lower but similar codon adaptability compared to MERS-CoV.