In silico evaluation of potential inhibitory activity of remdesivir, favipiravir, ribavirin and galidesivir active forms on SARS-CoV-2 RNA polymerase.

Since the outbreak emerged in November 2019, no effective drug has yet been found against SARS-CoV-2. Repositioning studies of existing drug molecules or candidates are gaining in overcoming COVID-19. Antiviral drugs such as remdesivir, favipiravir, ribavirin, and galidesivir act by inhibiting the vital RNA polymerase of SARS-CoV-2. The importance of in silico studies in repurposing drug research is gradually increasing during the COVID-19 process. The present study found that especially ribavirin triphosphate and galidesivir triphosphate active metabolites had a higher affinity for SARS-CoV-2 RNA polymerase than ATP by molecular docking. With the Molecular Dynamics simulation, we have observed that these compounds increase the complex's stability and validate the molecular docking results. We also explained that the interaction of RNA polymerase inhibitors with Mg++, which is in the structure of NSP12, is essential and necessary to interact with the RNA strand. In vitro and clinical studies on these two molecules need to be increased.

Introduction

Coronavirus disease 2019 (COVID-19) is a viral infection with high pathogenicity and contagiousness caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). First, cases of pneumonia of unknown etiology were reported on November 17, 2019, in Wuhan, Hubei Province, China [1]. On January 7, 2020, it was identified that the disease agent was a new coronavirus (2019-nCoV) that had not been detected in humans before. The was later named coronavirus disease-19 (COVID-19). In one study examining the 2019 new CoV genome, this new virus was found to have 89% nucleotide homology with bat SARS-like CoVZXC21 and 82% with human SARS-CoV BJ01 2003. Hence, it was named SARS-CoV-2 because of its close resemblance to SARS-CoV [2]. The co-crystal structure of 2019-nCoV protease, 6LU7, has ∼99% identity with SARS-CoV protease [3]. SARS-CoV-2 is a single-stranded, positive polarity enveloped RNA virus. COVID-19 spread rapidly to many countries and was officially declared a pandemic by the World Health Organization on March 11, 2020, having caused the deaths of more than 4000 people [4, 5]. The disease is highly contagious, and at the onset of the disease, the main symptoms were fatigue, fever, dry cough, myalgia, and dyspnea. Less common symptoms were nasal congestion, headache, runny nose, sore throat, vomiting, and diarrhea. The virus can be found in patients' respiratory secretions 1–2 days before the onset of clinical symptoms and two weeks after the disease symptoms. The incubation period of the disease is between 4 and 6 days. Severe cases usually have dyspnea and/or hypoxemia one week after onset and then go on to septic shock, followed by Acute Respiratory Distress Syndrome (ARDS) [6-8]. It has been shown that this disease is transmitted from person to person, especially adults, who are susceptible to COVID-19, and the severity of the disease is related to age. It has also been shown that the disease is more severe in people with comorbidities such as hypertension, diabetes, and cardiovascular disease [9]. The fact that infected people can infect others without showing any symptoms also makes it difficult to control COVID-19.

Coronaviruses are 80–220 nm in diameter and contain a single-stranded RNA of approximately 30 kbs in length as nucleic acid [10]. A study using the genetic sequencing method reported a 96% similarity between the gene sequence of SARS-CoV-2 and bat coronavirus [11]. Coronaviruses attach to host cells via the spike (S) protein on the outer surface and enter the cell; S protein recognizes the receptor in the target cell and regulates the entry of the virus into the host cell. The virus's life cycle begins with binding the S protein to the angiotensin-converting enzyme 2 (ACE2) receptor on the host cell surface. It has been shown that SARS-CoV-2 also enters the host cell using similar mechanisms after binding to the ACE2 receptor [12, 13]. Structural studies revealed that SARS-CoV-2 spike (S) glycoprotein binds ACE2 with higher affinity (~ 15–40 nM) than SARS-CoV spike protein [14, 15]. Details of the pathological mechanisms leading to multisystemic organ dysfunction in SARS-CoV-2 infection are as yet unknown [12, 16]. The treatments currently used in COVID-19 are the treatments applied to prevent the virus's entry into the cell—viz. inhibit or reduce its replication, and suppress the increased inflammation response in line with this physiopathogenesis [17]. Also, convalescent plasma treatment, which includes infected and recovered patients' antibodies, is among the options [18].

