Alkaloids from Cryptolepis sanguinolenta as Potential Inhibitors of SARS-CoV-2 Viral Proteins: An In Silico Study

The ongoing global pandemic caused by the human coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has infected millions of people and claimed hundreds of thousands of lives. The absence of approved therapeutics to combat this disease threatens the health of all persons on earth and could cause catastrophic damage to society. New drugs are therefore urgently required to bring relief to people everywhere. In addition to repurposing existing drugs, natural products provide an interesting alternative due to their widespread use in all cultures of the world. In this study, alkaloids from Cryptolepis sanguinolenta have been investigated for their ability to inhibit two of the main proteins in SARS-CoV-2, the main protease and the RNA-dependent RNA polymerase, using in silico methods. Molecular docking was used to assess binding potential of the alkaloids to the viral proteins whereas molecular dynamics was used to evaluate stability of the binding event. The results of the study indicate that all 13 alkaloids bind strongly to the main protease and RNA-dependent RNA polymerase with binding energies ranging from -6.7 to -10.6 kcal/mol. In particular, cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine exhibited very strong inhibitory potential towards both proteins. Results from the molecular dynamics study revealed that a stable protein-ligand complex is formed upon binding. Alkaloids from Cryptolepis sanguinolenta therefore represent a promising class of compounds that could serve as lead compounds in the search for a cure for the corona virus disease.


