Theoretical investigation on the interactions of microplastics with a SARS-CoV-2 RNA fragment and their potential impacts on viral transport and exposure

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes the coronavirus disease-19 (COVID-19) pandemic spread across the world and remains difficult to control. Environmental pollution and habitat conditions do facilitate SARS-CoV-2 transmission as well as increase the risk of exposure to SARS-CoV-2. The coexistence of microplastics (MPs) with SARS-CoV-2 affects the viral behavior in the indoor and outdoor environment, and it is essential to study the interactions between MPs and SARS-CoV-2 because they both are ubiquitously present in our environment. To determine the mechanisms underlying the impact of MPs on SARS-CoV-2, we used molecular dynamic simulations to investigate the molecular interactions between five MPs and a SARS-CoV-2 RNA fragment at temperatures ranging from 223 to 310 K in vacuum and in water. We furthermore compared the interactions of MPs and SARS-CoV-2 RNA fragment to the performance of SARS-CoV-1 and Hepatitis B virus (HBV) RNA fragments in interacting with the MPs. The interaction affinity between the MPs and the SARS-CoV-2 RNA fragment was found to be greater than the affinity between the MPs and the SARS-CoV-1 or HBV RNA fragments, independent of the environmental media, temperature, and type of MPs. The mechanisms of the interaction between the MPs and the SARS-CoV-2 RNA fragment involved electrostatic and hydrophobic processes, and the interaction affinity was associated with the inherent structural parameters (i.e., molecular volume, polar surface area, and molecular topological index) of the MPs monomers. Although the evidence on the infectious potential of SARS-CoV-2 RNA is not fully understood, humans are exposed to MPs via their lungs, and the strong interaction with the gene materials of SARS-CoV-2 likely affects the exposure of humans to SARS-CoV-2.