Coronaviruses are one of the few RNA viruses that have a genomic regulation mechanism. The development of nucleoside analogs effective against coronaviruses has been particularly difficult due to a unique exoribonuclease (ExoN). ExoN functions as a corrective enzyme that corrects errors in the growing RNA chain [19]. The RdRp (RNA-dependent RNA polymerase) nonstructural protein 12 (NSP12) is highly conserved among coronaviruses, making it an attractive target for broad-spectrum antiviral drugs. The RdRp enzyme allows the viral genome to be copied into new RNA copies using the host cell's mechanism. Those in adenine or guanine analogs target the RdRp and block viral RNA synthesis across a broad spectrum of RNA viruses, including human coronaviruses. Generally, the rate-limiting step for activation of nucleoside analogs is the production of the nucleoside monophosphate. Nucleoside phosphoramidites such as remdesivir, favipiravir, ribavirin are bioisosteres of monophosphates and bypass this rate-limiting step. However, nucleoside phosphoramidites should be administered as prodrugs to mask the charged phosphonate group and allow faster cell entry. The nucleoside analog inhibits the viral RNA-dependent RNA polymerase (RdRp) by competing with the usual counterpart ATP, GTP, UTP, or CTP [20, 21]. Therefore, favipiravir, remdesivir, ribavirin, and galidesivir are thought to be potential drugs against SARS-CoV-2 [22-25].

Prodrugs are inactive drug forms designed by chemical changes in the original drug to increase the utilization of drugs in the human body. It is called a prodrug because after it is taken into the body by any means, it metabolizes kidney, liver, stomach acid, intestinal alkalinity (basicity) and changes its chemical structure. It creates a drug with a different structure, reaching the target tissue and having therapeutic properties. Even if prodrug drugs enter the blood for any reason, if they do not change, they do not show a pharmacological effect and have no therapeutic effect [26]. By way of example, Oseltamivir phosphate is a prodrug that is rapidly and extensively hydrolyzed in vivo to its active metabolite, oseltamivir carboxylate, which is a potent and selective inhibitor of influenza A and B virus neuraminidase [27]. Likewise, ribavirin, galidesivir, favipiravir, and remdesivir are in prodrug form and are converted into active metabolites after being taken into the body.

Ribavirin (tribavirin) is a broad-spectrum guanosine analog antiviral prodrug widely used to treat hepatitis C, RSV (respiratory syncytial virus) infection, Crimean-Congo hemorrhagic fever, and hantavirus. For the effect of SARS-CoV-2, high doses of 1.2 to 1.4 g orally every 8 h are required, which causes serious side effects in the patient in terms of hematological and hepatotoxic side effects potential. Inhalation administration also showed no advantage over oral or intravenous administration. It is also not recommended for pregnant women due to its teratogenic activity [28]. Clinical studies have shown that ribavirin can work synergistically with interferon β and inhibit SARS-CoV replication (Fig. 1). The poor in vitro activity of ribavirin against coronaviruses has been attributed to its removal by ExoN [29].

Fig. 1

The entry into the cell of the antiviral drugs ribavirin and galidesivir and the formation reaction of the active triphosphate forms

Galidesivir (BCX4430) is a broad-spectrum antiviral prodrug by BioCryst, a pharmaceutical company based in Durham that has potential in the treatment of COVID-19 and is safe and well-tolerated in Phase 1 trials in healthy subjects [30, 31]. Phase 1 clinical safety and pharmacokinetic studies of galidesivir by both intravenous and intramuscular administration have been completed in healthy volunteers. Currently, Galidesivir is under phase 2 human trial for Coronavirus in Brazil and around the world. Galidesivir has demonstrated broad-spectrum activity in vitro with EC50 ranging from ~ 3 to ~ 68 μM in preclinical studies against more than 20 RNA viruses in nine different families—filoviruses, togaviruses, bunyaviruses, arenaviruses, paramyxoviruses, flaviviruses, coronaviruses [31, 32]. Initially developed for hepatitis C virus, it phosphorylates Galidesivir, an adenosine nucleoside analog that inhibits viral RNA polymerase into a triphosphate that mimics the cellular kinases ATP. Viral RNA polymerases incorporate the drug's monophosphate nucleotide into the growing RNA chain, resulting in early chain termination, thereby disrupting the activity of viral RNA-dependent RNA polymerase [33, 34] (Fig. 1).