Introduction
The novel human coronavirus was reported in China in late 2019. Ever since, the virus has spread throughout the world and has been designated a pandemic [1]. This novel human coronavirus has been named as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and the disease it causes is also called the corona virus disease 2019   [2,3]. As at May 21st, 2020, over 5 million cases have been recorded globally and the number of confirmed deaths is in excess of 300,000. The impact of the disease has been felt in all spheres of life, with devastating effects on the healthcare, social, and economic fabric of many countries.
As a ribonucleic acid (RNA) virus, SARS-CoV-2 is similar to the viruses that caused the Middle East respiratory syndrome (MERS) and the severe acute respiratory syndrome (SARS). In terms of biological classification, coronaviruses belong to the order Nidovirales and the family Coronaviridae [4]. There are four genera of the coronaviruses, namely, alpha-, beta-, gamma-, and delta-coronavirus [5]. SARS and MERS are both in the genus beta-coronavirus. These coronaviruses possess single strands of RNA and are enveloped in a protein with spike-like structures that project out of the surface of the envelope. SARS-CoV-2 uses angiotensinconverting enzyme 2 (ACE2) as a receptor during infection [6]. Infection, persistence, pathogenesis, and reemergence after viral maturation in hosts are mostly regulated by essential viral proteins.
Due to the crucial role they play in replication of a virus, proteases are usually earmarked as key targets in antiviral drug development [7]. The genome of SARS-CoV-2 is made up of about 30,000 nucleotides, and this RNA encodes polyproteins required for replication and transcription of the virus. When translated, the polyproteins are extensively processed via proteolysis. The main protease (M pro ) of SARS-CoV-2 processes the polypeptides by cleaving them at about 11 conserved locations to yield functional proteins. The functional 33.8 kDa M pro is a symmetrical dimer, with each unit made up of 3 domains. The first 2 domains (I and II) possess an antiparallel β-barrel structure whereas the third domain (III) is made up of 5 α-helices arranged into an antiparallel globular knot. Domain II is joined to domain III via an extended loop region. A Cys-His catalytic dyad is at the heart of the catalytic activities of the SARS-CoV-2 M pro , and the substrate binding pocket in the protease can be found in a cleft situated between domains I and II [8,9]. Humans lack closely homologous proteins to M pro , and this makes the SARS-CoV-2 M pro a potentially important target in antiviral therapy development [10].
The SARS-CoV-2 RNA-dependent RNA polymerase is also one attractive drug target. The RNA-dependent RNA polymerase utilizes metal cofactors to catalyze the formation of phosphodiester bonds between 2 ribonucleotides during RNA replication. The structure of the polymerase is characterized by the presence of 2 important domains: a polymerase or RdRp domain and a nidovirus RdRp nucleotidyltransferase (NiRAN) domain. Connecting these 2 is an interface domain. There is an N-terminal β-hairpin structure that is found as an insert in a groove clamped on either side by the NiRAN domain and the RdRp domain. The conserved polymerase motifs A-G in the palm domain of the RdRp constitute the active site pocket of the SARS-CoV-2 RNAdependent RNA polymerase [11]. The points of entry of the RNA template or primer and the nucleoside triphosphates, as well as the exit point of the emerging RNA strand, are all positively charged and solvent accessible. These points of entry, composed of motifs A and C through motifs F and G supported by E and the thumb, and exit points (the RNA exit tunnel), converge in a central cavity where template-directed RNA synthesis is carried out [12].
The rapid global spread of COVID-19 has underscored the need for the development of potent anti-COVID-19 therapeutics to combat this pandemic. Lessons from SARS and MERS, clinical drug repurposing, and in vitro high throughput screening have inspired novel insights towards the discovery of anti-COVID-19 drug candidates [13]. While these have suggested a few old drugs (such as remdesivir, lopinavir, hydroxychloroquine, and their azithromycin combinations) with new tricks against COVID-19, experimental techniques have also proposed macromolecular targets for attenuating viral replication [9,12]. The presence of high resolution structures of important viral proteins provides an avenue for their use in silico techniques such as molecular docking and molecular dynamics simulations to screen and evaluate potential inhibitors [7].
According to the World Health Organization (WHO), 65-80% of the world's population depends on herbal medicine in treating various diseases [14]. Herbal preparations and medicinal plants represent a potential source of therapeutics in this time of great need for antiviral agents that can help in fighting COVID-19. Cryptolepis sanguinolenta is a widely used plant in West African herbal medical practice. Extracts from the plant are used in treating ailments like diabetes, hypertension, malaria, respiratory diseases, and diarrhea [15][16][17][18][19]. Several alkaloids have been isolated from the plant, and these compounds and the plant extracts possess broad spectrum antipathogenic activity [17,18,20,21]. The plant extract is also used in managing hepatitis B viral infection and liver damage [18,22]. Available data from the literature suggests the effectiveness of Cryptolepis sanguinolenta extracts in interfering with viral replication of the herpes simplex virus type 1 [22]. The extensive use of the plant in folkloric viral therapy and the indication that it interferes with viral replication motivated us to evaluate compounds isolated from the plant as potential inhibitors of SARS-CoV-2 viral proteins.
This work examined alkaloids from Cryptolepis sanguinolenta as potential inhibitors of the SARS-CoV-2 main protease and RNA-dependent RNA polymerase using in silico techniques. The isolated alkaloids examined are quindoline, cryptospirolepine, cryptolepine, hydroxycryptolepine, neocryptolepine, cryptomisrine, cryptolepicarboline, 11-isopropylcryptolepine, cryptolepinone, biscryptolepine, isocryptolepine, cryptoheptine, and cryptoquindoline [17]. Molecular docking was used to estimate binding affinities of the alkaloids towards the proteins and determine important interactions that mediate binding whereas molecular dynamics simulations were used to assess stability of protein-ligand complexes. We herein report that the alkaloids of Cryptolepis sanguinolenta showed strong inhibitory potentials towards both the main protease and the RNA-dependent RNA polymerase and the association exhibited remarkable stability.

Molecular Docking and Visualization
2.2.1. Protein Preparation. The proteins were prepared by removing all complexed ligands and maintaining water molecules, especially those in the active sites where applicable. Replacement of incomplete side chains was done using the default Dunbrack rotamer library. Protonation states were assigned for histidine while addition of polar hydrogen atoms to correct the calculation of partial charges for standard residues was also done. Rotamer and backbone secondary structure preferences were improved with AMBER ff14SB, and then Gasteiger charges were computed for each atom using Antechamber implemented in Chimera [28,29].

