key: cord-0683620-z4nqzt74
authors: Junior, Nilson Nicolau; Santos, Igor Andrade; Meireles, Bruno Amaral; Nicolau, Mariana Sant’Anna Pereira; Lapa, Igor Rodrigues; Aguiar, Renato Santana; Jardim, Ana Carolina Gomes; José, Diego Pandeló
title: In silico evaluation of lapachol derivatives binding to the Nsp9 of SARS-CoV-2
date: 2021-01-22
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2021.1875050
sha: 8aeb65fa5931e1f1486b849318e6ba960b1e8758
doc_id: 683620
cord_uid: z4nqzt74

SARS-CoV-2 is the etiological agent of COVID-19, which represents a global health emergency that was rapidly declared a pandemic by the World Health Organization. Currently, there is a dearth of effective targeted therapies against viruses. Natural products isolated from traditional herbal plants have had a huge impact on drug development aimed at various diseases. Lapachol is a 1,4- naphthoquinone compound that has been demonstrated to have therapeutic effects against several diseases. SARS-CoV-2 non-structural proteins (nsps) play an important role in the viral replication cycle. Nsp9 seems to play a key role in transcription of the RNA genome of SARS-CoV-2. Virtual screening by docking and molecular dynamics suggests that lapachol derivatives can interact with Nsp9 from SARS-CoV-2. Complexes of lapachol derivatives V, VI, VIII, IX, and XI with the Nsp9 RNA binding site were subjected to molecular dynamics assays, to assess the stability of the complexes via RMSD. All complexes were stable over the course of 100 ns dynamics assays. Analyses of the hydrogen bonds in the complexes showed that lapachol derivatives VI and IX demonstrated strongest binding, with a stable or increasing number of hydrogen bonds over time. Our results demonstrate that Nsp9 from SARS-CoV-2 could be an important target in prospecting for ligands with antiviral potential. In addition, we showed that lapachol derivatives are potential ligands for SARS-CoV-2 Nsp9. Communicated by Ramaswamy H. Sarma

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged recently as the causative agent of the current widespread coronavirus disease . It was quickly declared a pandemic and represents a global health emergency (Chhikara et al., 2020) . Since its identification in December 2019, SARS-CoV-2 has infected more than 84 million people, and led to 1,848,704 deaths (WHO, 2021), posing a worldwide challenge. It is most common modes of human-human transmission include direct transmission via droplets produced during coughing, sneezing, and talking (Doremalen et al., 2020) , and indirect transmission via contaminated inanimate objects (Lai et al., 2020) . Some infected people can be asymptomatic (Lai et al., 2020) . Most symptomatic cases present with symptoms including fever, cough, and acute respiratory disease syndrome (SARS). Severe cases can progress to pneumonia, kidney failure, and death .

SARS-CoV-2 is a Coronaviridae family member, and is an enveloped, single-stranded RNA virus (Hsu et al., 2020, p. 2) . The viral structure is primarily determined by structural proteins, including spike (S), membrane (M), envelope (E), and nucleocapsid (N) proteins (Kim et al., 2020, p. 2) . SARS-CoV-2 replication in host cells is not well understood, and the complete mechanism remains to be elucidated. However, knowledge of another coronavirus may suggest how SARS-CoV-2 replicates in host cells.

The SARS-CoV-2 genome has two open reading frames (ORFs), ORF1a and ORF1b, which encode the nonstructural proteins (nsps) (Andersen et al., 2020) . ORF1a encodes the polypeptide pp1a, which is cleaved into 11 nsps, and ORF1b encodes the polypeptide pp1ab, which is cleaved into 16 nsps (Kim et al., 2020; Lai & Stohlman, 1981; Yogo et al., 1977) . Proteolytic cleavage is performed by viral proteases Nsp3 and Nsp5 (Chan et al., 2020) , enabling complex assembly from other nsps to form the replicase-transcriptase complex (RTC) (Gordon et al., 2020) . The RTC consists of several enzymes, including RNA-dependent RNA polymerase (Nsp12), the main protease (Nsp5), an endonuclease (Nsp15), and a papain-like protease (Nsp3) (Chan et al., 2020) .