Favipiravir is a broad-spectrum prodrug that is phosphoribosylated inside the cell to form the active metabolite favipiravir-ribofuranosyl-5′-triphosphate. T-705-RTP competes with purine nucleosides and interferes with viral replication by incorporating into the viral RNA, potentially inhibiting RNA-dependent RNA polymerase (RdRp) of RNA viruses [20, 22, 35, 36]. It affects the RdRp of the influenza virus with a semi-maximal inhibitory concentration of 0.022 μg/mL. Still, it does not affect the human DNA polymerase α, β, and γ subunits up to 100 μg/mL. In some clinical studies conducted in the patient group receiving favipiravir, shorter viral clearance time, a higher recovery rate in chest imaging significantly reduced fever reduction and cough relief time in ordinary COVID-19 patients and patients with hypertension/or diabetes. Immediate approval of favipiravir (formulation: tablet, 0.2 g) for a clinical trial in adult patients with novel coronavirus pneumonia was announced by the National Medicines Administration (NMDA) in China in March 2020 because it has been observed to exhibit good overall efficacy and safety with minimal side effects [37].

Remdesivir (RDV, formerly GS-5734) is a monophosphoramidate nucleoside analog prodrug with broad-spectrum action against various DNA and RNA viruses initially developed in response to the 2014–2016 Ebola epidemic in West Africa [3, 19, 24, 28]. RDV changes from its pro-form to the active triphosphate form that inhibits RdRp by non-essential RNA chain termination. The newly formed structure is an adenosine analog attached to viral RNA chains and leads to premature termination [23, 38–40]. The exoribonuclease of the virus, which usually reads and corrects replication errors, cannot work against the active form of remdesivir. Remdesivir is more active against viruses lacking ExoN but may partially escape correction and maintain strong antiviral activity in the presence of ExoN [19]. Recently, the FDA granted an emergency use permit for Remdesivir for potential COVID-19 treatment [41]. Because the drug is highly selective for the virus's RNA polymerase, it is expected that humans' toxic side effects will be low. Thanks to its long half-life, it provides ease of use once a day. The therapeutic dose used for COVID-19; in the form of intravenous administration of 200 mg/day on the first day and 100 mg/day on the following days for 10 days. When the remdesivir treatment results are evaluated with clinical studies, 53 cases (patients with 94% oxygen saturation or receiving oxygen support) were followed for at least 28 days using remdesivir at the dose mentioned above. While 68% of the cases showed improvement in terms of oxygen support, 15% got worse. Of 30 intubated patients, 57% were extubated, and ECMO (extracorporeal membrane oxygenation) application could be terminated in three of four patients connected to ECMO (extracorporeal membrane oxygenation). After 28 days of follow up, the cumulative incidence of clinical improvement was evaluated as 84%. Mortality was 13% after treatment was completed. Serious side effects were encountered in 23% of the cases. The most common side effects; ALT elevation was noted as diarrhea, rash, renal dysfunction, and hypotension [42].

In this study, we aimed to analyze the interaction pathway of the active forms of ribavirin, galidesivir, favipiravir, and remdesivir prodrug antiviral drugs in the SARS-Cov-2 RNA polymerase active site by molecular docking and dynamic studies. Unlike many other computational studies, we carried out our study by considering the interaction of active forms of drug molecule candidates with both NSP12 and RNA strand. Unfortunately, many in silico studies have focused on inactive prodrug forms of drug candidates [43-47], and these studies are unlikely to be successful. Likewise, in many computational studies, classical discourse includes the statement that other structures other than NSP12 (NSP7, NSP8, and RNA strand) were deleted [46, 48]. Still, RNA polymerase inhibitors must interact with both the NSP12 and the RNA strand for activity.

Material method

Homology modeling

In the NSP12-NSP7-NSP8 complex (PDB: 7BV2), missing residues 51–83, 101–117, and 896–910 were detected in the nsp12 chain and reconstructed according to homology modeling. The missing residues were added with the homology modeling tool included in the Chimera software regarding the protein sequence in the 7BV2 PDB file [49, 50].

Molecular docking