Ligand Preparation.
A total of thirteen (13) alkaloids isolated from Cryptolepis sanguinolenta ( Figure 1) were used for the docking studies. The ligands were sketched in Spartan'14 (Wavefunction Inc., Irvine California, USA) Chem-Draw interface, modelled in 3D followed by SYBYL force field minimization as well as geometry optimization by 3 BioMed Research International equilibrium geometry estimation with the density functional theory (DFT) B3LYP/6-31G * basis set in vacuum. They were converted into pdb formats and prepared for docking by addition of polar hydrogens and Gasteiger charges. Additionally, curcumin, luteolin-7-glycoside, hydroxychloroquine, adenosine triphosphate (ATP), remdesivir, lopinavir, and nelfinavir were treated as above for validation purposes.
2.2.3. Molecular Docking. Precision docking protocol was developed by docking the bound ligand into the prepared active sites. In order to verify that the pose resulting from in silico docking represent correctly bound conformations, it was visually inspected and compared to the experimentally determined binding modes and conformations of N3 in 6LU7. The DINC 2.0 web server was used to dock N3 with M pro using estimated grid dimensions provided by the software since N3 is peptidic in nature [30]. The presence of conserved binding pocket interactions between N3 and the active site residues of M pro validated the docking protocol used in the study. In all cases, docking was performed in three technical runs. Free energies were obtained from runs with the most consistent poses. These free energies were used for the estimation of the ligand binding constant (K d ) and ligand efficiency (Δg) [31]. Protein-ligand complexes were analyzed in Discovery Studio for interactions with pocket atoms based on default settings in the software. Discovery Studio was used to visualize docked structures as well as the generation of 3D and 2D schematic representations of protein-ligand interactions.

Molecular Dynamics Simulation.
Ligands with binding affinities between -8.5 and -11.0 kcal/mol were considered for molecular dynamics (MD) simulations. In this regard, cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine were used for MD simulation studies of M pro whereas cryptomisrine, cryptospirolepine, and cryptoquindoline were used in MD simulation studies of RdRp and RdRpol. The Ligand and Receptor Molecular Dynamics (LARMD) web server was utilized for all MD simulations. The default MD simulation method on the LARMD web server was modified to suit the purpose of this study. The explicit water model was used in the simulation. Where necessary, Na + or Clwas added to obtain a neutral system. Trajectory analysis and free energy calculation protocols were used in the default mode as found on the LARMD web server [32].
Protein-ligand complexes were written from the best pose obtained from Vina. Ligand atoms were assigned AM1-BCC charges using the Antechamber module whereas the coordinate and topology files of the complex were con-structed with the teLeap module, all in the AMBER16 package [33,34]. The AMBER ff14SB force field [28] and GAFF (General AMBER Force Field) [29,35] were used for amino acid residues and ligands, respectively. The molecules were solvated in an octahedral box of TIP3P water [36] extended at least 10 Å in each direction from the solute [36,37]. Initially, all atoms were fixed except for water, ions, and hydrogens. Thereafter, only the backbone atoms of the protein were fixed. Finally, the residues around the ligand within 6 Å were minimized, and all the atoms were relaxed. The SANDER (Simulated Annealing with NMR Derived Energy Restraints) module in the AMBER16 program was utilized to perform the four-step minimization before the MD simulation. In all minimization processes, 2000 steps for steepest descent method and 3000 steps for the conjugated gradient method were used, followed by application of the Particle Mesh Ewald Molecular Dynamics module in the MD simulation [33]. The system was heated from 10 to 300 K in 30 ps. The subsequent release process was similar to the minimization. Finally, all the atoms were relaxed at 300 K and 1 atm by applying periodic boundary conditions and equilibration for 50 ps before the production run of 3.99 ns.