Introduction
The global pandemic of the coronavirus disease-19 (COVID-19) has suddenly made us realize that viruses have become important biological pollutants. The outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) not only seriously threatens human health (Topol, 2020;Turner et al., 2021), but also greatly increases environmental stress (Adelodun et al., 2021a(Adelodun et al., , 2021bBedrosian et al., 2021). It is thus essential to understand the environmental fate and the behavioral dynamics of the coronavirus. The SARS-CoV-2 can travel in all environmental compartments like water (Navarro et al., 2021;Sala-Comorera et al., 2021), air (Dubey et al., 2021;Razzini et al., 2020), and soil (Anand et al., 2021;Steffan et al., 2020). A nucleic acid material (DNA or RNA) enclosed in a nucleocapsid protein is referred to as the non-enveloped structure of a virus particle (Müller et al., 2019). This is in contrast to the enveloped structure of a virus particle which contains a biological membrane. An envelope increases viral sensitivity to external physical stressors (pH, heat, dryness, etc.) as biological membranes are relatively fragile structures. Consequently, the SARS-CoV-2 as an enveloped virus is more sensitive to environmental factors than non-enveloped viruses (Achak et al., 2021). Thus, it is reasonable to believe that the non-enveloped structural materials of the SARS-CoV-2 could be more resistant to these inactivation factors and are likely to maintain their stability for a long time. Furthermore, studies on the nucleic acid material of SARS-CoV-2 are used for its detection and control in the environment and even for the implementation of personal health prevention measures.
The plethora of sources that can contribute to the release of MPs into air have been summarized in (Catarino et al., 2018;UNEP, 2016). In addition, the main sources of indoor and outdoor plastic debris released into the air and subject to human inhalation are illustrated by Amato-Lourenço et al. (2020). The indoor concentrations ranged between 1.0 and 60.0 fibers/ m 3 whereas outdoor concentrations were significantly lower as they range between 0.3 and 1.5 fibers/m 3 (Dris et al., 2017). This is important to quantify and realize, because MPs have been reported as carriers or vectors for concurrent pollutants, e.g. metals , organic pollutants , and they exhibit diverse interactive effects (Bhagat et al., 2021;Kim et al., 2017;Sun et al., 2021). In addition, MPs are becoming a novel ecological habitat termed the plastisphere (Zettler et al., 2013), and could facilitate the survival and dissemination of bacterial and fungal pathogens (Moresco et al., 2021), and antibiotic resistance genes . Importantly, plastic pollution could be a secondary pathway for the transmission of human pathogenic viruses (Moresco et al., 2021) via the respiratory exposure route. We focus here on the MPs-SARS-CoV-2 interactions because both the virus as well as the sources of MPs (like fibers from clothes, building materials, household objects, polymer fragments in urban dust) are closely correlated to the presence of human. It was also reported that SARS-CoV-2 remains more stable on plastic surfaces than on stainless steel, glass, and ceramics (Gidari et al., 2021) which has its consequences for the oral and hand contact exposure routes for humans. Amato-Lourenço et al. (2022) found that SARS-CoV-2 aerosols may bind to total suspended particles, such as MPs, and facilitate virus entry into the human body. Moreover, SARS-CoV-2 virus particles have the ability to sorb to the surface of MPs released during washing processes (Belišová et al., 2022). Hence, there is an urgent need to further explore the interactions and mechanisms of MPs and SARS-CoV-2.
Virus stability in the environment is strongly influenced by the size and structure of the virus particle (including the presence or absence of an envelope), the type of genome (DNA or RNA), a transmission route such as faecal-oral and air droplets, the presence of vectors or carriers like the MPs, and the viral concentration of the contamination source. As known, the intrinsic properties such as polymer type of MPs dictate their interaction affinity with other co-contaminants (Fred-Ahmadu et al., 2020;Menéndez-Pedriza and Jaumot, 2020). Besides, many environmental factors can affect the stability of viruses in the environment (Aboubakr et al., 2021;Achak et al., 2021;Paul et al., 2021), in humans (Matson et al., 2020), and on common touch surfaces (Aboubakr et al., 2021). Notably, temperature (Paul et al., 2021) and relative humidity (Zhao et al., 2020) are the two critical factors that determine the fate and transport of coronaviruses given certain environmental conditions. Therefore, searching for some key characteristics that may affect the interaction of the MPs and SARS-CoV-2 is a noteworthy issue.
In silico methods are a promising approach and play a significant role in elucidating the mechanisms of the interactions of pollutants and biomacromolecules Ge et al., 2011). In particular, the molecular simulation method such as molecular dynamics (MD) simulation is a practical in silico method in environmental applications Feng et al., 2022;Sun et al., 2013). In addition, the molecular simulation method has shown to be an effective tool in exploring the interactions between MPs and SARS-CoV-2, and offered theoretical insights into the adsorption/separation and inactivation of carbon nanoparticles with a SARS-CoV-2 RNA fragments (Zhang et al., 2021b). This way in silico methods can not only contribute to minimizing the challenge of time-consuming and labor-intensive virus experiments under high risks of infection, but also to meeting our precautionary demand for options to handle any new versions of the coronavirus that might emerge in the future.
In light of the demands from the exploration of the interaction and mechanism between MPs and SARS-CoV-2, this knowledge gap needs to be addressed. Hence, in this work for the first time MPs were studied theoretically by MD simulation to characterize their interactions with the nonenveloped structural materials of SARS-CoV-2 including a nucleocapsid protein and a SARS-CoV-2 RNA fragment in the water phase and in the vacuum phase (as a reference for the water phase and as an approximation to the gas phase). Two reference viruses, namely SARS-CoV-1 (homologous coronavirus similar to SARS-CoV-2) and Hepatitis B virus (HBV, noncoronavirus dissimilar to SARS-CoV-2) were selected to compare the performance in interacting with MPs. The influence of five different MP types and the temperature as an environmental factor is considered. The objectives of this study were divided in several parts: 1) Comparison of the interactions of the MPs with the nucleocapsid protein and with the viral RNA fragments; 2) Interaction mechanisms between the MPs and viral RNA fragments; and 3) Correlation of the interaction affinity and molecular parameters of MP monomers.

