Molecular Modeling and Simulation Analysis of Antimicrobial Photodynamic Therapy Potential for Control of COVID-19

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can enter the host cells by binding the viral surface spike glycoprotein (SG) to angiotensin-converting enzyme 2. Since antiviral photodynamic therapy (aPDT) has been described as a new method for inhibiting viral infections, it is important to evaluate whether it can be used as a photoactivated disinfectant to control COVID-19. In this in silico study, SARS-CoV-2-SG was selected as a novel target for curcumin as a photosensitizer during aPDT to exploit its physicochemical properties, molecular modeling, hierarchical nature of protein structure, and functional analysis using several bioinformatics tools and biological databases. The results of a detailed computational investigation revealed that SARS-CoV-2-SG is most similar to 6VXX_A, with 100% query cover and identity. The predicted structure of SARS-CoV-2-SG displayed that it is a protein with a positive charge and random coil dominates other secondary structures located outside the viral cell. The protein-protein interaction network showed that SARS-CoV-2-SG interacted with ten potential interacting partners. In addition, primary screening of binding modes through molecular docking showed that curcumin desires to bind and interact with residues of SARS-CoV-2-SG as the main site to enhance the yield of aPDT. Overall, the computer simulation reveals that SARS-CoV-2-SG can be a suitable target site for interaction with curcumin during aPDT.