QSAR, ADME, and Toxicity Prediction.
Quantitative structure-to-activity relationship (QSAR) properties were estimated using Molsoft and the SwissADME web servers [38]. Toxicity predictions were obtained from the ADME-Tlab web server [39]. For the best four performing ligands (cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine), the parameters investigated included absorption, distribution, metabolism, elimination (ADME), and toxicity (T). ADME/T properties of lopinavir, omeprazole, and ibuprofen were computed for control purposes.

Results and Discussion
The search for therapeutic agents to combat the current COVID-19 pandemic is of utmost importance to all nations of the world. The absence of vaccines and potent antiviral agents against SARS-CoV-2 and effective treatment regimens amongst other factors has led to a significant increase in global mortality due to COVID-19. To date, the major strategy for identifying new COVID-19 drug candidates has been via drug repurposing, aided by high throughput screening of existing antiviral agents. Reliance on plant-derived natural products and formulations for the treatment of various diseases has proven beneficial globally. Extracts and compounds isolated from many of these plants have been shown to possess various pharmacological activities. Cryptolepis sanguinolenta is one such plant with antiviral, antiplasmodial, antiinflammatory, and hypotensive capabilities [17]. Due to its antiviral potential where it has been suggested that it interferes with the viral replication machinery of herpes simplex virus [22], and the fact that this plant is present in many approved herbal preparations in Ghana, we sought to evaluate the potential of the alkaloids isolated from the plant as potential inhibitors of proteins expressed by SARS-CoV-2 in in silico studies.

BioMed Research International
3.1. Protein Targets. The protein targets selected were the main protease and the RNA-dependent RNA polymerase primarily due to the key roles these proteins play in viral replication and sustainability as well as the availability of experimental structures in the protein database. Two different groups have recently solved the structures of the main protease of SARS-CoV-2 [7,8]. These 2 structures have been given the PDB codes 6LU7 and 6M0K. 6LU7 was chosen for this study because of its better resolution (1.5 Å) and the presence of a known in vitro inhibitor in the crystal structure. The presence of the ligand provides a simple route for the identification of binding pockets, selection of grid coordinates, and validation of docking protocols. The active site residues of this protease with a Cys-His catalytic dyad were identified. Two different structures were used for the RNA-dependent RNA polymerase. Initially, an experimental 3D structure of this protein was unavailable and was only obtained via homology modeling. A valid, high quality homology model of the RNA-dependent RNA polymerase was constructed by the SWISS-MODEL web server. The sequence similarity between RNA-dependent RNA polymerase and the template, 6NUR, on which the homology model was based was 96.35% [24,25]. The validity of the model is evidenced by the fact that for the Ramachandran plot, 100% of the residues were in the allowed regions and 97.5% were in the most favored region. The utility of the homology model of the RNA-dependent RNA polymerase (RdRp) has been presented in a number of molecular docking literature [40,41]. Recently, an electron microscopy structure of the RNA-dependent RNA polymerase (RdRpol) has been published [26]. Comparison of RdRpol to RdRp indicates high similarities within chain A (nsp12), and the root mean square deviation (RMSD) between 791 pruned atom pairs was 0.515 Å and that of 803 atom pairs was 0.933 Å (Fig.  S1 A) indicating over 96% similarity and about 3% dissimilarity. Despite the similarities in the binding pockets of both RdRp and RdRpol, both structures were chosen for docking studies for comparison.