MD simulation
The selected three-dimensional structure models of the SARS-CoV-2 RNA fragment determined by Zhang et al. (2021c), the SARS-CoV-1 RNA fragment determined by Robertson et al. (2005), and the HBV RNA fragment determined by LeBlanc et al. (2021) were used as model compounds for the simulation of the interactions between MPs and the viral RNA fragments. It should be noted that the SARS-CoV-2 RNA fragment is a model molecule of a frameshift stimulation element (FSE) from the SARS-CoV-2 RNA genome (Zhang et al., 2021c). The FSE plays an important role in the virus replication cycle and has emerged as a major drug target (Lan et al., 2022). The selected three-dimensional structure models of the SARS-CoV-2 nucleocapsid protein determined by Kang et al. (2020), the SARS-CoV-1 nucleocapsid protein determined by Huang et al. (2004), and the HBV nucleocapsid protein determined by Böttcher and Nassal (2018) were used as model compounds for the simulation of the interactions between MPs and the viral nucleocapsid protein. The structures of the RNA fragments [PDB ID: 6XRZ (SARS-CoV-2), 1XJR (SARS-CoV-1), 6VAR (HBV)] and nucleocapsid proteins [(PDB ID: 6M3M (SARS-CoV-2), 1SSK (SARS-CoV-1), 6HU7 (HBV)] were obtained from the RCSB Protein Data Bank (Burley et al., 2019).
The polymer chains derived from five plastic monomers were built as model compounds for MPs including polybutene (PB), polyethylene (PE), polypropylene (PP), polystyrene (PS), and polyvinyl chloride (PVC) within the simulation. All the simulations were carried out in a box with three-dimensional boundary conditions. The dimensions of the simulation boxes were a = b = c = 85 Å, a = b = c = 90°. The length of the simulation box in each direction was large enough to enable the interactions between the MP polymer chain and the materials of the viruses. The process of building the MP models refers to the simulation methods developed by Guo et al. (2019) with slight modifications. The MP polymer chains were built and energy minimized using the smart geometry optimization algorithm, which is the combination of steepest descent, conjugate gradient, and quasi-Newton geometry optimization algorithms. Then the optimized polymer chain was randomly packed in rectangular boxes with three-dimensional periodic boundary conditions by Amorphous Cell Construction. For each box, only one polymer chain was added. The amount of PB, PE, PP, PS, and PVC monomer molecules were 200, 600, 500, 200, and 600. The MP-virus systems included one polymer chain, one RNA fragment or one nucleocapsid protein, and either a vacuum layer (83 Å) or a water layer (83 Å). For the water system, 1000 water molecules were incorporated in each unit cell. The smart geometry optimization algorithm was used to minimize the energy of the simulation systems. Then the MD calculations were performed in the canonical ensemble NVT system in which the number of molecules [N], volume [V], and temperature [T] of the system are kept constant at 223, 263, 273, 298, and 310 K. These temperatures represent the range from a low-temperature environment to the temperature of the human body. The universal force field was adopted in the simulation framework. The van der Waals interaction cut-off was 12.5 Å, and the Ewald method (accuracy 0.001 kcal/mol) was used. The simulation was performed for 100 ps which allowed the studied system to reach equilibrium, and each step was 1.0 fs. A Nose thermostat was adopted. All the simulations were performed with the Materials Studio software package (ver. 8.0).

Interaction energy
For the interaction systems, the magnitude of the interaction energy (E int ) is an indication of the magnitude of the driving force towards complexation. A negative value reflects stable adsorption on the plastisphere. E int was calculated by where E MP-virus , E MP , and E virus represent the energies of the complex, the isolated MPs, and the viral RNA fragment or nucleocapsid protein, respectively.

Molecular parameters and linear correlation models
The MP monomers' molecular parameters (Table S1) such as volume of molecule (V M ), polar surface area (PSA), and molecular topological index (MTI) were selected to correlate with E int so as to develop a quantitative relationship between the inherent properties of MPs and E int . The molecular parameters were calculated using Multiwfn 3.8 software Chen, 2012a, 2012b). Correlation of interaction affinity and molecular parameters of the MP monomers was described using a polynomial relationship by performing linear regression models in Sigma Plot, ver. 14.0 (Systat Software Inc., San Jose, CA).

Statistical analysis
Statistically significant differences between test groups were determined by independent t-test and one-way analysis of variance with the Waller-Duncan test post hoc, at a significance level of p < 0.05 (IBM SPSS Statistics for Windows, ver. 23.0, IBM Corp., Armonk, NY). Linear regression analysis at the significant level of p < 0.05 was carried out using the SPSS.

Comparison of interactions of MPs with viral nucleocapsid protein and RNA fragments
To fully understand the interactions between the MPs and the nonenveloped structures of the virus, the interactions of the MPs with the nucleocapsid protein and with the viral RNA fragments were compared after geometry optimization (Fig. 1). As shown in Fig. 1A, for the SARS-CoV-2, the absolute E int values between the MPs and the nucleocapsid protein were significantly lower (p < 0.05) than those between the MPs and the RNA fragment. In contrast, for the HBV (Fig. 1C), the absolute E int values between the MPs and the nucleocapsid proteins were significantly higher (p < 0.05) than those between the MPs and the RNA fragments. For the SARS-CoV-1 (Fig. 1B), the absolute E int values between the MPs and the nucleocapsid proteins were higher than the corresponding values between the MPs and the RNA fragments, but the two groups showed no significant difference (p > 0.05). Moreover, there was no significant difference in the absolute E int values between the interactions of the MPs with the nucleocapsid proteins of the SARS-CoV-2 and the SARS-CoV-1. However, the absolute E int values between the MPs and the nucleocapsid proteins of the HBV were significantly higher than those between the MPs and the nucleocapsid proteins of the SARS-CoV-2 or the SARS-CoV-1 (p < 0.05). In addition, the absolute E int values between the MPs and the RNA fragments of the SARS-CoV-2 were significantly higher than those between the MPs and the nucleocapsid proteins of the SARS-CoV-1 or the HBV (p < 0.05). Moreover, no significant difference in the absolute E int values between the interactions of the MPs with the RNA fragments of the SARS-CoV-1 and the HBV was found.
Generalizing, when comparing the nucleocapsid protein and the RNA fragment, then the MPs exhibited a stronger interaction with the RNA fragment for the SARS-CoV-2, while the MPs exhibited a stronger interaction with the nucleocapsid protein for the HBV. Furthermore, this difference in the interactions was not affected by the type of MP. The plastic types were a bit more discriminative for SARS-CoV-1 and HBV compared to the SARS-CoV-2 that had interactions energies all similar for each type of MPs.

Interaction mechanisms between MPs and viral RNA fragments
To reveal the mechanisms of the interactions of the MPs with the viral RNA fragments, the values of E int as derived from the total energy (E t ), the potential energy (E p ), the 'van der Waals' energy (E v ), and the electrostatic energy (E e ) are summarized in Figs. 2, S1, and S2. As shown in Fig. 2A  temperature range in vacuum and the full temperature range in water. This indicates that the MPs can form stable complexes with the SARS-CoV-2 RNA fragment. Furthermore, the computed E int derived from the E e between the MPs and the SARS-CoV-2 RNA fragment were generally closer to the E int values derived from the E t /E p than the E int values derived from the E v in both vacuum and water. Moreover, there were no significant differences between the E int values derived from the E e and E t /E p (p > 0.05) in vacuum, but significant differences between the E int values derived from the E v and E t /E p (p < 0.05). This implies that the electrostatic interaction contributed mainly to the mechanism of interaction between the MPs and SARS-CoV-2 RNA fragment. The genetic material of the SARS-CoV-2 is positive single-stranded RNA (Zhang et al., 2021c), whereas the studied MPs are neutral and the electrostatic interactions are mainly ion-induced dipole interactions. Moreover, the absolute E int values derived from the E t , E p , or E e for the interactions between the MPs and the SARS-CoV-2 RNA fragment ( Fig. 2B) in water were significantly greater than those in vacuum (p < 0.05) ( Fig. 2A), implying that the interaction affinity of the MPs with the SARS-CoV-2 RNA fragment in water was stronger compared with the affinity in vacuum. This may be caused by the hydrophobicity of MPs (Ding et al., 2020;Zhang et al., 2020), which can provide stronger interactions with the viral RNA fragment in water.  As depicted in Figs. S1 and S2, the E int values derived from the E t and E p for the interaction between the MPs and the SARS-CoV-1 RNA or HBV RNA fragments in vacuum and water phases were significantly lower than those for the interaction between the MPs and SARS-CoV-2 RNA fragment (p < 0.05). This means that the MPs exhibited stronger interaction with the SARS-CoV-2 RNA fragment than with the SARS-CoV-1 RNA and the HBV RNA fragments. Moreover, most of the E int values for the interaction between the MPs and the SARS-CoV-1 RNA fragment or the HBV RNA fragment tended to be positive. This implies that the complexes of the MPs with SARS-CoV-1 RNA fragment or HBV RNA fragment were instable. As a result, it is difficult to analysis the interaction mechanisms of the MPs and the SARS-CoV-1 RNA fragment or the HBV RNA fragment.

Correlation of interaction affinity and temperatures
To test the impact of the studied temperature on the interactions of the MPs with the viral RNA fragments, the variation of the interaction affinity with the temperatures was plotted (Figs. 3, S3, and S4). In general, for each of the MPs, the E int values derived from the total energies fluctuated with the temperature. In particular, the E int values between the MPs and SARS-CoV-2 RNA fragment tended to reach the highest value at 298 K in vacuum (Fig. 3A), implying that the interaction affinity between the MPs and SARS-CoV-2 RNA fragment was lowest at 298 K. In water, the E int values between the PS MPs and SARS-CoV-2 RNA fragment decreased with an increase of the temperature (Fig. 3B). A similar phenomenon occurs in the interaction between the PS MPs and SARS-CoV-1 RNA fragment in water (Fig. S3B). In terms of considering the various types of MPs as a whole, the E int values were not significantly different between the temperatures (Figs. 3, S3, and S4). This also means that temperature was not a determinative factor affecting the interaction affinity between the MPs and viral RNA fragments in the present simulation study.

Correlation of interaction affinity and molecular parameters of MP monomers
To explore the impact of the inherent properties of MPs on their interactions with the viral RNA fragments, a correlation was conducted between the interaction affinity and molecular parameters of MP monomers (Tables 1, S2, S3, and Fig. 4). As shown in Table 1, the E int values derived from the total energies for the interaction of each of the MPs with the SARS-CoV-2 RNA fragment in vacuum and water phases correlated with the molecular parameters V M , PSA, and MTI of the MP monomers to varying degree. The degree of correlation tended to be higher in vacuum and at 310 and 298 K in water except for the PS MPs with aromatic hydrocarbons. In particular, the E int values correlated highly (Fig. 4A-C) and significantly ( Fig. 4D-F) with the molecular parameters except for the PS MPs. On the whole, the greater the V M , PSA, and MTI values, the stronger the interactions between the MPs and the SARS-CoV-2 RNA fragment (Fig. 4).
Generally, the E int values derived from the total energies for the interaction of the SARS-CoV-1 (Table S2) or the HBV (Table S3) RNA fragment with the MPs in vacuum and water phases correlated moderately or weakly with the molecular parameters V M , PSA, and MTI of the MP monomers. It can be also found that there was a higher correlation between the E int values for the interaction of the SARS-CoV-1 RNA fragment with the MPs and the molecular parameters of the MP monomers at 310 K in vacuum except for the PS MPs. In addition, for the interaction of the MPs with the SARS-CoV-1 and the HBV RNA fragment, no significant correlation was found between the E int values and the molecular parameters of the MP monomers.

Discussion
Owing to the high prevalence of both enteric and respiratory viruses in the population and the environment, there is significant potential for human viruses to become associated with the plastisphere (Moresco et al., 2021). There are many sources of MPs in the environment and potential pathways for the interaction, colonisation, and dissemination of viruses. We have studied the interaction between three different viruses and five different MPs in water and vacuum air. For these exposure routes we have taken different conditions; being different temperatures, and different coating of the virus. These coatings have been modelled theoretically how  Fig. 3. Variation of the interaction energies derived from the total energies of the five types of MPs with the SARS-CoV-2 RNA fragment in vacuum (A) and in water (B) with the studied temperatures (223, 263, 273, 298, and 310 K). Different letters represent statistically significant differences between the treatments (p < 0.05).

Table 1
Correlation coefficients between the E int values derived from the total energies between the MPs and SARS-CoV-2 RNA fragment and the molecular parameters of the MP monomers. a Correlation model Temperature  the genetic material such as the RNA of a virus is released into cells after the virus undergoes fusion. There, the RNA segments are covered with the nucleocapsid protein enabling to travel to specific organelles such as the ribosome.
The first pathway described is via the respiratory path: MPs can enter the human body through breathing, mainly due to the presence of MP pollution in the air (Amato-Lourenço et al., 2020); indoor dust as well as air in cities were shown to be large contributors. So not only the virus and MPs dose will be higher indoors, also interaction effectivity is large. It makes it a large potential exposure route for humans. It has been proven that face masks can release large numbers of MPs, which were detected in nasal mucus of mask wearers and can be inhaled by human beings (Ma et al., 2021). In a way face masks are preventing inhalation of virus for human not infected, but those infected may even breathe the virus out. It is speculated that the virus can bind to MPs from the mask and human beings inhale them again as an agglomerate.
Second, the SARS-CoV-2 is transmitted primarily through respiratory droplets (Stadnytskyi et al., 2020) and/or aerosols (Liu et al., 2020a). Airborne dust is another transmission route linked to infectious diseases (Maestre et al., 2021;Moreno et al., 2021). More severe weather phenomena such as sandstorms may exacerbate the migration of the virus (Meo et al., 2021). The adsorption of the SARS-CoV-2 on these airborne media can contribute to the long-range transport of the virus. Note that the airborne transmission route refers to the presence of particles with diameter < 5 μm, who can remain in the air for long periods (Morawska and Cao, 2020). The particle sizes of the MPs are also in this scale range. Thus, MPs dispersed in air can be inhaled by humans (Amato-Lourenço et al., 2020). The MPs can be released into the atmospheric air via several sources, e.g., synthetic textiles , tire wear particles (Lee et al., 2020), domestic laundry dryers (O'Brien et al., 2020), etc. Hence, there is a high probability that the MPs and the SARS-CoV-2 will meet in the atmospheric environment. It has been reported that SARS-CoV-2 aerosols may bind to MPs and facilitate virus entry into the human body (Amato-Lourenço et al., 2022). Our results show that the MPs stabilized the SARS-CoV-2 RNA fragment in both vacuum and water. This also means that the MPs could act as a carrier capable of carrying the gene materials of the SARS-CoV-2 and become a new airborne media for the transport of the virus.
Third, a non-droplet transmission is also possible, as the infectious SARS-CoV-2 particles are also present in human excretions (Wiktorczyk-Kapischke et al., 2021). The fragment of the SARS-CoV-2 RNA has been frequently detected in various countries in wastewater (Kumar et al., 2020;La Rosa et al., 2020;Randazzo et al., 2020), particularly hospital effluent (Gonçalves et al., 2021). The transmission of SARS-CoV-2 via the fecaloral route highlights the presence and persistence of SARS-CoV-2 in the aquatic environment (Arslan et al., 2020). Moreover, the SARS-CoV-2 RNA is relatively stable in sewage and non-chlorinated drinking water (Ahmed et al., 2020). The viral RNA was also found to be relatively stable in contrast to the rapid inactivation of infectious SARS-CoV-2 in river and in sea water (Sala-Comorera et al., 2021). The COVID-19 pandemic has a huge impact on the plastic waste management in many countries, in large due to the sudden surge of medical waste which has led to a potential significant release of MPs (Khoo et al., 2021). Recent studies indicated that MPs have a significant abundance in sewage. Therefore, the sewage treatment system may be an important site for the interaction between the MPs and the gene materials of SARS-CoV-2. Belišová et al. (2022) also confirmed the ability of SARS-CoV-2 virus particles to sorb to the surface of MPs, specifically microfibers in wastewater. The present results implied that the MPs stabilized the SARS-CoV-2 RNA fragment in the water phase, regardless of temperature and MP types. Additionally, the persistence of the SARS-CoV-2 RNA fragment when present on the MPs was different from the persistence of the SARS-CoV-1 and HBV RNA fragments. In comparison, the SARS-CoV-2 RNA fragment preferred to maintain on the MPs, which may cause the gene materials of the SARS-CoV-2 to be long lasting on the MPs.
The fourth path is via the oral route such as food and water. The results in Fig. 2 also indicated that the interaction affinity of the MPs with the SARS-CoV-2 RNA fragment in water was stronger compared with the affinity in vacuum by a factor of 10 at least. This means the MPs and viral genetic material may co-present in dairy products we eat. If the MPs are entered through the food chain (Bouwmeester et al., 2015;Mercogliano et al., 2020), the MPs enter cells via endocytosis and then are released into the cytoplasm. Particularly, the intestinal tract is the main place where MPs exist and is the channel into the circulatory system (Fournier et al., 2021;Visalli et al., 2021). In the meanwhile, it is confirmed that the SARS-CoV-2 can effectively infect intestinal epithelial cells and their precursors (Lamers et al., 2020), which reveals the fact that the intestinal tract is the potential infection site of the SARS-CoV-2 in humans. Taken together, an intercellular environment provides an opportunity for interaction between the MPs and the viral RNA segments/nucleocapsid protein. In our study, we revealed that the MPs showed stronger interaction with the SARS-CoV-2 RNA fragment than with its nucleocapsid protein. Comparison and analysis on the E int also supported the finding that the MPs interacted with the SARS-CoV-2 RNA fragment more strongly than with the SARS-CoV-1 or HBV RNA fragments. This also means that the MPs are more apt to stabilize the genetic materials of the SARS-CoV-2 in the intercellular environment, whereas this interaction may limit the transcription and replication of the viral RNA genomes.
The fifth potential route is via inanimate surfaces such as plastic, stainless steel, and glass has been established (Corpet, 2021;Gidari et al., 2021) on which the persistence of the SARS-CoV-2 is detected. For instance, Gidari et al. (2021) showed the ability of SARS-CoV-2 to persist on most common materials such as glass, stainless steel, and plastic with half-lives of 4.2, 4.4, and 5.3 h respectively. The SARS-CoV-2 is thus more stable on plastics than on steel or on glass. With the global outbreak and spread of COVID-19, disposable surgical masks as effective and cheap protective medical equipment have been widely used by the public. The random disposal of masks may result in new and greater MP pollution, because masks made of polymer materials would release MPs after entering the environment. More importantly, potential co-release of the MPs and the SARS-CoV-2 into the environment will be ineluctable. This might be expected as the result of the unreasonable disposal of the masks, especially the masks contaminated with the virus. MPs have been detected in the air. Thus, MPs can deposit upon the surface of various materials. Thus, there may be an opportunity for the interaction of MPs and the virus RNA. There is evidence that the SARS-CoV-2 RNA fragment has been detected on frozen food packaging Liu et al., 2020b), and aquatic products can be a route of transmission of COVID-19. Positive detection of COVID-19 nucleic acid in the samples of frozen food packaging is still occurring. Our theoretical investigation also indicated that the MPs stabilized the SARS-CoV-2 RNA fragment at very low temperatures ranging from 273 to 223 K. The presence of the genetic material of SARS-CoV-2 on the surfaces is not the same as the presence of the infectious virus, but indicates the transit and contact of infected individuals (Casabianca et al., 2022). Therefore, theoretical evidence of interactions between the MPs and the SARS-CoV-2 RNA fragment could support practices (e.g. strict sanitization of medical equipment, supplies, fabrics, environmental surfaces, and air contaminated with pathogens) that reduce the risk of SARS-CoV-2 infection and cut off its transmission route.
The plastisphere is a diverse microbial community of heterotrophs, autotrophs, predators, and symbionts (Zettler et al., 2013). Several studies demonstrated that the gene materials of microorganisms can be extracted from MPs and subsequently identified (Debeljak et al., 2017;Zettler et al., 2013). Regardless of environmental media and temperature, a stable binding between the MPs and the SARS-CoV-2 fragment was proven theoretically. After such a binding, the SARS-CoV-2 fragment is more difficult to degrade in the natural environment. This also means that entering the plastisphere appears to be an important process that significantly affects the global environmental fate of SARS-CoV-2.
SARS-CoV-2 belongs to the family of enveloped, single-strand RNA viruses (Mei and Tan, 2021). The viral membrane of SARS-CoV-2 surrounds a helical nucleocapsid in which the viral genome is encapsulated by the nucleocapsid protein (Savastano et al., 2020). The biological membrane, known as an envelope, contains lipids and proteins. An envelope may increase the viral sensitivity to physical influencing factors (pH, heat, dryness, etc.) as biological membranes are relatively fragile structures. The nucleocapsid protein of SARS-CoV-2 is produced at high levels within infected cells, enhances the efficiency of viral RNA transcription, and is essential for viral replication (Savastano et al., 2020). It is reported that the SARS-CoV-2 RNA is likely to persist for a long time in untreated wastewater (Ahmed et al., 2020). Consequently, it is essential to elucidate the interactions of the MPs with the nucleocapsid protein and SARS-CoV-2 RNA fragment. Further studies are warranted to evaluate the interaction of the MPs with other structural proteins of SARS-CoV-2, e.g. spike, membrane, and envelope. Furthermore, the interactions as addressed in this study are the first stepping stone to meet our precautionary demand for options to handle any new versions of the coronavirus that might emerge in the future.
It was also found that there are differences in the interaction affinity between the MPs with different compositions and SARS-CoV-2 RNA fragment ( Figs. 1 and 2). Notably, the molecular parameters of the PS monomer performed very different in affecting the interaction affinity as compared to the other MP monomers (Fig. 4). The benzene ring contained in PS allowed it to form π-π interactions with the SARS-CoV-2 RNA fragment that might modulate the interaction affinity. The differences in the composition of MPs are most directly reflected in the functional groups contained in their polymeric structural units. The properties of the MP monomer compounds can determine the mechanism of interaction of MPs with organic pollutants, which in turn exhibit a different interaction affinity for organic pollutants (Lee et al., 2014). In addition, changes in environmental conditions such as temperature can modulate the interaction between the MPs and SARS-CoV-2 RNA fragment (Fig. 3). Other factors such as pH, salinity, and dissolved organic matter which may result in differences in the interaction can also not be neglected. Accordingly, the single and combined effects of different environmental factors on the interaction of the MPs and SARS-CoV-2 will need to be considered in subsequent studies.
It is undeniable that in silico methods still have limitations in both space and time scales, which weakens their correlation with experimental observations and available experimental data. Moreover, quality assurance is required to minimize uncertainty in the calculation of toxicological data. In spite of this, in the face of the urgency of the COVID-19 pandemic, in silico methods are a useful tool to investigate the interaction of environmental pollutants such as MPs with the novel coronavirus, particularly the proposed methodologies that rely upon alternatives to biological testing with high risk of infection. Furthermore, in silico methods have the advantages of preliminary screening of high-risk combinations of multiple co-existing pollutants (e.g. SARS-CoV-2 and MPs) in the environment, and it will save valuable research time and efforts (e.g. model validation) as well as prevent infection during experimental testing.

Conclusions
In this work, we carried out MD simulations to investigate the interactions between five MPs and RNA fragments of three viruses including, SARS-CoV-2, SARS-CoV-1, and HBV at temperatures ranging from 223 to 310 K, in vacuum and in water phases. The estimated E int implied that the interactions of the MPs with the SARS-CoV-2 RNA fragment were stronger than those with the SARS-CoV-1 and HBV RNA fragments, regardless of the environmental media, temperature, and MP types. Furthermore, the electrostatic and hydrophobic processes were the predominant mechanisms for the interactions between the MPs and the SARS-CoV-2 RNA fragment, and the interaction affinity was associated with the inherent structural parameters (i.e., V M , PSA, and MTI) of the MP monomers. Our theoretical results suggest that MPs are capable of regulating the behavior and fate of the SARS-CoV-2 RNA fragment in the environment. While MPs are within air, food and water, this plastic pollution could be a secondary pathway for the transmission of human pathogenic virus and hence have consequences for the exposure of humans to SARS-CoV-2, both by the respiratory pathway (enhancing potential exposure) and the touch pathway where the plastic surface binds the SARS-CoV-2 RNA fragment and thus lowers potential exposure and infectious risks for human. It should be noted that the SARS-CoV-2 RNA fragment can be immobilized by MPs which are ubiquitous in the human environments and thus their persistence and circulation would prolong the presence of virus RNA in the environment. This in silico work serves to minimize the challenges of conducting time-consuming and labor-intensive virus experiments with a high risk of infection, while meeting our precautionary need for options to deal with any new versions of coronaviruses that may emerge in the future.

Declaration of competing interest
The authors declare no competing financial interest.