skip to main content

Title: Mathematical Modeling of Remdesivir to Treat COVID-19: Can Dosing Be Optimized?
The antiviral remdesivir has been approved by regulatory bodies such as the European Medicines Agency (EMA) and the US Food and Drug administration (FDA) for the treatment of COVID-19. However, its efficacy is debated and toxicity concerns might limit the therapeutic range of this drug. Computational models that aid in balancing efficacy and toxicity would be of great help. Parametrizing models is difficult because the prodrug remdesivir is metabolized to its active form (RDV-TP) upon cell entry, which complicates dose–activity relationships. Here, we employ a computational model that allows drug efficacy predictions based on the binding affinity of RDV-TP for its target polymerase in SARS-CoV-2. We identify an optimal infusion rate to maximize remdesivir efficacy. We also assess drug efficacy in suppressing both wild-type and resistant strains, and thereby describe a drug regimen that may select for resistance. Our results differ from predictions using prodrug dose–response curves (pseudo-EC50s). We expect that reaching 90% inhibition (EC90) is insufficient to suppress SARS-CoV-2 in the lungs. While standard dosing mildly inhibits viral polymerase and therefore likely reduces morbidity, we also expect selection for resistant mutants for most realistic parameter ranges. To increase efficacy and safeguard against resistance, we recommend more clinical trials with more » dosing regimens that substantially increase the levels of RDV-TP and/or pair remdesivir with companion antivirals. « less
Award ID(s):
Publication Date:
Journal Name:
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. We extend our established agent-based multiscale computational model of infection of lung tissue by SARS-CoV-2 to include pharmacokinetic and pharmacodynamic models of remdesivir. We model remdesivir treatment for COVID-19; however, our methods are general to other viral infections and antiviral therapies. We investigate the effects of drug potency, drug dosing frequency, treatment initiation delay, antiviral half-life, and variability in cellular uptake and metabolism of remdesivir and its active metabolite on treatment outcomes in a simulated patch of infected epithelial tissue. Non-spatial deterministic population models which treat all cells of a given class as identical can clarify how treatment dosage and timing influence treatment efficacy. However, they do not reveal how cell-to-cell variability affects treatment outcomes. Our simulations suggest that for a given treatment regime, including cell-to-cell variation in drug uptake, permeability and metabolism increase the likelihood of uncontrolled infection as the cells with the lowest internal levels of antiviral act as super-spreaders within the tissue. The model predicts substantial variability in infection outcomes between similar tissue patches for different treatment options. In models with cellular metabolic variability, antiviral doses have to be increased significantly (>50% depending on simulation parameters) to achieve the same treatment results as with the homogeneous cellularmore »metabolism.« less
  2. Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. In this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (∼1.5 kcal mol −1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick basemore »pairing with the template nt, and inserts into the active site with a medium low free energy barrier (∼2.6 kcal mol −1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.« less
  3. Abstract

    RNA dependent RNA polymerase (RdRp), is an essential in the RNA replication within the life cycle of the severely acute respiratory coronavirus-2 (SARS-CoV-2), causing the deadly respiratory induced sickness COVID-19. Remdesivir is a prodrug that has seen some success in inhibiting this enzyme, however there is still the pressing need for effective alternatives. In this study, we present the discovery of four non-nucleoside small molecules that bind favorably to SARS-CoV-2 RdRp over the active form of the popular drug remdesivir (RTP) and adenosine triphosphate (ATP) by utilizing high-throughput virtual screening (HTVS) against the vast ZINC compound database coupled with extensive molecular dynamics (MD) simulations. After post-trajectory analysis, we found that the simulations of complexes containing both ATP and RTP remained stable for the duration of their trajectories. Additionally, it was revealed that the phosphate tail of RTP was stabilized by both the positive amino acid pocket and magnesium ions near the entry channel of RdRp which includes residues K551, R553, R555 and K621. It was also found that residues D623, D760, and N691 further stabilized the ribose portion of RTP with U10 on the template RNA strand forming hydrogen pairs with the adenosine motif. Using these models of RdRp,more »we employed them to screen the ZINC database of ~ 17 million molecules. Using docking and drug properties scoring, we narrowed down our selection to fourteen candidates. These were subjected to 200 ns simulations each underwent free energy calculations. We identified four hit compounds from the ZINC database that have similar binding poses to RTP while possessing lower overall binding free energies, with ZINC097971592 having a binding free energy two times lower than RTP.

    « less
  4. The novel coronavirus disease 19 (Covid-19) which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a pandemic across the world, which necessitate the need for the antiviral drug discovery. One of the potential protein targets for coronavirus treatment is RNA-dependent RNA polymerase. It is the key enzyme in the viral replication machinery, and it does not exist in human beings, therefore its targeting has been considered as a strategic approach. Here we describe the identification of potential hits from Indonesian Herbal and ZINC databases. The pharmacophore modeling was employed followed by molecular docking and dynamics simulation for 40 ns. 151 and 14480 hit molecules were retrieved from Indonesian herbal and ZINC databases, respectively. Three hits that were selected based on the structural analysis were stable during 40 ns, while binding energy prediction further implied that ZINC1529045114, ZINC169730811, and 9-Ribosyl-trans-zeatin had tighter binding affinities compared to Remdesivir. The ZINC169730811 had the strongest affinity toward RdRp compared to the other two hits including Remdesivir and its binding was corroborated by electrostatic, van der Waals, and nonpolar contribution for solvation energies. The present study offers three hits showing tighter binding to RdRp based on MM-PBSA binding energy prediction formore »further experimental verification.« less
  5. Abstract

    Structure-based drug design targeting the SARS-CoV-2 virus has been greatly facilitated by available virus-related protein structures. However, there is an urgent need for effective, safe small-molecule drugs to control the spread of the virus and variants. While many efforts are devoted to searching for compounds that selectively target individual proteins, we investigated the potential interactions between eight proteins related to SARS-CoV-2 and more than 600 compounds from a traditional Chinese medicine which has proven effective at treating the viral infection. Our original ensemble docking and cooperative docking approaches, followed by a total of over 16-micorsecond molecular simulations, have identified at least 9 compounds that may generally bind to key SARS-CoV-2 proteins. Further, we found evidence that some of these compounds can simultaneously bind to the same target, potentially leading to cooperative inhibition to SARS-CoV-2 proteins like the Spike protein and the RNA-dependent RNA polymerase. These results not only present a useful computational methodology to systematically assess the anti-viral potential of small molecules, but also point out a new avenue to seek cooperative compounds toward cocktail therapeutics to target more SARS-CoV-2-related proteins.