Nsp9, encoded by ORF1a forms a six b strand core that projects extended loops outward, and includes an N-terminal b-strand, and a C-terminal a-helix (Figure 1a ) (Littler et al., 2020) . The N-and C-terminal regions are essential for Nsp9 dimer's replicase function (Figure 1b) (Littler et al., 2020) and are important for replication of CoV family viruses, including SARS-CoV (Frieman et al., 2012) . SARS-CoV Nsp9 has been demonstrated to bind RNA (Egloff et al., 2004) , stabilizing the nascent nucleic acid as it emerges from the nsp7-nsp8 complex during viral replication, and protecting it from nucleases Zhai et al., 2005) . In addition, Nsp9 may interact with a secondary structure element at the 3 0 end of the viral genome in several coronaviruses (Robertson et al., 2004) . The loops connecting the b2-3 (L23) and b3-4 (L34) strands of Nsp9 are positively charged and glycine rich, suggesting that these regions participate in RNA binding and interaction (Littler et al., 2020) . Soluble SARS-CoV Nsp9 dimerizes via interaction of an a-helical GxxxG motif, and any disruption in key conserved residues can impair RNA binding, consequently suppressing SARS-CoV viral replication (Frieman et al., 2012) . Mutations in the GxxxG motif of Nsp9 from porcine delta coronaviruses disrupts nucleotide binding capacity (Zeng et al., 2018) . Nsp9 replicase is considered an important SARS-CoV virulence factor (Miknis et al., 2009) . Since SARS-CoV-2 Nsp9 is 97% homologous to SARS-CoV Nsp9 (Littler et al., 2020; Ponnusamy et al., 2008) , this protein may play a crucial role in SARS-CoV-2 replication. It therefore represents a potential target for design and development of antiviral drugs.

Natural products isolated from plants used in traditional herbal remedies have had a huge impact on drug development against various diseases (Dias et al., 2012) . Naphthoquinones are natural metabolites found in plants, bacteria, and fungi that show extraordinary synthetic versatility, and are of great interest in medicinal chemistry (da Silva J unior et al., 2019) . Some naphthoquinones are currently used as medicinal drugs (Aminin & Polonika, 2020) . Lapachol is a 1,4-naphthoquinone compound originally extracted from the Bignoniaceae family, first isolated in Tabebuia impetiginosa, but also present in other plants (Hussain & Green, 2017 (Pinto et al., 1987; Sacau et al., 2003) . In addition, lapachol and its derivatives have been used in pre-clinical and clinical phases of drug development, demonstrating significant promise for cancer treatment (Block et al., 1974) .

Owing to widespread use of naphthoquinones, especially lapachol (I), b-lapachone (II), and a-lapachone (III) (Figure 2) , and their derivatives, there is great interest in studying their biological applications and their use as templates for synthesis of new derivatives. The chemistry of lapachol derivatives has been widely and extensively studied, and derivatives such as nor-b-lapachone (XII), hydroxylapachol (XI), and oxime derived from b-lapachol (IV) have been chemically synthesized Santos et al., 2016; Nasiri et al., 2013) . Therefore, lapachol derivative chemistry is a source of inspiration for the design and synthesis of new molecules, such as lapachol derivatives V, VI, VIII, IX and X ( Figure 2 ), which to our knowledge are novel, and currently under study in our laboratory. There is a dearth of knowledge on the antiviral activity of lapachol against coronaviruses, especially SARS-CoV-2. Therefore, our in silico data provide evidence that lapachol and its derivatives may serve as SARS-CoV-2 Nsp9 ligands, providing leads for development of candidate drugs against COVID-19.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) non-structural protein (nsp) 9 replicase protein structure (PDB:6W4B) (Figure 1 ) was used as a target to search for protein-ligand interactions with lapachol and its derivatives ( Figure 2 ). Virtual screening comprised prospecting for potential SARS-CoV-2 Nsp9 ligands from a library of 12 lapachol-related compounds. The GOLD 5.8.1 program (Jones et al., 1995) was used to calculate flexible anchoring between proteins and ligands. This program uses specific parameters for its virtual search of compounds in the selected databases. The GOLD program uses a genetic algorithm, propagating multiple copies of flexible ligand models in receptor active sites, randomly recombining segments of these copies until a set of structures converges. The position of the Nsp9 RNA-binding site was obtained from a study of SARS-CoV conducted by Sutton et al. (2004) , which defined the RNA binding site as composed of the L23, L45, and L7H1 loops ( Figure 1a ). These data were corroborated by virtual prediction using the metaPocket 2.0 server (Huang, 2009) .