Molecular docking studies were carried out with the Schrödinger Maestro module [51-53]. Firstly, the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) enzyme (PDB: 7BV2) crystal structure was imported from the protein data bank (http://www.rcsb.org/pdb/) into the 'Protein Preparation Wizard' module. Hydrogens were added, non-bonding command with metals, the formation of disulfide bonds, deleting water at 5 Å distance from het groups, and preprocess by creating pH: 7.00 ± 2.00 het states using Epik. Subsequently, water molecules and Zn++ and Mg++ metals contained in protein crystal structures were retained, and remdesivir monophosphate and POP molecules other than the SARS-CoV-2 NSP7, NSP8, NSP12 RNA polymerase protein, and RNA structure were deleted. Finally, H bond determination was optimized using PROKA pH: 7.00 with water sample orientation and the complex structure was prepared by minimizing OPLS3 force fields.

Galidisevir, ribavirin, remdesivir, favipiravir drug molecules, and ATP were downloaded in 3D SDF file format from PubChem (https://pubchem.ncbi.nlm.nih.gov/). Those prodrugs' metabolite was drawn, minimized, and saved as an SDF file using Chem3D 17.01. Formed compounds were entered into the 'Ligprep' module. The OPLS3 force field was preferred, and ionization was carried out using Epik in the range of possible pH: 7.00 ± 2.00, desalinated, and tautomer-formed. Chirality from the 3D structure of ligands was determined and prepared. The complex structure's active site coordinates were prepared based on the remdesivir monophosphate in the 7BV2 structure in the size x: 91.517, y: 92.382, z: 103.734, and 25 × 25 × 25 Å3 using the 'Receptor Grid Generation' module. Finally, the calculation of theoretical ligand–protein-RNA interactions was performed with Extra Precision (XP) using flexible ligand options of the 'Glide Ligand Docking' module. The Glide score and Glide emodel were evaluated, and 2D/3D interactions of ligand and protein were determined and exhibited.

Molecular dynamics simulation

The Glide docking tool's conformation was obtained in the Schrodinger suite to be located in the protein's active center. Compounds were extracted with Chimera software (Pettersen et al. 2004). Ligand topology parameters were obtained in accordance with the Amber force field using the AM1-BCC semi-empirical quantum calculation method with ACPYPE software [54, 55].

All simulations were run with GROMACS 2020.3 software [56]. The "Leap Frog" integration was used, and 2 fs time steps were applied. The system was prepared in periodic boundary conditions (PBC) as "Rhombic Dodecahedron." Cube sizes were calculated at a distance of 1.2 nm from the corner of the protein–ligand complex. Amber99SB-ildn force field and "TIP3P" water model were preferred [57]. The system was neutralized with 0.15 mM Na–Cl. Energy minimization was performed by applying the 50,000 step Steepest Descent algorithm. To bring the system to equilibrium, 300 ps NVT, and 1 ns NPT simulations were performed. According to the NVT stage, all bonds and atomic positions were restricted at 1000 kJ/mol force in the LINCS algorithm. In NPT simulation, only bonds were fixed with LINCS Restriction Algorithm. All simulations were performed at 2 fs time intervals. Verlet was used as the cut-off scheme during the production phase, and the Nose–Hoover algorithm was used as the temperature coupling algorithm. The temperature was set to 300 K. The Perinello-Rahman algorithm with isothermal compressibility was preferred as the pressure coupling algorithm, and the pressure was set to one atmosphere. All restrictions were removed, and a 50 ns long MD simulation was performed for each protein compound complex. "Particle-mesh-Ewald" was used as a long-range electrostatic algorithm. Cut-off of Van Der Waals and short-term electrostatic interactions were set at 10 Å.

Results

The crystal structure of SARS-CoV-2 RNA polymerase consists of NSP12 (A), NSP8 (B), NSP7 (C) chains, and RNA strand. There are two-piece Zn++ ions and two-piece Mg++ metal ions in the NSP12 structure. ATP, which activates RNA polymerase, and interacts with Mg++ (1004 numbered) by forming a salt bridge [58]. As the RNA-NSP12 complex and remdesivir interaction were examined, it was observed that the first phosphate group of remdesivir created an attractive charge with Mg++ 1004 and 1005. While it was seen that the 20th uracil of RNA formed a covalent bond with the same phosphate group, it was observed that the same uracil base formed a pi–pi interaction with the bicyclic ring. It was observed that while Arg555 interacted with the bicyclic structure with an attractive charge, the 10th uracil of RNA formed hydrogen bonds with the nitrogen atoms of the hexane ring of the bicyclic structure. Also, it was observed that the remdesivir and the Val557, Cys622, Asp623, Ser682, Ala688, Thr687, Asn691, and Asp760 residues contacted with Van Der Waals interactions. According to molecular docking results, galidesivir, remdesivir, and ribavirin (except favipiravir) prodrugs form a salt bridge with Mg++ 1004. When binding energies were examined, galidesivir (− 6.187 kcal/mol), galidesivir triphosphate (− 8.994 kcal/mol), ribavirin (− 6.128 kcal/mol), ribavirin triphosphate (− 9.280 kcal/mol), remdesivir (− 7.555 kcal/mol), remdesivir triphosphate (− 9.278 kcal/mol), favipiravir (− 5.336 kcal/mol), favipiravir triphosphate (− 8.951 kcal/mol), ATP (− 8.552 kcal/mol) formed the glide gscore energies (Table 1). Also, glide emodel scores were compatible with the glide gscore. The monophosphate and diphosphate interactions of favipiravir, remdesivir, ribavirin, and galidesivir were lower than triphosphate active forms, and the binding energy of triphosphate forms was lower than ATP. Based on these data, triphosphate forms of the prodrugs showed higher interactions. The interactions obtained as a result of molecular-docking studies of other intermediate metabolites are presented in Table 1. Interaction poses of GLT, RBT, FTP, RTP, and ATP at the NSP12-RNA strand active site are shown in Fig. 2.

Table 1

RNA polymerase and RNA residues (PDB: 7BV2) interacting with galidesivir, ribavirin, remdesivir, favipiravir, and their metabolites through H-bonding, salt bridge, Pi-cation, and other interactions

	Compound	Glide gscore	Glide emodel	Protein-RNA-ligand interactions
GLS		− 6.187	− 58.201	ARG836, ARG858, PHE441, ASP845, LYS849, LYS545, ALA547, ILE548, and P: A19, P: U18, P: U17
GLM		− 7.855	− 87.796	Mg1004, ARG55, ARG555, VAL557, LYS545, ASP618, ASP761, ASP760, SER682 and P: U20, T: U10
GLD		− 8.150	− 96.831	Mg1004, ILE548, ALA547, SER549, ASP618, ASP760, ASP76, SER814, ARG55, ARG836 and P: U20
GLT		− 8.994	− 107.152	Mg1004, LYS551, ILE548, ARG555, ALA547, SER549, ARG836, SER814, ASP760, ASP761 and P: U20
RBV		− 6.128	− 62.459	HIE439, SER814, LYS551, ALA550, ARG836, GLU811, ASP760, ASP761 and P: U20
RBM		− 7.983	− 88.177	Mg1004, ARG555, VAL557, SER682, LYS545, ASP760, ASP761
RBD		− 7.949	− 92.564	Mg1004, Mg1005, ARG555, ARG555, LYS545, VAL557, SER682, ASP760, ASP761, GLU811, SER814 and P: U20, T: U10
RBT		− 9.820	− 101.863	Mg1004, SER553, LYS621, SER549, ILE548, LYS551, HIE439, ALA550, ARG555, CYS622, PRO620, TYR619, ASP760 and P: U20
RMD		− 7.755	− 109.151	ARG553, ARG555, ARG555, ALA554, VAL557, CYS622, ASP623, ALA547, ILE548, SER549, SER682, SER759, ASP760, ARG836, and P: A19, P: A19
RMP		− 7.607	− 77.977	Mg1004, ARG555, VAL557, LYS45, SER682, ASP760, ASP761 and P: U20, T: U10
RDP		− 7.312	− 90.081	Mg1004, ARG553, ARG555, ALA554, LYS551, ALA550, SER549, ILE548, LYS545, ALA547, SER814 and P: U20
RTP		− 9.278	− 108.979	Mg1004, LYS621, LYS551, ARG555, ARG553, LYS545, VAL557, SER682, LYS545, ASP760, CYS622, PRO620 and P: U20, T: U10
FPR		− 5.366	− 30.956	Mg1004, CYS622, ASP623, LYS621, TYR619, ASP760
FRR		− 5.693	− 59.769	SER549, ILE548, ALA550, HIE439, ASP760, ASP761, GLU811, SER814 and P: U20
FMP		− 7.609	− 96.463	Mg1004, LYS545, ARG555, VAL557, SER682, ASP760, ASP761 and P: U20, T: U10
FDP		− 8.223	− 97.563	Mg1004, SER682, ARG555, VAL557, LYS545, ALA547, ASP760, ASP761, SER814 and P: U20
FTP		− 8.951	− 102.778	Mg1004, LYS551, LYS621, LYS551, LYS621, ARG555, SER682, LYS545, ARG553, CYS622, PRO620, TYR619 and P: U20
ATP		− 8.552	− 104.854	Mg1004, ARG55, ARG555, LYS545, VAL557, SER682, THR680, THR687, ASN691, CYS622, ASP623, ASP760, ASP761, SER759 and P: U20, T: U10, T: A11

Italic: salt bridge, bold: Hbond, red: Pi-cation İnteraction, normal: hydrophobic and other interactions

GLS galidesivir, GLM galidesivir monophosphate, GLD galidesivir diphosphate, GLT galidesivir triphosphate, RBV ribavirin, RBM ribavirin monophosphate, RBD ribavirin diphosphate, RBT ribavirin triphosphate, RMD remdesivir, RMP remdesivir monophosphate, RDP remdesivir diphosphate, RTP remdesivir triphosphate, FPR favipiravir, FRR favipiravir ribofuranose, FMP favipiravir monophosphate, FDP favipiravir diphosphate, FTP favipiravir triphosphate, U uracil, T thymine, P and T: chain

Fig. 2

RBT (purple), GLT (yellow), RTP (green), FTP (blue), and ATP (gray) interactions in the active site of SARS-CoV-2 RNA polymerase

Simulation of compound-structure complexes and apo-structure was carried out for 50 ns. As the protein compound complexes' structural stability was examined, it was observed that the ATP-protein complex remained very stable under 0.1 nm, and all the protein–ligand complexes were generally stable under 0.2 nm. It was observed that the RMSD value of the ligand-free structure gradually increased and could not remain stable (Fig. 3). As the stability of protein–ligand complexes was examined based on residues, it was found that the fluctuation in the residues between 107 and 109 of the NSP12 chain in non-ligand apoprotein structure was higher than the protein-compound complexes. No significant change was observed in other chains and residues (Fig. 4). As the Radius of the gyration graph, which is the measure of protein compactness, was examined, it was seen that the ATP-bound structure was the most stable, followed by the FTP and RBT-bound structures (Fig. 5). As the protein compound interactions were examined at the end of 50 ns, it was observed that ATP formed one hydrogen bond with Adenine and two hydrogen bonds with Uracil in RNA structure, and the compound also made one hydrogen bond with the Thr650. One of the Mg++ atoms was observed to act as a metal acceptor, while the other Mg++ atom was observed to participate in the ATP interaction with the attractive charge. In GLT residue interactions, it was observed that the adenine in the RNA structure formed a double hydrogen bond with GLT, and Arg555 formed a single hydrogen bond.

Fig. 3

Root-mean-square deviation (RMSD) analysis of the conformational stability of protein–ligand interactions throughout 50 ns compared to the non-liganded structure

Fig. 4

Root-mean-square fluctuation (RMSF) analysis of structures with and without ligand throughout 50 ns

Fig. 5

The radius of gyration graph of structures with and without ligands over 50 ns

On the other hand, the compound created many negative interactions. As the RBT residue interaction was examined, it was observed that uracil in the RNA structure formed a hydrogen bond with the compound and formed a pi-pi-stacked interaction. Many hydrogens and salt bridge interactions were observed. In RTP residue interactions, it was observed that uracil in the RNA structure made double hydrogen bonding. The uracil in the structure of the other RNA chain made hydrogen bond. Many salt bridges with phosphate groups were observed (Figs. 6, 7).

Fig. 6

Protein-compound interactions after 50 ns MD Simulations

Fig. 7.

2D Diagram of Protein-compound interactions after 50 ns MD Simulations

Conclusion

Since the RNA polymerase of SARS-CoV-2 is not present in humans and is important for the virus's survival, this enzyme's inhibition comes to the fore in drug research and development studies. Antiviral prodrugs such as favipiravir, remdesivir, galidesivir, and ribavirin inhibit the RNA-dependent RNA polymerase of the virus, preventing the virus's replication. In this study, molecular docking studies of the prodrugs mentioned above and all their metabolites and molecular dynamics studies of active triphosphate metabolites were carried out. According to molecular docking results, triphosphate active metabolite forms showed higher interaction than prodrug and other intermediate metabolites. We have also observed that these compounds increase the complex's stability and validate the molecular docking results with the RMSD, RMSF, and radius of gyration analyses we obtained with MD simulation. For RNA polymerase inhibition to occur, they must form a salt bridge of Mg++ with the first ribofuranose linked phosphate group and interact with Hie439, Lys545, Ile 548, Ser549, Arg553, Arg555, Val557, Lys621, Cys622, Ser682, Asn691, and Asp760 required. It must also form a hydrogen bond with U20 in the RNA helix. In computer-aided drug design studies, it should not be overlooked that it is indispensable for both the inhibitors' active metabolites and the inhibition of interacting with the RNA strand.