Validation of Docking Protocols.
To ensure that the docking methods used in this study are suitable, a number of checks were employed. The SARS-CoV-2 main protease, M pro , was docked against N3. N3 is the ligand cocrystallized with M pro in the crystal structure. Comparison of its docked output with the native revealed comparable pocket interactions. Interactions of N3 with Met49, Thr190, Gln189, Glu166, Ala191, Leu141, His41, Met165, His41, His163, and Gly143 were identified in both conformations (Fig. S2) covering significant catalytic subdomains in the binding pocket. Ten (10) conventional hydrogen bonds (distances between 2.3 and 3.7 Å) and eight (8) hydrophobic interactions were observed in the docking output. The hydrophobic interactions included amide pi-stacked, alkyl, and pi-alkyl interactions. Hydrogen bonding interactions existed predominantly between amino acids and amide moieties in N3 with active site residues of M pro . However, the orientation of the ligand was inverted in the docking output when compared to the crystal structure. A binding score of -7.3 kcal/mol was recorded which was consistent to that obtained by Huynh and coworkers [42] (Table 1).
Both RNA-dependent RNA polymerase structures used, RdRp and RdRpol, had similar binding pockets. Since no bound ligand was present in the RdRpol structure, ATP and a nucleotide analog, remdesivir, were chosen as ligands for validation purposes. The binding affinity of ATP towards RdRp was computed to be -7.4 kcal/mol whereas remdesivir towards RdRp was -6.9 kcal/mol. Both remdesivir (-7.3 kcal/mol) and   Figure 2: (a) View of 3D interaction of cryptomisirine with M pro pocket residues (left with black labels) and 2D interactions colored by interaction type (right). (b) View of 3D interaction of cryptospirolepine with M pro pocket residues (left with black labels) and 2D interactions colored by interaction type explained in legend (right). 6 BioMed Research International ATP (-7.4 kcal/mol) had similar binding affinities towards RdRpol. Interactions observed in the docking of the two ligands with RdRp and RdRpol were also similar (Fig. S3). Finally, docking of some ligands reported in the literature against the main protease or the RNA-dependent RNA polymerase was also carried out. The binding affinities obtained from the docking of these ligands with M pro and RdRp and those reported were similar. Slight differences in results may be attributed to variation in software and computing power [43]. Ligand efficiencies were also found to be similar (Table 1). Hydroxychloroquine exhibited very low binding energies towards the main protease and the two RNAdependent RNA polymerases used in this study (-6.2 kcal/mol for M pro , -5.4 kcal/mol for RdRp, and -5.9 kcal/mol for RdRpol), as shown in Table 1.

