key: cord-0981771-l5bwwbqf authors: Oso, Babatunde Joseph; Ogidi, Clement Olusola title: Repositioning of Fungal-based Peptides as Modulators of Angiotensin-converting Enzyme-related Carboxypeptidase, SARS-coronavirus HR2 Domain, and Coronavirus Disease 2019 Main Protease date: 2021-09-28 journal: J Transl Int Med DOI: 10.2478/jtim-2021-0038 sha: 8323ef7a1278b46b2263cc0e3f28aba37c589638 doc_id: 981771 cord_uid: l5bwwbqf BACKGROUND AND OBJECTIVES: Angiotensin-converting enzyme-related carboxypeptidase, SARS-Coronavirus HR2 Domain, and COVID-19 main protease are essential for the cellular entry and replication of coronavirus in the host. This study investigated the putative inhibitory action of peptides form medicinal mushrooms, namely Pseudoplectania nigrella, Russula paludosa, and Clitocybe sinopica, towards selected proteins through computational studies. MATERIALS AND METHODS: The respective physicochemical properties of selected peptides were predicted using ProtParam tool, while the binding modes and binding free energy of selected peptides toward proteins were computed through HawkDock server. The structural flexibility and stability of docked protein–peptide complexes were assessed through iMODS server. RESULTS: The peptides showed an optimum binding afinity with the molecular targets; plectasin from P. nigrella showed the highest binding free energy compared to peptides from R. paludosa and C. sinopica. Besides, molecular dynamic simulations showed all fungal-based peptides could influence the flexibility and stability of selected proteins. CONCLUSION: The study revealed fungal-based peptides could be explored as functional modulators of essential proteins that are involved in the cellular entry of coronavirus. The coronavirus disease 2019 (COVID- 19) , an outbreak caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), which are enveloped, nonsegmented, single-stranded, positive-sense RNA viruses of the family Coronaviridae, subfamily Orthocoronaviridae, and genera Betacoronavirus, has been regarded as a global pandemic disease. [1, 2] The COVID-19 outbreak emerged in Wuhan, Hubei Province, China in December 2019, where it was initially named as the 2019-novel coronavirus (2019-nCoV) on January 12, 2020 by the World Health Organization (WHO). The WHO officially named the disease as COVID-19, and the Coronavirus Study Group (CSG) of the International Committee proposed to name the new coronavirus as SARS-CoV-2, both issued on February 11, 2020. [3, 4] It uncontrollably spread to North America, Oceania, Europe, Africa, and South America during January-February 2020. [5] Since then, the COVID-19 outbreak has created a challenging period in many geographic locations. It has suspended socioeconomic activities, caused JOURNAL OF TRANSLATIONAL INTERNAL MEDICINE / JUL-SEP 2021 / VOL 9 | ISSUE 3 a setback in the world economy, and created public fear and panic, with psychological instability across individuals. Active surveillance, contact tracing, quarantine, and early strong social distancing efforts were reported as the nonpharmacological strategies to limit transmission of the virus. [6] Most countries affected by COVID-19 adopted temporary measures like isolation of suspected individuals, quarantine of people from risk countries, lockdown or curfew in cities or countries, use of hand sanitizer, wearing of face mask in the public, and social distancing as the basic measures to curtail the spread of the virus. With these unconstitutional preventive measures, COVID-19 has spread to 216 countries with 4,619,477 confirmed cases and a death toll of 311,847 as of May 18, 2020. [7] Yet, no vaccine or drug has been substantiated to cure the disease, but there are several ongoing clinical trials of different drugs (vaccines). Computational investigation of natural bioactive compounds could be a potential pathway to discover novel drugs (vaccines) to end the global pandemic infection. Bioactive proteins and peptides from different sources have various bioactivity attributes like antiviral, antibacterial, antifungal, antiparasitic, antihypertensive, antisporadical, and anticancer. These biological activities are exhibited by peptide molecules with different specificities and potencies depending on their structure and amino acid composition. [8] The use of peptides, especially short peptide fragments, to engineer the induction-targeted immune responses has been adopted as an alternative approach in treating various infectious diseases, and thus has proven to be more effective in reducing mortality and morbidity. [9] As peptides alone are weakly immunogenic due to their relatively small size, they require carrier molecules to stabilize and act as an adjuvant in order to achieve a strong immune response as well as toward appropriate health delivery systems. [10] To increase the complexity of peptide-peptide/protein with regard to drug design, repositioning of active ingredients is performed to produce a cost-effective drug capable of activating receptors and acting indirectly or directly to regulate a series of enzymatic activities including ion transportation, transcription of nucleic acids, and various post-translational modifications of translated protein. [11, 12] Drug repositioning is an emergent approach to new drug development from compounds that have already gone through some level of clinical testing and thus examined to have efficacy against diseases divergent from their original application. [13] This can be accelerated by using molecular software to screen existing chemical databases to identify suitable drug candidates. The development of identified small-molecule drugs that target a specific receptor by exploring the conformational binding space of peptide ligands could be effective efficacy against diseases divergent. [14] Computational methods enhance drug discovery when appropriate molecular applications such as cheminformatics, bioinformatics, network biology, and systems biology are adopted. Hence, this serves as a pathway to establish novel methods towards achieving crucial clinical trials against recent diseases. [15, 16] Isolation and identification of drug molecules of natural origin (fungi) with potentially tolerable side effects may, therefore, help to create modulators of SARS-CoV-2 enzymes. Medicinal mushrooms have a high content of proteins and other bioactive peptides that have immense therapeutic applications. [17] The fruiting bodies of Pseudoplectania nigrella contain plectasin, a peptide (fungal defensin) with antiviral potential and strong antibacterial action against various Gram-positive and negative bacteria. [18] Peptide from Russula paludosa (PRP) has a length of 10 amino acids and a molecular mass of 4.5 kDa. [19] PRP has shown effective inhibition against HIV-1 reverse transcriptase. [19] Moreover, a peptide from Clitocybe sinopica (PCS) with a molecular weight of 44 kDa, composed of two subunits of protein with N-terminal amino acid sequence, has been found to possess strong pharmacological effects against microbial agents. [20] This study investigated the putative influence of peptides of wild medicinal mushrooms on angiotensin-converting enzyme-related carboxypeptidase (ACE2), SARS-Coronavirus HR2 Domain (HRD), and COVID-19 main protease (CMP). Peptides from edible mushrooms, namely P. nigrella, R. paludosa, and C. sinopica, were selected for this hypothetical study based on their reported antiviral and other medicinal activities. [17] The 3D structures of selected peptides were predicted by submitting the respective FASTA amino acid sequence of peptides into I-TASSER modeling server (https://zhanglab.ccmb.med.umich.edu/I-TASSER/). [21, 22] The best model for each peptide was selected based on its C-scores and was further validated using Structure Analysis and Verification server SAVES v5.0 (https://servicesn.mbi. ucla.edu/SAVES/). [23] Physicochemical characterization Physicochemical characterization of the peptides was predicted using ProtParam tool on the SIB Expasy server ProtParam tool (https://www.expasy.org/tools/) using mammalian as defined organism. [24] The structures of ACE2, HRD, and CMP were retrieved from the Protein Data Bank (https://www.rcsb.org/) with their PDB IDs being 1R42, 2FXP, and 6LU7, respectively, and presented as receptors. Structural modeling analysis of selected peptides with the receptors was performed using the HawkDock server (http://cadd.zju.edu.cn/ hawkdock/) integrated with the ATTRACT docking algorithm, the HawkRank scoring function, and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) free energy decomposition analysis. [25, 26] The interactions were visualized using PyMOL ver. 1.1eval (De Lano Scientific LLC, CA, USA). All peptides and the protein were converted into Auto Dock Pdbqt format. The binding energy between the peptides and the protein was computed. The conformational stability of protein-ligand interactions was evaluated using molecular dynamics simulation analysis performed through iMODS server (http:// imods.chaconlab.org) by normal mode analysis (NMA), predicting properties such as deformability, mobility profiles, eigenvalues, variance, covariance map, and elastic network of the protein-ligand interactions. [27] In this study, a total of three fungal-based peptide sequences were selected from edible and medicinal mushrooms. The selected peptides from P. nigrella (plectasin), R. paludosa (PRP), and C. sinopica (PCS) with their corresponding FASTA sequence, amino acid length, and sources are presented in Table 1 . Table 2 shows the physicochemical properties of fungalbased peptides, which include molecular weight, theoretical isoelectric point, instability index, aliphatic index, and grand average hydropathy of selected peptides predicted through Expasy's ProtParam server. The results indicated that plectasin had the highest isoelectric point and estimated half-life, while PRP had the highest instability index. The computational study also revealed PCS had the highest aliphatic index and grand average of hydropathicity compared to other selected peptides. The outcome of the molecular docking is presented in Table 3 . The results reveal the energy scores for different models. The first model of each computation with the lowest energy score was selected as the score for each interaction. Plectasin had the highest affinity with all the selected proteins, as indicated by the least value of each binding energy and assessment by MM/GBSA free energy computation through an approximation of solvent electrostatic contribution in the Poisson-Boltzmann model. The molecular interactions of peptides with the selected proteins are shown in Figure 1 , while the interacting residues of the protein and peptide interface are presented in Table 4 . The interactions of the selected peptides with the proteins evaluated via molecular dynamics simulations analysis to examine the flexibility and dynamics of each of the proteins are presented in Figure 2 and Table 5 . Significant changes were noticed in the structural flexibility of all the selected proteins based on main chain deformability, B-factor, variance associated with each normal mode, covariance matrix, elastic network model, and Eigenvalues upon binding of the peptides comparable to wild types, which were without bound peptides. Therapeutic use of bioactive peptides is growing significantly as they are a popular candidate for production of various drugs. The state-of-the-art technology for using peptides as drugs for different clinical applications with advanced activity in therapeutic development for cancer and infection has a great prospect that has attracted billion The interactions between the selected peptides and the coronavirus-associated proteins were assessed through molecular docking analysis. This study predicted the putative molecular mechanism of action of these compounds through computational analyses. The first model of each molecular interaction was selected as the best docked complex based on the c-scores. [26] The binding affinity obtained through the semi-empirical docking scoring function for the complex of plectasin and HRD was found to be the lowest among other complexes, indicating an interaction with the highest binding affinity. Plectasin showed a relatively high affinity toward all the selected proteins. This putative computation positively correlates with the free energy of binding of the complexes computed through molecular mechanics combined with MM/GBSA, which involved minimization of the docked complexes. [34] The docking studies showed that the peptides could have reasonable biological effect that could interfere with the target proteins and limit the activities of the ACE2, HRD, and CMP. Further investigations were carried out on the interactions between the selected peptides and proteins during a timeframe through molecular dynamics simulation studies. These studies are helpful for predicting the effect of binding of ligands on the biological processes involving proteins. [35] Predominantly, a notable increase in deformity of all complexes of ACE2 and HRD with the peptides is found in simulation studies (Figure 2A) , when compared to the unbound ACE2. This shows that binding of selected peptides could influence the biological function of these proteins and their interactions with their respective targets. [36] However, the binding of all the peptides affected the experimental B-factor, which is an atomic displacement parameter computed from the PDB field and the NMA ( Figure 2B ). This substantiated the effect of peptides on the flexibility and rigidity of proteins, which are critical due to their role in the cellular entry and subsequent replication of coronavirus. [37] PCS increases the variance associated with each normal mode of the amino acid residues in all the studied proteins ( Figure 2C ). The covariance matrix shows the connection between the pairs of amino acid residues in the proteins; it showed that plectasin induces correlated motions among the amino acids, which are indicated by an increase in the intensity of red coloration ( Figure 2D ). Correspondingly, plectasin increased the elasticity of HRD and rigidity in CMP, as revealed by the pattern map of the elastic network ( Figure 2E ). The associated eigenvalue, which is directly correlated to the energy required to deform proteins, is enhanced by the binding of PCS to all the proteins (Table 5 ). This implies that PCS could increase the motion stiffness of the proteins and be useful as a better therapeutic agent over small molecules, based on its high target specificity and potency. This observation could dollar markets. [28] Mushrooms produce a large number of pharmaceutically active proteins (lectins, laccases, ribosome-inactivating proteins [RIPs], nucleases, and fungal immunomodulatory proteins [FIPs]) as a natural source of antiviral, antitumor, antimicrobial, immune-enhancing agents. [29] Many medicinal mushrooms have attracted increasing attention as functional foods or medicines for the prevention and treatment of several chronic diseases. [30] The bioactive peptides in mushrooms could be considered as functional products with substantial health benefits beyond basic nutrients. The diversity of bioactive compounds in medicinal mushrooms with new action of mechanisms has attracted the attention of scientists to use them as potential options of treatment against current life-threatening diseases. [31] Physicochemical characterization of selected peptides predicted through the Expasy's ProtParam server showed the molecular weight, theoretical isoelectric point, and halflife of each of the peptides corresponding to the amino acid length of selected peptides. These factors were deduced from the corresponding peptide sequence and could affect the stability and optimal actions of the peptides. The theoretical isoelectric points and half-life revealed charge stability and estimated half-life of each of the peptides. The half-life predicts the time it takes for half of the amount of each peptide to degrade in a cell. Plectasin was predicted to have an optimal in vivo half-life, which indicates its bioavailability could be essential for optimal function. [32] The corresponding estimated half-life in mammalian reticulocytes through in vitro evaluation correlated with the length of their amino acid composition. The instability index, which is an estimate of stability of the peptides in a test tube, revealed that the peptides were hypothetically stable when smaller than 40. [33] The respective aliphatic index of peptides, which was determined by the volume occupied by aliphatic side chains, showed that plectasin and PCS were thermally stable over a wide temperature range. The low grand average of hydropathicity indicates the provable hydrophilicity of peptides. [24] These properties demonstrate that the putative therapeutic application of the selected peptides would depend on their respective sequence-derived physicochemical characteristics. substantiate the use of peptides as better therapeutic agents over small molecules based on their high target specificity and potency. [38] Furthermore, this observation is inversely related to the hypothetical influence of peptides on the associated variance to each normal mode. This study has shown that bioactive peptides from medicinal mushrooms could be developed and explored for their use in peptide-based drugs against coronavirus infection through their modulatory effects on ACE2, HRD, and CMP and probable effect in preventing the cellular entrance and replication of coronavirus. Further experimental studies are, therefore, recommended to validate the putative theoretical potential of the selected peptides. Severe Acute Respiratory Syndrome Coronavirus 2 (SARSCoV-2): an overview of viral structure and host response Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak -an update on the status Coronavirus disease 2019 (COVID-19): current status and future perspectives The COVID-19 Pandemic: A Summary High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2 Coronavirus disease (COVID-19) WHO; 2020 Bioactive peptides: A review Peptide vaccine: progress and challenges More than one reason to rethink the use of peptides in vaccine design Specificity quantification of biomolecular recognition and its implication for drug discovery A comprehensive review on current advances in peptide drug development and design Review of drug repositioning approaches and resources SPIDR: small-molecule peptide-influenced drug repurposing Informatics and computational methods in natural product drug discovery: a review and perspectives Systems biology approaches for advancing the discovery of effective drug combinations Therapeutic proteins and peptides from edible and medicinal mushrooms-review Plectasin is a peptide antibiotic with therapeutic potential from a saprophytic fungus A peptide with HIV-1 reverse transcriptase inhibitory activity from the medicinal mushroom Russula paludosa Purification and characterization of an antibacterial protein from dried fruiting bodies of the wild mushroom Clitocybe sinopica I-TASSER server: new development for protein structure and function predictions improved protein function prediction by combining structure, sequence and proteinprotein interaction information Deviations from Standard Atomic Volumes as a Quality Measure for Protein Crystal Structures SIB bioinformatics resource portal HawkRank: A new scoring function for protein-protein docking based on weighted energy terms HawkDock: a web server to predict and analyze the structures of protein-protein complexes based on computational docking and MM/GBSA Internal coordinates normal mode analysis server Peptide Drugs in: Hughes AB. Amino acids, peptides and proteins in organic chemistry: protection reactions, medicinal chemistry, combinatorial synthesis Research progress of bioactive proteins from the edible and medicinal mushrooms Bioactivities and health benefits of mushrooms mainly from China Medicinal mushrooms as a new source of natural therapeutic bioactive compounds Serum Stability of Peptides Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence Assessing the performance of the MM/PBSA and MM/GBSA methods: I. The accuracy of binding free energy calculations based on molecular dynamics simulations The role of dynamic conformational ensembles in biomolecular recognition Predictions of Protein Flexibility: First Order Measures Structural basis for FEN-1 substrate specificity and PCNA-mediated activation in DNA replication and repair The future of peptide-based drugs Repositioning of fungalbased peptides as modulators of angiotensin-converting enzymerelated carboxypeptidase, SARS-coronavirus HR2 domain, and coronavirus disease 2019 main protease The authors declare no competing interests.