Introduction
Coronavirus disease (COVID- 19), known as a new form of pneumonia disease that emerged in Wuhan, Hubei province (China), in winter 2019, is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) which has confronted the world with the dire state of global emergency health [1]. SARS-CoV-2 is an enveloped virus with a positive-sense, single-stranded RNA genome with sizes ranging from 26 to 32 kilobases in length [2]. e SARS-CoV-2 genome encodes four major structural glycoproteins to produce a structurally complete viral particle which includes the nucleocapsid (N), membrane (M), envelope (E), and spike (S) proteins [3]. It has been described that SARS-CoV-2 can enter the host cells by binding the viral surface spike glycoprotein to angiotensin-converting enzyme 2 [4]. e spike glycoprotein forms homo-trimers protruding from the viral surface; especially, the spike 1 surface unit allows the attachment of the virus to cellular receptors [3,4]. Despite extensive research conducted by scientists and researchers, there is still no effective antiviral therapy to eradicate the virus. However, reports are available on the application of antimicrobial photodynamic therapy (aPDT) to combat viral infections [5][6][7]. e initial clinical applications of aPDT as a potentially effective and safe treatment modality for infectious diseases were directed against superficial viral diseases. PDT is associated with the administration of a photosensitizing agent as a photosensitizer, which selectively accumulates in the target cells and activates following local irradiation with visible light.
Curcumin, 1,7-bis(4-hydroxy-3-methoxyphenyl)-1,6heptadien-3,5-dione, as a nontoxic lipophilic molecule, can rapidly permeate host cell membrane [8]. As mentioned previously, curcumin has been proven to be effective in a wide range of activities against various viruses [9][10][11]. As commented on by Chen et al. [12], curcumin may display antiviral activity against influenza A virus propagation due to being a strong inhibitor of NF-lB signaling and inhibits hemagglutinin (HA) as a viral protein in the influenza A virus. In addition, curcumin significantly decreased the viral titer at subcytotoxic doses. As studies indicate, curcumin targets the early stages of virus infection, especially the attachment stage in the virus life cycle [12]. erefore, in this in silico study, curcumin was considered a photosensitizer in aPDT.
Recently, several clinical trials of systemic and topical aPDT as a treatment approach and adjuvant therapy with surgery, cryotherapy, or chemotherapy have been conducted against viral infections [5,7,13]. Currently, there are two main clinical and medical aspects for the aPDT: (1) the treatment of local viral infections and (2) the area of blood product decontamination. Many reports have also emerged over the years of the use of aPDTas an effective methodology for sterilization of blood or blood products [13][14][15]. According to the literature, SARS-CoV-2 not only causes respiratory tract infections but also can be detected in blood plasma [16].
Several studies have revealed that lipid-enveloped viruses are more susceptible to inactivation using aPDT than nonenveloped viruses [17]. aPDTdamages to the structure of the viral envelope can inhibit the binding of the virus to the cells. Nevertheless, it is important to identify the viral relevant molecular targets with appropriate sensitivity to aPDT. Based on the results of our previous studies [18][19][20], in silico analysis of proteins enhances our knowledge of the target site to modify the photosensitizer for increasing the efficiency of aPDT.
In the current study, we used a number of bioinformatics tools and biological databases to predict the physicochemical properties, molecular modeling, and structure validation of spike glycoprotein of SARS-CoV-2 (SARS-CoV-2-SG) as one of the important factors to virus entry into the host cells in order to evaluate the effectiveness of aPDT with curcumin as a photosensitizer. and Expasy ProtParam server (http://www.expasy.org/cgi-bin/protpraram) were used to predict the physicochemical properties and functional characterization of SARS-CoV-2-SG such as such isoelectric point (pI), molecular weight (Mw), extinction coefficient (EC-quantitative study of protein-protein and protein-ligand interactions), instability index (II-stability of proteins), aliphatic index (AI-relative volume of protein occupied by aliphatic side chains), and grand average of hydropathicity index (GRAVY-sum of all hydropathicity values of all amino acids divided by the number of residues in a sequence).

Protein
ree-Dimensional Model Prediction. e query sequence taken was the amino acid sequence of spike glycoprotein (severe acute respiratory syndrome coronavirus 2). So the comparative homology protein model was predicted through the SWISS-Model Workspace (https:// swissmodel.expasy.org/interactive). e quality of the SARS-CoV-2-SG structure was validated by means of programs PROCHECK, ERRAT, VERIFY3D, and SAVES that are available at the structure analysis and verification server (http://nihserver.mbi.ucla.edu/SAVES/). SARS-CoV-2-SG was analyzed as the target based on the Global Model Quality Estimation and Qualitative Model Energy Analysis (QMEAN) scores. e interaction energy per residue and confirmation of the local and overall model quality of SARS-CoV-2-SG were calculated by means of the ProSA-web server (https://prosa.services.came.sbg.ac.at/ prosa.php). As well as, the signal peptide was predicted using the SignalP-5.0 server (http://www.cbs.dtu.dk/ services/SignalP/).

Protein-Protein Interactions Network Analysis.
STRING (http://string-db.org) was utilized as an approach to assess the protein-protein interactions network and rank their significance or validity as targets. Moreover, the query sequence was also analyzed to determine the family of proteins using the motif finder server (http://www.genome. jp/tools/motif/). e Protter server (http://wlab.ethz.ch/ protter/) was then used to determine the interactive integration and embodiment of annotated and predicted protein sequence features together with experimental proteomic evidence.

Protein-Photosensitizer Docking.
As previously reported, docking studies were performed to identify the preferred orientation and molecular interactions of natural compounds with targeted proteins. In this current study, the binding patterns and affinity estimations for the interaction between SARS-CoV-2-SG and curcumin were performed using molecular docking (SYBYL-X 2.1 program). e topranked conformation with the lowest binding energy was selected.
From the BLAST results, PDB ID: 6VXX_A was found with maximum similarity. e total score for protein alignment between SARS-CoV-2-SG and 6VXX_A with 1281 amino acids was 2657, with 100% and 100% query cover and identity, respectively.

Physicochemical Characterization.
Physicochemical properties have a key role in the characterization of a specific protein. Basic information from the initial structure showed that SARS-CoV-2-SG had 103 positively charged residues (Arg + Lys) and 110 negatively charged residues (Asp + Glu). e aliphatic index was 84.67 which suggested that the protein is thermostable. Its theoretical isoelectric point (PI), extinction coefficient at 280 nm measured in water, and the grand average of hydropathicity (GRAVY) values were 6.24, 148960 M −1 cm −1 , and −0.079, respectively. e lower range of GRAVY indicated that a better interaction was established between protein and water. e estimated half-life of SARS-CoV-2-SG was 30 h (mammalian reticulocytes, in vitro), >20 h (yeast, in vivo), and >10 h (Escherichia coli, in vivo). As a result, it is a stable protein that can be a therapeutic target for aPDT.

Secondary Structure Prediction.
Secondary structure prediction using PSIPRED and GOR IV showed the location and spatial arrangement of each 6VXX_A amino acid separately (Figures 1(a) and 1(b)). As the figure shows, the random coil (Cc: 57.69%, 739 residues) has dominance over other secondary structures, whereas the alpha helices (Hh: 23.33%, 286 residues) and the extended strand (Ee: 19.98%, 256 residues) were the least frequent.

Protein ree-Dimensional Modeling and Functional Analysis.
A three-dimensional structure of SARS-CoV-2-SG is shown in Figure 2(a). e selected template sequence was 6VXX_A with a query sequence of 100% and a resolution of 2.80 angstroms (Ǻ). e oligo state of the predicted protein model was a homo-trimer. e QMEAN and global score were −2.59 and 0.76 ± 0.05, respectively ( Based on an analysis of 118 structures of resolution of at least 2.0 Ǻ and an R-factor no greater than 20%, a good quality model would be expected to have over 90% in the most favored regions. According to the Ramachandran plot analysis, 90.0% (2277/2529) of residues fall in the most favored region [A, B, L]; among them, 243 residues (9.6%) were present in the additional allowed region [a, b, l, p]; and 3 (0.1%) and 6 (0.2%) residues in generously allowed [∼a,∼b,∼l,∼p] and disallowed regions [xx], respectively ( Figure 3). e numbers of end-residues (excl. Gly and Pro), glycine residues, and proline residues were 69, 174, and 144, respectively.
e mean overall quality factor of the predicted model by ERRAT was 92.20, which implicated that any bonded atomic interactions of the generated model were not within the normal range ( Figure 4).
As shown in Figure 5, the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. Also, expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. Good high-resolution structures generally produce values around 95% or higher. For lower resolutions (2.5 to 3A) the average overall quality factor is around 91%. e results of SAVES in the current study show that the overall quality factor is 95.05.
VERIFY3D a table computed from the atomic coordinates of the structure was used to score the compatibility of the three-dimensional structure model with each amino acid sequence. e finding in Figure 6 verified the correctness of a protein model by the compatibility of its three-dimensional profile with its sequence.
e ProSA-web results of local and overall model quality of the predicted model exhibited that the energy profile of the SARS-CoV-2-SG is consistent with the reliable conformation with the Z-Score of −10.65 (Figure 7). Hence, results gained from PROCHECK, ERRAT, VERIFY3D, SAVES, and ProSA indicate that the modeled structure of SARS-CoV-2-SG is reasonable and reliable as a therapeutic target for aPDT.

Protein-Protein Interactions Network Analysis.
Functional analysis by STRING revealed ten potential interacting partners of transmembrane protease serine 2 (TMPRSS2) in the protein interaction network ( Figure 8). Serine protease proteolytically cleaves and activates the viral SP. As shown in Figure 6, the protein-protein interaction network, their length, and their type of relationship with TMPRSS2 are displayed. e closest interacting proteins having the shortest node were found to be AR (an androgen receptor) and FKBP5 (as a peptidyl-prolyl cis-trans isomerase FKBP5), while the distant interacting proteins were found to be PTEN (phosphatase and tensin homolog) and ETV4 (ETS translocation variant 4). e result of the signal peptide identification in SARS-CoV-2-SG is shown in Figure 9. e signal peptide is located at amino acids 1-13 (MFUFLLFLTLTSG). As well, the positions of posttranslational modification (PTM), variants, and disulfide bonds are shown in Figure 9.
In addition, nine functional motifs were also detected from the functional study in the case of SARS-CoV-2-SG ( Figure 10).

S F E L L H A P A T V C G P K K S T C E E E C E E E E E E E C C E E E E E E E E E C C C C C C C C C E E E E C C C C C E E E E E C C C C N L V K N K C V N F N F N G L T G T G V L T E S N K K F L P F Q Q F G R D I A D T T D A V R D P Q T H H H C C C C C C C C C C C C E E E E E C C C C C C C E E E E E E E E E C C C C E E E E E E E E E E C V L G Q S K R V D F C G K G Y H L M S F P Q S A P H G V V F L H V T Y V P A Q E K N F T T A P A I E E C C E E E C C C C C E E E E E C C E E E E C C C C C C C C C C C C C C C E E E E C C C C E E E E C H D G K A H F P R E G V F V S N G T H W F V T Q R N F Y E P Q I I T T D N T F V S G N C D V V I G E E C C C C C C C C C C C C C C H H H H H H H H H C C C C C C C C C C C C H H H H H H E E C C H H H I V N N T V Y D P L Q P E L D S F K E E L D K Y F K N H T S P D V D L G D I S G I N A S V V N I Q K H H H H H H H H H H H H H H C C C C H H H H C E E E E E E C C C C C C C C C E E E E C C E E E E E E E H H C C C C C C C C C
A P Coil Helix Strand  e molecular interaction of curcumin with SARS-CoV-2-SG was identified by performing a molecular docking study (Figure 11). e best docking orientation was chosen based on the lowest energy value (DG bind ) and negative. Docking analysis showed a strong bond between curcumin and SARS-CoV-2-SG with the lowest energy value (−6.2 kcal/mol). Key hydrogen bond interactions play a significant role for all ligands involving

Discussion
e World Health Organization (WHO) has already declared COVID-19 infection as a global pandemic problem and that there is no antiviral treatment against SARS-CoV-2 [16]. e effect of aPDTas adjunct therapy on viral infections such as the Middle East respiratory syndrome (MERS), and influenza caused by Middle East respiratory syndrome coronavirus (MERS-CoV) and orthomyxoviruses, respectively, can be generalized to COVID-19 due to the similarity of viral pathogenesis to cause primarily mild to severe respiratory infections. In humans, these viruses can be detected with higher viral load and longer duration in the respiratory tract and have been detected also in faces, serum, urine, and blood samples, causing millions of infections each year [21][22][23][24][25][26][27].
Although there is no evidence specifically investigating the effect of aPDT on viral lipids and/ or proteins, there are investigations about the effect of ROS on viral lipids [29][30][31][32]. Since there are viral lipids in the envelope of COVID-19, it should probably be sensitive to the effects of aPDT. e disparate effectiveness of aPDT can be due to variations in the target sites of the photosensitizer. On the other hand, it is unknown which structure of the virus can bind to the photosensitizer. In this study, we introduced the SG as a target site against SARS-CoV-2 using molecular modeling and simulation analysis through numerous data banks. As we know, prediction of the three-dimensional structure of a protein by bioinformatics research because of the nonavailability of the crystal structures is a highly challenging aspect to confirm the data obtained from the NMR or X-ray crystallographic based methods [33].
In concordance with the results obtained, SARS-CoV-2-SG was a homo-trimer, thermostable protein, having an average molecular weight of 139125.14 Da with nine functional motifs. As reported [34,35], molecular docking is very useful and reasonably reliable for the prediction of putative binding modes and affinities of macromolecules (receptors)  e Scientific World Journal and small molecules (ligands). e more negative docking score reveals stronger binding between ligand and receptor. e results of the molecular docking study by Sohilait et al. [34] revealed the binding orientations of curcumin analogues in the active sites of cyclooxygenase-2 towards the design of potent inhibitors with anti-inflammatory activities. In another study, Hobani et al. [35] studied the interaction of curcumin with the enzyme dihydrofolate reductase (DHFR) by molecular docking. ey reported curcumin can be considered an appropriate compound in the development of new inhibitors of DHFR, which is a potential target of anticancer drugs. e findings of this study revealed that curcumin and SARS-CoV-2 had a strong bond with a binding energy value    is in silico study might be a valuable contribution in the field of bioinformatics research in the control of COVID-19 and may help other researchers to get an idea about the protein structure, its physicochemical properties, and protein-protein interaction to conduct extensive examinations on the design and manufacture of drugs and vaccines against COVID-19.

Conclusion
Selected protein structural properties were visualized to select an appropriate photosensitizer and a specific wavelength of light to show the effectiveness and interpretability of the aPDT against COVID-19. Taken together, computer simulation reveals that SARS-CoV-2-SG can be a suitable target site for interaction with curcumin. Also, it suggests the aPDT can be effective against COVID-19 by targeting the spike glycoprotein of SARS-CoV-2 as well as the virus entry step. Further studies to assess the current outcome are needed to achieve better performance of aPDT against COVID-19.

Data Availability
All datasets supporting the conclusions of this article are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.