Docking of Alkaloids from Cryptolepis sanguinolenta.
Thirteen alkaloids isolated from Cryptolepis sanguinolenta were docked against M pro , RdRp, and RdRpol to identify possible binding interactions between proteins and the alkaloids. Against M pro , the binding energies observed for all 13 alkaloids were between -6.9 and -10.6 kcal/mol (Table 2). Cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine possessed the best binding affinities, with binding energies less than -8.50 kcal/mol. In contrast, the binding energies of hydroxycryptolepine, cryptolepinone, neocryptolepine, isocryptolepine, quindoline, and cryptolepine were greater than -7.50 kcal/mol. The extent of association between the ligands and M pro was estimated using the dissociation constant, K d . Cryptomisrine and cryptospirolepine had dissociation constants less than 100 nM, implying  very strong association. Ligand efficiency, which represents the binding energy of the ligand per atom, was also calculated. Ligand efficiency of the ligands ranged between -0.26 and -0.41. All ligands interacted favorably and strongly between domains I and II of M pro where the cysteinehistidine catalytic residues exist [44]. This pocket was observed to be hydrophobic with solvent accessible regions existing at the periphery. Of the best four ligands that interacted well with M pro , cryptomisrine and cryptospirolepine interacted with both Cys145 and His41. The nature of the interactions included pi-alkyl, pi-pi T-shaped pi-sulfur, and pi-pi stacked hydrophobic interactions with distances between 4 and 5 Å (Figure 2). Hydrogen bonding was observed for the interaction between the quinoline nitrogen of cryptomisrine and Met165 of M pro at a distance of 2.88 Å and at angles 116.352°and 102.64°(D-H-A and H-A-Y, respectively) leading to a resultant positively charged nitrogen (Figure 2). For cryptospirolepine, a conventional hydrogen bonding was found between the carbonyl oxygen and Gly143 at a distance of 2.41 Å ( Figure 2) and angles 142.895°and 124.549°(DHA and HAY, respectively) whereas that with Asn142 was that of a carbon hydrogen bond at a distance of 2.17 Å and angles 131.604°and 138.885°(D-H-A and H-A-Y, respectively). Cryptoquindoline and biscryptolepine interacted with Met165 and Met49 via pi-alkyl hydrophobic forms at a distance between 4 and 5 Å whereas interactions with Glu166 were of a pi-sulfur and pi-anion kind. A conventional hydrogen bond was also found between Gln189 and biscryptolepine indole NH, at a distance of 1.91 Å (Fig. S4) Figure 4: (a) View of 3D interaction of cryptoquindoline with RdRpol pocket residues (left with black labels) and 2D interactions colored by interaction type (right). (b) View of 3D interaction of cryptomisrine with RdRpol pocket residues (left with black labels) and 2D interactions colored by interaction type explained in legend (right). 8 BioMed Research International and Thr190 at distances between 2.13 and 2.8 Å (Fig. S4). All remaining ligands exhibited strong interactions with Cys145 and/or His41 as well as other residues in the various protease subdomains (Fig. S2). When compared to other antiviral compounds suggested as possible lead compounds for COVID-19, the alkaloids of Cryptolepis sanguinolenta had noteworthy binding energies. For example, the binding energies of lopinavir, nelfinavir, and hydroxychloroquine when docked against M pro were -8.7, -8.3, and -6.2 kcal/mol, respectively (Table 1). Low binding energies were associated with the binding of all the Cryptolepis sanguinolenta alkaloids against M pro . Of the 13 compounds tested, only cryptolepine had a binding energy greater than -7.0 kcal/mol ( Table 2). The interaction of all alkaloids with at least one of the catalytic residues of M pro suggests a probable mode of inhibition where these active site residues required for proteolysis are unavailable in the presence of the substrate. Recently, Gyebi and coworkers docked various alkaloids against the main protease of SARS-CoV-2, SARS-CoV, and MERS-CoV [45]. Interestingly, some alkaloids from Cryptolepis sanguinolenta exhibited strong binding to SARS-CoV-2 M pro , similar to the results obtained in this work. Estimated binding affinities were similar to that obtained in this work. Cryptospirolepine, with a binding affinity of -9.2 kcal/mol, was found to interact with Cys145. This is in agreement to the results of this work where cryptospirolepine interacted with both amino acids in the catalytic dyad-Cys145 and His41-with a binding affinity of -9.5 kcal/mol.
In general, all the alkaloids bind to essential domains in the active sites of the two structures of the RNA-dependent RNA polymerase (Fig. S1 B and C). The differences in binding energy of alkaloids with both RdRp and RdRpol were about 0.5 kcal/mol. Of note, the binding energies of cryptolepicarboline, 11-isopropylcryptolepine, cryptoheptine, hydroxycryptolepine, cryptolepinone, neocryptolepine, and isocryptolepine with both RdRp and RdRpol were exactly the same (Table 2). Cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine alkaloids exhibited strong binding to RdRp and RdRpol, with binding energies less than -8.50 kcal/mol. The binding constants (K d ) for cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine were well below 1 μM for both RdRp and RdRpol, indicating strong association between ligands and protein targets. Cryptomisrine and cryptospirolepine exhibited much stronger binding to M pro in contrast to both RdRp and RdRpol. Ligand efficiencies for all 13 alkaloids were very similar for M pro , RdRp, and RdRpol (Table 2). Electrostatic and hydrophobic interactions were the driving force of binding, in general. A positive charge on an aryl nitrogen of cryptomisrine interacted with the negative charge of an oxygen in Asp623 at a distance of 4.5 Å (Figure 3). The guanidino groups of Arg553 and Arg624 were involved in pi-cation interactions with pi-orbitals in cryptomisrine. Other common interactions observed were pi-anion, pi-pi T-shaped, and pi-alkyl interactions. Overall, 13 nonbonding interactions were obtained for cryptospirolepine and cryptoquindoline towards RdRp whereas 11 were found for biscryptolepine. Most of the interactions between alkaloids and RdRpol ( Figure 4) were similar to those of RdRp, with residues in the palm, interface, and thumb domains (Fig. S1).
Elfiky recently described ribavirin, remdesivir, sofosbuvir, galidesivir, and tenofovir as potentially potent inhibitors of the RNA-dependent RNA polymerase of SARS-CoV-2 as they exhibited binding energies of -7.0 to -8.0 kcal/mol and suggested that the tight binding may inhibit polymerase function [40]. In a similar molecular docking study by Shah and coworkers, binding energies between -6.0 and -9.0 kcal/mol were obtained by docking a series of antiviral agents  [41]. In comparison to both works, the alkaloids studied in this work showed much tighter binding as evidenced by the low binding energies. In particular, the high binding affinities of cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine towards both RdRp and RdRpol makes them potential inhibitors that can be further explored.