The scoring function used for ranking was ChemPLP, which is defined by its use of a hydrogen bonding term, and multiple linear potentials to model van der Waals and repulsive terms. The parameter for searching an interaction's degrees of freedom was defined as "Very Flexible," corresponding to a search efficiency of 200%. All other parameters used options pre-defined by the program. All selected ligands were subjected to 10 iterations using this genetic algorithm. When the best three iterations of each molecule converged to the same docking pose, the five best compounds, according to PLP ranking and convergence, were selected for the next stage. Docking results were also analyzed by generation of 2 D diagrams using LigPlot þ (Laskowski & Swindells, 2011) . To evaluate potential toxic features, selected molecules were submitted to in silico toxicity prediction using the pkCSM webtool (Pires et al., 2015) .

The complexes selected from the results of molecular docking, and the Nsp9 apo protein were subjected to molecular dynamics (MD) simulations using GROMACS 2016 (Abraham et al., 2015) . Parameterization of each ligand was performed on the SwissParam server (Zoete et al., 2011) , and the MD of each protein-ligand complex was simulated using the CHARMM27 force field with the TIP3P water model. The MD unit cell was triclinic in shape, and water and ions were added to the system. System energy was minimized using the steepest descent algorithm, and minimization was stopped when the maximum force was less than 10.0 kJ/mol. NVT and NPT equilibrium phases were performed at 100 ps each. Solutions were subjected to 100 ns simulations, comprising 10 million simulation steps during the production phase under NTP conditions. Electrostatic interactions were calculated using particle mesh Ewald (PME) for long-range electrostatics (Essmann et al., 1995) . Temperature was maintained constant at 300 K using a modified Berendsen thermostat (Berendsen et al., 1984) , and pressure was maintained at 1 atm using the Parrinello-Rahman method (Parrinello & Rahman, 1980 ). Shortrange van der Waals and electrostatic forces were set to 1.2.

To assess ligand stability, molecular dynamics results were evaluated using tools provided by the GROMACS program. Analyses performed were as follows: root mean square deviation (RMSD), root mean square fluctuation (RMSF), and number of hydrogen bonds between protein and ligand.

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method (Massova & Kollman, 2000) with energy decomposition was used to calculate binding free energy between Nsp9 and each ligand. According to (Genheden & Ryde, 2015) crude approximations, it succeeds in rationalizing observed differences between docking studies. The MM-PBSA method was performed using the g_mmpbsa tool (g_mmpbsa-A GROMACS tool for high-throughput MM-PBSA calculations), which calculates three energetic terms: potential energy in vacuum, polarsolvation energy, and non-polar solvation energy.

Principal component analysis (PCA) or covariance analysis can be used in molecular dynamics to calculate overall protein molecular motion by analyzing a protein's dynamic nature and its atomic positional fluctuations (Aarthy et al., 2020) . The covar Gromacs tool was used to perform covariance analysis to create and diagonalize a covariance matrix using protein backbone atoms. The two most dominant eigenvectors were used to calculate the overlap between principal components and coordinates of the trajectory.

To prospect for ligands to non-structural protein (nsp) 9 of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) (Figure 1 ), we selected lapachol and 11 lapachol derivatives ( Figure 2 ) for in silico analysis. These compounds were subjected to molecular docking targeting the RNA binding site of Nsp9 composed of loops L23, L45, and L7H1 ( Figure 1a ) (Sutton et al., 2004) . The top five compounds (V, VI, VIII, IX, and XI) were chosen based on best ChemPLP scores, and on their capacity to converge to a single pose within three iterations after docking. Comparison of 3 D docking with an analysis of 2 D diagrammatic docking demonstrated that all five molecules formed a minimum of one hydrogen bond within the RNA binding site from Nsp9, together with hydrophobic interactions (Figures 3 and 4) . Molecule V formed two hydrogen bonds to Ser35 in the Nsp9 RNA binding site (2,86 and 2,56 Ð) and seven hydrophobic interactions (to Met8, Arg35, Val37, Phe52, Thr63, Ile61, and Lys54), producing a docking score of 54.52. Compound VI formed one hydrogen bond to Val37 (2,88 Å), and ten hydrophobic interactions (to Gly33, Gly34, Arg35, Ser55, Met8, Pro53, Phe52, Ile61, Phe36, and Thr63), with a docking score of 56.34. Compounds VIII and IX formed two hydrogen bonds to Arg35 (2.65 and 2.91 Ð) or one hydrogen bond (3.09 Ð), respectively. Compounds VIII and IX also established six hydrophobic interactions (to Ser55, Phe52, Ile61, Val37, Thr63, and Gly34), producing docking scores of 51.05 and 64.23, respectively. Molecule XI formed two hydrogen bonds to Arg35 (3,00 and 2,96 Ð) and one to Pro53 (2,65 Ð), and formed six hydrophobic interactions with residues Gly33, Gly34, Ile61, Val37, Phe52 and Ser55, with a docking score of 52.25 (Figure 4 ).

