ORIGINAL RESEARCH article
Computational Simulations Identified Marine-Derived Natural Bioactive Compounds as Replication Inhibitors of SARS-CoV-2
- 1Division of Life Sciences, Department of Bio & Medical Big Data (BK4 Program), Research Institute of Natural Science, Gyeongsang National University, Jinju, South Korea
- 2Division of Applied Life Science, Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Gyeongsang National University (GNU), Jinju, South Korea
The rapid spread of COVID-19, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a worldwide health emergency. Unfortunately, to date, a very small number of remedies have been to be found effective against SARS-CoV-2 infection. Therefore, further research is required to achieve a lasting solution against this deadly disease. Repurposing available drugs and evaluating natural product inhibitors against target proteins of SARS-CoV-2 could be an effective approach to accelerate drug discovery and development. With this strategy in mind, we derived Marine Natural Products (MNP)-based drug-like small molecules and evaluated them against three major target proteins of the SARS-CoV-2 virus replication cycle. A drug-like database from MNP library was generated using Lipinski’s rule of five and ADMET descriptors. A total of 2,033 compounds were obtained and were subsequently subjected to molecular docking with 3CLpro, PLpro, and RdRp. The docking analyses revealed that a total of 14 compounds displayed better docking scores than the reference compounds and have significant molecular interactions with the active site residues of SARS-CoV-2 virus targeted proteins. Furthermore, the stability of docking-derived complexes was analyzed using molecular dynamics simulations and binding free energy calculations. The analyses revealed two hit compounds against each targeted protein displaying stable behavior, binding affinity, and molecular interactions. Our investigation identified two hit compounds against each targeted proteins displaying stable behavior, higher binding affinity and key residual molecular interactions, with good in silico pharmacokinetic properties, therefore can be considered for further in vitro studies.
The world is experiencing a global public health emergency due to the outbreak of coronavirus disease 2019 (COVID-19), which is caused by a novel virus known as severe acute respiratory syndrome coronavirus-2 (SARS- CoV-2; Gorbalenya et al., 2020; Sohrabi et al., 2020). The outbreak was first identified in Wuhan, China (Huang et al., 2020; She et al., 2020; Zhu et al., 2020), but has since spread all around the world. As of December 1, 2020, more than 61.8 million confirmed cases of COVID-19, including 1.4 million deaths, have been reported to the WHO. According to the WHO and recent studies, the most common symptoms of SARS-CoV-2 infection are fever, a dry cough, and fatigue. The infection can also lead to pneumonia, SARS, and death in severe cases (Chen et al., 2020a; Wang et al., 2020a). The novel virus belongs to the Coronaviridae family, a type of positive-sense, single-stranded enveloped RNA virus. There are four classes of coronaviruses designated as alpha, beta, gamma, and delta. The betacoronavirus class includes SARS coronavirus (SARS-CoV), Middle East Respiratory Syndrome (MERS) coronavirus (MERS-CoV), and the COVID-19 causative agent SARS- CoV-2 (Pal et al., 2020). In the past, before the discovery of SARS-CoV-2, six types of coronaviruses were known to cause diseases in humans, including SARS-CoV and MERS-CoV (Mathpal et al., 2020; Zhu et al., 2020). The genome analysis of SARS-CoV-2 sequences revealed that it shares a 79% sequence identity with SARS-CoV and 89% with bat SARS-CoV (Lu et al., 2020). This implies that SARS-CoV-2 might have originated in bats (Boni et al., 2020). The genome of the virus is approximately 30 kB long, consisting of two large polyproteins; ORF1a and ORF1b (that proteolytically cleaves to form 16 nonstructural proteins), four structural proteins, and six accessory proteins (Snijder et al., 2003; Khailany et al., 2020; Yoshimoto, 2020). The structural proteins of SARS-CoV-2 are the enveloped protein (E protein) which helps in morphogenesis during the assemblage of the viral genome, the membrane protein (M protein) which maintains the shape of the virus envelope, and the nucleocapsid (N protein) which is responsible for RNA binding and dimerization (Fehr and Perlman, 2015; Schoeman and Fielding, 2019; Astuti and Ysrafil, 2020). The viral glycosylated spike (S) protein is responsible for the entry of both SARS-CoV and SARS-CoV-2 into the host via interacting with a receptor protein called angiotensin-converting enzyme 2 (ACE2), found on the surface membrane of the host cells (Tai et al., 2020). After the entry of the SARS-CoV-2 into the cell, the replication of the coronavirus starts with the translation of two overlapping open reading frames (ORF1a and ORF1b) to produce pp1a and pp1ab, generated by a ribosomal frameshift mechanism (Astuti and Ysrafil, 2020; Gul et al., 2020; Kim et al., 2020). These polyproteins are then processed proteolytically by two essential viral cysteine proteases-papain-like protease or PLpro (non-structural protein 3 or nsp 3) and 3C-like protease (3CLpro, also known as the main protease, Mpro, or nsp5) (Chen et al., 2020b). Cleavage of pp1a and pp1ab results in the production of 16 non structural protein (NSP)s, with PLpro and 3CLpro (Prajapat et al., 2020). Nsp12 is the first ORF1b encoded protein and has viral RNA-dependent RNA polymerase (RdRp) activity, it has been reported that it requires nsp7 and nsp8 co-factors for maximal activity of its replication and transcription of positive and negative strands of RNA (Kirchdoerfer and Ward, 2019; Hillen et al., 2020). The role of other NSPs were also reported to play specific functions in viral replication and transcription and has been reviewed in detail (Astuti and Ysrafil, 2020; Chen et al., 2020b). Since the outbreak, researchers have repurposed several drugs that could have potential effectiveness against COVID-19. Currently there are no clinically proven specific antiviral agents available for SARS-CoV-2 infection and researchers have been battling to discover promising cures to save lives. Recently, researchers have targeted viral structural and non-structural proteins using computational drug repurposing as well as experimental approaches (Beck et al., 2020; Caly et al., 2020; Chang et al., 2020; Elfiky and Ibrahim, 2020; Gao et al., 2020; Jin et al., 2020; Li and De Clercq, 2020; Li et al., 2020; Ma et al., 2020; Sekhar, 2020; Ton et al., 2020; Wang et al., 2020b). Identified drugs were not able to treat infections to a significant extent and further research is therefore required for effective drugs or vaccine candidates against SARS-CoV-2. It has been reported that drugs like Lopinavir, Ritonavir, Remdesivir, Umifenovir, and Favipiravir are under clinical trials against SARS- CoV-2 target proteins (Drożdżal et al., 2020; Pawar, 2020), however, these drugs have been reported to show severe side effects in COVID-19 patients (Fan et al., 2020a; Khan et al., 2020b; Türsen et al., 2020). Simultaneously, researchers have focused on the identification of natural compounds via computational as well as experimental methods for SARS-CoV-2 inhibition and have identified potent molecules (Alagu Lakshmi et al., 2020; Alexpandi et al., 2020; Antonio et al., 2020; Bhuiyan et al., 2020; Chikhale et al., 2020; Pandey et al., 2020; Wang and Yang, 2020; Kumar Verma et al., 2021). Recently published reports summarized the role of natural compounds and revealed that a few of them are under clinical trials against COVID-19 (Khalifa et al., 2020; Kumari et al., 2020). Alternative sources of natural compounds against SARS-CoV-2 are of marine origin, and have not been explored in great detail and this can serve as an invaluable approach in the ongoing search for novel COVID-19 inhibitors. The marine environment is a rich source of chemically varied, biologically active natural compounds. Marine chemical compounds include polyketides, terpenes, nitrogen containing compounds, and polysaccharides (Uzair et al., 2011). It has been well-documented that marine compounds show antibacterial, antifungal, anti-malarial, and antiviral activities (Choudhary et al., 2017; Riccio et al., 2020). Antiviral activity of marine compounds has been tested against human immunodeficiency virus-1 (HIV-1), herpes simplex virus-2 (HSV-2), polio virus, dengue virus, measles virus, influenza virus, and SARS virus (Yasuhara-Bell and Lu, 2010; Uzair et al., 2011; Gogineni et al., 2015). A polyphenol from marine algae, namely dieckol, was described to be the most active agent against SARS-CoV Mpro through in vitro studies (Park et al., 2013). Similarly, coumarin derivatives from a marine sponge were also observed to display in vitro activities against SARS-CoV Mpro (De Lira et al., 2007). In the same way, molecular docking-based studies recently identified potential leads of marine origin against SARS-CoV-2 Mpro (Gentile et al., 2020; Khan et al., 2020a; Vijayaraj et al., 2020). This literature survey and the recent studies on marine compounds have therefore inspired us to intensively search for additional potent inhibitors. We therefore designed a virtual screening study for the search of novel drug-like molecules against NSPs 3CLpro, PLpro, and RdRp—deemed as key proteins required for the virus replication cycle of SARS-CoV-2—hence considered as attractive drug targets (Chellapandi and Saranya, 2020; Li and De Clercq, 2020; Muhammed, 2020).
Materials and Methods
The Marine Natural Product (MNP) library was downloaded from Prof. Encinar’s website.1 Compounds were filtered for their drug-likeness properties such as the number of hydrogen bond donors, hydrogen bond acceptors, and AlogP values using Lipinski’s rule of five (Lipinski et al., 2001) and then for absorption, distribution, metabolism, and excretion properties using ADMET descriptors protocol in Discovery Studio (DS) v18 (www.accelrys.com Accelrys Inc. San Diego, CA, United States). The filtered compounds were energy minimized using the steepest descent algorithm with CHARMm forcefield implemented in DS. The optimized drug-like compounds were saved in Sybyl Mol 2 format for further processing of molecular docking.
Preparation of Protein Structures
The crystal structure of the selected proteins 3CLpro (PDB ID: 6Y2F) (Zhang et al., 2020), PLpro (PDB ID: 6WX4) (Rut et al., 2020) and RdRp (PDB ID: 6M71) (Gao et al., 2020) were downloaded from RCSB Protein Data Bank (PDB).2 The 3CLpro was in complex with its substrate α-ketoamide inhibitor 13b, which was reported to inhibit 3CLpro with IC50 0.67 μM, while PLpro was complexed with a peptide inhibitor VIR251, which was identified through a substrate hydrolysis assay. Additionally, two anti-retroviral drugs, Lopinavir and Ritonavir, which were found effective against SARS-CoV-2 proteases in in silico and in vitro studies were considered for comparison with our results against 3CLpro and PLpro (Fan et al., 2020b; Joshi et al., 2020; Mohamed et al., 2021; Servidio and Stellacci, 2021). During the course of this study, the crystal structure of RdRp complexed with its inhibitor substrate was unavailable. However, as per recent reports, Remdesivir was identified as a potent inhibitor of RdRp with EC50 0.77 μM (Wang et al., 2020b). Therefore, Remdesivir was considered as a reference compound against RdRp in our study. The three protein structures were then prepared in DS using Clean Protein tool, water molecules were removed, and hydrogens were added. Protein structures were then energy minimized with Steepest Descent algorithm in DS.
The prepared protein structures (3CLpro, PLpro, and RdRp) were subjected to molecular docking with energy minimized drug-like molecules generated from the MNP library. The bound inhibitors of 3CLpro and PLpro were used to define the active site for docking, whereas the catalytic site for RdRp was defined by considering Remdesivir binding residues in DS (Gao et al., 2020). The docking calculation was performed using Genetic Optimization of Ligand Docking (GOLD v5.2.2) (Jones et al., 1997). The Goldscores (high) and Chemscores (low) of drug-like molecules were used for scoring and rescoring, respectively (Kumar et al., 2017, 2019). The Goldscore is used as a major scoring function as it measures the binding affinity of the ligand, whereas Chemscore calculates the total free energy change when a ligand binds to the protein (Verdonk et al., 2003; Bavi et al., 2017). Ten poses were generated for each drug-like ligand and the optimal pose was selected on the basis of scoring functions and key intermolecular interactions with active site residues of the protein.
Molecular Dynamics Simulations
Selected compounds were further subjected to molecular dynamics (MD) simulations using Groningen Machine for Chemical Simulations (GROMACS v5.1.5) under physiological conditions (Pronk et al., 2013). The parameters and coordinate files for 3CLpro, PLpro, RdRp, and selected potential hit compounds were generated by CHARMMm27 forcefield (Sapay and Tieleman, 2011) in GROMACS and SwissParam (Zoete et al., 2011), respectively. The TIP3P water model was used for each simulation system which was neutralized by addition of 2Na+ ions (3CLpro), 1Cl– ion (PLpro), and 12 Na+ ions (RdRp) in a dodecahedron periodic box. Energy minimization was performed for 50,000 nsteps using the steepest descent algorithm to avoid steric clashes. Equilibration of each system was performed in two stages; the first phase was carried out under a constant number of particles, volume, and temperature (NVT) ensemble for 500 ps at 300 K, using the V-rescale thermostat (Bussi et al., 2007) and in the second phase, the pressure of each system was equilibrated for 500 ps at a constant number of particles, pressure, and temperature (NPT) at 1 bar using a Parrinello–Rahman barostat (Parrinello and Rahman, 1981). Each equilibrated system was simulated for 30 ns under periodic boundary conditions to avoid edge effects. During simulation, the LINCS algorithm (Hess et al., 1997) was used to restrain the bond of heavy atoms and Particle Mesh Ewald (PME) was used for electrostatic interactions (Darden et al., 1993). The MD trajectories were analyzed using DS, GROMACS, and Visual Molecular Dynamics (VMD; Hess et al., 1997).
Binding Free Energy Calculations
Prediction of the binding affinities of small-molecule inhibitors with their biological targets plays an important role in prioritizing compounds to be evaluated experimentally (Rifai et al., 2019). In this study, we have exploited Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) to predict the protein-ligand binding free energies (BFE), by implementing a GROMACS compatible tool “g_mmpbsa” (Kumari et al., 2014). To calculate the BFE, 40 snapshots of protein–ligand complexes were selected evenly from the last 10ns of the MD trajectories. The detailed methodology of the BFE calculation was followed as described in previous reports (Bavi et al., 2017; Kumar et al., 2019; Zeb et al., 2019). The protein–ligand BFE ΔGbind is calculated as:
The final ΔGbind values for protein–ligand complexes were the average values from the 20 to 30 ns of the MD simulation trajectories.
Evaluation of Pharmacokinetic Properties
Investigation of key pharmacokinetic properties (PK) of the final compounds is crucial before considering the identified compounds for further drug development (van de Waterbeemd and Gifford, 2003). Therefore, pharmacokinetic properties such as the absorption, distribution, metabolism, excretion, and toxicity (ADMET) of the final compounds were predicted by a machine-learning platform pkCSM3 (Pires et al., 2015). This computational information could affect the decision to begin with synthesis either via a traditional medicinal chemistry or via a combinatorial chemistry strategy (van de Waterbeemd and Gifford, 2003).
In the present study, a series of computational methods have been applied for the identification of potential compounds against COVID-19 replication proteins. The detailed schematic representation of the approach is shown in Figure 1.
Figure 1. Schematic representation of methodologies used in the current study for the identification of potential COVID-19 replication inhibitors.
Selection of Marine Compounds
In the past, various marine natural compounds have been reported to display antibacterial, antifungal, anti- malarial, and antiviral activities (Choudhary et al., 2017; Riccio et al., 2020). Therefore, in the present study, the MNP library of 14,492 compounds was utilized for virtual screening against SARS-CoV-2 replication proteins. A drug-like database was generated as previously mentioned, and the compounds were then minimized for further screening. A total of 2,033 marine drug-like compounds were selected for further processing of molecular docking against three important SARS-CoV-2 replication proteins (3CLpro, PLpro, and RdRp). The three proteins were prepared by removing all the heteroatoms, adding hydrogens, and consequently minimized prior to docking.
Molecular docking studies were performed to study the initial binding interactions of ligands with proteins using the GOLD program (Jones et al., 1997). The docking parameters were validated using REF inhibitors and the pose showing acceptable RMSD <2 Å for respective proteins were selected (Supplementary Figure 1). Ten conformers were generated for each drug-like compound and were subsequently analyzed on the basis of high Goldscore, low Chemscore, and key intermolecular interactions with target proteins (Kumar et al., 2017). The three-dimensional crystal structure of 3CLpro with PDB ID: 6Y2F was used for molecular docking studies. The active site docking sphere was defined using the co-crystallized α- ketoamide inhibitor as a reference (REF) compound with a radius of 10.18 Å. The molecular docking analysis revealed that the REF compound displayed a Goldscore of 61.23 and a Chemscore of −17.76, and therefore, these scores were used as cut-off criteria to filter the drug-like compounds. Accordingly, 12 compounds exhibited a higher Goldscore and lower Chemscore than the REF compound, while only five of them demonstrated desirable hydrogen bond interactions with either one or more key residues of 3CLpro reported earlier; –His41, Cys145, and Glu166 (Supplementary Table 1; Iftikhar et al., 2020; Wang, 2020; Zhang et al., 2020; Kumar et al., 2021). The molecular docking study for SARS-CoV-2 PLpro was performed with PDB ID: 6WX4. The co-crystallized bound peptide inhibitor VIR251 was used as a REF for defining the active site with a 12 Å radius. The docking analyses revealed that the REF compound displayed a Goldscore of 56.09 and a Chemscore of -23.41. Accordingly, 10 compounds exhibited a higher Goldscore and lower Chemscore than the REF compound, while only four of them demonstrated hydrogen bond interactions with either of the key residues; Cys111, Asp164, Gly271, and Tyr278 (Supplementary Table 1; Naidoo et al., 2020; Rut et al., 2020; Shin et al., 2020). Molecular docking analyses revealed that both the drugs Lopinavir and Ritonavir displayed a lower docking score than the selected potential compounds against 3CLpro and PLpro (Table 1). In case of RdRp, the crystal structure with PDB ID: 6M71 was used as a receptor for the molecular docking study. The catalytic site was defined with a 13 Å radius using the only known inhibitor, Remdesivir, as a REF. The REF compound displayed a Goldscore of 60.95 and a Chemscore of -22.41. Correspondingly, 14 compounds were observed to display better scores than REF, while five compounds demonstrated favorable interactions with key residues Arg553, Arg555, Asp618, Asp760, Asp761, and Arg836 which are deemed essential for inhibition (Alexpandi et al., 2020; Elfiky, 2020; Iftikhar et al., 2020; Koulgi et al., 2020; Yin et al., 2020). The 2D structures, Goldscore, and Chemscore values of the selected compounds are shown (Supplementary Table 1).
Molecular Dynamics Simulation Analysis
Molecular dynamics simulation studies were performed to understand the stability, binding mode, and molecular interactions of potential compounds selected from molecular docking (Supplementary Table 1). A 30ns production run was performed using the GROMACS program. MD simulation analyses was performed for important properties like root mean square deviation (RMSD), BFE and identification of binding mode interaction. Only those compounds which displayed stable behavior throughout the simulation, the interaction with key residues, and the better binding affinity in terms of BFE than REF inhibitor, were considered for detailed analysis. The average structure was calculated from the last 10 ns of the simulation run for each protein–ligand complex to scrutinize the binding mode of the interaction of selected hit compounds with its respective protein structures. The compounds were further ranked according to BFE, calculated using the MM/PBSA tool (Table 1).
Analysis of Stability and Binding Free Energy
The stability of the simulated complexes was analyzed through backbone RMSD analysis and the BFE calculations. The RMSD analyses of Hit1, Hit2, and REF inhibitors for 3CLpro, PLpro, and RdRp (Figure 2) demonstrated that all hits remained stable throughout the 30ns of the production run and displayed convergent RMSD below 0.3 nm. Our results also confirms that Lopinavir and Ritonavir showed the previously mentioned behavior against 3CLpro and PLpro during the MD-simulation run, confirming that the Hit compounds can interact with the target proteins under simulated conditions (Supplementary Figures 2A,B). Simultaneously, binding free energies of the hits were calculated using MM/PBSA methodology. For this purpose, 40 frames were generated from the last 10 ns of trajectories displaying stable RMSD. The binding energies obtained were averaged and displayed (Table 1). The SARS-CoV-2 3CLpro bound REF inhibitor α-ketoamide displayed a BFE of −96.85 kJ/mol, whereas Hit1 and Hit2 exhibited −104.37 and −102.90 kJ/mol of BFE, respectively. This analysis revealed that our hits have greater binding affinity towards 3CLpro and can be considered for further evaluation (Table 1 and Figure 3A). Interestingly, the binding affinity of Lopinavir and Ritonavir for 3CLpro was also found to be less than the selected Hits but the binding affinity of Lopinavir and Ritonavir was found to be better than the co-crystallized ligand 13b, used as REF in the study (Table 1 and Supplementary Figure 2C). In case of SARS-CoV-2 PLpro, the REF inhibitor VIR251 was found to display a BFE value of -51.82 kJ/mol, whereas Hit1 and Hit2 displayed a lower BFE value of -92.57 and -82.17 kJ/mol, respectively (Table 1 and Figure 3B). The binding affinity of Lopinavir and Ritonavir for PLpro was found to be less than the selected Hits but interestingly, both drugs displayed better binding affinity than the Co-crystal inhibitor VIR251 (Table 1 and Supplementary Figure 2D). Similarly, SARS-CoV-2 RdRp REF inhibitor Remdesivir displayed a BFE value of -55.50 kJ/mol whereas, Hit1 and Hit2 was observed to demonstrate a BFE value of -111.52 and -63.12 kJ/mol, respectively (Table 1 and Figure 3C).
Binding Mode and Intermolecular Interactions
The binding mode of interaction for hits was analyzed by calculating the average structure from the last 10ns simulation trajectories. SARS-CoV-2 3CLpro bound α-ketoamide REF inhibitor was observed to display three hydrogen bonds with catalytic site residues His41, Glu166, and Gln189. On the other hand, Hit1 formed four hydrogen bonds with catalytic site residues His41, Val186, Thr190, and Gln192, while Hit2 displayed three hydrogen bonds with residues His41, Glu166, and Gln192. The hit compounds were observed to demonstrate a significant number of non-covalent interactions including hydrogen bonds, van der Waals, and π-alkyl bonds with 3CLpro (Figure 4 and Supplementary Table 2). Most significantly, our analysis revealed that both the hit compounds were found to target the catalytic triad (His41, Cys145, and Glu166) of 3CLpro through hydrogen bonds or hydrophobic interactions. The binding mode of Lopinavir and Ritonavir for 3CLpro was also analyzed and it was observed that both drugs target the catalytic triad residues, but the total number of hydrogen bonds observed are less than the Hit compounds (Supplementary Figure 3). The average structure of SARS-CoV-2 PLpro bound REF inhibitor VIR251 was calculated, and the analysis revealed that the REF inhibitor was observed to demonstrate hydrogen bond interactions with Gly163, Gly271, Tyr264, Thr301, and Asp302. The identified Hit1 formed hydrogen bonds with key residues Asp164 and Tyr268 whereas, Hit2 exhibited hydrogen bonds with Asp164, Arg166, and Tyr268 (Figure 5 and Supplementary Table 2). It is worth noting that our hit compounds were found to target catalytically essential residues (Asp164 and Tyr268) of PLpro through hydrogen bonds. The binding mode of Lopinavir and Ritonavir for PLpro revealed that both drugs target the previously mentioned key residues (Supplementary Figure 4).With respect to RdRp, it was observed that the REF inhibitor Remdesivir, formed hydrogen bonds with Arg555 and Arg836 of RdRp. The binding mode of Hit1 was observed to be slightly different than the REF, and was found to target key residues Arg553, Arg555, Asp623, and Ala762 of the RdRp catalytic pocket. Moreover, Hit1 demonstrated significantly better binding affinity in terms of BFE than the REF inhibitor and was therefore considered as a potential hit (Table 1). The binding mode of Hit2 was observed to be similar to the REF and formed hydrogen bonds with key residues Ser549, Ala550, Arg555, and Asp761 (Figure 5 and Supplementary Table 2).
Figure 4. Binding mode of interaction in 3D (upper panel) for hits with the active site of SARS-CoV-2 3CLpro, (A) REF inhibitor, (B) Hit1 and (C) Hit2. The protein in background is shown as wire representation with gray color, interacting residues are displayed as cyan sticks while REF, Hit1 and Hit2 are displayed as red, purple and orange color sticks, respectively. The lower panel displaying 2D interactions, (D) REF, (E) Hit1 and (F) Hit2 with active site interacting residues of SARS-CoV-2 3CLpro. The hydrogen bonds are shown as green dashed lines while the π-alkyl and van der Waals interactions are displayed as pink and light green spheres, respectively.
Figure 5. Binding mode of interaction in 3D (upper panel) for hits with the active site of SARS-CoV-2 PLpro, (A) REF inhibitor, (B) Hit1 and (C) Hit2. The protein in background is shown as wire representation with gray color, interacting residues are displayed as green sticks while REF, Hit1 and Hit2 are displayed as red, purple and orange color sticks, respectively. The lower panel displaying 2D interactions, (D) REF, (E) Hit1 and (F) Hit2 with active site interacting residues of SARS-CoV-2 PLpro. The hydrogen bonds are shown as green dashed lines, while the π-alkyl and van der Waals interactions are displayed as pink and light green spheres, respectively.
The overall MD simulation analyses revealed that our identified hits displayed stabilities throughout the production run, and demonstrated better binding affinities in terms of their energies with their respective protein partners and interactions with key residues of the active site. We therefore propose that our hits can be considered fit for further analysis against COVID-19.
In silico Prediction of Pharmacokinetic Properties
The competence of therapeutic drugs mainly depends on their pharmacokinetic (PK) properties once within the system. The major parameters considered for prediction of these PK properties against 3CLpro, PLpro, and RdRp were calculated using an online machine-learning platform pkCSM (Table 2). Canonical SMILES IDs of ligands were generated using BIOVIA Draw and were then used as input in pkCSM. The server predicts absorption property of compounds considering Caco-2 permeability, which is widely used as an in vitro model to check the absorption level of orally administered drugs. The compounds showing values greater than 0.90 display a high Caco-2 permeability. Our hits presented with better Caco-2 permeability properties than REF compounds. Hits identified against 3CLpro and PLpro can easily cross the Caco-2 cell lines, whereas hits against RdRp displayed properties similar to the REF inhibitor Remdesivir, exhibiting less potential to cross the cell lines. Intestinal absorption (IA) below 30% indicates poor absorbance. The six identified hits displayed better IA properties than REF compounds. Interestingly, REF compound against PLpro displayed poor IA. The blood brain barrier (BBB) acts as an extra border amidst the circulating blood and the extracellular space of the brain. It has been endorsed that the standard value to cross the BBB is >0.3 high to <-1 low logBB. The detailed BBB permeability (BBBP) analysis revealed that potential hits against three proteins cannot easily cross BBB, thereby demonstrating a lower chance of brain associated toxicity. Cytochrome P450s is considered as an essential enzyme for the metabolism of drugs in the liver. The two main subtypes of cytochrome P450 are CYP2D6 and CYP3A4. The inhibitor of these two enzymes can be metabolized in liver and henceforth, cannot be considered as an idyllic drug-like compound. Metabolism analysis revealed that the six identified hits against SARS-CoV-2 3CLpro, PLpro, and RdRp were not involved in CYP2D6 inhibition. Importantly, the hits against SARS-CoV-2 3CLpro and PLpro were not found to be CYP3A4 inhibitors, whereas Hit1 of RdRp including the REF compound can lead to its inhibition, as they were predicted to be CYP3A4 substrate inhibitors. The total clearance (TC) parameter helps to determine the dosing rate, indicating the rate of drug elimination by the body. The drug approaching the zero value is indicative of a lower clearance rate. The six identified hits were found to have a good TC property than REF compounds. As a final criterion, it is essential to test the toxicity of identified hits as this plays a pivotal role in further drug selection. In the present study, three important toxicity parameters namely AMES mutagenicity, hepatotoxicity (HP Tox), and skin sanitization (SS) were considered and studied in detail. It was observed that Hit2 of both 3CLpro and PLpro may be mutagenic, while other hits were demonstrated to be non-mutagenic. With respect to HP Tox PK property, Hit1 of RdRp and Hit2 of both 3CLpro and PLpro displayed a toxicity similar to the REF compounds of the three proteins, whereas other hits were observed to be non-toxic. Another parameter, namely SS, to test the toxicity of our hits was analyzed and our hits were perceived to be non-sensitive to the skin. From the above overall analyses, we can conclude that Hit1 of 3CLpro, Hit1 of PLpro, and Hit2 of RdRp have excellent in silico PK properties among our tested hits, whereas the remaining hits also display comparable properties when compared with standard values set by pkCSM and values of REF inhibitors.
Novelty and Source Identification
As a final assessment, IUPAC names of our hits were generated using BIOVIA draw, and their chemical name and source of origin was identified (Table 3). Identified hits against 3CLpro were found to originate from sponge and fungus, while PLpro hits originated from red algae, and the RdRp identified hits were found to originate from sponge and starfish (John and Murray, 2008). Furthermore, the PubChem search engine4 was used for the purpose of a novelty search for our identified marine hits and it was found that these compounds were not reported against SARS-CoV-2 in any in silico and in vitro studies. Our identified hits can therefore be recommended as novel therapeutic SARS-CoV-2 replication inhibitors.
Table 3. Molecular structures and properties of identified hits against SARS-CoV-2 3CLpro, PLpro, and RdRp.
Identification of inhibitors that can bind and inhibit the proteins involved in viral replication has been the most effective strategy in antiviral drug discovery (Magden et al., 2005; Iftikhar et al., 2020; Wu et al., 2020). The viral proteases 3CLpro and PLpro are the key proteins for the proteolytic processing of the polyproteins that are helpful for initiating replication, whereas RdRp is a vital replication machinery of the virus, making multiple copies of the RNA genome. Therefore, these proteins have been considered as promising drug targets in the treatment of viral diseases (Fehr and Perlman, 2015; Yadav et al., 2020). Recently, several drugs including Lopinavir, Ritonavir, Remdesivir, and Favipiravir were repurposed and are currently under clinical trials to assess their anti-viral efficacy and safety for the treatment of COVID-19 (Drożdżal et al., 2020). Similarly, the recent study reviewed various natural compounds such as Artemisinin, Berbrine, Cholcgicine, Glycyrrhizin, Hanfangchin A, Lactoferrin, Quercetin, Resveratrol Vitamin C, and Vitamin D were found effective against SARS-CoV-2 druggable targets and are currently in clinical trials (Romeo et al., 2021). The previously mentioned study confirms that natural compounds can act as lead scaffolds for drug design or help in improving immunity. Islam et al., also reviewed the potential of natural products and their derivatives such as lycorine, homoharringtonine, silvestrol, ouabain, tylophorine, and 7-methoxycryptopleurine which were reported to inhibit SARS-CoV and its variant key proteins in nanomolar concentrations (Islam et al., 2020). Unfortunately, to date, none of the drugs have been identified as specific inhibitors for COVID-19 viral infection and research is ongoing to explore novel specific compounds. Correspondingly, the emerging source of natural products that have not been explored against COVID-19, is the marine habitat, home to more than 20,000 natural products discovered thus far (Blunt et al., 2016). It has been reported that marine compounds have antimicrobial, antiviral, antidiabetic, anticancer, and anti- inflammatory properties (Mayer et al., 2010, 2011; Pangestuti and Kim, 2017). Intriguingly, coumarin derived compounds from marine sponge Axinella sp. and phlorotannins isolated from marine algae Eckolina cava were able to target SARS 3CLpro through in vitro studies (De Lira et al., 2007; Park et al., 2013). A molecular docking-based study recently reported the antiviral potential of marine-derived compounds including oleic acid, saringosterol, β-Sitosterol, Caulerpin, Glycoglycerolipids, Kjellmanianone, and Loliolide isolated from red, green, and brown macroalgae as candidate inhibitors against 3CLpro, the Spike protein, and the ACE-2 receptor of SARS-CoV-2 (El-Mageed et al., 2021). Another study also reviewed the therapeutic role of Spirulina algae-derived nutraceuticals in boosting the adaptive and innate immunity against coronavirus as well as other viral diseases (Ratha et al., 2021). Similarly, Ray et al. (2020) reviewed the potential role of marine-derived bioactive polysaccharides including fucan sulfates, ulvan, alginate, agarans, carrageenans, and galactans from red seaweeds against diverse. Correspondingly, the authors have reported on the key role of phlorotannins derived from Sargassum spinuligerum brown algae for COVID-19 therapeutics (Ray et al., 2020). We therefore designed our study to target viral key replication proteins using marine derived natural compounds, by using a series of computational methods. A drug-like database was prepared employing Lipinski’s rule of five and ADMET descriptors from a MNP library of 14,492 compounds. A total of 2,033 compounds were retrieved. Molecular docking based virtual screening was performed using GOLD software for selected proteins. Before the docking experiments, the docking method was validated with co- crystallized ligands (REF) and an acceptable binding pose with lower RMSD was selected (Supplementary Figure 1). The docking results confirm the identification of five, four, and five compounds binding to 3CLpro, PLpro, and RdRp, respectively, on the basis of a cluster analysis, higher docking scores than the reference compounds (Supplementary Table 1), and key residue interactions. Since molecular docking does not involve physiological conditions and real-time behavior of the protein-ligand interactions, selected compounds were subjected to MD simulations to check their stability and binding affinity towards respective target proteins (Blunt et al., 2016; Mathpal et al., 2020).
The structural analysis of the SARS-CoV-2 3CLpro crystal structure bound with α-ketoamide revealed the inhibitor targets of His41, Phe140, Gly143, Cys145, His164, and Glu166 through hydrogen bonds. Furthermore, it has been elucidated that Glu166 is indispensable in maintaining the shape of the catalytic pocket and also for maintaining the enzyme conformation in an active state (Zhang et al., 2020). In a drug repurposing study, interaction energy of Glu166 with hit compounds was reported to be significantly low, thus confirming its imperative role in drug interaction (Wang, 2020). Recently published studies targeting 3CLpro revealed that hydrogen bonds with catalytic residues His41 and Cys145 is essential for the inhibition of 3CLpro (Alexpandi et al., 2020; Iftikhar et al., 2020; Kumar et al., 2020; Tahir ul Qamar et al., 2020; Vardhan and Sahoo, 2020; Zhang et al., 2020). In the present study, it was observed that our hits occupied a catalytic site similar to the REF compound. As explained earlier, the interaction with catalytic triad His41, Cys145, and Glu166 is significant for the inhibition of 3CLpro. Interestingly, it was noticed that our hits formed hydrogen bonds with either His41 or Glu166 or both (Figure 4). Additionally, Hit1 formed stable hydrogen bonds with other active site residues Val186, Thr190, and Gln192, whereas Hit2 formed a hydrogen bond with Gln192. In a recent repurposing study on COVID-19, residues His41, Met49, Asn142, His164, Met165, Glu166, and Gln189 were reported to be common hotspots for inhibiting 3CLpro (Wang, 2020). Likewise, inspection for these residues revealed that Hit1 targeted hotspot residues His41, Met49, Cys145, Met165, and Gln189, while Hit2 targeted residues His41, Met49, Cys145, His164, Met165, Glu166, and Gln189 characterized by various hydrogen, hydrophobic, and π-alkyl bonds (Figure 4 and Supplementary Table 2). Similarly, a previously established binding mode for Lopinavir and Ritonavir against 3CLpro was obtained during our analysis which provides support for the binding mode of identified Hit compounds (Bolcato et al., 2020; Kumar et al., 2021). Recently, a terpenoid (T3) from marine sponge Cacospongia mycofijiensis was also identified through in silico studies against SARS CoV-2 3CLpro showing similar binding interactions like our hit molecules (Sepay et al., 2021). A recent pharmacophore-based study was published where authors screened the same MNP library against SARS-CoV-2 3CLpro (Gentile et al., 2020) and identified 17 different hits with different scaffolds. It is worth noting that their identified hits are different from our proposed ones (Gentile et al., 2020). The BFE calculations also revealed significantly better binding energies for our hits with Hit1 demonstrating a BFE of -104.37 kJ/mol and Hit2 with a BFE of -102.90 kJ/mol than REF compound α-ketoamide with a BFE of -96.64 kJ/mol (Table 1). Additionally, BFE values for Lopinavir and Ritonavir was found to be -101.13 and -97.40 kJ/mol, respectively. It was observed that both drugs have better binding affinity than the co-crystallized REF inhibitor but interestingly, our Hit compounds displayed slightly better affinity than both the drugs (Table 1). The above overall analyses clearly indicates that our hits have better binding affinity than α-ketoamide against 3CLpro. Another SARS-CoV-2 protease PLpro performs the hydrolysis of the peptide bond on the carboxyl side of glycine residue which then results in the release of nsp1, nsp2, and nsp3 proteins essential for virus replication. The peptide inhibitor VIR251 bound crystal structure of PLpro was solved recently and revealed that the inhibitor occupies the active site of protease and makes numerous interactions that are important for inhibition (Rut et al., 2020). Several studies involved in the inhibition of SARS-CoV-2 PLpro published recently reveal that Cys111, Leu162, Asp164, Met208, Pro247, Pro248, Tyr264, and Tyr268 are the crucial residues for its inhibition (Amin et al., 2020; Arya et al., 2020; Naidoo et al., 2020; Rut et al., 2020; Sen Gupta et al., 2020; Shin et al., 2020). Our identified Hit1 and Hit2 was found to target key residues Asp164 and Tyr268 through the hydrogen bond, while Hit2 also formed additional hydrogen bonds with Arg166 (Figure 5). An interesting finding by Shin et al. (2020) reported that when a mutant of Tyr268 is treated with inhibitor GRL-0617, its activity reduces significantly, thus revealing the role of Tyr268 in PLpro inhibition. Our identified hits were observed to bind with Tyr268 characterized by hydrogen bonds. This analysis revealed that our hits can inhibit PLpro in a similar way as GRL-0617. The similar residues also targeted by the REF drug Ritonavir provides support for Hit compound binding behavior (Supplementary Figure 4). Moreover, our hits targeted additional residues essential for inhibition, like Cys111, Leu162, Met208, Pro247, Pro248, and Tyr264 as stated earlier, via π-alkyl and van der Waals interactions (Figure 5 and Supplementary Table 2). BFE calculations for PLpro revealed substantial differences in BFE values for identified Hits and REF inhibitors. The Hit1 and Hit 2 displayed a BFE value of -92.57 kJ/mol and -82.17 kJ/mol, whereas the BFE value for the REF peptide inhibitor VIR251, Lopinavir, and Ritonavir was observed as -51.82, -60.84, and -70.74 kJ/mol respectively (Table 1). It is worth noting that Lopinavir and Ritonavir displayed better binding affinity towards PLpro than the co-crystalized ligand but not better than identified Hit compounds. The overall analysis provides adequate support to propose our hits as potential leads to be targeted against PLpro.
The RdRp, which is a key component of the replication machinery of the virus to make multiple copies of the RNA genome, is a potential therapeutic target (Jiang et al., 2020). The cryo-EM structure of RdRp nsp12 in complex with nsp7 and nsp8 has been recently solved, revealing that RdRp consists of three domains- fingers, palm, and thumb (Gao et al., 2020). Recently, it has been shown that Remdesivir can bind in the active site of RdRp which is mainly in the palm domain (Yin et al., 2020). A literature survey further revealed Ala550, Arg553, Arg555, Asp618, Ser759, Asp760, and Asp761 of SARS-CoV-2 RdRp to be vital residues for its inhibition (Alexpandi et al., 2020; Barage et al., 2020; Elfiky, 2020; Gao et al., 2020; Iftikhar et al., 2020; Koulgi et al., 2020; Singh et al., 2020; Vardhan and Sahoo, 2020; Yin et al., 2020). Our identified Hit1 formed hydrogen bonds with key residues Arg553 and Arg555 and additionally targeted other catalytic residues Asp623 and Ala762. Hit2 was also observed to bind with Ala550 and Arg555 along with residues Ser549 and Asp761 via hydrogen bonds essential for inhibition of RdRp. More importantly, Remdesivir was perceived to form a smaller number of interactions than our identified hits (Figure 6). Similar interactions of Remdesivir with RdRp was shown recently through MD simulations thus providing support for our findings (Koulgi et al., 2020). Additionally, our hits displayed numerous interactions with residues Asp618, Cys622, Phe694, Trp800, His810, Cys813, Asp760, Ser814, and Arg836 characterized by π-alkyl and van der Waals interactions (Figure 6 and Supplementary Table 2). Similarly, BFE of identified RdRp hits through MM/PBSA was calculated with a Hit1 BFE value of -111.52 kJ/mol, followed by a Hit2 BFE value of -63.12 kJ/mol. The identified RdRp hits displayed a significantly better binding energy than the Remdesivir BFE value of -55.50 kJ/mol (Table 1). In addition to the above analysis, identified hits were also scrutinized for their pharmacokinetic properties such as intestinal absorption, CYP P450 inhibition, and their hepatoxicity, by an online platform pkCSM and were compared with reference compounds. The in silico PK results demonstrated that identified MNPs are non-mutagenic, are easily absorbed in the human intestine, display Caco-2 permeability, do not inhibit CYP enzymes, and are essentially non-toxic (Table 2). Overall, we believe that these MNP hits may be efficient and effective drug candidates as novel therapeutics against COVID-19 and can be recommended for further in vitro studies.
Figure 6. Binding mode of interaction in 3D (upper panel) for hits with the active site of SARS-CoV-2 RdRp, (A) REF inhibitor, (B) Hit1 and (C) Hit2. The protein in background is shown as wire representation with gray color, interacting residues are displayed as pink sticks while REF, Hit1 and Hit2 are displayed as red, purple and orange color sticks, respectively. The lower panel displaying 2D interactions, (D) REF, (E) Hit1 and (F) Hit2 with active site interacting residues of SARS-CoV-2 RdRp. The hydrogen bonds are shown as green dashed lines, while the π-alkyl and van der Waals interactions are displayed as pink and light green spheres, respectively.
In this study, marine derived drug-like natural products were used for the inhibition of SARS-CoV-2 replication proteins 3CLpro, PLpro, and RdRp. A drug-like database was generated by applying Lipinski’s rule of five and ADMET descriptors on a database of 14,492 compounds generating a total of 2,033 molecules. Detailed docking analysis revealed that 14 compounds were able to display desirable interactions with the protein active site. Selected compounds were then subjected to MD simulations to check their stability and binding affinity against respective target proteins, supplemented with a rigorous MM/PBSA analysis. Six marine compounds displayed better molecular characteristics and binding free energies than the reference inhibitors, characterizing them as potential leads. Detailed analyses of in silico pharmacokinetic properties of the six identified hits demonstrated exceptional properties indicative of their non-toxic and non-mutagenic nature. We anticipate that our identified hits could substantially influence the development of novel therapeutics against COVID-19 in the future.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
VK, SP, and KL conceived the idea of the project and analyzed the results. VK and SP performed the experiments and compiled the manuscript. All the authors have read and approved the manuscript for submission.
This research was supported by the Bio & Medical Technology Development Program of the National Research Foundation (NRF) and funded by the Korean government (MSIT) (No. NRF-2018M3A9A70-57263).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.647295/full#supplementary-material
Supplementary Figure 1 | Validation of docking parameters using co-crystallized ligands (red) and docked pose of (A) 3CLpro (cyan), (B) PLpro (green), and (C) RdRp (pink).
Supplementary Figure 2 | MD simulation analyses displaying backbone RMSD (A) 3CLpro, (B) PLpro and MM/PBSA binding free energy analyses (C) 3CLpro, (D) PLpro for Lopinavir and Ritonavir.
Supplementary Figure 3 | Binding mode of interaction in 3D (upper panel) and 2D (lower panel) for Lopinavir and Ritonavir with the active site of SARS-CoV-2 3CLpro. The protein in background is shown as wire representation with gray color, interacting residues are displayed as cyan sticks while Lopinavir and Ritonavir are displayed as black and brown color sticks, respectively. The hydrogen bonds are shown as green dashed lines, while the π-alkyl and van der Waals interactions are displayed as pink and light green spheres, respectively.
Supplementary Figure 4 | Binding mode of interaction in 3D (upper panel) and 2D (lower panel) for Lopinavir and Ritonavir with the active site of SARS-CoV-2 PLpro. The protein in background is shown as wire representation with gray color, interacting residues are displayed as cyan sticks while Lopinavir and Ritonavir are displayed as black and brown color sticks, respectively. The hydrogen bonds are shown as green dashed lines, while the π-alkyl and van der Waals interactions are displayed as pink and light green spheres, respectively.
Supplementary Table 1 | Molecular structures and docking scores of drug-like and reference compounds.
Supplementary Table 2 | Intermolecular interactions formed by identified hits with 3CLpro, PLpro, and RdRp.
- ^ http://docking.umh.es/downloaddb
- ^ www.rcsb.org
- ^ http://biosig.unimelb.edu.au/pkcsm/
- ^ https://pubchem.ncbi.nlm.nih.gov/
Alagu Lakshmi, S., Shafreen, R. M. B., Priya, A., and Shunmugiah, K. P. (2020). Ethnomedicines of Indian origin for combating COVID-19 infection by hampering the viral replication: using structure-based drug discovery approach. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1778537 [Epub ahead of print].
Alexpandi, R., De Mesquita, J. F., Pandian, S. K., and Ravi, A. V. (2020). Quinolines-based SARS- CoV-2 3 CLpro and RdRp inhibitors and Spike-RBD-ACE2 inhibitor for drug-repurposing against COVID-19: an in silico analysis. Front. Microbiol. 11:1796. doi: 10.3389/fmicb.2020.01796
Amin, S. A., Ghosh, K., Gayen, S., and Jha, T. (2020). Chemical-informatics approach to COVID-19 drug discovery: monte carlo based QSAR, virtual screening and molecular docking study of some in-house molecules as papain-like protease (PLpro) inhibitors. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1780946 [Epub ahead of print].
Arya, R., Das, A., Prashar, V., and Kumar, M. (2020). Potential inhibitors against papain-like protease of novel coronavirus (COVID-19) from FDA approved drugs. ChemRxiv[Preprint] doi: 10.26434/chemrxiv.11860011
Astuti, I., and Ysrafil (2020). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2): an overview of viral structure and host response. Diabetes Metab. Syndr. Clin. Res. Rev. 14, 407–412. doi: 10.1016/j.dsx.2020.04.020
Barage, S., Karthic, A., Bavi, R., Desai, N., Kumar, R., Kumar, V., et al. (2020). Identification and characterization of novel RdRp and Nsp15 inhibitors for SARS-COV2 using computational approach. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1841026 [Epub ahead of print].
Bavi, R., Kumar, R., Rampogu, S., Kim, Y., Kwon, Y. J., Park, S. J., et al. (2017). Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation. J. Recept. Signal Transduct. 37, 224–238. doi: 10.1080/10799893.2016.1212376
Beck, B. R., Shin, B., Choi, Y., Park, S., and Kang, K. (2020). Predicting commercially available antiviral drugs that may act on the novel coronavirus (SARS-CoV-2) through a drug-target interaction deep learning model. Comput. Struct. Biotechnol. J. 18, 784–790. doi: 10.1016/j.csbj.2020.03.025
Bolcato, G., Bissaro, M., Pavan, M., Sturlese, M., and Moro, S. (2020). Targeting the coronavirus SARS-CoV-2: computational insights into the mechanism of action of the protease inhibitors lopinavir, ritonavir and nelfinavir. Sci. Rep. 10:20927. doi: 10.1038/s41598-020-77700-z
Boni, M. F., Lemey, P., Jiang, X., Lam, T. T. Y., Perry, B. W., Castoe, T. A., et al. (2020). Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic. Nat. Microbiol. 5, 1408–1417. doi: 10.1038/s41564-020-0771-4
Caly, L., Druce, J. D., Catton, M. G., Jans, D. A., and Wagstaff, K. M. (2020). The FDA-approved drug ivermectin inhibits the replication of SARS-CoV-2 in vitro. Antiviral Res. 178:104787. doi: 10.1016/j.antiviral.2020.104787
Chen, N., Zhou, M., Dong, X., Qu, J., Gong, F., Han, Y., et al. (2020a). Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan. China Descriptive Study. Lancet 395, 507–513. doi: 10.1016/S0140-6736(20)30211-7
Chikhale, R. V., Sinha, S. K., Patil, R. B., Prasad, S. K., Shakya, A., Gurav, N., et al. (2020). In- silico investigation of phytochemicals from Asparagus racemosus as plausible antiviral agent in COVID-19. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1784289 [Epub ahead of print].
Choudhary, A., Naughton, L. M., Montánchez, I., Dobson, A. D. W., and Rai, D. K. (2017). Current status and future prospects of marine natural products (MNPs) as antimicrobials. Mar. Drugs 15:272. doi: 10.3390/md15090272
De Lira, S. P., Seleghim, M. H. R., Williams, D. E., Marion, F., Hamill, P., Jean, F., et al. (2007). A SARS-coronovirus 3CL protease inhibitor isolated from the marine sponge Axinella cf. corrugata: structure elucidation and synthesis. J. Braz. Chem. Soc. 18, 440–443. doi: 10.1590/s0103-50532007000200030
John, W. B., and Murray, H. G. M. (2008). Dictionary of Marine Natural ProductsWith CD-ROM - Google Books. Available online at: https://books.google.co.kr/books?id=w1bLBQAAQBAJ&pg=PA912&dq=153660-23-2&hl=en&sa=X&ved=2ahUKEwikkp7c4YTtAhWZxYsBHbXEAzYQ6AEwAHoECAAQAg#v=onepage&q&f=false (accessed November 15, 2020).
Drożdżal, S., Rosik, J., Lechowicz, K., Machaj, F., Kotfis, K., Ghavami, S., et al. (2020). FDA approved drugs with pharmacotherapeutic potential for SARS-CoV-2 (COVID-19) therapy. Drug Resist. Updat. 53:100719. doi: 10.1016/j.drup.2020.100719
El-Mageed, H. R. A., Abdelrheem, D. A., Ahmed, S. A., Rahman, A. A., Elsayed, K. N. M., Ahmed, S. A., et al. (2021). Combination and tricombination therapy to destabilize the structural integrity of COVID-19 by some bioactive compounds with antiviral drugs: insights from molecular docking study. Struct. Chem. doi: 10.1007/s11224-020-01723-5 [Epub ahead of print].
Elfiky, A., and Ibrahim, N. (2020). Anti-SARS and anti-HCV drugs repurposing against the papain- like protease of the newly emerged coronavirus (2019-nCoV). Res. Square [Preprint] doi: 10.21203/rs.2.23280/v1
Fan, S., Xiao, D., Wang, Y., Liu, L., Zhou, X., and Zhong, W. (2020b). Research progress on repositioning drugs and specific therapeutic drugs for SARS-CoV-2. Future Med. Chem. 12, 1565–1578. doi: 10.4155/fmc-2020-0158
Fehr, A. R., and Perlman, S. (2015). “Coronaviruses: an overview of their replication and pathogenesis,” in Coronaviruses: Methods and Protocols, eds H. Maier,! E. Bickerton, and P. Britton (New York, NY: Springer), 1–23. doi: 10.1007/978-1-4939-2438-7_1
Gentile, D., Patamia, V., Scala, A., Sciortino, M. T., Piperno, A., and Rescifina, A. (2020). Putative inhibitors of SARS-COV-2 main protease from a library of marine natural products: a virtual screening and molecular modeling study. Mar. Drugs 18:225. doi: 10.3390/md18040225
Gorbalenya, A. E., Baker, S. C., Baric, R. S., de Groot, R. J., Drosten, C., Gulyaeva, A. A., et al. (2020). The species severe acute respiratory syndrome-related coronavirus: classifying 2019- nCoV and naming it SARS-CoV-2. Nat. Microbiol. 5, 536–544. doi: 10.1038/s41564-020-0695-
Gul, S., Ozcan, O., Asar, S., Okyar, A., Barıs, I., and Kavakli, I. H. (2020). In silico identification of widely used and well-tolerated drugs as potential SARS-CoV-2 3C-like protease and viral RNA- dependent RNA polymerase inhibitors for direct use in clinical trials. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1802346 [Epub ahead of print].
Hess, B., Bekker, H., Berendsen, H. J. C., and Fraaije, J. G. E. M. (1997). LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18, 1463–1472. doi: 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4<3.0.CO;2-H
Huang, C., Wang, Y., Li, X., Ren, L., Zhao, J., Hu, Y., et al. (2020). Clinical features of patients infected with 2019 novel coronavirus in Wuhan. China. Lancet 395, 497–506. doi: 10.1016/S0140-6736(20)30183-5
Iftikhar, H., Ali, H. N., Farooq, S., Naveed, H., and Shahzad-ul-Hussan, S. (2020). Identification of potential inhibitors of three key enzymes of SARS-CoV2 using computational approach. Comput. Biol. Med. 122:103848. doi: 10.1016/j.compbiomed.2020.103848
Islam, M. T., Sarkar, C., El-Kersh, D. M., Jamaddar, S., Uddin, S. J., Shilpi, J. A., et al. (2020). Natural products and their derivatives against coronavirus: a review of the non-clinical and pre- clinical data. Phyther. Res. 34, 2471–2492. doi: 10.1002/ptr.6700
Jiang, Y., Yin, W., and Xu, H. E. (2020). RNA-dependent RNA polymerase: structure, mechanism, and drug discovery for COVID-19. Biochem. Biophys. Res. Commun. 538, 47–53. doi: 10.1016/j.bbrc.2020.08.116
Jones, G., Willett, P., Glen, R. C., Leach, A. R., and Taylor, R. (1997). Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 267, 727–748. doi: 10.1006/jmbi.1996.0897
Khalifa, S. A. M., Yosri, N., El-Mallah, M. F., Guo, R. G. Z., Musharraf, S. G., Du, M., et al. (2020). Screening for natural and derived bio-active compounds in preclinical and clinical studies: one of the frontlines of fighting the coronaviruses pandemic. Phytomedicine doi: 10.1016/j.phymed.2020.153311 [Epub ahead of print].
Khan, M. T., Ali, A., Wang, Q., Irfan, M., Khan, A., Zeb, M. T., et al. (2020a). Marine natural compounds as potents inhibitors against the main protease of SARS-CoV-2—a molecular dynamic study. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1769733 [Epub ahead of print].
Khan, Z., Karataş, Y., Ceylan, A. F., and Rahman, H. (2020b). COVID-19 and therapeutic drugs repurposing in hand: the need for collaborative efforts. Pharm. Hosp. Clin. 56, 3–11. doi: 10.1016/j.phclin.2020.06.003
Koulgi, S., Jani, V., Uppuladinne, M. V. N., Sonavane, U., and Joshi, R. (2020). Remdesivir-bound and ligand-free simulations reveal the probable mechanism of inhibiting the RNA dependent RNA polymerase of severe acute respiratory syndrome coronavirus 2. RSC Adv. 10, 26792–26803. doi: 10.1039/d0ra04743k
Kumar, R., Bavi, R., Jo, M. G., Arulalapperumal, V., Baek, A., Rampogu, S., et al. (2017). New compounds identified through in silico approaches reduce the α-synuclein expression by inhibiting prolyl oligopeptidase in vitro. Sci. Rep. 7, 1–14. doi: 10.1038/s41598-017-11302-0
Kumar, R., Kumar, V., and Lee, K. W. (2021). A computational drug repurposing approach in identifying the cephalosporin antibiotic and anti-hepatitis C drug derivatives for COVID-19 treatment. Comput. Biol. Med. 130:104186. doi: 10.1016/j.compbiomed.2020.104186
Kumar, R., Parameswaran, S., Bavi, R., Baek, A., Son, M., Rampogu, S., et al. (2019). Investigation of novel chemical scaffolds targeting prolyl oligopeptidase for neurological therapeutics. J. Mol. Graph. Model. 88, 92–103. doi: 10.1016/j.jmgm.2018.12.006
Kumar Verma, A., Kumar, V., Singh, S., Goswami, B. C., Camps, I., Sekar, A., et al. (2021). Repurposing potential of Ayurvedic medicinal plants derived active principles against SARS- CoV-2 associated target proteins revealed by molecular docking, molecular dynamics and MM–PBSA studies. Biomed. Pharmacother. 137, 111356. doi: 10.1016/j.biopha.2021.111356
Kumar, Y., Singh, H., and Patel, C. N. (2020). In silico prediction of potential inhibitors for the main protease of SARS-CoV-2 using molecular docking and dynamics simulation based drug- repurposing. J. Infect. Public Health 13, 1210–1223. doi: 10.1016/j.jiph.2020.06.016
Kumari, P., Rawat, K., and Saha, L. (2020). Pipeline pharmacological therapies in clinical trial for COVID-19 pandemic: a recent update. Curr. Pharmacol. Rep. 6, 228–240. doi: 10.1007/s40495-020-00226-5
Li, J. Y., You, Z., Wang, Q., Zhou, Z. J., Qiu, Y., Luo, R., et al. (2020). The epidemic of 2019-novel- coronavirus (2019-nCoV) pneumonia and insights for emerging infectious diseases in the future. Microbes Infect. 22, 80–85. doi: 10.1016/j.micinf.2020.02.002
Lipinski, C. A., Lombardo, F., Dominy, B. W., and Feeney, P. J. (2001). Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 46, 3–26. doi: 10.1016/S0169-409X(00)00129-0
Lu, R., Zhao, X., Li, J., Niu, P., Yang, B., Wu, H., et al. (2020). Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395, 565–574. doi: 10.1016/S0140-6736(20)30251-8
Ma, C., Sacco, M. D., Hurst, B., Townsend, J. A., Hu, Y., Szeto, T., et al. (2020). Boceprevir, GC- 376, and calpain inhibitors II, XII inhibit SARS-CoV-2 viral replication by targeting the viral main protease. Cell Res. 30, 678–692. doi: 10.1038/s41422-020-0356-z
Mathpal, S., Joshi, T., Sharma, P., Joshi, T., Pundir, H., Pande, V., et al. (2020). A dynamic simulation study of FDA drug from zinc database against COVID-19 main protease receptor. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1821785 [Epub ahead of print].
Mayer, A. M. S., Glaser, K. B., Cuevas, C., Jacobs, R. S., Kem, W., Little, R. D., et al. (2010). The odyssey of marine pharmaceuticals: a current pipeline perspective. Trends Pharmacol. Sci. 31, 255–265. doi: 10.1016/j.tips.2010.02.005
Mayer, A. M. S., Rodríguez, A. D., Berlinck, R. G. S., and Fusetani, N. (2011). Marine pharmacology in 2007-8: marine compounds with antibacterial, anticoagulant, antifungal, anti- inflammatory, antimalarial, antiprotozoal, antituberculosis, and antiviral activities; affecting the immune and nervous system, and other miscellaneous mechanisms of action. Comp. Biochem. Physiol. C Toxicol. Pharmacol. 153, 191–222. doi: 10.1016/j.cbpc.2010.08.008
Mohamed, K., Yazdanpanah, N., Saghazadeh, A., and Rezaei, N. (2021). Computational drug discovery and repurposing for the treatment of COVID-19: a systematic review. Bioorg. Chem. 106:104490. doi: 10.1016/j.bioorg.2020.104490
Naidoo, D., Roy, A., Kar, P., Mutanda, T., and Anandraj, A. (2020). Cyanobacterial metabolites as promising drug leads against the Mpro and PLpro of SARS-CoV-2: an in silico analysis. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1794972 [Epub ahead of print].
Pandey, P., Rane, J. S., Chatterjee, A., Kumar, A., Khan, R., Prakash, A., et al. (2020). Targeting SARS-CoV-2 spike protein of COVID-19 with naturally occurring phytochemicals: an in silico study for drug development. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1796811 [Epub ahead of print].
Park, J. Y., Kim, J. H., Kwon, J. M., Kwon, H. J., Jeong, H. J., Kim, Y. M., et al. (2013). Dieckol, a SARS-CoV 3CLpro inhibitor, isolated from the edible brown algae Ecklonia cava. Bioorganic Med. Chem. 21, 3730–3737. doi: 10.1016/j.bmc.2013.04.026
Pires, D. E. V., Blundell, T. L., and Ascher, D. B. (2015). pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. Med. Chem. 58, 4066–4072. doi: 10.1021/acs.jmedchem.5b00104
Pronk, S., Páll, S., Schulz, R., Larsson, P., Bjelkmar, P., Apostolov, R., et al. (2013). GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics 29, 845–854. doi: 10.1093/bioinformatics/btt055
Ratha, S. K., Renuka, N., Rawat, I., and Bux, F. (2021). Prospective options of algae-derived nutraceuticals as supplements to combat COVID-19 and human coronavirus diseases. Nutrition 83:111089. doi: 10.1016/j.nut.2020.111089
Ray, B., Schütz, M., Mukherjee, S., Jana, S., Ray, S., and Marschall, M. (2020). Exploiting the amazing diversity of natural source-derived polysaccharides: modern procedures of isolation, engineering, and optimization of antiviral activities. Polymers 13:136. doi: 10.3390/polym13010136
Riccio, G., Ruocco, N., Mutalipassi, M., Costantini, M., Zupo, V., Coppola, D., et al. (2020). Ten- year research update review: antiviral activities from marine organisms. Biomolecules 10:1007. doi: 10.3390/biom10071007
Rifai, E. A., Van Dijk, M., Vermeulen, N. P. E., Yanuar, A., and Geerke, D. P. (2019). A comparative linear interaction energy and MM/PBSA study on SIRT1-ligand binding free energy calculation. J. Chem. Inf. Model 59, 4018–4033. doi: 10.1021/acs.jcim.9b00609
Rut, W., Lv, Z., Zmudzinski, M., Patchett, S., Nayak, D., Snipas, S. J., et al. (2020). Activity profiling and crystal structures of inhibitor-bound SARS-CoV-2 papain-like protease: a framework for anti–COVID-19 drug design. Sci. Adv. 6, 4596–4612. doi: 10.1126/sciadv.abd4596
Sen Gupta, P. S., Biswal, S., Singha, D., and Rana, M. K. (2020). Binding insight of clinically oriented drug famotidine with the identified potential target of SARS-CoV-2. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2020.1784795 [Epub ahead of print].
Sepay, N., Sekar, A., Halder, U. C., Alarifi, A., and Afzal, M. (2021). Anti-COVID-19 terpenoid from marine sources: a docking, admet and molecular dynamics study. J. Mol. Struct. 1228:129433. doi: 10.1016/j.molstruc.2020.129433
She, J., Jiang, J., Ye, L., Hu, L., Bai, C., and Song, Y. (2020). 2019 novel coronavirus of pneumonia in Wuhan, China: emerging attack and management strategies. Clin. Transl. Med. 9:19. doi: 10.1186/s40169-020-00271-z
Shin, D., Mukherjee, R., Grewe, D., Bojkova, D., Baek, K., Bhattacharya, A., et al. (2020). Papain- like protease regulates SARS-CoV-2 viral spread and innate immunity. Nature 587, 657–662. doi: 10.1038/s41586-020-2601-5
Singh, S., Sk, M. F., Sonawane, A., Kar, P., and Sadhukhan, S. (2020). Plant-derived natural polyphenols as potential antiviral drugs against SARS-CoV-2 via RNA-dependent RNA polymerase (RdRp) inhibition: an in-silico analysis. J. Biomol. Struct. Dyn.. doi: 10.1080/07391102.2020.1796810 [Epub ahead of print].
Snijder, E. J., Bredenbeek, P. J., Dobbe, J. C., Thiel, V., Ziebuhr, J., Poon, L. L. M., et al. (2003). Unique and conserved features of genome and proteome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage. J. Mol. Biol. 331, 991–1004.
Sohrabi, C., Alsafi, Z., O’Neill, N., Khan, M., Kerwan, A., Al-Jabir, A., et al. (2020). World health organization declares global emergency: a review of the 2019 novel coronavirus (COVID-19). Int. J. Surg. 76, 71–76. doi: 10.1016/j.ijsu.2020.02.034
Tahir ul Qamar, M., Alqahtani, S. M., Alamri, M. A., and Chen, L. L. (2020). Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants. J. Pharm. Anal. 10, 313–319. doi: 10.1016/j.jpha.2020.03.009
Tai, W., He, L., Zhang, X., Pu, J., Voronin, D., Jiang, S., et al. (2020). Characterization of the receptor-binding domain (RBD) of 2019 novel coronavirus: implication for development of RBD protein as a viral attachment inhibitor and vaccine. Cell. Mol. Immunol. 17, 613–620. doi: 10.1038/s41423-020-0400-4
Ton, A., Gentile, F., Hsing, M., Ban, F., and Cherkasov, A. (2020). Rapid identification of potential inhibitors of SARS-CoV-2 main protease by deep docking of 1.3 billion compounds. Mol. Inform. 39:2000028. doi: 10.1002/minf.202000028
Vardhan, S., and Sahoo, S. K. (2020). In silico ADMET and molecular docking study on searching potential inhibitors from limonoids and triterpenoids for COVID-19. Comput. Biol. Med. 124:103936. doi: 10.1016/j.compbiomed.2020.103936
Vijayaraj, R., Altaff, K., Rosita, A. S., Ramadevi, S., and Revathy, J. (2020). Bioactive compounds from marine resources against novel corona virus (2019-nCoV): in silico study for corona viral drug. Nat. Prod. Res. doi: 10.1080/14786419.2020.1791115 [Epub ahead of print].
Wang, J. (2020). Fast Identification of Possible Drug Treatment of Coronavirus Disease-19 (COVID- 19) through Computational Drug Repurposing Study. J. Chem. Inf. Model. 60, 3277–3286. doi: 10.1021/acs.jcim.0c00179
Wang, L., Wang, Y., Ye, D., and Liu, Q. (2020a). Review of the 2019 novel coronavirus (SARS- CoV-2) based on current evidence. Int. J. Antimicrob. Agents 55:105948. doi: 10.1016/j.ijantimicag.2020.105948
Wang, M., Cao, R., Zhang, L., Yang, X., Liu, J., Xu, M., et al. (2020b). Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro. Cell Res. 30, 269–271. doi: 10.1038/s41422-020-0282-0
Wang, Z., and Yang, L. (2020). Turning the tide: natural products and natural-product-inspired chemicals as potential counters to SARS-CoV-2 infection. Front. Pharmacol. 11:1013. doi: 10.3389/fphar.2020.01013
Wu, C., Liu, Y., Yang, Y., Zhang, P., Zhong, W., Wang, Y., et al. (2020). Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta Pharm. Sin. B 10, 766–788. doi: 10.1016/j.apsb.2020.02.008
Chang, Y., Tung, Y., Lee, K., Chen, T., Hsiao, Y., and Chang, H. (2020). Potential therapeutic agents for COVID-19 based on the analysis of protease and RNA polymerase docking. Preprint. doi: 10.20944/PREPRINTS202002.0242.V1
Yadav, M., Dhagat, S., and Eswari, J. S. (2020). Emerging strategies on in silico drug development against COVID-19: challenges and opportunities. Eur. J. Pharm. Sci. 155:105522. doi: 10.1016/j.ejps.2020.105522
Yin, W., Mao, C., Luan, X., Shen, D. D., Shen, Q., Su, H., et al. (2020). Structural basis for inhibition of the RNA-dependent RNA polymerase from SARS-CoV-2 by remdesivir. Science 368, 1499–1504. doi: 10.1126/science.abc1560
Zeb, A., Kim, D., Alam, S. I., Son, M., Kumar, R., Rampogu, S., et al. (2019). Computational Simulations Identify Pyrrolidine-2,3-Dione Derivatives as Novel Inhibitors of Cdk5/p25 Complex to Attenuate Alzheimer’s Pathology. J. Clin. Med. 8, 746. doi: 10.3390/jcm8050746
Zhang, L., Lin, D., Sun, X., Curth, U., Drosten, C., Sauerhering, L., et al. (2020). Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved a-ketoamide inhibitors. Science 368, 409–412. doi: 10.1126/science.abb3405
Keywords: COVID-19, 3CLcpsdummypro, PLcpsdummypro, RdRp, molecular dynamics simulations (MD), MM/PBSA binding free energy, pharmacokinetic properties
Citation: Kumar V, Parate S, Yoon S, Lee G and Lee KW (2021) Computational Simulations Identified Marine-Derived Natural Bioactive Compounds as Replication Inhibitors of SARS-CoV-2. Front. Microbiol. 12:647295. doi: 10.3389/fmicb.2021.647295
Received: 29 December 2020; Accepted: 25 February 2021;
Published: 21 April 2021.
Edited by:Mohamed L. Ashour, Ain Shams University, Egypt
Reviewed by:Muhammad Akram, Paracelsus Medical University, Austria
Ruchi Tiwari, U.P. Pandit Deen Dayal Upadhyaya Veterinary University, India
Moustafa E. El-Araby, King Abdulaziz University, Saudi Arabia
Copyright © 2021 Kumar, Parate, Yoon, Lee and Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Keun Woo Lee, email@example.com
†These authors have contributed equally to this work