In December 2019, the World Health Organization was alerted to the occurrence of a new strain of coronavirus which originated in Wuhan, China (Elfiky, 2020; Enayatkhani et al., 2020; Gupta et al., 2020; Huang et al., 2020; Pant et al., 2020; Xu et al., 2020) . Since then, the virus has cataclysmically spread across the world and has infected more than 11 million patients effecting a death toll of around five hundred thousand (updated on 03 July, 2020) (https://www.worldometers.info/coronavirus/). The coronavirus disease , caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been reported to be more fatal to elderly individuals and patients with comorbidities (Al-Khafaji et al., 2020; Elfiky, 2020) . Though bats are believed to be the primary reservoirs for this novel coronavirus, intermediate hosts are suggested to be involved before its final spillover event into humans (Lu et al., 2020) . Of the handful of coronaviruses that are known to infect humans, causing either mild illness or severe disease, SARS-CoV-2 is the seventh and is associated with symptoms including cough, fever and difficulty breathing which may progress to pneumonia and death if unattended (Andersen et al., 2020; Boopathi et al., 2020; Das et al., 2020; Elfiky, 2020; Lobo-Galo et al., 2020) . Significant variability exists between the genome sequence of SARS-CoV-2 in comparison to SARS-CoV (approximately 21%), and MERS-CoV (approximately 50%), the causative agents of the SARS and MERS outbreaks that occurred in recent past (Lu et al., 2020) . The SARS-CoV-2 genome encodes several structural and non-structural proteins (nsp) that play pivotal roles in attaching the virus to host cellular receptors, viral replication and subsequent infection Tai et al., 2020) . The viral spike (S) protein is an important structural protein that facilitates viral attachment, fusion and subsequent entry into the human host via its interaction with the human (host) cellular receptor angiotensin converting enzyme II (ACE 2) (Abdelli et al., 2020; Sinha et al., 2020; Tai et al., 2020) . After attachment to the ACE 2 receptor, the genome of SARS-CoV-2 is released into the cytoplasm and takes control of cellular apparatus to implement its own translation process (Hasan et al., 2020) .
The viral genome encodes several other vital structural proteins that include the envelop protein (a hydrophobic protein responsible for morphogenesis during assemblage of the viral genome), the membrane (M) protein (responsible for the maintenance of the shape of the viral envelope) and the nucleocapsid (N) protein (responsible for RNA binding and dimerization) (Andersen et al., 2020; Sarma et al., 2020; Wahedi et al., 2020) . These structural proteins are produced via a complex proteolytic process involving two proteases namely the chymotrypsin-like cysteine protease (main protease [M pro ] or 3 C-like protease [3CL pro ]) and the papain-like cysteine protease (PL pro ) (Adeoye et al., 2020; Enmozhi et al., 2020; Islam et al., 2020; Joshi et al., 2020; Umesh et al., 2020) . The vital role of the SARS-CoV M pro and PL pro in proteolytic cleavage of the large viral polyprotein orf1ab and viral replication has established them as promising drug and vaccine targets in the areas of therapeutic research against COVID-19 (Aanouz et al., 2020; Elmezayen et al., 2020; Gyebi et al., 2020; 2020) . The sequence of the M pro /3CL pro has been conserved to up to 96.08% in SARS-CoV-2 Ul Qamar et al., 2020) . As such, several protease inhibitors including lopinavir and ritonavir and combinations thereof have been tried to inhibit SARS-CoV-2 protease function and arrest the infection and spread of the virus. However, lack of specific antiviral drugs targeted against SARS-CoV-2 proteases has complicated and impeded treatment efforts.
In recent decades, the global pharmaceutical industry has relied heavily on natural products from microbial sources for the development of novel and effective therapeutics. The present research employed a library of 23 chemically diverse biologically active metabolites (Figure 1(A,B) ) from cyanobacteria to target the main protease (M pro /3CL pro ) as well as the papain-like protease of SARS-CoV-2 using a molecular docking approach. We undertook efforts to provide a comprehensive analysis of the physicochemical and drug-like features of the cyanobacterial metabolites that displayed encouraging inhibitory potential. Our results were further validated by robust molecular dynamics simulations and assessment of molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) binding free energy with a motive to profile the most promising inhibitory candidates against the SARS-CoV-2 M pro and PL pro , targeted towards the development of effective therapeutics.
The high-resolution X-ray diffraction crystal structure of the M pro (PDB ID: 6M03 Chain A; 2.0 Å resolution) and the PL pro (PDB ID: 6W9C Chain A; 2.7 Å resolution) of SARS-CoV-2 were downloaded from the PDB. The respective protein structures were prepared for docking analysis by removing the associated water and ligand molecules. Polar hydrogen atoms and Gasteiger charges were added to the structure using AutoDock tools. Ligands of interest were biologically active compounds isolated from cyanobacteria and were sourced from the NCBI PubChem and ZINC databases. Three dimensional structures of the respective compounds were prepared using the Open Babel server (O'Boyle et al., 2011) and the smiles structures were energy optimized using the PRODRG server (Sch€ uttelkopf & Van Aalten, 2004) . Gromos 96 forcefield was applied to the structures for energy minimization. The refined protein and ligand structures were subjected to molecular docking using AutoDock Vina software (Trott & Olson, 2010) in light of the existing information regarding the active binding pockets of the respective viral proteins (B aez- Santos et al., 2015; Jin et al., 2020) . Molecular docking was performed employing a grid-based docking method opting a rigid protein receptor and flexible ligand docking protocol (Naidoo et al., 2020) . Grid boxes of coordinates 30, 30 and 30 in x, y and z dimensions encompassing the active site residues were employed to carry out molecular docking at the binding pockets of the viral proteins (Naidoo et al., 2020) . Protein-ligand complexes with the lowest binding energy scores and RMSD < 2.0 Å (Naidoo et al., 2020) were used for further interaction analysis using the PyMol software (version 1.7.4) and the Protein-Ligand Interaction Profiler (PLIP) web-server (Salentin et al., 2015) . Binding affinities (DG) derived from AutoDock Vina were used to determine the inhibition constant (K i ) using the following equation:
where, DG 1 represents the binding affinity of the ligand exposed to a conversion factor (4.184) to render kJ À1 mol À1 , R is the gas constant (8.314 J mol À1 K À1 ) and T is the temperature of the system held at 298 K À1 . In order to corroborate the predicted potency of the concerned compounds, we determined the ligand efficiency (LE; Equation 2) and ligand lipophilic efficiency (LLE; Equation 3) as measures to quantify the effectiveness of the structural components of each compound to bind the concerned viral proteases (Hopkins et al., 2004; , Shultz, 2013 .
where, HA is the number of heavy atoms (non-hydrogen atoms) and pK i is the negative logarithm of the inhibition constant.
where, pK i is the negative logarithm of the inhibition constant and cLogP is the calculated partial coefficient.
Cyanobacterial metabolites that displayed binding energy scores À7.5 kcal/mol with the respective viral proteins in AutoDock Vina were chosen for further estimation of physicochemical properties and drug-like features (Elfiky, 2020; Naidoo et al., 2020) . The physicochemical properties were estimated employing the SwissADME server and were validated with the pkCSM server (Daina et al., 2017; Pires et al., 2015) . SwissADME and FAF-Drugs4 (Lagorce et al., 2017) servers were used to detect the presence of pan-assay interference structures (PAINS) associated with structural alerts among the selected compounds. ProTox-II (Banerjee et al., 2018) and vNN-ADMET web servers (Schyman et al., 2017) were used to assess the toxicity indices like mutagenicity (AMES mutagenesis) and cytotoxicity of the compounds.
Molecular dynamics (MD) simulations offer a thorough validation of molecular docking results and allow reliable assessment of potential stability of a protein-ligand complex Kumar et al., 2020) . The ligands displaying the best interaction scores (in terms of binding energy values) with each SARS-CoV-2 M pro and PL pro were subjected to MD simulation for a timescale of 120 nanoseconds (ns) using GROMACS software (version 2019) (Abraham et al., 2015) . The MD simulations were performed as per the protocols opted by . The equilibration steps were set with constant pressure and temperature (NPT) ensemble (Umesh et al., 2020) The MD simulations were carried out at a standard temperature of 300 K and pressure level of 1.013 bar (Umesh et al., 2020) . The structural parameters pertaining to the MD simulations like the root mean square deviations (RMSD), the root mean square fluctuations (RMSF), the radius of gyration (RoG) and the hydrogen bonds of the respective protein-ligand complexes were estimated as a function of time (120 ns) to assess their conformational and structural stability . complexes were estimated employing the 'g_mmpbsa' script of GROMACS (Kumari et al., 2014) . The total free energy of binding of a protein-ligand complex involves the potential energy, polar and non-polar solvation energies and was calculated as:
where DG bind signifies the total energy of binding of the receptor-ligand complex, DG receptor and DG ligand represent the energies of the free receptor and unbound ligand respectively (Bhardwaj et al., 2020) . 
The M pro (3CL pro ) is encoded by the SARS-CoV-2 genome as an integral protease responsible for the cleavage of two polyproteins to form functional components (Andersen et al., 2020) . The M pro has been suggested as an important target candidate for antiviral drug development (Ul Qamar et al., 2020) . It has thus become essential to discover and develop new drugs that specifically target M pro and provide therapeutic relief against SARS-CoV-2. Cyanobacterial metabolites have been identified as promising therapeutics against a range of diseases and disorders owing to an array of biological activities (Tan & Phyo, 2020) . In the present work, a collection of cyanobacterial metabolites that have displayed significant anticancer, antimicrobial and antiviral activities were assessed to explore their inhibitory prospects against the M pro of SARS-CoV-2. Recent reports have suggested that the active site of the M pro displays a stereoselective recognition towards non-covalent inhibitors (Jacobs et al., 2013; Karypidou et al., 2018) . Our analyses thus focused on the identification of compounds that could potentially form strong non-covalent interactions with the catalytic dyad of the M pro . The detailed information regarding the binding energy scores of the 23 compounds with SARS-CoV-2 M pro have been provided in Table 1 (compounds displaying binding energy scores À7.5 kcal/mol) and Table S1 [Supplementary material] (compounds displaying binding energy scores > À7.5 kcal/mol). It was evident from our analysis that out of the 23 compounds, cylindrospermopsin, deoxycylindrospermopsin, carrageenan, cryptophycin 52, eucapsitrione, tjipanazole, tolyporphin and apratoxin A displayed encouraging results with binding energy scores À7.5 kcal/mol (Table 1 ) (Elfiky, 2020; Naidoo et al., 2020) . The cyanotoxin, deoxycylindrospermopsin originally isolated from Cylindrospermopsis sp. displayed the most promising binding efficiency (-8.6 ± 0.02 kcal/mol) and interacted with the important residues Thr26, Cys44, Phe140, Ser144, Cys145, His163 and Glu166 at the active binding pocket of M pro (Table 1 and Figure 2 ). The compound cylindrospermopsin displayed a high binding energy value of À8.1 ± 0.03 kcal/ mol and exhibited hydrophobic interactions with the residues Phe8 and Phe305 and formed hydrogen bonds with Gln110, Thr111, Gln127, Asn151, Asp153, Phe294, Asp295 and Arg298 at the active binding site of M pro (Table 1) .
Similarly, carrageenan displayed an encouraging binding energy score (-8.4 ± 0.02 kcal/mol) and interacted with the two crucial residues of the catalytic dyad of M pro namely, Cys145 through hydrogen bonding and His41 through a salt bridge with the carboxylate group (Table 1) (Jin et al., 2020) . Interestingly, carrageenan did not form hydrophobic interactions with any of the residues at the binding site. The compound is not typically isolated from cyanobacteria but was included in this study after considering the multitude of medicinal uses that include antiviral activity against HIV-1 (Schaeffer & Krylov, 2000) and might be a promising avenue for the development of anti-coronavirus drugs.
The cryptophycins are a large class of peptolides with potent anticancer activity isolated from Nostoc sp. that has undergone clinical trials in cancer treatments. Derivatives of Hydrophobic interactions are represented in italics, Hydrogen bonds are represented in bold and salt bridges are displayed with #. cryptophycin were used as conjugates for antibodies to bind to cysteine mutants that target HER2 and CD22 receptors in human breast cancer (Weiss et al., 2017) . Of the two analogs investigated, cryptophycin 52 exhibited a better binding energy score (-8.3 ± 0.01 kcal/mol) (Table 1 ) relative to cryptophycin 1 (-7.0 ± 0.01 kcal/mol) (Table S1 [Supplementary material]). The compound cryptophycin 52 was noted to interact with the residues Thr25, Leu27, Glu166, Pro168 and Gln189 at the active binding pocket of M pro (Table 1) . The tjipanazoles, isolated from Tolypothrix tjipanasensis are characteristic of carboxazoles exhibiting three benzenoid rings, a hexopyranose ring attached to the aromatic ring system and an indole type NH group (Bonjouklian et al., 1991) . The compound formed hydrogen bonds with essential amino acids Asn142, Leu144, His163, His164 and Glu166 of M pro and returned a binding energy score of À8.0 ± 0.03 kcal/mol (Table 1) . Similarly, tolyporphin isolated as a resistance reversing agent from T. nodosa exhibited a binding energy score of À8.2 ± 0.02 kcal/mol and involved the residues Val104, Ile106, Gln110, Asp153, Phe294, Asp295, Arg298 and Phe305 in its interaction with M pro ( Table 1) .
The apparent biosynthetic potential of cyanobacteria as a producer of diverse and biologically valuable chemicals is further highlighted by the apratoxins, a group of cyclodepsipeptides that contain polypeptide and polyketide fragments (Chen & Forsyth, 2003) . The cytotoxicity of this group of compounds has established them as interesting targets for the development of anticancer drug leads (Luesch et al., 2001 ). In the current investigation, apratoxin A revealed a binding affinity of À7.5 ± 0.02 kcal/mol with the M pro of SARS-CoV-2 forming interactions with the amino acid residues Met49, Asn142, Met165, Glu166 and Gln189 (Table 1) .
Eucapsitrione is an anthraquinone derivative that was isolated as a potent anti-Mycobacterium tuberculosis agent from the cyanobacterium Eucapsis sp. (Sturdy et al., 2010 ). In the current study, the compound was observed to interact with the residues Lys102, Val104, Gln110, Asp295 and Arg298 at the active site of the M pro and produced a binding energy score of À7.8 ± 0.01 kcal/mol (Table 1) .
Thus, it was evident from the molecular docking results that the metabolites cylindrospermopsin, deoxycylindrospermopsin, carrageenan, cryptophycin 52, eucapsitrione, tjipanazole, tolyporphin and apratoxin A (Table 1 ) exhibited promising inhibitory prospects (binding energy scores À7.5 kcal/mol) (Elfiky, 2020; Naidoo et al., 2020) against SARS-CoV-2 M pro and were accordingly selected for further analyses.
Similar to the M pro , the PL pro is responsible for the cleavage of several polypeptides to produce functional viral components (B aez- Santos et al., 2015) . The PL pro also plays an integral role in negating the effects of the host's innate immune system response against viral infection (B aez- Santos et al., 2015) . The human immune response is regulated by ubiquitin, a small protein and the ubiquitin-like modifier interferon stimulated gene 15 (ISG15) (Hu & Sun, 2016) . The protease is capable of removing ubiquitin and ISG15 from cellular conjugates thereby, downregulating the differentiation of T-cells which essentially limits the response of the immune system to infections (B ek es et al., 2015; Hu & Sun, 2016) . The PL pro is thus an essential target for antiviral drugs (B aez- Santos et al., 2015) . Present work involved a thorough molecular docking of the 23 selected biologically active cyanobacterial metabolites at the active binding pocket of the PL pro (B aez- Santos et al., 2015) . The binding energy scores of the 23 compounds with SARS-CoV-2 PL pro have been detailed in Table 2 (compounds displaying binding energy scores À7.5 kcal/mol) and Table S2 [Supplementary material] (compounds displaying binding energy scores > À7.5 kcal/mol). Of the 23 compounds, cryptophycin 1, cryptophycin 52 and deoxycylindrospermopsin showed promising inhibitory potential with binding energy scores À7.5 kcal/ mol (Table 2) (Elfiky, 2020; Naidoo et al., 2020) . The metabolites cryptophycin 1, cryptophycin 52 and deoxycylindrospermopsin displayed binding energy scores of À7.7 ± 0.05 kcal/mol, À7.6 ± 0.02 kcal/mol and À7.9 ± 0.04 kcal/ mol respectively (Table 2) . Deoxycylindrospermopsin seemed to be the most efficient and interacted with the amino acid residues Lys157, Gly163, Asp164, Arg166, Ala246, Tyr264, Thr301 and Asp302 at the active binding pocket of PL pro (Table 2 and Figure 3 ) (Adeoye et al., 2020) . Based on the decent binding energy scores with SARS-CoV-2 PL pro these metabolites were selected for further analyses.
The affinity of a proposed therapeutic candidate to a target receptor is a measure of the strength to which the concerned compound binds with the receptor. Binding affinity is an essential factor that determines the potency of a drug and is inversely proportional to the common measures of potency including the dissociation constant or the inhibition constant (K i ) (Hughes et al., 2011) . Thus, the stronger the affinity of a compound, the lower the concentration required to produce the desired activity. A compound that exhibits a k i within micromolar range (1-40 lM) is considered a hit or lead compound and should be considered for further drug development (Hughes et al., 2011) . The current study identified three compounds namely, carrageenan, cryptophycin 52 and deoxycylindrospermopsin that were predicted to be active within sub-nanomolar range against the viral proteins (Tables 1 and 2 ). However, the potency of a compound is not the only factor to consider when selecting a compound as a drug candidate. Several studies have resolved that the physicochemical features of a compound including molecular weight, lipophilicity, number of hydrogen bond donors and acceptors and non-hydrogen atom properties affect the solubility, membrane permeability, stability and toxicity of a compound (Daina et al., 2017; Hopkins et al., 2014) . Ligand metrics were therefore designed to provide a measure of the biological potency of a compound after considering its physicochemical properties (Hopkins et al., 2014) . For instance, the ligand efficiency (LE) metric deployed in the current study describes the potency of the compound per heavy atom (Hopkins et al., 2004) . The compounds deoxycylindrospermopsin, carrageenan and cryptophycin 52 displayed LE values that fell below the threshold of 0.3 (Hopkins et al., 2004) . In addition to the predicted biological activity of deoxycylindrospermopsin, the compound was predicted to be an efficient ligand. The lipophilic ligand efficiency (LLE) metric on the other hand, is a measure of the ability of the compound to bind specifically to the target receptor relative to its polar or hydrophobic state (log P) (Hopkins et al., 2014) . The compounds carrageenan, eucapsitrione, cylindrospermopsin and deoxycylindrospermopsin exhibited LLE values greater than the proposed acceptable threshold (Leeson & Springthorpe, 2007; Ul Qamar et al., 2020) .
The physicochemical properties of the most promising inhibitory candidates (binding energy score À7.5 kcal/mol) among the cyanobacterial metabolites were assessed in accordance with Lipinski's rule of five which predicts the oral bioavailability and membrane permeability of a drug candidate in terms of the pharmacokinetic properties of a concerned compound (Table 3 ) (Lipinski et al., 2001 ). Lipinski's rule of five proposes that a drug candidate with logarithm of partial coefficient (log P) 5, molecular weight 500 Dalton, number of hydrogen bond acceptors 10 and number of hydrogen bond donors 5 exhibits decent gastrointestinal absorption and enhanced membrane permeability (Lipinski et al., 2001) . Despite its promising inhibitory potential (Tables 1 and 2) , carrageenan was noted to violate the Lipinski's rule of five (molecular weight > 500 Dalton) ( Table  3 ). The compounds apratoxin A, cryptophycin 1, cryptophycin 52 and tolyporphin also showed violations owing to their higher molecular weights than the proposed threshold (Table 3) . Interestingly, deoxycylindrospermopsin displayed decent physicochemical features which were observed to be in accordance with the Lipinski's rule of five (Table 3 ). The compounds apratoxin A, cryptophycin 1, cryptophycin 52 and tolyporphin were predicted to display poor solubility (Table 3) , whereas, carrageenan, eucapsitrione and tjipanazole exhibited moderate solubility. Furthermore, apratoxin A, carrageenan and tolyporphin exhibited low bioavailability scores of 0.17 (Table 3 ) (Naidoo et al., 2020) . However, deoxycylindrospermopsin was predicted to be highly soluble and showed decent bioavailability score of 0.55 (Table 3) (Naidoo et al., 2020) . The solubility of a drug candidate is an essential pharmacokinetic index that defines how efficiently a drug will be metabolized (Naidoo et al., 2020) . The cyanobacterial metabolites were predicted not to possess any structural alerts associated with toxic consequences on PAINS toxicophore estimation (Table 3) (Naidoo et al., 2020; Pouliot & Jeanmart, 2016) . The compound tjipanazole A1 was predicted to be potentially mutagenic (Table 3 ). The cryptophycins have been reported to be associated with severe neurotoxicity in clinical trials as a result of their superior cytotoxic activity (Weiss et al., 2017) . Although antibody conjugates of cryptophycin were deemed safer for anticancer therapy, in the interest of producing a rapid anti-coronavirus treatment, this might not be an appropriate choice (Weiss et al., 2017) . Cylindrospermopsins are cyanotoxins produced by a variety of freshwater cyanobacteria (Berry et al., 2012) . Despite the relatively non-toxic nature of deoxycylindrospermopsin, variants cylindrospermopsin and 7-epicylindrospermopsin have been found to accumulate in the liver of rats causing DNA fragmentation and display toxic effects (Solter & Beasley, 2013) . Thus, a thorough analysis of the physicochemical features and toxicity assessment established deoxycylindrospermopsin as the most effective and safe candidate for further analyses as it fulfilled the rules of drug-likeness and was predicted to be non-toxic. Hydrophobic interactions are represented in italics, Hydrogen bonds are represented in bold and salt bridges are displayed with #. 
The compound deoxycylindrospermopsin exhibited better binding potential with SARS-CoV-2 M pro and PL pro in terms of binding energy in comparison to the other cyanobacterial metabolites (Tables 1 and 2 and Tables S1 and S2 [Supplementary material]). Furthermore, the physicochemical features and drug-like attributes of the compound was noted to be in accordance with the Lipinski's rule of five and it was predicted to be non-toxic (Table 3) . Accordingly, the complexes of deoxycylindrospermopsin with SARS-CoV-2 M pro and PL pro proteins were selected for MD simulations for a timescale of 120 ns to validate the results of molecular docking and assess their conformational stability. A detailed RMSD analysis of the deoxycylindrospermopsin-SARS-CoV-2 M pro complex revealed that it displayed fluctuations in RMSD values of Ca atoms until 94 ns and attained stability thereafter (Figure 4(A) ). The average RMSD value of the complex was observed to be 1.4 Å. On the contrary, the free protein was found to display more fluctuations throughout the course of MD simulation displaying a higher average RMSD value of 1.9 Å (Figure 4(A) ). Similar trend was noted from the RMSD analysis of the deoxycylindrospermopsin-SARS-CoV-2 PL pro complex which experienced fluctuations in RMSD values of Ca atoms until 92 ns and attained equilibrium after that and remained stable (Figure 4(B) ). The average RMSD value of the complex was found to be 1.3 Å. However, the free protein was noted to display fluctuations throughout the timescale of 120 ns and returned a higher average RMSD value of 1.8 Å (Figure 4(B) ). It has been suggested that lower RMSD values signify higher stability of a protein-ligand complex (Muralidharan et al., 2020) . Thus, it was evident from our analysis that the binding of deoxycylindrospermopsin with SARS-CoV-2 M pro and PL pro rendered considerable stability to the complexes. Our observations indicated towards the conformational stability and decent binding potential of deoxycylindrospermopsin with SARS-CoV-2 M pro and PL pro . Analysis of root mean square fluctuation (RMSF) of the residues of a protein-ligand complex is effective in addressing local changes along a protein chain and deciphering its conformational stability (Muralidharan et al., 2020) . Extensive analysis of RMSF of the free proteins and the complexes revealed comparatively lower fluctuations among the complexes compared to the free proteins ( Figures 5(A,B) ), thus indicating higher stability of the protein-ligand complexes. The average RMSF values were 1.6 Å and 1.9 Å for the complexes deoxycylindrospermopsin-SARS-CoV-2 M pro ( Figure 5(A) ) and deoxycylindrospermopsin-SARS-CoV-2 PL pro ( Figure 5(B) ) respectively. On the contrary the average RMSF values of the free receptors SARS-CoV-2 M pro and PL pro were observed to be 2.2 Å ( Figure 5(A) ) and 2.6 Å ( Figure 5(B) ) respectively. The peaks in the RMSF plot indicate the residues which experience major fluctuations during the MD simulations (Muralidharan et al., 2020) . It was evident from the RMSF analysis of the concerned protein-ligand complexes that the residues in the active binding pockets of the SARS-CoV-2 proteins that interacted with deoxycylindrospermopsin (Tables 1 and 2) were relatively stable ( Figures 5(A,B) ) throughout the course of the MD simulations of 120 ns and thus, reflected their conformation stability. Intermolecular hydrogen bonds (H-bonds) impart stability to a protein-ligand complex and play a crucial role in molecular interaction Kumar et al., 2020) . Proper assessment of the formation and deformation of H-bonds is a vital step in analyzing MD simulations and inferring about the stability of a receptor-ligand complex . A detailed analysis of the hydrogen bonding trajectories associated with the MD simulations of the deoxycylindrospermopsin-CoV-2 M pro and deoxycylindrospermopsin-CoV-2 PL pro complexes for a timescale of 120 ns revealed that both the complexes displayed a maximum number of 6 intermolecular hydrogen bonds (Figures 6(A,B) ). The minimum number of hydrogen bonds for the complexes deoxycylindrospermopsin-CoV-2 M pro and deoxycylindrospermopsin-CoV-2 PL pro were observed to be 2 and 1 respectively (Figures 6(A,B) ). The average number of H-bonds was found to be 4 in both the cases (Figures 6(A,B) ). It was evident that in the initial phases of MD simulations the number of H-bonds in both the complexes decreased, however, in the final stages the number of H-bonds increased to 6 in both the cases which implied the conformational stability of the complexes. It has been suggested that intermolecular H-bonds confer stability to a receptor-ligand complex and it was evident from our observations that the complexes of deoxycylindrospermopsin with SARS-CoV-2 M pro and PL pro were considerably stable with 6 hydrogen bonds each at the end of MD simulations of 120 ns (Figure 6(A,B) ).
To further validate our observations, a detailed analysis of the interaction profiles of the concerned protein-ligand complexes after the MD simulations of 120 ns was performed. It was evident that the complexes deoxycylindrospermopsin-SARS-CoV-2 M pro and deoxycylindrospermopsin-SARS-CoV-2 PL pro did not experience any change in interacting residues and pertaining hydrophobic and hydrogen bonds prior to and after the MD simulations of 120 ns (Figures 7(A,B) and Tables 1 and 2) . Thus, the complexes were inferred to be conformationally stable and correlated with the results of molecular docking.
3.6. Analysis of MM-PBSA free energies of binding of deoxycylindrospermopsin with M pro and PL pro and radius of gyration of the complexes
The deoxycylindrospermopsin-SARS-CoV-2 M pro and deoxycylindrospermopsin-SARS-CoV-2 PL pro complexes were selected for further estimation of the free energies of binding using the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. It has been suggested that MM-PBSA provides accurate estimates of free energies of binding of protein-ligand complexes and a more negative value indicates stronger binding . The detailed information regarding the electrostatic energy, SASA (Solvent Accessible Surface Areas) energy, van der Waals energy, polar solvation energy and final binding energy from 20 ns to 120 ns with 20 ns interval has been detailed in Table 4 . It was evident that all forms of energy, except the polar solvation energy, favourably contributed to the interaction deoxycylindrospermopsin with SARS-CoV-2 M pro and SARS-CoV-2 PL pro ( Table 4 ). The complexes deoxycylindrospermopsin-SARS-CoV-2 M pro and deoxycylindrospermopsin-SARS-CoV-2 PL pro displayed final binding energy scores of À45.44 ± 1.66 kcal/mol and À41.39 ± 1.59 kcal/mol respectively (Table 4 ) which implied considerable conformational stability of the complexes after MD simulations of 120 ns. Radius of gyration (RoG) confers meaningful information about the compactness associated with folding and unfolding of protein structures on binding with the ligands . Accordingly, RoG was calculated to assess the compactness and structural stability of the deoxycylindrospermopsin-SARS-CoV-2 M pro and deoxycylindrospermopsin-SARS-CoV-2 PL pro complexes. It has been suggested that higher RoG values signify lower conformational stability and compactness due to unfolding while lower RoG values indicate higher compactness and stability due to folded state of a protein-ligand complex . The average RoG values of the deoxycylindrospermopsin-SARS-CoV-2 M pro and deoxycylindrospermopsin-SARS-CoV-2 PL pro complexes were found to be 20.31 ± 0.2 Å and 21.43 ± 0.3 Å. Our findings revealed that the respective complexes were compact with significant conformational stability.
COVID-19 pandemic has shocked the world with its devastating impact on global public health. The numbers of infections and deaths have been escalating regularly with the rapid spread and transmission of the causative virus SARS-CoV-2. There has been a dire need to develop specific effective therapeutics and prophylactics against the menacing virus and combat COVID-19. The present research effort aimed at exploring the inhibitory prospects of 23 cyanobacterial metabolites against the SARS-CoV-2 M pro and PL pro proteases that have been considered as promising drug targets. Our analysis revealed that the metabolites cylindrospermopsin, deoxycylindrospermopsin, carrageenan, cryptophycin 52, eucapsitrione, tjipanazole, tolyporphin and apratoxin A exhibited promising inhibitory potential against SARS-CoV-2 M pro based on the binding energy scores of the interactions. The compounds cryptophycin 1, cryptophycin 52 and deoxycylindrospermopsin displayed significant inhibitory prospects against the PL pro of SARS-CoV-2, as evident from the binding energies. Subsequent analysis of physicochemical features, potential toxicity, molecular dynamics simulations and MM-PBSA energy scoring function established deoxycylindrospermopsin as the most promising common inhibitory candidate against both the SARS-CoV-2 proteases. Present findings promise to establish foundations to further validate the inhibitory potential of the compound deoxycylindrospermopsin in vitro and in vivo and deploy it as an effective inhibitor of SARS-CoV-2 M pro and PL pro in the quest for effective drug development against COVID-19.