Toxicity predictions evaluated eight features related to potential toxic effects of the selected molecules (Table 1) . Except for molecule V, all compounds yielded negative results for AMES toxicity, which assesses potential mutagenicity. The maximum tolerated dose in humans gave the highest values for compounds VIII and IX and the lowest for compound V. Values below or equal to 0.477 log(mg/ kg/day) were considered low according to the PKCSM webserver. All molecules showed negative results for hERG inhibition. This feature evaluated the potential for potassium channel inhibition that could lead to fatal ventricular arrhythmia. All molecules except compound VIII displayed negative hepatotoxicity results. All molecules yielded negative skin sensitization results. T. pyriformis is used as a model in toxicity studies. All molecules displayed log(ug/L) values above the limit expected for low toxicity (-0.5 log(mg/L)). Minnow toxicity tests, using fathead fish minnows as a model, gave positive results in silico for all molecules. Expected high acute toxicity occurs at low compound concentrations, with log LC 50 scores below 0.3. In general, the molecules displayed little toxicity, except for hepatotoxicity and T. pyriformis toxicity. Molecule V gave the broadest range of toxicities, and molecule VIII displayed the lowest predicted toxicities. It is important to note that as yet, no experimental data corroborate or disagree with these predicted results, as these lapachol derivatives have not been tested in vitro or in vivo.

Lapachol derivatives form a highly stable complex with the SARS-CoV-2 Nsp9 RNA binding site

Complexes between the selected lapachol derivatives V, VI, VIII, IX, and XI and the Nsp9 RNA binding site were subjected to molecular dynamics assays to assess the stability of each complex via RMSD. During the 100 ns dynamic simulation, all complexes showed a stable profile ( Figure 5) , except for the complex with derivative XI. This was confirmed by mean RMSD and low standard deviation values (Table 2) . Complexes with lapachol derivatives V, VI, VIII, and IX displayed similar stabilities compared with the apo enzyme RMSD ( Figure 5 ) and by the mean ± standard deviation (Table 2 ). However, XI showed a slight variation in the last 25 ns compared with the apoenzyme. The ligand-binding site is mainly comprised of residues 25 to 65 in Nsp9 chain A ( Figure 6 ). Visual inspection of RMSF analysis ( Figure 6 ) reveals a similar profile of residue fluctuation in all analyzed structures. However, the RMSF mean and standard deviation fluctuate less in residues 25 to 65 in complexes with molecules V, VI, and VIII, compared to apo enzyme (Table 2) . These results indicate increased residue restraints, probably due to ligand interactions with these residues.

Analysis of hydrogen bonds formed within each complex showed that lapachol derivative VI (Figure 7b ) had the most stable hydrogen bond profile with a mean of 1.51 hydrogen bonds plotted along the MD (Table 3) . By contrast, complexes with molecules V and XI showed lower means and standard deviations. Similarly, the Nsp9-VIII and IX complexes did not demonstrate a tendency to stabilize, showing high standard deviations (Table 3 ) and sporadic increased values for hydrogen bonding (Figure 7c and d) .

Binding energy and stability of lapachol derivative/Nsp9 RNA binding site complexes according to MD simulations

The MM/PBSA method is an approach used to estimate free energies of small ligand complexes with macromolecules that requires a small number of computational processes. It has been useful for analyzing binding energies of docked structures and rationalizing observed differences in biological macromolecules' behaviors (Genheden & Ryde, 2015) .

