Abstract

Besides being a common threat to farm animals and poultry, coronavirus (CoV) was responsible for the human severe acute respiratory syndrome (SARS) epidemic in 2002–4. However, many aspects of CoV behavior, including modes of its transmission, are yet to be fully understood. We show that the amount and the peculiarities of distribution of the protein intrinsic disorder in the viral shell can be used for the efficient analysis of the behavior and transmission modes of CoV. The proposed model allows categorization of the various CoVs by the peculiarities of disorder distribution in their membrane (M) and nucleocapsid (N). This categorization enables quick identification of viruses with similar behaviors in transmission, regardless of genetic proximity. Based on this analysis, an empirical model for predicting the viral transmission behavior is developed. This model is able to explain some behavioral aspects of important coronaviruses that previously were not fully understood. The new predictor can be a useful tool for better epidemiological, clinical, and structural understanding of behavior of both newly emerging viruses and viruses that have been known for a long time. A potentially new vaccine strategy could involve searches for viral strains that are characterized by the evolutionary misfit between the peculiarities of the disorder distribution in their shells and their behavior.

1. Introduction

1.1. Protein Intrinsic Disorder and Viral Behavior

Previously, we provided evidence that the behavior of viruses can be predicted from the analysis of their predicted intrinsic disorder in their protein shells, more specifically, by looking at the peculiarities of disorder distribution in their matrix and capsid proteins [13]. For example, the predicted disorder in the matrix of retroviruses was shown to vary with the mode of the viral transmission. The HIV and EIAV viruses, that are related but have distinctly different modes of transmission, were used to illustrate this point since. HIV is largely transmitted sexually, whereas EIAV is transmitted by a blood-sucking horsefly. It has been observed that the abundance of predicted intrinsic disorder (PID) in the HIV and EIAV matrix proteins was very different, with the HIV proteins being highly disordered, especially HIV-1. An explanation for this has to do with the need for a more rigid encasement in viruses that are not sexually transmitted, so as to protect the virion from harsher environmental factors [1].

1.2. Goals

Further development of a model that could predict how a virus will behave in terms of transmission would be extremely useful for both clinical and fundamental research. Such a model will also provide a tool to assist the implementation of public health policies for handling old and newly emerging pathogenic viruses. This paper extends the line of research on protein intrinsic disorder in viral proteins to coronaviruses (CoV), which have caught the attention of the scientific community because of the sudden appearance of the lethal virus causing severe acute respiratory syndrome, the SARS-CoV [4, 5]. Clinical, structural, and epidemiological data are available for SARS-CoV and its animal cousins, which remain to be a serious threat to farming communities [4, 610].

One goal of this research is to use of the concept of protein intrinsic disorder to shed light on behaviors of coronaviruses by creating a predictive model, which could provide insight into the differences between the transmission behavior of animal and human coronaviruses and also classify the various animal coronaviruses by their spread behavior. In this way, greater understanding of the viral evolution based on its hosts and its environment can be achieved along with the better understanding of the structural mechanisms involved in such adaptive evolution.

1.3. Coronaviruses

For a long time, coronaviruses have been known to cause respiratory diseases and gastroenteritis [4]. Greater attention was given to this family of viruses when the SARS-CoV moved into human hosts beginning in 2002 and inflicted 1091 deaths [11, 12]. During the outbreak of SARS-CoV, questions regarding its modes of transmission were raised [7, 13]. Greater understanding of the peculiarities of the transmission mode(s) would have enabled better decisions to control the spread of this ominous virus. However, as of today, complete understanding of the molecular mechanisms determining the transmission behavior of the virus remains elusive [6]. For example, there is no satisfactory explanation for the observation that many of those infected with SARS-CoV did not have any contact with the infected individuals [14].