Molecular Dynamics.
Ligand-driven molecular dynamics simulations were carried out to explore the stability of the interactions between ligands and proteins. The explicit water solvation model was applied to the four best ligands of M pro and three best ligands of RdRp and RdRpol.The ligand binding free energies were estimated using both the Poisson-Boltzmann and the Generalized Born surface area continuum solvation, and the observations were summarized in Table 3 (detailed energetic contributions for binding free energy computations can be found in Table S1). For M pro , ligands investigated were cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine. The root mean square deviation (RMSD) of ligands in the M pro -ligand complexes was between 0.39 and 0.62 Å and that of the protein (M Pro ) was in the range of 1.43-1.95 Å. The radius of gyration for all M pro -ligand complexes was about 22 Å. For both RdRp and RdRpol, the ligands involved in the protein-ligand complexes investigated were cryptomisrine, cryptospirolepine, and cryptoquindoline. The ligand and protein RMSD for the RdRp complexes ranged between 0.30 and 0.69 Å and 1.80 and 1.87 Å, respectively. For the RdRpol complexes, the ligand and protein RMSD ranged between 0.42-0.51 Å and 2.94-3.96 Å, respectively. It is interesting to note that protein RMSD was much greater in RdRpol complexes (>2 Å) than in both RdRp and M pro complexes ( Table 3). The binding pocket of M pro is characterized by tight structural packing and is probably responsible for observation of a largely stable protein during the entire simulation timescale. The presence of a metal cofactor in RdRp (which was absent in RdRpol) may also have contributed to its conformational stability shown in minimal overall configurational entropy. The radius of gyration for RdRp and RdRpol complexes was about 29 Å  D-score: drug-likeness model score.
All alkaloids violated none or only one of Lipinski's rules. Cryptomisirine donates the most hydrogen bonds followed by biscryptolepine and hydroxycryptolepine. Cryptoheptine had the ability to accept the most hydrogen bonds followed by cryptoquindoline and quindoline. With a similar molecular landscape, all 13 alkaloids were largely hydrophobic, with most being moderately soluble to poorly soluble. Only hydroxycryptolepine was completely soluble. All ligands have the ability to cross the blood-brain barrier with most having high gastrointestinal absorption indices. 11 BioMed Research International and 32 Å, respectively. However, the radius of gyration for the cryptospirolepine-RdRpol complex was much higher, at 38.65 Å. Entropic estimation for experimental structures were positive in M pro (13-22 kcal/mol) and RdRpol (15 to 24 kcal/mol), except for those of RdRp (Table S1>). More tightly bound ligands are expected to have high entropic costs than loosely bound ligands [46]. All alkaloids were tightly bound to respective protein targets with overall high entropic costs as a result of high hydrophobic free energy contributions to total binding free energies (ΔG/PBSA and GBSA). Free energy contribution arising from electrostatic interactions with pocket residues was less in contrast to that of remdesivir triphosphate (RemTP) affecting the end point free energies (89.22 and 32.7 kcal/mol, Table S1). The nature and kinds of interactions that mediated protein ligand binding were largely similar to those obtained in the docking experiments (Table 3). Overall, the protein ligand complexes probed in the molecular dynamics simulations were stable with enhanced ligand binding efficiencies as well as observed minimal protein conformational fluctuations over the experimental timescale ( Figure 5).
3.5. QSAR, ADME, and Toxicity. The QSAR and ADMET property and probability predictions were made to assess the key drug-like qualities of some of the alkaloids used. Cryptomisrine had the most hydrogen bond donors followed by biscryptolepine and hydroxycryptolepine. Cryptoheptine had the ability to accept the most hydrogen bonds followed by cryptoquindoline and quindoline. With a similar molecular landscape, all 13 alkaloids were largely hydrophobic, with most being moderately soluble to poorly soluble. Only hyroxycryptolepine was predicted to be completely soluble. All ligands have the ability to cross the blood-brain barrier with most having high gastrointestinal absorption indices (Table 4). In the QSAR analyses, all alkaloids investigated violated either one or none of Lipinski's rules and hence were drug-like. Lopinavir, omeprazole, and ibuprofen were used for comparison in ADMET probabilities and predictions (Table 4 and Tables S2 and S3)). Despite the low aqueous solubility of the four alkaloids studied, their permeability as well as metabolism was relatively high compared to that of lopinavir. The volume distribution of the four (0.291 L/kg, 0.59 L/kg, 0.78 L/kg, and 0.625 L/kg) indicates a better distribution. The predictions show low half-life times indicating good bioavailability with relatively higher clearance rates (see Table S1). Drug-likeness model score predicted values between -1 and -0.2, on a scale of -6 (largely for non-drugs) to +6 (highest score for drugs) [47]. The majority, however, had scores from -0.20 to -0.5, indicating that these compounds are fairly drug-like [47,48] confirming Lipinski's predictions. LD 50 's greater than 500 mg/kg with the exception of cryptospirolepine (LD 50 ; 184.902 mg/kg) were obtained for the best four tight binding alkaloids assessed. These three alkaloids were less toxic than lopinavir (LD 50 ; 570.85 mg/kg). Overall, these predictions indicate that these molecules taken together would be safe. Cryptoquindoline was predicted to be strongly mutagenic (probability 0.8), cryptomisrine was moderate (probability, 0.6) while in contrast, cryptospirolepine and biscrypto-lepine were not potential mutagens or carcinogens (Table S2). Toxicity studies of the extracts and alkaloids of Cryptolepis sanguinolenta are well documented in the literature [49][50][51][52], and have led to investigations into their use as potential anticancer agents [52,53]. It has been shown that the plant extract when administered at very high concentrations of 2000 mg/kg may result in some degree of toxicity but is safe at concentrations at or below 500 mg/kg [17,51]. Investigations into tolerable dosages in animal models could be used to confirm our findings on toxicity, mutagenicity, and carcinogenicity.

Conclusion
The use of natural plant-based extracts as antiviral agents is well established in the literature. One advantage of using plant extracts for therapeutic purposes is the potential for synergism by the compounds of that extract. As shown in this molecular docking study, the different alkaloids of Cryptolepis sanguinolenta have exhibited high binding affinity and hence potential inhibitory activity towards two of the major prions of SARS-CoV-2, the main protease and the RNAdependent RNA polymerase. The high binding affinity of cryptomisrine, cryptospirolepine, cryptoquindoline, and biscryptolepine for both protein targets, drug-like characteristics, tight binding, and residence times within protein binding pockets makes them ideal candidates for further in vitro and in vivo validation.

Data Availability
All data generated or analyzed during this study are available upon reasonable request.

Conflicts of Interest
All authors declare no competing financial, professional, or personal interests that might have influenced the performance or presentation of the work described in this manuscript.