Within each 100 ns simulation, 100 snapshots of the trajectory were obtained at constant intervals of 1 ns for analyses. All complexes showed negative binding energy profiles, with the exception of molecule V (Figure 8) . Molecule VI had the most stable profile (Figure 8) , and together with molecule IX, showed the lowest mean binding energy among the ligands (Table 4 ). Molecule IX also showed the lowest standard deviation. These results indicate that van der Waals interactions contribute the most energy to these protein-ligand interactions, and molecules VI and IX display the highest binding affinity.

A simple comparison of the dynamics over time of each ligand Nsp9 complex was performed. The protein-ligand interactions are illustrated by the complexes at 0, 50, and 100 ns in each simulation (Figure 9 ). Lapachol derivative VI displayed the most constant number of hydrophobic interactions and hydrogen bonds over the course of the simulation compared to the other ligands ( Figure 9b ). This result corroborates prior analysis defining VI as the most stable ligand.

PCA eigenvector projections represent the overlap between the first two principal components and coordinates of each trajectory. They plot the essential subspace occupied by each protein during molecular dynamics. The results show that complexes V, VI, and VIII filled conformational space similar to the apo protein ( Figure 10) . Thus, complex compactness was maintained during the simulation. By contrast, molecules IX and XI showed exploration of a different conformational space, and a less compact profile compared to the apo structure ( Figure 10 ).

To compare results with another coronavirus target, the severe acute respiratory syndrome-coronavirus (SARS-CoV) Nsp9 replicase protein (PDB:3EE7) was subjected to docking and dynamics with ligand VI, using the same parameters described in Methods. Lapachol derivative VI interacted with Ser59 from the Nsp9 RNA-binding site, and formed nine other hydrophobic interactions (to Met12, Pro57, Ser59, Arg39. Val41, Phe56, Thr67, Ile65, and Lys58), producing a docking score of 57.12 (Figure 11a ). Although the RMSD profile of SARS-COV Nsp9 does not seem to reach a stable plateau (Figure 11b ), the number of hydrogen bonds seems to increase from one bond to stabilize at four bonds over 100 ns of molecular dynamics (Figure 10c ). In addition, binding energy values were under À50 kJ/mol (Figure 11d ), similar to SARS-CoV-2 Nsp9-lapachol VI derivative complex values. 

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has threatened the global health system (Hsu et al., 2020) and accentuated the urgent need for drugs to combat viral outbreaks. The lack of antivirals and/or vaccines directly impacts our ability to bring pandemics under control, increasing morbidity/mortality, and dissemination of SARS-CoV-2/COVID-19. This situation also impacts populations' quality of life, overloading healthcare systems worldwide, and resulting in economic hardship. Therefore, drugs that abrogate SARS-CoV-2 replication are needed.

In this study, we docked ligands to the SARS-CoV-2 nonstructural protein (nsp) 9 RNA binding site (a promising drug target) to identify candidate antiviral agents. This methodology has been used to design and identify potential ligands for viruses including Chikungunya virus, dengue virus, and SARS-CoV. For emergency development of antiviral agents, studies have focused on identifying drug targets among essential proteins in virus replication, primarily targeting conserved structures and enzymes such as RNA-dependent-RNA-polymerase (RdRp), main 3-chymotrypsin-like cysteine protease (M PRO ), or the receptor-binding domain (RBD) of SARS-CoV-2 spike protein (Bharadwaj et al., 2020; Buonaguro et al., 2020; Rout et al., 2020) . However, many other targets are predicted to impact viral activity (Gordon et al., 2020; Wu et al., 2020) . Nsp9 is an essential enzyme in the replicasetranscriptase complex (RTC) required for coronavirus replication, that binds to genomic RNA, ensuring production of new RNA genomes, and protecting them from nucleases (Egloff et al., 2004, p. 9) . Additionally, it is highly conserved among coronaviruses (Littler et al., 2020; Ponnusamy et al., 2008) . Therefore, Nsp9 may represent a good target for direct-acting antivirals, given that inhibiting the RNA binding site of Nsp9 prevents it from binding to the nascent RNA from the nsp7-8 complex during SARS-CoV-2 replication. Once released from the nsp7-8 complex, the nascent RNA does not form a stable secondary structure Zhai et al., 2005) , and can be degraded by nucleases if Nsp9 action is blocked. Nsp9 has three loops (L23, L45, and L7H1) that have been associated with the RNA binding site and forms an OB-fold that is predicted to interact with oligonucleotides (Flynn & Zou, 2010; Sutton et al., 2004) . Moreover, the b-barrel in which the L23, L45, and L7H1 loops are found allows RNA binding capacity through nonspecific interactions (Sutton et al., 2004) . In this context, considering that SARS-CoV-2 Nsp9 is a conserved structure among coronaviruses (Littler et al., 2020) , we explored the RNA binding site formed by the cited loops as a possible target for lapachol and its derivates to assess their antiviral potential.

