Understanding COVID-19 via comparative analysis of dark proteomes of SARS-CoV-2, human SARS and bat SARS-like coronaviruses

The recently emerged coronavirus designated as SARS-CoV-2 (also known as 2019 novel coronavirus (2019-nCoV) or Wuhan coronavirus) is a causative agent of coronavirus disease 2019 (COVID-19), which is rapidly spreading throughout the world now. More than 1.21 million cases of SARS-CoV-2 infection and more than 67,000 COVID-19-associated mortalities have been reported worldwide till the writing of this article, and these numbers are increasing every passing hour. The World Health Organization (WHO) has declared the SARS-CoV-2 spread as a global public health emergency and admitted COVID-19 as a pandemic now. Multiple sequence alignment data correlated with the already published reports on SARS-CoV-2 evolution indicated that this virus is closely related to the bat severe acute respiratory syndrome-like coronavirus (bat SARS-like CoV) and the well-studied human SARS coronavirus (SARS-CoV). The disordered regions in viral proteins are associated with the viral infectivity and pathogenicity. Therefore, in this study, we have exploited a set of complementary computational approaches to examine the dark proteomes of SARS-CoV-2, bat SARS-like, and human SARS CoVs by analysing the prevalence of intrinsic disorder in their proteins. According to our findings, SARS-CoV-2 proteome contains very significant levels of structural order. In fact, except for nucleocapsid, Nsp8, and ORF6, the vast majority of SARS-CoV-2 proteins are mostly ordered proteins containing less intrinsically disordered protein regions (IDPRs). However, IDPRs found in SARS-CoV-2 proteins are functionally important. For example, cleavage sites in its replicase 1ab polyprotein are found to be highly disordered, and almost all SARS-CoV-2 proteins contains molecular recognition features (MoRFs), which are intrinsic disorder-based protein–protein interaction sites that are commonly utilized by proteins for interaction with specific partners. The results of our extensive investigation of the dark side of SARS-CoV-2 proteome will have important implications in understanding the structural and non-structural biology of SARS or SARS-like coronaviruses. Electronic supplementary material The online version of this article (10.1007/s00018-020-03603-x) contains supplementary material, which is available to authorized users.


Introduction
The emerging coronavirus disease 2019 (COVID-19) is a recent pandemic which has been recently declared as a public health emergency by the World Health Organization (WHO). Since its first appearance in visitors of the Wuhan's seafood and meat market, China, reported in December 2019, COVID-19 has now a large-scale socioeconomic impact [1]. According to WHO, till 6 April 2020, the infection has spread over to at least 170 countries and territories, where there have been more than 1.21 million confirmed cases, with more than 67,000 deaths due to COVID-19. One should also keep in mind that these data on the COVID-19 spread and related casualties are rapidly becoming outdated, almost with the speed of typing of these sentences [2].
According to the International Committee on Taxonomy of Viruses (ICTV), SARS-CoV-2 comes under the coronavirinae sub-family of coronaviridae family of order nidovirales. Viruses of the nidovirales order are enveloped, nonsegmented positive-sense, single-stranded RNA viruses [3]. The family coronaviridae comprises of vertebrate infecting viruses that transmit horizontally, mainly through the oral/ fecal route and cause gastrointestinal and respiratory problems to the host [4]. Sub-family coronavirinae consists of four genera, namely: alpha, beta, gamma, and delta coronaviruses based on the phylogenetic clustering of viruses [5,6]. Coronavirinae having the largest genomes among the RNA viruses incorporate their ~ 30 kb genomes inside an enveloped capsid [7].
SARS-coronavirus genomic RNA includes a 5′ cap, leader sequence, UTR, a replicase gene, genes for structural and accessory proteins, 3′ UTR, and a poly-A tail (Fig. 1). Two-third of the genome codes for the replicase polyproteins (~ 20 kb) containing all non-structural viral proteins, while the remaining part of the genome (~ 10 kb) contains genes for accessory proteins interspersed between the genes responsible for coding structural proteins [7,8]. The ~ 20 kb (replicase gene) ssRNA is translated first into two long polyproteins: replicase polyprotein 1a and 1ab inside host cells. The newly formed polyproteins, after cleavage by two viral proteases, result in 16 non-structural proteins (Nsps) that perform a wide range of functions for viruses inside the host cell [9,10]. The genomic sequence of SARS-CoV-2 is reported to have 29,903 nucleotides with GenBank accession number NC_045512 [11].
In this study, we analysed the dark side of SARS-CoV-2 proteome (i.e. a part of a proteome that includes proteins or protein regions, which are not amenable to experimental structure determination by existing means and inaccessible to homology modeling), to better understand an interplay between the ordered and disordered components of the proteome. According to the "heretic" viewpoint of the "presence of functional intrinsic disorder in proteins", a noticeable amount of biologically active proteins (of protein regions) fail to fold into the well-defined structures and instead remain disordered, existing as highly dynamic ensembles of rapidly interconverting conformations under the physiological conditions. These proteins and protein regions are known now as intrinsically disordered proteins (IDPs) and intrinsically disordered protein regions (IDPRs), respectively. The propensity of being functional intrinsically disordered proteins (similar to the propensity of forming unique biologically active structures of ordered proteins) is determined by the amino acid sequences [12][13][14]. IDPs exhibit their biological functions in numerous biological processes commonly associated with cellular signalling, gene regulation, and control by interacting with their physiological partners [15][16][17][18][19]. These functions of IDPs and IDPRs are regulated by their protein-protein, protein-RNA, and protein-DNA interactions [20,21]. Molecular recognition features (MoRFs) are the regions in IDPs implicated in the regulation of IDP function by protein-protein interactions and serve as the primary stage in molecular recognition.
It is known that the IDPs/IDPRs are present in all three kingdoms of life, and viral proteins often contain unstructured regions that have been strongly correlated with their virulence [22][23][24][25]. In this report, we have investigated the disordered side of SARS-CoV-2 proteome using a complementary set of computational approaches to check the prevalence of IDPRs in its proteins and to shed some light on their disorder-related functions. We also have comprehensively analysed IDPRs among the closely related viruses, human SARS and bat SARS-like CoVs. Furthermore, we have also identified protein functions related to protein-protein interactions, RNA binding, and DNA binding from all three viruses. Since these three viruses are closely related, our study provides an important means for a better understanding of the sequence and structural peculiarities of their evolution.

Sequence retrieval and multiple sequence alignment
The protein sequences of bat CoV (SARS-like) and human SARS CoV were retrieved from UniProt (UniProt IDs for individual proteins are listed in Table 1). The translated sequences of SARS-CoV-2 proteins [GenBank database [26] (Accession ID: NC_045512.2)] were obtained from GenBank. We used these sequences for performing multiple sequence alignment (MSA) and predicting the IDPRs. We have used Clustal Omega [27] for protein sequence alignment and Esprit 3.0 [28] for constructing the aligned images.

Per-residue predictions of intrinsic disorder predisposition
For the prediction of the intrinsic disorder predisposition of CoV proteomes, we used multiple predictors, such as members of the PONDR ® (Predictor of Natural Disordered Regions) family including PONDR ® VLS2 [29], PONDR ® VL3 [30], PONDR ® FIT [31], and PONDR ® VLXT [32], as well as the IUPred platform for predicting long (≥ 30 residues) and short IDPRs (< 30 residues) [33]. These computational tools predict residues/regions, which do not have the tendency to form an ordered structure. Residues with disorder scores exceeding the threshold value of 0.5 are considered as intrinsically disordered residues, whereas residues with the predicted disorder scores between 0.2 and 0.5 are considered flexible. Complete predicted percent of intrinsic disorder (PPID) in a query protein was calculated for every protein of all the three viruses from outputs of six predictors. The detailed methodology has been given in our previous reports [34,35].