While SARS-CoV had been observed to spread among human most easily by respiratory means, the SARS-CoV was also observed to be not as infectious as Influenza (http://www.who.int/csr/sars/sarsfaq/en/) [5]. It is likely that SARS-CoV, like most animal coronaviruses, spreads most efficiently by contact or by oral-fecal routes among animals.

Several human coronaviruses (HCoV) such as 229E, OC43, HKU1, and NL63 have been known for some time [15, 16], with HCoV-229E being the most studied HCoV. Comparison between the spread mechanisms of HCoV and animal coronaviruses such as SARS-CoV has been made. Using HCoV-229E as a representative virus, it has been shown that HCoV tends to spread more efficiently by a respiratory route, whereas the animal CoVs tend to spread most efficiently by direct contact [15]. Furthermore, gastroenteritis is a common disease provoked by coronaviruses, and oral-fecal route is also common mode of transmission, especially among animal hosts. Therefore, in contrast to human, contact and oral-fecal routes are generally the most efficient mode of CoV spread in animals [15, 17, 18].

1.4. Structural Proteins

In order to introduce the two structural proteins investigated in this paper, a quick overview of virus structure needs to be put forth. Proteins near or at the surface of the coronavirus include the spike protein (S, 150 kD), the hemagglutinin-esterase protein (HE, 65 kD), membrane/matrix glycoprotein (M, 25 kDa), and the small envelope glycoprotein (E, 9–12 kD). A protein that is closer to the RNA is the nucleocapsid protein (N, 60 kD). They are all likely to play roles in protecting the virions [1, 2, 19, 20]. Since the M glycoproteins are the most abundant structural proteins in the coronavirus, they are likely to play a greater role in protecting the virion from damage.

The M-protein is a triple-spanning transmembrane protein and has a short aminoterminal ectodomain. Other than its protective function, M-protein also plays a role in the capsid self-assembly and serves as a major determinant of the virion morphogenesis via selecting the S protein for incorporation into virions during viral assembly. Another protein of interest to this paper is the N-protein [5, 12]. The N-protein is an RNA binding protein. While many of its functions remain unknown, it is likely that the N-protein is involved in packaging and protecting the viral RNA and also plays a role in viral replication participating in transcription [7, 8].

1.5. Grouping of Coronaviruses

Although human coronaviruses, with the exception of SARS-CoV, are generally mild, their animal cousins (e.g., the Avian, Porcine, and Bovine coronaviruses) are often devastating to the farming industry. In fact, an outbreak could be costly in terms of the livestock loss [10, 21]. Coronaviruses are generally classified into three groups based on their genetic and antigenic makeup [12, 16]. A summary of this classification is shown in Table 1. The SAR-CoV did not fit into any of the current three groups but did have some resemblance to Group 2. For this reason, SARS-CoV fell into a new category, group 2b.

Previously, we showed that the similarity in protein intrinsic disorder prediction does not necessarily reflect the genetic proximity of proteins studied, but rather could be related to the evolution, the modes of transmission, and the environment that the virus lives in [1, 2]. Therefore, the categorization of coronaviruses based on the intrinsic disorder propensities of their proteins might provide important clues for useful hypotheses, especially when intrinsic disorder of the similar proteins is considered across coronaviruses, serotypes, or subtypes.

1.6. Protein Intrinsic Disorder

The main tools used here are the predictors of protein intrinsic disorder. Intrinsically disordered proteins have also been described by various other names such as “intrinsically unstructured” [22] and “natively unfolded” [23, 24]. The generality of the intrinsic disorder concept emerged from the occasional observations of exceptions to the structure-to-function paradigm according to which unique protein structures are necessary for specific protein functions [25]. Comparison of ordered and disordered proteins revealed that there is a noticeable difference between the amino acid sequences of ordered and intrinsically disordered proteins and that disordered proteins/regions share at least some common sequence features over many proteins [24, 2628]. In fact, the disordered proteins/regions were shown to be significantly depleted in bulky hydrophobic (Ile, Leu, and Val) and aromatic amino acid residues (Trp, Tyr, and Phe), which would normally form the hydrophobic core of a folded globular protein, and also possess low content of Cys and Asn residues. The depletion of disordered protein in Cys is also crucial as this amino acid residue is known to have a significant contribution to the protein conformation stability via the disulfide bond formation or being involved in coordination of different prosthetic groups. These depleted residues, Trp, Tyr, Phe, Ile, Leu, Val, Cys, and Asn, were proposed to be called order-promoting amino acids. On the other hand, ID proteins were shown to be substantially enriched in Ala, polar, disorder-promoting amino acids: Arg, Gly, Gln, Ser, Glu, and Lys and also in the hydrophobic, but structure-breaking Pro [2731]. Development of such predictors provided a direct support for the hypothesis that intrinsic disorder is encoded in protein amino acid sequences. These features made intrinsically disordered proteins/regions recognizable and were used to develop specific predictors of intrinsic disorder. Currently, design of algorithms for finding regions lacking ordered structure is a very active area of research, and more than 50 predictors of disorder have been developed [32]. The predictor used in this paper is PONDR VLXT (Predictors of Naturally Disordered Regions), which is a set of neural network predictors of disordered regions on the basis of local amino acid composition, flexibility, hydropathy, and other factors [3336]. These predictors classify each residue within a sequence as either ordered or disordered. Since PONDR VLXT is sensitive to local sequence peculiarities, it is frequently used for identifying functionally important sites within the disordered regions.

2. Results

2.1. General Trends of Viral Predicted Disorder

Figure 1 represents a graphical comparison of the percent of intrinsic disorder (PID) of the M- and N-proteins of coronaviruses with that of the influenza A virus and the matrix and nucleocapsid proteins of RNA viruses in general. The mean PIDs of matrix and nucleocapsid proteins of the influenza A virus and coronaviruses provide an interesting contrast since both viruses have similarities and dissimilarities. For example, the amounts of disorder in the N-proteins of both viruses are rather similar, whereas the M-proteins of coronaviruses are predicted to be more ordered. Since the major function of M-proteins is to protect the virion, it is tempting to hypothesize that these differences in the overall disorder of M-proteins can be related to the need to protect viruses from different environments, and therefore can reflect differences in the viral transmission mode. We know, for example, that both viruses are often spread via droplet transmission, but animal coronaviruses are often more associated with oral-fecal transmission routes. The differences in the PIDs may reflect this trend.

A more detailed comparison between the PID in M- and N-proteins of the HCoVs and those of animal coronaviruses in terms of predicted disorder is shown in Figure 2. Avian and SARS coronaviruses are also included for comparison. Figure 2 shows that the N-protein of the avian coronavirus, IBV, is characterized by the highest PID level. The high PID variance (standard deviation) of the HCoV M- and N-proteins should be also noted. Both features likely hint to a higher respiratory spread component.

2.2. Disorder of Nucleocapsid and Matrix Proteins of Coronaviruses Varies with the Type of the Human and Animal Hosts

While Figure 2 allows us to compare human coronaviruses (HCoV) with nonhuman ones in general, Figure 3 breaks down the data of animal coronaviruses further by animal hosts. We are able to see differences in PID of both matrix and nucleocapsid proteins of viruses by the various animal host species.

Figure 4 provides an overview of the mean PIDs in the M- and N-proteins of various animal coronaviruses by species of the hosts and by the type of virus. The differences in PIDs between host types should be noted. For instance, in comparison with its porcine counterparts, bovine coronaviruses tend to have higher PIDs in their N-proteins. The important note here is that the PIDs of the N-proteins for all these animal coronaviruses are below 50%, unlike human or avian coronaviruses.

Table 2 provides numerical PID values for the N- and M-proteins of various porcine strains shown in Figure 4. A close examination of these data revealed that the particular strains may be quite different from each other by the amount of intrinsic disorder they possess in their N- and M-proteins. An illustrative example is the remarkable differences in the PIDs of the N-proteins from TEGV and PEDV. Here, the CV777 strain of the PEDV has lower disorder levels in its matrix protein than that of the Br1/87 strain (8% versus 13%). Using SARS-CoV as a reference point, we could clearly see that each strain of the virus is characterized by a unique PID signature. For example, TGEV is characterized by the relatively low PID score for its N-protein (~43%), whereas its M-protein has a somewhat higher PID (~14%) especially when compared to certain strains of PEDV (~8%).

2.3. Human Coronaviruses (HCoV)

A breakdown of the HCoV PIDs by strains is shown in Figure 5 which shows the unusually high and low PIDs for the M- and N-proteins of the HCoVs 229E and HKU1. The noticeable differences in the PID values of HCoV 229E and SARS-CoV should also be noted. Some HCoV strains were studied clinically and the peculiarities of their transmission modes are listed in Table 1. The HCoV-229E has been studied more extensively than other human coronaviruses and has been shown to spread by respiratory modes more easily. This is also believed to be generally the case for all the HCoVs. The results shown in Figure 5 illustrate that, in general, the N- and M-proteins of human coronaviruses tend to be more disordered than those of the animal coronaviruses (cf. Figures 13).

Figure 5 also shows that both matrix and nucleocapsid proteins of the HCoV-229E are exceptionally disordered. Figure 5 also illustrates an important point that not all HCoV strains have the same characteristics as 299E and there is great variability in the PIDs of the matrix and nucleocapsid proteins of HCoVs. For example, the predicted disorder in HKU1 for some reason resembles the mean PIDs we would normally find in the animal coronaviruses.

2.4. Categorization of Coronaviruses by Shell Disorder and Transmission Behavior

The patterns of predicted disorder seen in Figures 15 allowed us to regroup the various coronaviruses by their PID levels. Table 3 summarizes the new grouping of coronaviruses. It could be seen that although the new grouping of CoVs shown in Table 3 is based on the disorder analysis of the viral shells, and in particular their N-proteins, the grouping seems to reflect how the viruses spread. In other words, a correlation between common transmission behavior and the disorder characteristics of the viruses in each group can be observed. We see, for example, that based on their PID analysis, the canine respiratory coronavirus falls in Category B, which contains viruses with moderate levels of respiratory and fecal-oral components. Its less disordered enteric counterpart falls into Category C, which has greater fecal-oral component and less respiratory component. This pattern is also seen for the porcine coronaviruses TGEV and PEDV.

Therefore, data shown in Table 3 can be used for the prediction of the levels of respiratory and oral-fecal transmission modes for each virus in a group based on the PIDs of their M- and N-proteins. For example, we can conclude that viruses in Group A will have greater respiratory component than those in Group B or C. MANOVA analysis reveals statistically significant influence of the PIDs of both N and M proteins on the transmission behaviour of coronaviruses ( for each variable, , ). Furthermore, because of the statistical significances of the model and its variables, the three groups should be easily identifiable via the analysis of the PID levels of the corresponding M and N proteins. However, any correlation between PIDs of the two proteins is not statistically significant (MANOVA, ). An interesting note is that not all the viruses in the Category C have low PID levels of their M-proteins, even if their N-proteins are characterized by low predicted disorder.

While the figures and tables above show the differences in mean PID levels between different CoV proteins, Figure 6 visualizes the differences in the predicted disorder levels by comparing the 3D structures of parts of the N-proteins of SARS-CoV and IBV. Large disordered regions (shown in red) can be seen in the IBV nucleocapsid, whereas the SARS-CoV nucleocapsid possesses much smaller amount of disordered regions. A rough disorder analysis revealed that the PID of the SARS-CoV nucleocapsid is 50%, whereas in the case of IBV nucleocapsid, the mean PID is of 56%. Figure 6 suggests that this difference may be reflected on their structural differences.

3. Discussion

3.1. Predicted Disorder Varies with the Evolution and Transmission Mode of the Virus
3.1.1. Protective Functions of the Viral Shells: Stepwise Disorder towards the Core

Data shown in Figure 1 suggest that the protective functions of the viral surface proteins can be seen from the results of protein disorder predictions. In fact, there is a correlation between the PID levels and the localization of surface proteins, with proteins located closer to the virion surface being generally more rigid than other surface proteins (e.g., in coronaviruses, nucleocapsid proteins possess noticeably higher PID than the matrix proteins). Generally, the stepwise decrease in the PIDs is seen for all viral nucleocapsid, capsid, and matrix proteins analyzed so far.

3.1.2. Disorder in the Nucleocapsid Associates with Infectivity Pathways Especially for Viruses Transmitted via the Respiratory Mode

Figure 1 shows that there are similarities and dissimilarities between the PID levels in the influenza and coronavirus shell proteins. Disorder peculiarities can be potentially mapped to the transmission behavior and the nature of the viruses. It can be seen that in comparison with other viruses, the coronaviruses possess a tendency to be more ordered at the level of the matrix proteins, but their nucleocapsid proteins are generally more disordered than nucleocapsid proteins RNA viruses in general. The fact that both influenza viruses and coronaviruses have high predicted disorder in their nucleocapsid proteins suggests that higher nucleocapsid disorder may be necessary for viruses with ability to spread via respiratory means. Further support for this hypothesis will be presented below as we inspect and categorize the various coronaviruses. Past analysis showed that many viruses that spread by respiratory modes are characterized by disordered shell proteins [3].

3.1.3. Higher Shell Rigidity, Higher Levels of the Oral-Fecal Transmission Mode

Viral shell proteins constitute the protective proteinaceous layer which is needed for virus to survive during its transmission between the hosts. Figure 1 shows that the mean PID levels of M-proteins in coronaviruses are strikingly lower than those of influenza A. On the other hand, coronaviruses are often transmitted via the oral-fecal and fecal-contact modes [10, 13, 15, 18], whereas the major transmission mechanism of the influenza viruses is the respiratory mode. Based on these and similar observations we hypothesize that there is a correlation between the PID levels in the viral shell proteins and the virus transmission modes. While higher levels of disorder in shell proteins are likely to be associated with greater probability of respiratory transmission, lower PID levels in the nucleocapsid and matrix proteins are conversely often associated with increased levels of fecal-oral transmission mode. These correlations can be due to the need of a given virus to adjust to the changes in the environmental conditions associated with the process of transmission between the hosts.

Obviously, the environmental conditions associated with the respiratory transmission are less harsh than those seeing in the fecal-oral transmission mode. However, the viruses that are spread via the respiratory mode might experience greater environmental pressure during transmission and they have to be able to adjust to a greater range of environmental changes to survive and to be successfully transmitted. This ability of the respiratory transmitted viruses to adjust and survive in a wide range of conditions can be due to the higher levels of intrinsic disorder in their shell proteins. In other words, shell proteins of these viruses are assembled into a flexible shield, the pliability of which helps viruses to rapidly adjust to the environmental changes. Table 3 supports this hypothesis by showing that for several viruses analyzed, there is a noticeable correlation between the level of disorder in the shell proteins and the transmission mode.

3.2. Patterns of Protein Intrinsic Disorder Might Determine the Peculiarities of the Animal Infection by Coronaviruses

Disorder patterns are linked to the behavior of the viruses and their hosts. Marked differences in disorder of shell proteins were predicted among human, avian, and other animal coronaviruses. There are also noticeable differences in the disorder propensities among the coronaviruses infecting different animal (Figures 26, Table 2). For example, the PID levels in the matrix and nucleocapsid proteins of the porcine and canine coronaviruses are generally lower than those of the bovine counterparts (Figure 3). A plausible explanation lies in the nature of the hosts of the various viruses and in the mode of the virus transmission. Furthermore, since the habits and diets of the various animals can also contribute to the tendency of them to be exposed to fecal material, that are crucial factors determining the peculiarities of the virus transmission. Both the porcine and canine coronaviruses have greater fecal-oral components, whereas the bovine coronaviruses are more efficiently transmitted via the respiratory routes. All these factors are apparently more definitive determinants of the intrinsic disorder levels in the viral shell proteins than the genetic proximities and the categorization of the viruses seen in Table 1.

3.2.1. Differences in Intrinsic Disorder Seen in Different Coronaviruses Infecting the Same Host Are Potentially Linked to the Types of the Infected Organs

Some variations in the levels of predicted intrinsic disorder might depend not only on the transmission mode but also on the nature of organ affected by the virus during infection. Figure 4 shows that the canine respiratory coronavirus (that affects the respiratory tract) has higher PID levels than the canine enteric coronavirus (preferentially affecting intestines). Similar patterns can also be seen among the porcine coronaviruses (see Table 3).

3.3. Porcine Coronaviruses: A Path towards Understanding of the Correlation between the Protein Disorder and Viral Behavior
3.3.1. Some Enigmatic Characteristics of Porcine Coronaviruses

The transmissible gastroenteritis coronavirus (TGEV) and porcine epidemic diarrhea virus (PEDV) have similar modes of transmission. Although they are preferentially passed on via the oral mode, both viruses can be transmitted through the lungs too [17, 18], with the PEDV being generally less efficiently transmitted. A puzzling aspect of the PEDV is that this virus has the ability to reemerge to infect hosts months after the breeding areas have been cleaned [17]. This suggests that PEDV is perhaps more persistent in the environment outside the host organism. As it will be seen from the next paragraphs, our analysis suggested that intrinsic disorder can produce useful information for better understanding of the viral behavior.

3.3.2. The Use of the Intrinsic Disorder Predictions to Address Some Viral Behavior Puzzles

Tables 2-3 show that the level of predicted disorder in the TGEV N-protein is noticeably lower than that of the PEDV analogue, and that at least one PEDV strain has a noticeably lower PID in its M-protein (Table 2). Based on the fact that the M-protein of PEDV possesses the greater rigidity we hypothesize that some strains of PEDV are likely to be more persistent outside the physiological environments when compared to TGEV. This hypothesis is also consistent with the observed PEDV ability to survive in the harsh nonphysiological environment (outside the host), to sustain the disinfection of the pig pens and infect new hosts many months after the first infection [17].

3.3.3. Protein Shell Disorder Can Explain Mechanisms of Infection and Persistence

Our analysis also suggests that the observed intrinsic disorder propensities of the PEDV shell proteins may be related to the way of how the virus spreads between hosts. Because of the higher levels of intrinsic disorder predicted in the N-protein by PONDR VLXT, it is expected that there is a relatively large “respiratory” component in the PEDV transmission (see Table 3). On the other hand, TGEV spreads more easily since contact and oral-fecal modes are the more advantageous forms of transmission among swine [17, 18].

3.3.4. SARS-CoV and PEDV Possess Similar Transmission Behaviors and Comparable Levels of Predicted Disorder

The comparison of the PID levels in the M- and N-proteins suggests that various coronaviruses can be grouped as shown in Table 3. An important result of this analysis is the conclusion that in terms of the intrinsic disorder in M- and N-proteins of various coronaviruses, the closest neighbor of the SARS-CoV is PEDV. The structural disorder similarity suggests that the behaviors of these two viruses might also be similar, insights that were not previously noticed.

3.3.5. Evolutionary Misfits Are Potential Targets for a New Vaccine

It should be noted, however, that there are some strains of the animal coronaviruses that could be used as counterexamples of the categorization shown in Table 3. One of such case is porcine respiratory coronavirus (PRCV), which is a mild variant of the transmissible gastroenteritis coronavirus (TGEV) that has been observed to transmit predominantly by respiratory means and can scarcely be detected in feces or small intestines [4, 10]. The PRCV transmission behavior was explained by the specific changes in the viral spike (S protein) which defined the ability of this protein to interact with specialized cells in the respiratory tract, and not with the receptors found in cells from the gastrointestinal region [10]. Since PRCV and TGEV have virtually identical levels of predicted disorder in their M- and N-proteins, both of them should be assigned to the Category C, that is, to the class of coronaviruses that are spread mostly via the fecal-oral and contact modes. Therefore, PRCV that relies on the respiratory transmission mode may be inept by evolution. This potential evolutionary misfit may also explain why the symptoms resulting from the PRCV infection are generally mild and why PRCV is a vaccine strain. The observation may also suggest a new vaccine search strategy based on the recognition of viruses or strains that are evolutionary misfits.

3.4. Disorder Comparison among HCoV, SARS-CoV, and Other Animal Coronaviruses Provides Clues to Viral Behavior
3.4.1. Large Variations Are Seen for Disorder in Human Coronaviruses: Complexity of Human Immune Systems

Data shown in Figures 15 and Table 3 can be used for the comparison of the human matrix and nucleocapsid proteins with the corresponding proteins from the nonhuman coronaviruses. A striking feature is the presence of large PID variability for both the M and N of human coronaviruses especially when compared to those of nonhuman animal coronaviruses. This difference is likely due to the greater complexity of the human immune system [37]. Furthermore, many of the HCoVs are characterized by the relatively high levels of predicted disorder in their N and M proteins (see Table 3 and Figure 5). Since our analysis revealed that higher PID levels in one or both of the capsid proteins typically correspond to the coronaviruses with greater respiratory transmission component, the high disorder levels seen for both M- and N-proteins may indicate that many HCoVs are more easily spread by respiratory means, an important hypothesis which is in agreement with the clinically and experimentally observed transmission behavior of HCoVs [15].

3.4.2. Similarities and Differences between the Avian IBV-CoV and Human HCoV-229E

Figures 26 and Table 3 show that nucleocapsid of the avian coronavirus (IBV) is characterized by the high predicted disorder levels. The only virus that possesses the comparable level of disorder is the human coronavirus 229E (Figures 35). However, HCoV-229E and IBV-CoV have some differences too. In HCoV-229E, both M- and N-proteins are disordered, whereas in IBV-CoV, the nucleocapsid is more disordered but the PID of the M is only slightly above average. This implies that although IBV-CoV and some HCoV are likely to be highly infectious mostly via the respiratory transmission, the IBV virus could also be present in the gastrointestinal region, and larger quantity of this virus can be found in stools.

3.4.3. Greater Disorder in Avian Coronavirus Nucleocapsid: Greater Respiratory Component

Previous studies revealed that predicted intrinsic disorder in the matrix and surface proteins of other avian viruses is often associated with the peculiarities of their evasion of immune system [1, 2, 38]. Here, the ability of the viruses to evade the immune systems may allow them to be spread between the different host species. Another reason for the high disorder in the avian influenza nucleocapsid is related to its predominant transmission via the respiratory modes. Comparison of the PID values in N- and M-proteins of these two viruses showed that despite the comparable disorder levels at the nucleocapsid, the M-protein of IBV-CiV is predicted to be relatively more ordered than the corresponding protein of the influenza virus (Figure 2). This is an indication of the greater likelihood of the presence of the virus in stools, despite its predominant respiratory transmission. In agreement with this hypothesis, IBV-CoV has been reported to be present in feces of infected birds, and respiratory mode is believed to be the main transmission means [21, 39]. This may point to the greater possibility of fecal-respiratory transmission in IBV. Similar interpolations can also be made using the analysis of the various porcine coronaviruses in terms of their shell disorder and their respective transmission behavior.

3.4.4. Characteristic Differences between HCoV and Other Animal Coronaviruses

Based on the experimental and clinical studies on the HCoV-229E and SARS-CoV, it has been concluded that HCoV-229E is easily transmissible by respiratory modes, whereas SARS-CoV is more easily transmitted via contact [15], and these observations were used to generalize the transmission differences between human and other coronaviruses [15]. Our data on protein intrinsic disorder support this conclusion and also provide some additional points for consideration.

3.4.5. Higher Respiratory Component in Many but Not All HCoVs

The fact that HCoVs are generally characterized by high PID levels of their M-proteins suggests that many human coronaviruses may have a higher “respiratory” component in their transmission mode. However, as we have seen above, the variability of predicted disorder in HCoV is very high, with HCoV-229E possessing the highest PID rates among HCoVs and with HCoV-HKU1 being characterized by lowest levels of predicted disorder in its N- and M-proteins. This suggests that it is difficult to expect that all the HCoVs would have the same transmission characteristics as HCoV-229E. In fact, some HCoVs resemble nonhuman coronaviruses. For example, HCoV-HKU1 is quite different from other HCoVs, even though several of them are found in the same group (see Table 1 and Figure 5). Table 3 further emphasizes the presence of a noticeable variability in the PID levels of HCoV surface proteins and also shows that different HCoVs can have noticeable variability in their transmission modes.

3.4.6. HCoV-HKU1 and MHV: An Enigmatic Relationship

It has been generally believed that HCoV-HKU1 is spread by respiratory means like the other HCoVs. However, more recent studies revealed that HCoV-HKU1 may be quite different from the other HCoVs. For example, there is a noticeable difference between the genetic makeups of HCoV-HKU1 and HCoV-OC43, and, genetically, HCoV-HKU1 seems to be close to the MHV. Furthermore, a case was described where HCoV-HKU1 was detected in a patient suffering from hepatitis [40]. Given this observation and the close genetic proximity of HCoV-HKU1 and MHV, it is suspected that HCoV-HKU1 may play a role in the pathogenesis of hepatitis. In agreement with these observations and hypotheses, Table 3 shows that there is a close similarity in the PID levels evaluated for the surface proteins of HCoV-HKU1 and MHV. This finding is rather surprising, even if the genetic similarity of these two CoVs were to be taken into account, since data shown in Tables 13 clearly illustrate that genetic proximity does not necessarily translate into the disorder pattern similarity. This also suggests that HCoV-HKU1 and MHV may also possess similar transmission behavior. This hypothesis is supported by the recent clinical studies which indicated the greater presence of HCoV-HKU1 in stools of patients and the greater association of gastrointestinal illness with HCoV-HKU1 infections when compared to cases of other HCoV infections [41]. Viruses, such as MHV, that infect the liver may need harder shell encasements to protect the virion during their exposure to bile, which is known to inflict damage on certain viruses [42].

3.4.7. Moderate Disorder in SARS-CoV and PEDV-CoV When Compared to HCoV-229E

Based on the results of the disorder analysis, we can hypothesize that there is a great likelihood that HCoV-229E is most easily spread via the respiratory mode. Based on the PID data shown in Table 3 one can expect that SARS-CoV should have a larger fecal-oral component than HCoV-229E and therefore should likely be spread by contact, oral-fecal, and respiratory means. In fact, more moderate PID levels in SARS-CoV and many animal CoVs imply that these viruses are likely to have reasonably high levels of both oral-fecal and respiratory components, which allow them to persist in the environment for a longer time and yet being able to infect via droplets or aerosols. These hypotheses are consistent with the previously described direct observations of the mostly respiratory transmission mode of HCoV-229E [15]. On the other hand, although many animal CoVs were shown to have a larger fecal-oral component, they also can be spread by touch, since a larger fecal-oral component will allow them to persist in the environment for a longer time, or via the fecal-respiratory mode.

3.5. SARS-CoV: Understanding the Transmission Behavior of Coronaviruses via Predicted Disorder
3.5.1. The SARS-CoV Puzzles

One of the difficulties that epidemiologists experienced during the SARS outbreak of 2003-4 was the way the virus spreads. From the way the virus passed between the hosts, it was known that the virus has some “respiratory” components already at the beginning of the outbreak. When it was then shown to be a coronavirus, a special attention was paid to possible oral-fecal routes, given the nature of coronaviruses. In fact, there were instances where the spread involved fecal routes, such as in the case of the outbreak at Amoy Gardens in Hong Kong, which involved a spread via the sewerage system [11, 43]. Despite of what was known from the outbreak at Amoy Gardens, the exact modes of SARS-CoV spread remain unclear as of yet [6]. For example, there were many cases where patients were infected without any known contact with an infected person [14]. Therefore, the precise spread mechanisms remain somewhat a mystery, and there are no appropriate animal models for the SARS-CoV transmission modes.

3.5.2. SARS-CoV and PEDV: Fecal-Respiratory Routes

The disorder analysis of SARS-CoV and other coronaviruses suggests that SARS-CoV falls in Category B (see Table 3), which includes viruses that are likely to have intermediate levels of both respiratory and fecal-oral components. These viruses are expected to have less respiratory component than HCoV-229E and IBV from the Category A (Table 3 and Figures 3, 5, and 6), but greater contribution of the respiratory mode than viruses of the Category C, such as TGEV and MHV. By the same token, the SARS-CoV has a higher oral-fecal component than the former but lower oral-fecal component than the latter. This suggests a greater likelihood of being spread via the fecal-respiratory mode, which affirms what has been long suspected clinically. The role of fecal-respiratory transmission, as suggested by the disorder analysis mentioned above, also explains many of features observed for PEDV, and not observed for its counterpart, TGEV. It also explains similarly in the spreading behaviors of SARS-CoV and PEDV, given that both have the closest shell disorder seen in our sample.

3.5.3. SARS-CoV: Signs of Persistence Seen in the Disorder Analysis Data

Another feature that was a puzzle during the SARS epidemic and remains unclear now is the question on the SARS-CoV persistence outside the physiological environment. Although the ability of SARS-CoV to survive on the surfaces and feces for extended periods is known [6, 43], the level of persistence and the corresponding evolutionary mechanisms have not been fully understood. The results shown in Tables 2-3 suggest a way to address the problem. In the case of PEDV and SARS-CoV, the PID levels of both N- and M-proteins are quite moderate. In PEDV and TGEV, the differences in the transmission behavior and the levels of shell disorder suggest that special attention should be paid to the PID peculiarities of the M-protein. A sign of persistence of SARS-CoV and some specific strains of PEDV (see Tables 2-3) can be found in the relatively low PID levels of their M-proteins (~8%), especially when compared to viruses such as TGEV (~14%) and HCoV-229E (~23%).

3.5.4. Infections without Contact with Infected Hosts

Another enigma is the ability of SARS-CoV to infect hosts without any apparent contact with an infected host [6, 44]. A comprehensive understanding of the mechanisms involved in this ability remains elusive. The differences between PEDV and TGEV in terms of shell disorder and observed spread characteristics offer some clues to this conundrum, since the ability of PEDV to infect pigs remains even after attempted cleaning of pens and isolation of infected pigs [17]. Both shell proteins of both PEDV and SARS-CoV are predicted to be moderately ordered. Therefore, the same reasoning used to predict the environmental persistence and the greater fecal-respiratory transmission component in PEDV can be extended to explain the behavior of SARS-CoV.

3.5.5. Spread Inefficiency of SARS-CoV and PEDV

Another similarity between the SARS-CoV and PEDV is their relative transmission inefficiency. The SARS-CoV did not exhibit higher spread efficiency especially when compared to influenza A, whereas PEDV has been observed to spread less quickly than TGEV. It would then be tempting to conclude that these inefficiencies are for the same reasons. A closer analysis of the corresponding intrinsic disorder data reveals that the reasons for the spread inefficiency are actually very different. Comparisons of the PIDs of N-proteins of SARS-CoV (Tables 2-3), IBV (Figure 6), and HCoV-229E (Figure 5, Table 2) revealed that the intrinsic disorder in the SARS-CoV N-protein is relatively moderate. Conversely, when the N-protein of SARS-CoV is compared to TGEV (Tables 2-3) and many other animal CoVs, the opposite is held true (Figure 2): the PID level in the SARS-CoV N-protein is higher than the corresponding levels in many animal coronaviruses, especially TGEV.

4. Conclusions

4.1. Predicting Viral Behaviors Using a Model Based on the Analysis of the Rigidity of the Shell Proteins

A model to relate disorder at the M- and N-proteins to the behavior of coronaviruses can be built. This model is based on the observation that viruses that are exposed to a less harsh environment, such as those that are transmitted via oral-fecal route, require a more rigid encasement so as to protect their virions from damage. One way to measure the rigidity of the encasement is by looking at the level of predicted intrinsic disorder in the capsid and matrix proteins. From the known behavior and shell disorder patterns of influenza A virus and various porcine coronaviruses, it can be deduced that higher disorder at M-protein can be linked to greater persistence, while higher (or lower) levels of predicted disorder in the nucleocapsid can be related to the levels of respiratory (or fecal-oral) transmission components, respectively. While the disorder levels in these two proteins tend to follow each other, this is not necessarily the case for all coronaviruses.

4.2. Observed and Predicted Evidence of the Respiratory Spread in HCoV-229E and IBV-CoV

Previous studies have already provided evidence that sexually transmitted viruses often have highly disordered M- and N-proteins since such viruses do not need a rigid protective shell but, instead, flexible layers of proteins that could help these viruses evade the host’s immune system. In this paper, a number of coronaviruses were sampled, and further support of the link was found. For instance, many CoVs that are known to spread predominantly by the respiratory mode are characterized by higher disorder levels in their shell proteins. This is especially the case for the HCoV-229E, for which our analysis revealed that both the M- and N-proteins are highly disordered. The data also suggest that many other human coronaviruses may rely on the respiratory route as the main mode of transmission since HCoVs generally have large PID levels in both matrix and nucleocapsid. An especially large respiratory component has been observed for HCoV-229E when compared with animal coronaviruses including the SARS-CoV [15]. Such high level of disorder was not generally detected in the animal CoVs, with the exception of the avian IBV. It is likely that IBV possesses higher ability to spread via fecal-respiratory route given the more ordered M-protein and more disordered nucleocapsid protein.

4.3. Behavioral Similarities of PEDV and SARS-CoV

The results of the disorder analysis of the shell proteins of the porcine coronaviruses not only gave more evidence to link predicted disorder to viral behavior but also showed how the model can be used to predict the CoV behavior. In fact, there is a close similarity between the PID levels in the matrix and nucleocapsid proteins of SARS-CoV and PEDV. Furthermore, the transmission behaviors of these two viruses are strikingly similar (both PEDV and SARS-CoV spread primarily by contact), despite their relatively distant genetic relationship (Table 1). On the other hand, based on the lower nucleocapsid PID found in TGEV, we could expect to see a larger oral-fecal component in the transmission of this virus. This suggests that since TGEV has much larger oral-fecal and contact components in transmission, it can spread much more efficiently than PEDV, the hypothesis that can be deduced from the presented PID analysis. This suggests that the levels of predicted disorder provide a definite way to evaluate the fitness of the virus.

4.4. PEDV and SARS-CoV: Higher Levels of the Aerosol Spread

The model also suggests a mechanism by which PEDV could reemerge to infect herds even after cleaning of the premises [17]. Our analysis suggests that this ability to reemerge arises from the existence of both oral-fecal and respiratory components. Therefore, the aerosol transmission arising from the viral particles left behind may have a crucial role in this reemerging. In its turn, the ability to have two transmission components can be determined by the peculiarities of the predicted intrinsic disorder distribution in the viral shell proteins. PID similarity between PEDV and SARS-CoV suggests that these two viruses might have similar transmission mechanisms and suggests that SARS-CoV might have propensity to reemerge, which can explain how people could be infected even though they had not been known to be in contact with the SARS-CoV-infected persons [14]. This model also suggests that reemerging behavior cannot be expected for TEGV. Overall, the disorder-based analysis constitutes the categorization algorithm shown in Table 3 that can be used to identify the transmission modes of various coronaviruses.

4.5. A New Tool for Understanding Behavior of Existing and Newly Emerging Viruses

Results of our analysis show that disorder predictors can be used to predict transmission behavior of newly emerging strains of viruses. In the cases of epidemic outbreaks, such a technique would be a useful tool for predicting the potential transmission behavior of a new virus. Furthermore, this analysis gives new clues for better understanding the spread behavior of viruses that have been known for years.

5. Materials and Methods

Searches for the appropriate samples were made using the Entrez website. Table 4 represents some of the proteins analyzed in this study. Proteins with known X-ray crystal structures were chosen to ensure consistency. Further details of implementation can be found in our previous papers [1, 2, 38]. The data extracted from PDB were placed into the database that enables further development of the JAVA codes to create 3D representation in Jmol [45]. The relational database was modified to allow import and storage of the data from UNIPROT (http://www.uniprot.org/) [46]. A graphical user interface program was implemented using JAVA to allow easy entry of protein details and accession codes.

The open-source statistical package, R (http://www.r-project.org/), was used to analyze data. Statistical methods used include MANOVA (Multivariate Analysis of Variance).

Conflict of Interests

Authors have no conflict of interests.

Authors’ Contribution

G. Goh proposed the idea of the study, implemented the experiments, carried out the analyses, and drafted the paper. A. Dunker helped to design experiments. V. Uversky helped to design experiments, coordinated the studies, and helped to draft the paper. All authors read and approved the final paper.

Acknowledgment

This work was supported in part by the “Molecular and Cellular Biology” Program of the Russian Academy of Sciences (V. Uversky).