Lapachol and its derivatives have been associated with several biological activities in preclinical and clinical trials (Block et al., 1974; Hussain & Green, 2017) . Here, we evaluated lapachol, and eleven lapachol derivatives as possible ligands for the RNA binding site of SARS-CoV-2 Nsp9. Among the analyzed compounds, lapachol derivatives V, VI, VIII, IX, and XI demonstrated strong interactions with the target site in molecular docking studies.

Lapachol derivatives have been shown to exhibit improved activity and solubility when associated with new radicals (Oliveira et al., 2002 (Oliveira et al., , 2017 . Our data corroborate this literature, since chemical modifications of lapachol were able to increase hydrophobic interactions with regions of the RNA binding site, and create new interactions with specific amino acids, as demonstrated by lapachol derivative VI.

To assess the overall stability of each complex (lapachol or a lapachol derivative bound in the SARS-CoV-2 Nsp9 RNA binding site), molecular dynamics analyses were performed to obtain RMSD profiles. Lapachol has previously been described as very stable in molecular dynamics assays against ZIKV nsp3 (Oguntade et al., 2017) . However, in our analyses, lapachol did not yield the best docking score. Lapachol derivative VI presented one of the best docking scores, showing stable interaction behavior in RMSD and RMSF. In addition, molecular dynamics assessed hydrogen bond formation in each complex over the course of a 100 ns simulation. Hydrogen bonds are essential for stabilizing ligands within a binding pocket (Kostal, 2016) , and are mediators of strong interactions (Mar echal, 2007) . Lapachol derivative VI showed a tendency to stabilize existing hydrogen bonds, while the other compounds, especially compound IX, showed an increasing number of hydrogen bonds as the simulation progressed. When the lapachol derivative VI/SARS-CoV Nsp9 complex was compared to that of compound VI/SARS-CoV-2 Nsp9, the results showed similar interactive behavior. The main discrepancy in molecular dynamics performance was that the protein RMSD profile of the SARS-CoV Nsp9 protein did not stabilize over time. The hydrogen bonds formed, and the binding energy profiles were highly similar between SARS-CoV and SARS-CoV-2 Nsp9-lapachol derivative VI complexes. The high sequence identity between these two Nsp9 homologs (90%), obtained by the EMBOSS Needle sequence alignment tool (Needleman and Wunsch, 1970) , and structural similarity (RMSD ¼ 3.89 Å), obtained by structural alignment using the Chimera program (Pettersen et al., 2004 ) may explain these functional similarities ( Figure S1 ). Moreover, derivative VI bonding with SARS-CoV Nsp9 displayed similar affinity to that of SARS-CoV-2 Nsp9, suggesting broad-spectrum activity (Littler et al., 2020; Ponnusamy et al., 2008) . The affinity of these molecules for other non-structural SARS- CoV-2 proteins is being investigated by our research group. In addition, potential toxic characteristics can be corrected with other rational modifications during in vitro and in vivo tests.

In summary, we have demonstrated that SARS-CoV-2 Nsp9 may be a viable target for antiviral development. We also showed that lapachol derivatives, mainly VI and IX, may function as potential ligands for the SARS-CoV-2 Nsp9 RNAbinding site. Our work shows the importance of studying protein-ligand complex stability via dynamics to corroborate docking studies. Although the results seem promising, it is important to validate this activity. Therefore, these molecules should be further investigated for in vitro and/or in vivo antiviral activity. They may also be used as templates for development of future drugs against SARS-CoV-2 and other coronaviruses. Thus, these data may provide relevant information to advance our ability to combat COVID-19.

No potential conflict of interest was reported by the authors. 

Structural dynamic studies on identification of EGCG analogues for the inhibition of Human Papillomavirus E7

Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX, 1-2, 19-25

1,4-naphthoquinones: Some biological properties and application

The proximal origin of SARS-CoV-2

Efficacy of lapachol on treatment of cutaneous and visceral leishmaniasis

Structural insights into SARS coronavirus proteins

Molecular dynamics with coupling to an external bath

SARS-CoV-2 Mpro inhibitors: Identification of anti-SARS-CoV-2 Mpro compounds from FDA approved drugs

Early clinical studies with lapachol (NSC-11905)

SARS-CoV-2 RNA polymerase as target for antiviral therapy

Genomic characterization of the 2019 novel humanpathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan

Corona virus SARS-CoV-2 disease COVID-19: Infection, prevention and clinical advances of the prospective chemical drug therapeutics

Controlled Release of Norb-lapachone by PLGA Microparticles: A Strategy for Improving Cytotoxicity against Prostate Cancer Cells

Synthesis of quinones with highlighted biological applications: A critical update on the strategies towards bioactive compounds with emphasis on lapachones

Thiosemicarbazone of lapachol acts on cell membrane in Paracoccidioides brasiliensis

A historical overview of natural products in drug discovery

Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1

The severe acute respiratory syndrome-coronavirus replicative protein nsp9 is a single-stranded RNA-binding subunit unique in the RNA virus world

A smooth particle mesh Ewald method

Oligonucleotide/oligosaccharide-binding fold proteins: A growing family of genome guardians

Molecular determinants of severe acute respiratory syndrome coronavirus pathogenesis and virulence in young and aged mouse models of human disease

The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing

The novel coronavirus (SARS-CoV-2) epidemic

MetaPocket: A meta approach to improve protein ligand binding site prediction

Lapachol and lapachone analogs: A journey of two decades of patent research(1997-2016)

Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation

The architecture of SARS-CoV-2 transcriptome

Computational Chemistry in Predictive Toxicology: Status quo et quo vadis

Asymptomatic carrier state, acute respiratory disease, and pneumonia due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2): Facts and myths

Comparative analysis of RNA genomes of mouse hepatitis viruses

LigPlotþ: Multiple ligand-protein interaction diagrams for drug discovery

The hydrogen bond: Formation, thermodynamic properties, classification

Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding

Severe acute respiratory syndrome coronavirus nsp9 dimerization is essential for efficient viral growth

Design, Synthesis, and Biological Testing of Novel Naphthoquinones as Substrate-Based Inhibitors of the Quinol/Fumarate Reductase from Wolinella succinogenes

A general method applicable to the search for similarities in the amino acid sequence of two proteins

Characterizing the ligand-binding landscape of Zika NS3 helicase-promising lead compounds as potential inhibitors

Ruthenium (II)/triphenylphosphine complexes: An effective way to improve the cytotoxicity of lapachol

New enamine derivatives of lapachol and biological activity

Crystal structure and pair potentials: A molecular-dynamics study

UCSF Chimera-a visualization system for exploratory research and analysis

Antiviral activity of naphthoquinones. I. Lapachol derivatives against enteroviruses

pkCSM: Predicting small-molecule pharmacokinetic and toxicity properties using graphbased signatures

Variable oligomerization modes in coronavirus non-structural protein 9

The Structure of a rigorously conserved RNA element within the SARS virus genome

In silico investigation of spice molecules as potent inhibitor of SARS-CoV-2

Inhibitory effects of lapachol derivatives on epstein-barr virus activation

Cytotoxicity Activity of Semisynthetic Naphthoquinone-1-oximes against Cancer Cell Lines

The nsp9 replicase protein of SARS-Coronavirus, structure and functional insights

Novel Coronavirus (2019-nCoV) situation reports

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

Polyadenylate in the virion RNA of mouse hepatitis virus

Dimerization of coronavirus nsp9 with diverse modes enhances its nucleic acid binding affinity

Insights into SARS-CoV transcription and replication from the structure of the nsp7-nsp8 hexadecamer

A novel coronavirus from patients with pneumonia in China

SwissParam: A fast force field generation tool for small organic molecules