Combined CH-CDF analysis to predict disorder predisposition of proteins
The charge hydropathy plot [36] and PONDR ® VLXT-based cumulative distribution function are two binary predictors of disorder (i.e. tool evaluating entire protein as mostly ordered or mostly disordered), which are available on the PONDR web server (https ://www.pondr .com). Combining the result from these binary predictors helps to classify the proteins into different groups, depending on their global disorder [37].

Molecular recognition feature (MoRF) determination in CoV proteomes
The authentic online bioinformatics predictors that use a different set of algorithms for the prediction of MoRFs were used. These include MoRFchibi_web [38], ANCHOR [39,40], MoRFPred [41], and DISOPRED3 [42]. The protein residues with ANCHOR, MoRFPred, and DISOPRED3 scores above the threshold value of 0.5 and MoRFchibi_ web score above the threshold value of 0.725 are considered MoRF regions.

Identification of DNA-and RNA-binding regions in CoV proteomes
Often, IDPs and IDPRs facilitate interactions with RNAs and DNAs and regulate many cellular functions [43]. Thus, for predicting the DNA-binding residues in CoV proteins, we used two online servers: DRNAPred [20] and DisoR-DPbind [43]. For RNA-binding residues, we used PPRInt (Prediction of Protein RNA-Interaction) [44] and DisoRD-Pbind servers [43].

Comprehensive computational analysis of intrinsic disorder in structural and accessory proteins of SARS-CoV-2, human SARS and bat CoV (SARS-like)
The mean values of the predicted percentage of intrinsic disorder scores (mean PPIDs) that were obtained by averaging the predicted disorder scores from six disorder predictors (Tables S1-S3) for structural and accessory proteins of SARS-CoV-2 as well as human SARS, and bat CoV are represented in Table 1.
Figure 2a-c are 2D-disordered plots generated for SARS-CoV-2, human SARS and bat CoV proteins, respectively, and represent the PPID PONDR-FIT vs. PPID Mean plots. Based on their predicted levels of intrinsic disorder, proteins can be classified as highly ordered (PPID < 10%), moderately disordered (10% ≤ PPID < 30%) and highly disordered (PPID ≥ 30%) [46]. From the data in Table 1, Fig. 2a-c, as well as the PPID based classification, we conclude that the nucleocapsid protein from all three CoVs possess the highest percentage of disorder and, therefore, is classified as a highly disordered protein. The ORF3b protein in bat CoV, ORF6 protein of all three CoVs, and ORF9b proteins of SARS-CoV-2 and bat CoV belongs to the class of moderately disordered proteins. While the structured proteins, namely, spike glycoprotein (S), envelope protein (E) and membrane protein (M) as well as accessory proteins ORF3a, ORF7a, ORF8 (ORF8a and ORF8b in case of human SARS) of all three strains of CoVs are ordered proteins. ORF14 and ORF10 proteins are also ordered proteins.
To further investigate the nature of the disorder in proteins of all three CoVs, we utilized the combined CH-CDF tool that uses the outputs of two binary classifiers of disorder: charge hydropathy (CH) plot and cumulative distribution function (CDF) plot. This helped in retrieving more detailed characterization of the global disorder predisposition of query proteins and their classification according to the disorder "favors". The CH plot is a linear classifier that differentiates between proteins that are predisposed to extended disordered conformations which includes random coils and pre-molten globules from proteins that have compact conformations (ordered proteins and molten globule-like proteins). The other binary predictor, CDF, is a nonlinear classifier that uses PONDR ® VLXT scores to discriminate ordered globular proteins from all disordered conformations, which include native molten globules, premolten globules, and random coils. The CH-CDF plot can be divided into four quadrants: Q1 (bottom right quadrant) containing ordered proteins; Q2 (bottom left quadrant) includes proteins predicted to be disordered by CDF and compact by CH (i.e. native molten globules and hybrid proteins containing high levels of both ordered and disordered regions); Q3 (top left quadrant) contains proteins that are predicted to be disordered by both CH and CDF analysis (i.e. highly disordered proteins with the extended disorder); and Q4 (top right quadrant) possesses proteins disordered according to CH but ordered according to CDF analysis [34]. Figure 2d-f represent the CH-CDF analysis of proteins of SARS-CoV-2, human SARS, and bat CoV and shows that all the proteins are located within the two quadrants Q1 and Q2. The CH-CDF analysis leads to the conclusion that all proteins of all three CoVs are ordered except nucleocapsid protein, which is predicted to be disordered by CDF as well as CH and hence lies in Q3. Molecular recognition features (MoRFs) are short interaction-prone disordered regions found within IDPs/IDPRs that commence a disorder-to-order transition upon binding to their partners [47,48]. In this study, we have analysed and compared MoRFs (protein-binding regions) in SARS-CoV-2 with human SARS and bat CoVs. The results of this analysis are summarized in Table 2, which clearly shows that most of the SARS-CoV-2 proteins contain at least one MoRF. This is indicative of an important role played by disorder in functionality of these viral proteins. All of the SARS-CoV-2 proteins have been predicted to contain MoRFs except ORF7b and Nsp13 proteins. MoRFs in human SARS and bat CoV proteomes are listed in Tables S7 and S8. Similar to SARS-CoV-2 proteome, bat CoV proteins ORF7b, and Nsp13 are not predicted to have any MoRF by any of the servers used. In human SARS proteome, proteins ORF7b, Nsp13, Nsp2, and Nsp15 do not show the presence of any MoRF. Interestingly, the N protein from SARS-CoV-2, human SARS, and bat CoV shows high number of variable MoRFs, signifying its central role in virus pathogenesis.

Nucleotide-binding propensity in proteins of coronaviruses
In addition to protein-protein interactions/protein-binding functions, IDPs and IDRs also mediate functions by facilitating their interactions with nucleotides (DNA and RNA) [21,49]. Therefore, we have used a combination of two different online servers for locating protein residues that show the propensity to bind with DNA as well as RNA. The nucleotide-binding residues in proteins of the three studied coronaviruses are listed in Tables S9-S11. Interestingly, all the viral proteins of SARS-CoV-2, human SARS, and bat CoV have shown the propensity to bind to nucleic acids. In particular, structural (S, M, and N proteins) and non-structural (Nsp 2, 3, 4, 5, 6, 12, 13, 14, 15, and 16) proteins of all three viruses display a e human SARS and f bat CoV, the Y coordinate of each protein spot signifies the distance of corresponding protein from the boundary in the CH plot and X coordinate value corresponds to the average distance of CDF curve for respective protein from the CDF boundary large number of RNA-binding residues. However, ORF3a, ORF3b, ORF6, ORF7a, ORF7b, ORF8, and ORF14 proteins show less RNA-binding and more DNA-binding residues.

Intrinsic disorder analysis of structural proteins of coronaviruses
Coronaviruses encode four structural proteins, namely, spike (S), envelope (E) glycoprotein, membrane (M), and nucleocapsid (N), which are translated from the last ~ 10 kb nucleotides and forms the outer cover of the CoVs, encapsulating their single-stranded genomic RNA.

Spike (S) glycoprotein
The S protein is a large multifunctional protein forming the exterior of the CoV particles [50,51]. It forms surface homotrimers and contains two distinct ectodomain known as S1 and S2. Subunit S1 initiates viral infection by binding to the host cell receptors, while S2 acts as a class I viral fusion protein that mediates the fusion of the virion and cellular membranes and thereby promotes the viral entry into the host cells [52,53]. It binds to specific surface receptor angiotensin-converting enzyme 2 (ACE2) on host cell plasma membrane through its N-terminal receptor-binding domain (RBD) [54].
S protein consists of an N-terminal signal peptide, a long extracellular domain, a single-pass transmembrane domain, and a short intracellular domain [55]. A 3.60 Å resolution structure (PDB ID: 6ACC) of human SARS S protein complexed with its host-binding partner ACE2 is obtained using cryo-electron microscopy (cryo-EM) (Fig. 3b). In this PDB structure, few residues (1-17, 240-243, 661-673, 812-831 and 1120-1203) are missing [56], suggesting their flexible nature. Also, the structure of S protein (3.5 Å) from SARS-CoV-2 has been recently deduced by Wrapp et al. using   [57]. MSA analysis among all three coronaviruses demonstrates that the S protein of SARS-CoV-2 has a 77.71% sequence identity with bat CoV and 77.14% identity with human SARS (Fig. S1). As observed, there is a significant sequence variation in RBD located at the N-terminal region which might affect its virulence, receptor-mediated binding and entry into the host cell.
According to our intrinsic disorder propensity analysis, the S protein from all three CoVs are found to be highly structured ( Table 1). The mean PPID scores of SARS-CoV-2, human SARS, and bat CoVs are calculated to be 1.41%, 1.12%, and 1.85%, respectively. Figure 3c-e represent the intrinsic disorder profiles of S proteins from SARS-CoV-2, human SARS and bat CoV obtained from six disorder predictors. Finally, Fig. 3f shows aligned disorder profiles of S proteins from these CoVs and illustrates remarkable similarity in their disorder propensity, especially in the C-terminal region.
It is of interest to map known functional regions of S proteins to their corresponding disorder profiles. The maturation of S protein requires specific post-translational modifications (PTM), proteolytic cleavage that happens at two stages. First, host cell furin or another cellular protease nicks the S precursor to generate S1 and S2 proteins, whereas the second cleavage takes place after the viral attachment to host cell receptors which leads to the release of a fusion peptide generating the S2′ subunit. In human SARS, the first and second cleavage site are located at residues R 667 and R 797 , respectively, whereas in bat CoV, the corresponding cleavage sites are residues R 654 and R 784 . As it follows from Fig. 3, these cleavage sites are located within the IDPRs. In human SARS S protein, fusion peptide (residues 770-788) located within a flexible region is characterized by a mean disorder score of 0.232 ± 0.053. Similarly, in bat CoV S protein, fusion peptide (residues 757-775) has a mean disorder score of 0.320 ± 0.046. It contains two heptad repeat regions that form coiled-coil structure during viral and target cell membrane fusion, assuming a trimer-of-hairpins structure needed for the functional positioning of fusion peptide. In human SARS, heptad repeat regions are formed by residues 902-952 and 1145-1184, which have mean disorder scores of 0.458 ± 0.067 and 0.353 ± 0.062, respectively. The analogous situation is observed for S protein of bat CoV, where these heptads repeat regions are positioned at residues 889-939 (0.44 ± 0.11) and 1132-1171 (0.353 ± 0.062).
Another functional region found in S proteins is the RBD (residues 306-527 and 310-514 in human SARS and bat CoVs, respectively) containing a receptor-binding motif responsible for interaction with human ACE2. In S protein of human SARS, this motif (residues 424-494) is not only characterized by structural flexibility, possessing a mean disorder score of 0.30 ± 0.16, but also contains a disordered region (residues 461-466). Since S protein is known as spike glycoprotein, it contains numerous glycosylation sites. Due to rather close similarity of disorder profiles of S proteins analysed here, we can assume that all the aforementioned indications of the functional importance of disorder and flexible regions in S proteins from human SARS and bat CoVs are also applicable to SARS-CoV-2 S protein. Finally, Table 2 shows that S protein from SARS-CoV-2 contains a MoRF region at its C-terminal (residues 1265-1272) as predicted by MoRFchibi_web, two MoRF regions ((residues 2-6, 819-823) by MoRFPred, and one MoRF region at the N-terminal (residues 1-10) by DISOPRED3. These results indicate that intrinsic disorder is important for its interaction with binding partners. Strikingly, the N-terminal region of S protein (residues 1-10) from all three viruses are predicted to be a MoRF by two servers (MoRFPred and DISOPRED3). This displays its role in viral interaction with host receptor, while the C-terminal MoRF is engaged in interaction with M protein for assembly of viral particles [58]. Moreover, MoRF regions lying in the N-and C-terminal regions suggest their possible role during cleavage as well. In addition to proteinbinding regions, S protein also shows many nucleotide-binding residues. Tables S9-S11 shows numerous RNA-binding residues predicted by PPRInt in all three viruses. Further, DRNApred and DisoRDPbind predicted the presence of many DNA binding residues in all three S proteins. These results signify the molecular recognition (protein-protein interaction, RNA binding, and DNA binding) and interactions with host cell membrane and further viral infection. Therefore, IDPs/IDPRs and residues/regions in S proteins that are crucial for molecular recognition can be targeted for disorder-based drug discovery.

Envelope (E) small membrane protein
Envelope (E) protein is a small, multifunctional membrane protein that plays an important role in the assembly and morphogenesis of virions in the cell [59][60][61]. It consists of two ectodomains associated with N-and C-terminal regions, and a middle transmembrane domain. It homo-oligomerizes into a pentameric membrane destabilizing transmembrane hairpins to form a pore necessary for its ion channel activity [62]. Figure 4a shows the NMR-structure (PDB ID: 2MM4) of human SARS envelope glycoprotein of 8-65 residues [63].
MSA results illustrate (Fig. 4b) that this protein is highly conserved, with only three amino acid substitutions in E protein of SARS-CoV-2 conferring its 96% sequence similarity with human SARS and bat CoV. Also, bat CoV shares 100% sequence identity with human SARS. Mean PPID calculated for SARS-CoV-2, human SARS, and bat CoV E proteins are 5.33%, 6.58%, and 6.58%, respectively ( Table 1). The E protein is found to have a reasonably well-predicted structure; however, residues of N-and C-terminals display a higher tendency for the disorder (disorder profiles in Fig. 4c-e). Evidences show that the last 18 hydrophilic residues (residues 59-76) adopt a random-coil conformation with and without the addition of lipid membranes [64]. Further, the last four amino acids of the C-terminal region containing a PZD-binding motif are involved in protein-protein interactions with a tight junction protein PALS1. PALS1 is involved in maintaining the polarity of epithelial cells in mammals [65]. Our results support the existing literature, as we identified a long N-terminal region of ~ 30 residues as a MoRF region in all three viruses (see Tables 2, S7, S8). We speculate that the disordered regions may facilitate interactions with other proteins as well. In agreement with this hypothesis, the C-terminal domain of SARS-CoV-2 E protein serves as a protein-binding region. We also found that residues from 45-75 is a long MoRF in E proteins of all three viruses as predicted by MoRFchibi_web. As aforementioned, these randomly coiled binding residues at the C-terminus may gain structure while assisting the protein-protein interaction mediated by E protein. One more MoRF region (residues [26][27][28][29][30] in the transmembrane domain was observed by DISOPRED3 in all three E proteins. Since these residues are part of the ion channel, they may be involved in guiding the specific function of ion channel activity. Few nucleotidebinding residues are predicted for all three E proteins (Tables S9-S11).

Membrane (M) glycoprotein
It plays an important role in virion assembly by interacting with the nucleocapsid (N) and E proteins [66][67][68]. Protein M interacts specifically with coronavirus RNA containing a short viral packaging signal in the absence of N protein, highlighting an important nucleocapsid-independent viral RNA packaging mechanism inside the host cells [69]. Cryo-EM and tomography data reveal its two distinct conformations, a compact structure having high flexibility and low spike density, and an elongated M protein having a rigid structure and narrow range of membrane curvature [70]. Although no structural information is available for fulllength M protein, a short peptide of the membrane glycoprotein (residues 88-96) from human SARS is co-crystallized with a complex of A-2 α chain of the HLA class I histocompatibility antigen and β2-microglobulin (PDB ID:3I6G) [71]. Figure 5a shows the extended conformation of M protein.
The M protein of SARS-CoV-2 has a sequence similarity of 90.1% with bat CoV and 89.6% with human SARS M proteins (Fig. 5b). Disorder profiles in Fig. 5c-e show a relatively low level of disorder in M proteins of SARS-COV-2 (2.70%), human SARS CoV (1.36%,), and bat CoV (1.36%). This is consistent with a previous publication by Goh et al. on human SARS HKU4, where they found the mean PPID of 4% using additional predictors such as TopIDP and FoldIndex along with the predictors used in our study [72]. The last 20 residues of MERS-CoV M protein are important for intracellular trafficking and contain a determinant that localizes it into the Golgi network [73]. MoRF analysis revealed that the disordered C-tail of M protein contains a MoRF region which can serve as a binding site for its partner required during localization inside the host cell. A long MoRF region (residues 186-220) at the C-terminal of M protein in all three viruses is located by MoRFchibi_web. Two MoRF regions [one at N-terminus (residues 1-16) and one at the C-terminus (residues 205-221)] are predicted by DISO-PRED3 in human SARS and bat CoV. However, only a single MoRF (residues 117-132) is observed in SARS-CoV-2 (by DISOPRED3) (Tables 2, S7, S8). Furthermore, the M protein from all three viruses displays strong tendency to bind with RNA (as predicted by PPRInt and DisoRDPbind) and DNA (as predicted by DRNApred and DisoRDPbind) (see Tables S9-S11). Our understanding of M protein of CoVs (IDPs and MoRF at C-terminus) elucidates its critical role in interaction with the N and E proteins for viral assembly.

Nucleocapsid (N) protein
It is one of the major viral proteins playing an essential role during transcription, and virion assembly of CoVs [74]. It binds to viral genomic RNA forming a ribonucleoprotein core required for RNA encapsidation during viral particle assembly [75]. It consists of two structural domains, the N-terminal RNA-binding domain (NTD: 45-181 residues) and the C-terminal dimerization domain (CTD: 248-365 residues) with a disordered patch between these domains. It is demonstrated to bind with viral RNA using both NTD and CTD [76]. Recently, residues 50-173 of the N protein of SARS-CoV-2 has been crystallized (PDB ID: 6VYO) (Fig. 6a). Figure 6b1 displays the NMR solution structure of NTD (45-181 residues) of human SARS N protein (PDB ID: 1SSK) [77]. Figure 6b2 shows an X-ray crystal structure of CTD of human SARS N protein (270-366 residues) (PDB ID: 2GIB) [78]. A model of domain organization of N-protein from SARS-CoV-2 is shown in Fig. 6c.
The 419 amino acid-long N protein of SARS-CoV-2 shows a percentage identity of 88.76% and 89.74% with N proteins of bat and human SARS CoVs (Fig. S2). Our analysis revealed the highest levels of intrinsic disorder in N proteins of all three CoVs (graphs in Fig. 6d-f), which is in accordance with the previously evaluated intrinsic disorder predisposition [72]. In fact, N proteins from SARS-CoV-2, human SARS, and bat CoVs are characterized by the mean PPIDs of 64.91%, 71.09%, and 65.80%, respectively. This is further supported by Fig. 6g, where PONDR ® VSL2-generated disorder profiles of these three proteins are overlapped to show almost complete coincidence of their major disorder-related features. In particular, SARS-CoV-2 N protein residues 1-57, 64-102, 145-162, 166-289, and 362-422 are found to be disordered (Fig. 6d). Many of these residues lie within the NTD and CTD regions, which due to their structural plasticity does not get crystallized in human SARS N protein crystal structure. Overall, all three N proteins are found to be highly disordered.
Tables 2, S7 and S8 show that the N protein is heavily decorated with MoRFs, suggesting that this protein is a promiscuous binder. Long disorder-based protein bonding regions at the N-and C-terminus of the N protein of all three viruses are observed by all four predictors. The N protein from human SARS has one phosphorylation site (residue S 177 ) and several regions with compositional biases, such as Ser-rich (residues 181-213), Poly-Leu, Poly-Gln, and Poly-Lys (residues 220-225, 240-245, and 370-376), all predicted to be disordered. Similarly, the N protein of bat CoV, S 176 is phosphorylated and has Ser-rich, Poly-Leu, and Poly-Lys regions (residues 176-206, 219-224, and 369-375, respectively), all of which are disordered. It has been reported to interact using the central disordered region with M protein, hnRNP A1, and self-N-N interaction [79][80][81]. The middle flexible region is also responsible for its RNA-binding activity [82]. Deletion of 184-196 residues, 169-308 residues, and 161-210 residues of N abolishes its multimerization, RNA-binding capacity, and hnRNP A1 interactions, respectively. The MoRFs present in the aforementioned regions may mediate these interactions of N proteins. Figure 6b2 represents another important disorder-related functional feature of N protein. CTD homodimer shown is characterized by highly intertwined morphology, which is typically a result of binding-induced folding [83][84][85], indicating that a very significant part of CTD gains structure during dimerization. We identified numerous RNA-binding residues in all three viruses using PPRInt server. This finding supports the function of N protein as it interacts with genomic RNA for a ribonucleoprotein core formation, which is a crucial step for RNA encapsidation. Additionally, DRNApred and DisoRDPbind predict multiple DNA-binding residues in N protein of all the studied CoVs. The flexible (IDPRs) regions at the N-and C-terminus of SARS-CoV-2 have long protein-binding as well as nucleotide-binding regions that may play a vital role in its interaction with viral RNA. These flexible regions can be targeted to inhibit the interaction of N protein with viral genomic RNA.

Intrinsic disorder analysis of accessory proteins of coronaviruses
Literature suggests that some viral proteins are translated from the genes interspersed in between the genes of structural proteins. These proteins are known as accessory proteins, and many of them are proposed to be involved in viral pathogenesis [86].

Proteins ORF3a and ORF3b
ORF3a is a multifunctional protein (of molecular weight ~ 31 kDa) that performs a major function during virion assembly by co-localizing with E, M, and S viral proteins [87][88][89][90][91]. The homo-tetrameric complex of ORF3a has been demonstrated to form a potassium-ion channel on the host cell plasma membrane [92]. ORF3b protein can be found in the cytoplasm, nucleolus, and outer membrane of mitochondria of the host cells [93,94]. In Huh 7 cells, its over-expression has been linked with the activation of AP-1 via the ERK and JNK pathways [95].

Protein ORF6
Also known as P6, this membrane-associated protein serves as an interferon (IFN) antagonist [96]. Using its C-terminal residues, ORF6 disrupts karyopherin import complex in cytosol and, therefore, hampers the movement of transcription factors like STAT1 into the nucleus resulting in downregulation of the IFN pathway [96,97]. It contains a YSEL motif near its C-terminal region that functions in protein   Fig. 3 internalization from the plasma membrane into the endosomal vesicles [98].
MSA results demonstrate that (Fig. 9d), SARS-CoV-2 ORF6 is closer to human SARS ORF6, having a sequence similarity of 68.85% than to bat CoV ORF6 (67.21%). Novel SARS-CoV-2 ORF6 is predicted to be the second most disordered structural protein with a PPID of 22.95%, containing a disordered C-terminal region.
The mean PPID of the other two ORF6 proteins are listed in Table 1. Graphs in Fig. 9a-c illustrate that all three ORF6 proteins are moderately disordered with the presence of high disorder near C-terminal residues. As aforementioned, this hydrophilic region contains lysosomal targeting motif (YSEL) and diacidic motif (DDEE) responsible for its binding and recognition during translocation [98], this region is important for the biological activities of ORF6. Moreover, the N-terminus does not contain any prominent disorder. First, 38 amino acids of human SARS ORF6 are described to form an α-helical structure spanning the membrane [99]. A long MoRF region [(residues 26-61 in SARS-CoV-2), (residues 31-63 in Human SARS), and (residues 30-60 in bat Cov)] is also present near the C-terminus. It represents very few RNA-and DNA-binding residues.

ORF7a and ORF7b proteins
ORF7a is a type I transmembrane protein [100,101]. It contributes to viral pathogenesis by activating the release of proinflammatory cytokines and chemokines, such as IL-8 and RANTES [102,103]. The presence of a KRKTE motif near the C-terminal region is needed for its import from ER to Golgi apparatus [100,101]. On the other hand, ORF7b is an integral membrane protein that has been shown to localize in the Golgi complex [104,105]. These reports also confirm the role of ORF7b as an accessory as well as a structural protein in SARS-CoV virion [104,105]. Figure 10d represents the 1.8 Å X-ray crystal structure of the 14-96 fragment of the ORF7a from human SARS (PDB ID: 1XAK) and demonstrates the compact seven β-stranded topology of this protein similar to the Ig-superfamily members [106]. Importantly, in this crystal structure, residues 82-96 constitute the region with missing electron density, indicating the highly dynamic nature of this segment. In line with this hypothesis, NMR solution structure of the 16-99 fragment of ORF7a of human SARS (PDB ID: 1YO4) showed that residues 81-99 are highly disordered [107].
As observed from Table 1 and Fig. 10a-c, our disorder predisposition analyses resulted in the overall PPID values for ORF7a proteins-1.65% (SARS-CoV-2), 0.82% (bat CoV), and 0.82% (Human SARS). The mean PPIDs estimated for ORF7b proteins are 9.30% for SARS-CoV-2, 4.55% for bat CoV and 4.55% human SARS. Table 2 shows  Fig. 3 the presence of several MoRFs in ORF7a, indicating its potential involvement in disorder-dependent protein-protein interactions. At the N-terminus, one MoRF region (residues 1-10) is predicted by DISOPRED3 in all three ORF7a proteins. In addition to protein-binding regions, ORF7a also contains several RNA-and DNA-binding residues. Analysis also reveals the low disorder content in all three ORF7b proteins ( Fig. S3A-C), and subsequently no MoRFs. Although ORF7b does not contain protein-binding regions, it has many nucleotide (both RNA and DNA)-binding residues. Figures S3A, 3B, and 3C depict the residues predisposed for disorder in ORF7b proteins of SARS-CoV-2, human SARS CoV, and bat CoV, respectively. In particular, both proteins in all three studied viruses have ordered structures.

Proteins ORF8a and ORF8b
In isolates from early human infections, the ORF8 gene codes for a single ORF8 protein. However, in late infections, more specifically, at middle and late stages, a 29 nucleotide deletion in the ORF8 gene led to the formation of two distinct proteins, ORF8a and ORF8b containing 39 and 84 residues, respectively [108,109]. Both proteins have conformations different from that of the longer ORF8 protein and interacts with different structural proteins [110]. The disorder-based protein-binding regions of this protein identified in this study may have an important role in interaction with other proteins.

ORF9b protein
This protein is expressed from an alternative ORF within the N gene through a leaky ribosome-binding process [111]. This protein is shown to interact with a nuclear export protein receptor Exportin 1 (Crm1), using which it's translocated out of the nucleus [112]. Our MoRFs analysis shows the presence of disorder-based protein-binding regions in ORF9b protein which may have a role in its interaction with Crm1 for translocation outside the nucleus. A 2.8 Å resolution crystal structure of ORF9b protein from human SARS CoV (PDB ID: 2CME) shows the presence of a dimeric tentlike β-structure along with the central hydrophobic amino acids (Fig. 11d) [113].
Based on the sequence availability (accession ID NC_045512.2), translated protein sequence of ORF9b is not reported for SARS-CoV-2. However, based on a report by Wu and colleagues [45], corresponding annotated sequence is used for intrinsic disorder analysis. According to the MSA (results shown in Fig. 11e), ORF9b protein from SARS-CoV-2 shares 73.2% identity with human SARS and 74.23% identity with bat CoV.
Our IDP analysis (Table 1) exposed moderate disorder content in ORF9b of human SARS having a mean PPID of 26.53%. As depicted in Fig. 11a-c, disorder in human SARS ORF9b protein mainly lies near the N-terminal end (residues 1-10) and near the central region (residues 28-40) with a well-ordered inner core. The X-ray crystal structure of ORF9b has a missing electron density of first 8 residues and 26-37 residues near the central region. This indicates that the corresponding regions are disordered, which are difficult to crystallize due to their highly dynamic structural organization. SARS-CoV-2 ORF9b with a mean PPID of 10.31% also has an N-terminal (1-10 residues) disordered segment. ORF9b of bat CoV is shown to have an intrinsic disorder content of 9.28%, comparatively lower than the other two ORF9b proteins. MoRFs lies in the N-terminal region of ORF9b proteins of all three viruses (Tables 2, S7, S8). In the absence of other viral proteins, its first 41 residues are demonstrated to induce membranous structures similar to DMVs [99]. The available crystal structure also has a missing electron density in the N-terminal region suggesting that these flexible amino acids are likely to interact with host lipids. The 3-29 amino acid segment of SARS-CoV-2 is identified as disorder-based protein-binding region that may mediate its interaction with host lipids for the formation of DMVs.

ORF10 protein
The newly emerged SARS-CoV-2 has an ORF10 protein of 38 amino acids. ORF10 of SARS-CoV-2 has a 100% sequence similarity with ORF10 of bat CoV strain  Fig. 3 bat-SL-CoVZC45 [11]. However, we did not conduct the disorder analysis for ORF10 from the bat-SL-CoVZC45 strain, since all our studies reported here are related to a different strain of bat CoV (reviewed strain HKU3-1). Therefore, we have only reported the results of disorder analysis for the ORF10 protein from SARS-CoV-2, according to which this protein has a mean PPID of 0.00% (see also Fig. S6 for disorder profile of ORF10). This protein contains a MoRF from three to seven residues at its N-terminus as predicted by MoRFchibi_web. Further, we predicted its binding tendency to nucleotides and found the presence of few RNA-binding sites; however, it does not contain DNAbinding residues.

Protein ORF14
This is a 70 amino acid long uncharacterized protein of unknown function. According to the MSA, ORF14 of SARS-CoV-2 has 77.1% identity with human-SARS and 72.9% identity with bat CoV as represented in Fig. S7D. Figure S7A-C shows the resulting disorder profiles of all three ORF14 proteins (mean PPIDs are listed in Table 1). Further, these proteins have calculated mean PPID values of 0.00%, 2.86%, and 0.00%, respectively. These proteins have flexible N-and C-terminal regions. It can use intrinsic disorder or structural flexibility for protein-protein interactions since it possesses MoRFs. It mainly contains MoRFs at the N-and C-terminal regions (Tables 2, S7, S8) and several RNA-and DNA-binding residues (Tables S9-S11). These regions indicate its vital role in protein function related to protein-RNA and protein-DNA interaction.

Intrinsic disorder analysis of non-structural proteins of coronaviruses
In coronaviruses, due to ribosomal leakage during translation, two-third of the RNA genome is processed into two polyproteins: (i) replicase polyprotein 1a and (ii) replicase polyprotein 1ab. Both contains non-structural proteins (Nsp1-10) in addition to different proteins required for viral replication and pathogenesis. Replicase polyprotein 1a contains an additional Nsp11 protein of 13 amino acids, the function of which has not been investigated yet. The longer replicase polyprotein 1ab of 7073 amino acids accommodates five other non-structural proteins (Nsp12-16) [114,115]. Table 3 represents the mean PPID scores of 15 Nsps derived from the replicase polyprotein 1ab in SARS-CoV-2, human SARS, and bat CoV. These values were obtained by combining the results from six disorder predictors (see Tables S4-S6). Figure 12a-c represents 2D-disordered plots of the Nsps coded by ORF1ab in SARS-CoV-2, Human SARS, and bat CoV, respectively. Based on the mean PPID scores in Table 3, Fig. 12a-c, and taking into PPID based classification [46], we conclude that none of the Nsps in SARS-CoV-2, human SARS, and bat CoV are highly disordered. Only Nsp1 and Nsp8 proteins are found to be moderately disordered (10% ≤ PPID ≤ 30%). We also observed that Nsp2, Nsp3, Nsp5, Nsp6, Nsp7, Nsp9, Nsp10, Nsp15, and Nsp16 have less than 10% disordered residues and hence, belong to the category of mostly ordered proteins. Other non-structural proteins, namely, Nsp4, Nsp12, Nsp13, and Nsp14 have negligible levels of disorder (PPID < 1%) and are concluded to be highly structured.

Global analysis of intrinsic disorder in the replicase polyprotein 1ab
The CH-CDF analysis of Nsps from SARS-CoV-2, human SARS and bat CoV is depicted in Fig. 12d-f respectively. It was observed that all Nsps of the three CoVs are located within the quadrant Q1 of the CH-CDF phase space, which is indicative of their ordered structure.

Replicase polyprotein 1ab
The longer replicase polyprotein contains 15 Nsps listed in Table 3. Nsp1, Nsp2, and Nsp3 are cleaved using a viral papain-like proteinase (Nsp3/PL-Pro), while the rest of the Nsps are cleaved by another viral 3C-like proteinase, Nsp5/3CL-Pro.We mapped the cleavage sites of the replicase 1ab polyprotein from human SARS CoV to the disorder profile of this polyprotein. Figure 13 represents the results of this analysis by showing zoomed-in regions surrounding all the cleavage sites with few residues spanning at both terminals. Interestingly, we observed that all the cleavage sites are largely disordered, suggesting that intrinsic disorder may have a crucial role in the maturation of individual non-structural proteins. As Nsps of human SARS are evolutionarily closer to Nsps of SARS-CoV-2, we hypothesize that cleavage sites in the SARS-CoV-2 replicase 1ab polyprotein are also intrinsically disordered or flexible. To shed more light on other implications of IDPRs, the structural and functional properties of Nsps and their predicted IDPRs are thoroughly described below.

Non-structural protein 1 (Nsp1)
This protein acts as a host translation inhibitor as it binds to the 40S subunit of ribosome and blocks the translation of cap-dependent mRNAs as well as mRNAs that uses the internal ribosome entry site (IRES) [116]. Figure 14a shows the NMR solution structure (PDB ID: 2GDT) of human SARS Nsp1 protein (13-128 residues) [117].
SARS-CoV-2 Nsp1 shares 84.44% and 83.80% sequence identity with Nsp1s of human SARS and bat CoV, respectively (Fig. 14b). The respective mean PPIDs of Nsp1s from SARS-CoV-2, Human SARS, and bat CoV are 12.78%, 14.44%, and 12.85% (disorder profiles in Fig. 14c-e). In particular, the following regions are predicted to be disordered: SARS-CoV-2 (residues 1-7 and 165-180), human SARS (residues 1-5 and 165-180), and bat CoV (residues 1-5 and 165-179). NMR solution structure of Nsp1 from human SARS revealed the presence of two unstructured segments near the N-terminal (1-12 residues) and C-terminal (129-179 residues) regions [117]. The disordered region   (128-180 residues) at the C-terminus is already mapped important for its expression [118]. Based on sequence homology with human SARS Nsp1, the predicted disordered C-terminal region of SARS-CoV-2 Nsp1 may play a critical role in its expression. Alanine mutants at K164 and H165 near the C-terminal region are reported to abolish its binding with the 40S subunit of the host ribosome [119]. In conjunction with this data, several MoRFs are present in the unstructured segments of Nsp1 proteins. These regions are shown in Tables 2, S7 and S8.

Non-structural protein 2 (Nsp2)
This protein functions by disrupting the host survival pathway via interaction with the host proteins prohibitin-1 and prohibitin-2 [120]. Reverse genetic deletion in the coding sequence of Nsp2 of SARS virus attenuated little viral growth as well as replication and allowed the recovery of mutant virulent viruses [121]. The sequence identity of Nsp2 protein of SARS-CoV-2 with Nsp2s of human SARS and bat CoV amounts to 68.34% and 68.97%, respectively (Fig. S8). We have estimated the mean PPIDs of Nsp2s of SARS-CoV-2, human SARS, and batbat CoV to be 5.17%, 2.04%, and 2.03% respectively (see Table 3) (per-residue predisposition of intrinsic disorder is depicted in Fig. S9A-C). According to the results, residues 570-595 (SARS-CoV-2), residues 110-115 (Human SARS), and residues 112-116 (bat CoV) are predicted to be disordered. As listed in Tables 2, S7 and S8, human SARS does not contain MoRF, while SARS-CoV-2 and bat CoV have a N-terminally located MoRF region predicted by MoRFchibi_web.

Non-structural protein 3 (Nsp3)
Nsp3 is a viral papain-like protease (PLP) that affects the phosphorylation and activation of IRF3 and, therefore, antagonizes the IFN pathway [122]. It's also reported to stabilize NF-κβ inhibitor which further blocks the NF-κβ pathway [122]. Figure 15d represents the 1.85 Å resolution X-ray crystal structure of the catalytic core of Nsp3 protein from human SARS CoV (PDB ID: 2FE8) [123]. The structure consisting of residues 723-1036 revealed folds similar to a deubiquitinating enzyme in vitro, the deubiquitinating activity of which was found to be efficiently high [123]. A 1.45 Å resolution structure (PDB ID: 6W6Y) of SARS-CoV-2 Nsp3 homodimer (chains A and B from 207-374 residues) is recently generated using X-ray diffraction (Fig. 15e) [124].
Nsp3 protein of SARS-CoV-2 contains several substituted residues throughout the protein. It is equally close to both Nsp3 proteins of human SARS and bat CoV, sharing 76.69% and 76.31% identity respectively (Fig. S10). According to our results, the mean PPIDs of Nsp3 proteins of SARS-CoV-2, human SARS, and bat CoV are 7.40%, 7.91%, and 7.78% respectively (Table 3). Disorder profiles in Fig. 15a-c shows that all three Nsp3 proteins are highly structured. This is further supported by Fig. 15f, where PONDR ® VSL2generated disorder profiles of these three proteins are overlapped to show almost complete coincidence of their major disorder-related features. According to the mean disorder analysis (see Fig. 15a-c), Nsp3 proteins are predicted to have the following IDPRs: SARS-CoV-2 (1-5, 105-199, 1221-1238), human SARS (102-189, 355-384, 1195-1223) and bat CoV (107-182, 352-376, 1191-1217). The first 112 residues represent a ubiquitin-like globular fold, while 113-183 residues form the flexible acidic domain rich in glutamic acid. It is thought to bind and ubiquitinate viral E protein using the N-terminal acidic domain [125,126]. This unstructured segment has many MoRFs predicted by ANCHOR and MoRFPred servers which may facilitate the protein-protein interaction (Table 2). Interestingly, Nsp3 of all three viruses is found to have the highest number of RNA-binding residues (Tables S9-S11).

Non-structural protein 4 (Nsp4)
Nsp4 is reported to induce the formation of DMVs for optimal replication inside host cells [127][128][129]. Although no crystal or NMR solution structure is reported, Nsp4 is demonstrated to contain a tetra-spanning transmembrane region with its N-and C-terminals present in cytosol [130].
Nsp4 protein of SARS-CoV-2 has multiple substitutions near the N-terminal region and has a quite conserved C-terminus (Fig. S11). It is found to be closer to Nsp4 of bat CoV (81.40% identity) than to human SARS Nsp4 (80%). The low level of intrinsic disorder illustrated in Fig. S12A-C and mean PPIDs of Nsp4 proteins (Table 3) classify it as a highly structured protein which, however, contains some flexible regions. Likewise, only N-and C-terminal MoRFs which possibly assist in its cleavage from long polyproteins 1a and 1ab are shown in Table 2.
Nsp5 protein is found to be highly conserved in all three studied CoVs. SARS-CoV-2 Nsp5 shares a 96.08% sequence of ADP ribose phosphatase of Nsp3 [residues 207-374 (orange colour)] of SARS CoV-2. f Aligned disorder profiles generated for all three Nsp3 is based on the outputs of the PONDR ® VSL2. Colour schemes are similar to those given in Fig. 3 identity with human SARS Nsp5 and 95.42% with bat CoV Nsp5 (Fig. S13). Therefore, it is not surprising that our analysis demonstrated the identical mean PPID values of 1.96% for all three Nsp5s (Table 3). As the graphs (Fig. 16a-c) depict, Nsp5s have several flexible regions and N-terminally IDPR of six residues. Due to the low flexibility of this protein, a single MoRF predicted by MoRFchibi_web is present in the N-terminal region (residues 3-8) in all Nsp5s (Tables 2, S7, S8). Further, the identified nucleotide-binding residues in Nsp5 proteins are tabulated in Tables S9-S11.

Non-structural protein 6 (Nsp6)
Nsp6 protein is involved in blocking ER-induced autophagosome/autolysosome vesicle formation that functions in restricting viral production inside host cells. It induces autophagy by activating the omegasome pathway, which is normally utilized by cells in response to starvation [135].
Nsp6 of SARS-CoV-2 is equally close to Nsp6s from both human SARS and bat CoV, having a sequence identity of 87.24% (Fig. S14D). Similarly, mean PPIDs for all three Nsp6 proteins is calculated to be 1.03%. The graphs in Fig.  S14A-C further illustrates its highly structured nature. As Nsp6 is a membrane protein, all three proteins are predicted to have a single MoRF near the N-terminal region (residues 1-19 in SARS-CoV-2, residues 1-22 in human SARS, and residues 1-21 in bat CoV) by the DISOPRED3 server.
The role of this protein-binding region for the induction of autophagy needs to be elucidated.

Non-structural protein 9 (Nsp9)
Nsp9 protein is a single-stranded RNA-binding protein [140]. It might protect RNA from nucleases by binding and stabilizing viral nucleic acids during replication or transcription [140]. Our results on nucleotide-binding tendency of Nsp9 shows the presence of several RNA-binding and few DNA-binding residues in Nsp9 of SARS-CoV-2, Human SARS, and bat CoV (Tables S9-S11). Presumed to evolve from a protease, Nsp9 forms a dimer using its GXXXG motif [141,142]. Figure 19d shows a 2.7 Å crystal structure of human SARS Nsp9 homodimer (PDB ID: 1QZ8) that identified a unique and previously unreported oligosaccharide/oligonucleotide fold-like fold [140]. Here, each monomer contains a cone-shaped β-barrel and a C-terminal α-helix arranged into a compact domain [140].
Nsp9 of SARS-CoV-2 is equally similar to other two Nsp9 proteins (with a percentage identity of 97.35%). The difference in the three amino acids at 34, 35 and 48 positions accounts for its similarity (Fig. 19e). Mean PPIDs of all Nsp9s are listed in Table 3. Graphs in Fig. 19a-c show that all three Nsp9s are rather structured, but contain flexible regions. It contains conserved residues (R10, K52, Y53, R55, R74, F75, K86, Y87, F90, K92, R99, and R111) of positively charged side chains suitable for binding with the negatively charged phosphate backbone of RNA and aromatic side-chain amino acids providing stacking interactions [140]. These residues are a part of multiple disorder-based binding sites predicted by MoRFchibi_webserver (Tables 2,  S7, S8).

Non-structural protein 10 (Nsp10)
Nsp10 forms a complex with Nsp14 for hydrolysing dsRNA in 3′-5′ direction [143]. In addition to activating the exonuclease activity of Nsp14, it also stimulates its methyltransferase (MTase) activity required during RNA-cap formation after replication [144]. Figure 20d represents the X-ray crystal structure of the Nsp10/Nsp14 complex (PDB ID: 5C8T) [145]. In agreement with the results of previous biochemical experimental studies, the structure identified important interactions with the ExoN (exonuclease domain) of Nsp14 without affecting its N7-MTase activity [143,144].
SARS-CoV-2 Nsp10 protein is quite conserved having a 97.12% sequence identity with Nsp10 of human SARS and 97.84% with Nsp10 of bat CoV (Fig. 20e). Mean PPIDs of all three studied Nsp10 proteins are found to be 5.04%. Figure 20a-c represents the disorder profiles of Nsp10s and signifies the lack of long IDPRs. Furthermore, Tables 2, S7 and S8 shows that all three Nps10 proteins have multiple MoRFs. For SARS-CoV-2, three MoRFs (residues 25-32, 91-99, and 133-138) were identified by MoRFchibi_web server and one MoRF (residues 11-18) was predicted by MoRFPred server. Interestingly, the SARS-CoV Nsp10 residues F16, F19, and V21 form van der Waals interactions with many of the Nsp14 amino acids [145] out of which one residue (F16) is located in the MoRF region identified in this study. Furthermore, many nucleotide-binding residues which are found in all three Nsp10s are listed in Tables S9-S11.

Non-structural protein 12 (Nsp12)
In coronaviruses, Nsp12 acts an RNA-dependent RNA polymerase (RDRP). It accomplishes both primer-independent and primer-dependent synthesis of viral RNA   In this structure, A and C chains (cornflower blue colour) signifies 1-131 residues of Nsp10, while B and D chains corresponds to residues 1-453 and 465-525 of Nsp14 (dim grey colour). e MSA profile of all three Nsp10 proteins. Colour schemes are similar to those given in Fig. 3 with Mn 2+ as its metallic co-factor and viral Nsp7 and 8 as protein co-factors [146]. As aforementioned, a 3.1 Å resolution structure of human SARS Nsp12 in association with Nsp7 and Nsp8 proteins (PDB ID: 6NUR) has been reported using electron microscopy (Fig. 17e). Nsp12 has a polymerase domain similar to "right hand" containing finger subdomain (398-581, 628-687 residues), palm subdomain (582-627, 688-815 residues) and a thumb subdomain (816-919) [139].
SARS-CoV-2 Nsp12 protein has a highly conserved C-terminal region (Fig. S16). It is found to share a 96.35% sequence identity with human SARS Nsp12 and 95.60% with bat CoV Nsp12. Mean PPID values for all three Nsp12s are estimated to be 0.43% (Table 3). Graphs in Fig.  S15A-C show that although Nsp12s are mostly ordered, they have multiple flexible regions. As RDRP protein is observed to be mostly structured, significant MoRFs in disordered regions are not found (Tables 2, S7, S8).
Nsp13 of SARS-CoV-2 is found to be almost conserved as it shares 99.83% with Nsp13 of human SARS and 98.84% with Nsp13 of bat CoV (Fig. S17). Accordingly, mean PPIDs of all three Nsp13 proteins are estimated to be 0.67%. Graphs in Fig. 21a-c show that Nsp13s contain multiple flexible regions but does not possess significant disorder. As expected, being ordered proteins, Nsp13s does not contain any MoRF (Tables 2, S7, S8), but has several nucleotidebinding residues (Tables S9-S11).
SARS-CoV-2 Nsp14 protein shares a 95.07% identity with human SARS Nsp14 and 94.69% with bat CoV Nsp14 (Fig. S18). Low mean PPID values for all three Nsp14s (Table 3) and disorder profiles depicted in Fig. S19A-C shows its highly structured nature. Likewise, all three Nsp14 proteins contains two protein binding regions (residues 8-13 and 441-445) predicted by the MoRFPred.
SARS-CoV-2 Nsp15 shares 88.73% sequence identity with human SARS and 88.15% with bat CoV (Fig. S20). The calculated mean PPIDs of Nsp15s from SARS-CoV-2, human SARS, and bat CoV are 1.73%, 2.60%, and 2.60%, respectively. Similar to many other Nsps, all three Nsp15 proteins are predicted to possess multiple flexible regions but contain virtually no IDPRs (see Fig. 22a-c). Also, no significant disorder-binding regions are predicted in Nsp15 proteins ( Table 2). SARS-CoV-2 and bat CoV Nsp15s possesses very short binding regions, while human SARS Nsp15 does not contain any MoRF (Tables S7, S8). Tables S9-S11 depict the presence of many RNA-binding residues and few DNA-binding residues in Nsp15 of all three viruses.
Nsp16 of SARS-CoV-2 is found to be identical with other two Nsp16 proteins (93.29%) (Fig. S21). Mean PPIDs for Nsp16s from SARS-CoV-2, human SARS, and bat CoV are 5.37%, 3.02%, and 3.02%, respectively. In line with these PPID values, graphs in Fig. 23a-c show that these proteins are mostly ordered having several flexible regions. Correspondingly, only a single MoRF (residues 151-156) is present in all three Nsp16s. Further, several RNA-binding and few DNA-binding residues are also identified (Tables  S9-S11).

Replicase polyprotein 1a
Since replicase polyprotein 1a contains non-structural proteins 1-10 identical to those found in replicase polyprotein 1ab, we did not perform their disorder analysis separately. However, replicase polyprotein 1a has one additional nonstructural protein designated as Nsp11.

Non-structural protein 11 (Nsp11)
Nsp11 is a small uncharacterized protein with unknown function and requires extensive experimental insights to reveal its structural indentity. The intrinsic disorder-predicting software used in this study requires amino acid sequences which are at least 30-residue long. Therefore, because of their short sequences (just 13 residues), Nsp11s from all three studied coronaviruses are not checked for the intrinsic disorder, disorder-based protein-binding regions, and nucleotide-binding residues. Based on the MSA outputs, Nsp11 from SARS-CoV-2 is found to have a sequence Fig. 23 Analysis of intrinsic disorder propensity of Nsp16. Graphs ac represent the intrinsic disorder profiles of Nsp9 protein of a SARS-CoV-2, b human SARS, and c bat CoV. Colour schemes are similar to those given in Fig. 3. d A 1.95 Å resolution crystal structure (PDB ID: 6W75) of the Nsp10-Nsp16 complex of SARS-CoV-2. Nsp16 of residues 2-298 is represented using pink colour, while Nsp10 of residues 18-139 is shown in cornflower blue colour. e A 2.60 Å crystal structure (PDB ID: 3R24) of human SARS Nsp10-Nsp16 complex. Chain A shown in turquoise colour corresponds to residues 3-294 of Nsp16 identity of 84.62% with Nsp11s from human SARS and bat CoV (Fig. S22).

Concluding remarks
Emergence of new viruses and associated deaths around the globe represent one of the major concerns of modern times. Despite its pandemic nature, there is very little information available in the public domain regarding the structures and functions of SARS-CoV-2 proteins. Based on its similarity with human SARS CoV and bat CoV, the published reports have suggested the functions of SARS-CoV-2 proteins. In this study, we utilized information available on SARS-CoV-2 genome as well as translated proteome from GenBank, and carried out a comprehensive computational analysis of the prevalence of intrinsic disorder in SARS-CoV-2 proteins. Additionally, a comparison is also made with proteins from close relatives of SARS-CoV-2 from the same group of beta coronaviruses, human SARS CoV and bat CoV. Our analysis revealed that in these three CoVs, the N proteins are highly disordered, possessing the PPID values of more than 60%. These viruses also have several moderately disordered proteins, such as Nsp8, ORF6, and ORF9b. Although other proteins have shown lower disorder content, almost all of them contain at least one IDPR. Importantly, our study provides novel information on the presence of intrinsic disorder at the cleavage sites of replicase polyprotein 1ab of SARS CoVs. This observation confirms the crucial role of IDPRs in maturation of individual proteins. We also established that many of these proteins contain disorder-based binding motifs. Since IDPs/IDPRs might undergo structural transition upon association with their physiological partners, our study generates important grounds for better understanding of the functionality of these proteins, their interactions with other viral proteins, as well as interaction with host proteins in different physiological conditions.

Future perspective
The periodical outbreaks of pathogens worldwide always pinpoint the lack of suitable drugs or vaccines for proper cure or treatment. In 2003, nearly 750 deaths were reported due to the SARS outbreak in more than 24 countries. But this time, the outbreak of Wuhan's novel coronavirus (SARS-CoV-2) has quickly surpassed this number, indicating more casualities soon. The lack of accurate information and ignorance of primary symptoms are major reasons, which cause many infection cases. Although efficient transmission from human to human has been confirmed, the actual reasons for fast SARS-CoV-2 spread are still unknown, but some assumptions are made by researchers and the Chinese authorities. The fast spread of SARS-CoV-2, COVID-19 pandemic, and associated introduction of quarantine also have made major impacts on economy and education worldwide due to several restrictions, such as limited transportation, restrained or frozen travel, halted attendance of mass events and the introduction of distant teaching and learning. Due to advancements in sequencing techniques, the full genomic sequence of SARS-CoV-2 was made available in a few days of the first infection report from Wuhan, China. However, massive subsequent research needs to be done to identify the actual cause of SARS-CoV-2 infectivity and to design suitable treatment in the future. Certain possibilities can be explored with the available information. The mutational pressure study on this virus will be very interesting to see if this virus transforms from bat SARS to human SARS to SARS-CoV-2. More in-depth experimental studies using molecular and cell biology techniques to establish structure-function relationships are required for a better understanding of the functioning of SARS-CoV-2 proteins. Additionally, based on the sequence homology and information on protein-protein interactions, the associated viral and host proteins should be explored, for finding means suitable for limiting replication, maturation, and ultimately pathogenesis of this virus. Although structural biology techniques (so-called rational drug design) can be used in drug development utilizing high-throughput screening of compounds virtually or experimentally, the applicability of these techniques is limited by the presence of intrinsic disorder in target proteins. Therefore, the thorough disorder analysis of three coronaviruses conducted in this study will help structural biologists to rationally design experiments keeping this information in mind.