Blog

May 22, 2023

Activation of the urotensin

Communications Biology volume

Communications Biology volume 6, Article number: 511 (2023) Cite this article

1070 Accesses

20 Altmetric

Metrics details

Remdesivir is an antiviral drug used for COVID-19 treatment worldwide. Cardiovascular side effects have been associated with remdesivir; however, the underlying molecular mechanism remains unknown. Here, we performed a large-scale G-protein-coupled receptor screening in combination with structural modeling and found that remdesivir is a selective, partial agonist for urotensin-II receptor (UTS2R) through the Gαi/o-dependent AKT/ERK axis. Functionally, remdesivir treatment induced prolonged field potential and APD90 in human induced pluripotent stem cell (iPS)-derived cardiomyocytes and impaired contractility in both neonatal and adult cardiomyocytes, all of which mirror the clinical pathology. Importantly, remdesivir-mediated cardiac malfunctions were effectively attenuated by antagonizing UTS2R signaling. Finally, we characterized the effect of 110 single-nucleotide variants in UTS2R gene reported in genome database and found four missense variants that show gain-of-function effects in the receptor sensitivity to remdesivir. Collectively, our study illuminates a previously unknown mechanism underlying remdesivir-related cardiovascular events and that genetic variations of UTS2R gene can be a potential risk factor for cardiovascular events during remdesivir treatment, which collectively paves the way for a therapeutic opportunity to prevent such events in the future.

Nucleoside analogs have a long history in the field of drug design for antiviral treatment because nucleosides are used as building blocks for both DNA and RNA synthesis during viral replication1. The primary mechanism of nucleoside analogs’ antiviral activities is attributed to the inhibition of the viral RNA-dependent RNA polymerase2. In response to the global pandemic of COVID-19, several nucleoside analogs, such as remdesivir, molnupiravir, and favipiravir, have been developed to treat the disease.

Remdesivir (GS-5734; Veklury) is a modified adenosine analog that contains a McGuigan prodrug moiety, including phenol and l-alanine ethylbutyl ester, which increases its lipophilicity and cell permeability3. Following intravenous administration, remdesivir is rapidly converted into the mono-nucleoside form (GS-441524) and is intracellularly metabolized by multiple host enzymes to its pharmacologically active triphosphate form, which, in turn, acts as a potent and selective inhibitor of the RNA-dependent RNA polymerase of multiple viruses4,5. Remdesivir was initially used for the treatment of Ebola virus6 and has been approved for the treatment of Coronavirus Disease 2019 (COVID-19) amid the global pandemic. Remdesivir shortened the time to recovery in adults who were hospitalized with COVID-19 and had evidence of lower respiratory tract infection7. More recent data demonstrated a significant reduction in hospitalizations with a 3-day course of intravenous remdesivir8. Although remdesivir is generally well tolerated by most individuals, common adverse events for remdesivir have been reported, including rash, headache, nausea, diarrhea, and elevated transaminases7. Current guidelines for remdesivir suggest careful monitoring of liver function during treatment and recommend against its use in patients with renal dysfunction9. Additionally, cardiovascular events, including hypotension, bradycardia, QT prolongation, and T-wave abnormality, have been reported10,11,12,13. When administrated intravenously, remdesivir shows broad tissue distribution, including to the heart14, but the precise molecular mechanism underlying the cardiovascular side-effects of remdesivir remains unclear.

Molnupiravir (EIDD-2801/MK-4482; Lagevrio) has been authorized for emergency use by the FDA under an Emergency Use Authorization for the treatment of mild-to-moderate COVID-19 in adults who are at high-risk for progression to severe COVID-19. Molnupiravir is an oral cytosine analog that contains an isopropylester prodrug of the β-d-N4-hydroxycytidine (NHC). The active form of NHC is a substrate of viral RNA-dependent RNA polymerase and impairs the fidelity of SARS-CoV-2 replication, provoking error catastrophe15. A clinical trial in nonhospitalized adults showed that early treatment with molnupiravir effectively reduced the risk of hospitalization or death in at-risk, unvaccinated adults with COVID-1916. In addition to molnupiravir, favipiravir (T-705; Avigan), which is an anti-influenza drug, has been under clinical trial for the treatment of COVID-19. Favipiravir s a nucleobase analog derived from pyrazine carboxamide (6-fluoro-3-hydroxy-2-pyrazinecarboxamide)17. The suggested modes of action of favipiravir comprise a mix of both chain termination and mutator events18. Several adverse events have been reported during the use of molnupiravir and favipiravir, including diarrhea, dizziness, and nausea for molnupiravir19, and hyperuricemia and increased alanine aminotransferase for favipiravir20. Importantly, unlike remdesivir, cardiovascular side effects have not been reported with the use of molnupiravir or favipiravir.

In addition to their function as building blocks for DNA/RNA synthesis, nucleotides/nucleosides can act as endogenous ligands for G-protein-coupled receptors (GPCR) and induce diverse pathophysiological responses21,22,23,24,25. Given the presence of nucleoside-mimicking structures in remdesivir, molnupiravir, and favipiravir, we hypothesized that these drugs might cause side effects by directly activating GPCR. Here, we report that remdesivir, but not molnupiravir and favipiravir, is a selective ligand for the urotensin-II receptor (UTS2R) and causes cardiac dysfunction.

We screened anti-COVID-19 drugs, including remdesivir, molnupiravir, and favipiravir, against 348 GPCRs using an alkaline phosphatase-tagged transforming growth factor-α (AP-TGFα) shedding assay26. We used chimeric Gα subunit proteins for the initial screening to efficiently detect the receptor activation regardless of the type of Gα subunit involved26. Among the three drugs, we discovered that remdesivir is a selective activator for the urotensin-II receptor (UTS2R) (Fig. 1a, Supplementary Fig. 1a, Supplementary Data 1). As remdesivir potently induced UTS2R response without any chimeric Gα protein, the subsequent UTS2R analysis was performed without the exogenous addition of the chimeric Gα protein. A concentration–response analysis revealed that the half-maximal effective concentration (pEC50) of remdesivir was 4.89 ± 0.03 (EC50 = 13 ± 0.9 μM, Fig. 1b left). It should be noted that both the potency and the efficacy of remdesivir (Emax = 47 ± 1.4 %AP-TGFα release) toward UTS2R are lower than those of the endogenous peptidic ligand urotensin-II (UT2, pEC50 = 10.72 ± 0.04; EC50 = 21 ± 2.1 fM, Emax = 59 ± 0.70 %AP-TGFα release, Fig. 1b right). Nevertheless, remdesivir administration at a clinical dosage is considered to be within the working range of agonistic effects because the maximum plasma concentration of remdesivir reaches 9.03 μM after intravenous injection in healthy adults27. Interestingly, unlike UT2, remdesivir failed to induce a β-arrestin recruitment response (Fig. 1c, Supplementary Fig. 1b).

a Chemical structures of tested compounds and heatmaps showing the relative levels of GPCR activation as measured by the TGFα-shedding assay. The color scale represents the % GPCR activation compared to TPA (12-O-tetradecanoylphorbol 13-acetate)-mediated receptor activation, which induces the maximum TGFα-shedding response independently of GPCR. The yellow cell represents activation of UTS2R by remdesivir. The tested compounds’ concentrations were 10 μM for remdesivir, molnupiravir, GS-441524, and sofosbuvir, and 100 μM for favipiravir (n = 3). b TGFα-shedding response curves for UTS2R by remdesivir (left) and urotensin-II (UT2, right) in the presence or absence of urantide, a UTS2R antagonist (Data are shown as means ± SEM, n = 3). The half-maximal effective concentration (pEC50) of remdesivir was 4.89 ± 0.03 (EC50 = 13 μM, Emax = 47 ± 1.4 %AP-TGFα release). The half-maximal effective concentration (pEC50) of UT2 was 10.72 ± 0.04 (EC50 = 21 fM, Emax = 59 ± 0.70 %AP-TGFα release). c Remdesivir (left)- or urotensin-II (UT2, right)-mediated β-arrestin 1 recruitment assay for UTS2R. Data are shown as means ± SEM (n = 3). d, e Competitive binding assay using biotin-UT2 peptide. The membrane fraction of HEK293 cells expressing FLAG-UTS2R (input) was incubated with biotin-UT2 in the presence or absence of remdesivir or urantide. FLAG-UTS2R was pulled down by streptavidin-coated magnetic beads (M280), and FLAG-UTS2R was detected by western blot. d Representative image of three independent trials was shown. e Band densitometry. ***p < 0.001 and **p < 0.01 versus vehicle + 10 μM UT2 group (one-way ANOVA followed by Dunnett's multiple comparisons test); means ± SEM.

To evaluate whether remdesivir can directly interact with UTS2R, we performed a competitive binding assay. A previous study has shown that a synthetic UT2 peptide with an N-terminal biotin tag (biotin–UT2) can form stable binding with UTS2R28. Importantly, the binding was stable enough for biochemical pulldown of UTS2R using streptavidin resin. Taking advantage of the unique property of biotin–UT2 peptide, we performed a pulldown assay using membrane fraction of HEK293 cells expressing UTS2R in the presence or absence of remdesivir (Supplementary Fig. 1c). Among various types of magnetic beads that we tested, magnetic beads with hydrophobic coating showed the maximum pulldown efficacy with low non-specific binding (Supplementary Fig. 1d). Using the optimized protocol, we found that remdesivir and uranide significantly impaired biotin–UT2-mediated UTS2R pulldown (Fig. 1d, e, Supplementary Fig. 1e). These results suggest that remdesivir can directly interact with UTS2R, thereby interfering the biotin–UT2–UTS2R binding.

We further investigated whether the metabolites of remdesivir activate UTS2R. GS-441524 and GS-704277, the major and minor metabolites of remdesivir6,27, respectively, showed no effect on UTS2R activation (Fig. 1a, Supplementary Fig. 1a). To investigate whether the McGuigan prodrug moiety is responsible for remdesivir's UTS2R activation, we tested sofosbuvir (GS-7977), an FDA-approved McGuigan-class prodrug3, against UTS2R. However, sofosbuvir did not activate UTS2R (Fig. 1a, Supplementary Fig. 1a). These results suggest that both the McGuigan prodrug moiety and nucleoside base of remdesivir are required to activate the receptor.

UTS2R belongs to the class A GPCR family29 and consists of the canonical 7 transmembrane helices (TM), the amphipathic helix 8 at the C-terminus (H8), two antiparallel β-strands in the extracellular loop 2, a relatively short N-terminal domain with two N-glycosylation sites, and a palmitoylation anchor at the C-terminal tail (Fig. 2a upper). Meanwhile, remdesivir is an analog of adenosine and a phosphoramidate prodrug of the McGuigan class (phenol and l-alanine ethylbutyl ester), which masks the anionic phosphate moiety on remdesivir, thus improving drug delivery. To further elucidate the molecular basis for ligand-receptor binding, we performed in silico structural docking of UTS2R in the presence of remdesivir (Fig. 2a lower). Our analysis showed multiple amino acid residues in the orthosteric pocket that potentially stabilize binding to remdesivir. First, the cyano group of the nucleosugar of remdesivir forms a hydrogen bond to the UTS2R residue T3047.42×41 (superscripts denote the generic GPCR numbering system30). Second, the phenyl group of the McGuigan prodrug moiety of remdesivir is capped with N2977.35×34. Finally, the amino group of the nucleobase of remdesivir forms a hydrogen bond with M1343.36 at the bottom of the pocket (NH···S hydrogen bond31).

a Upper panel, Schematic of UTS2R; Lower panel, Docking model of UTS2R with remdesivir. AlphaFold structure of human UTS2R is shown as green ribbons, omitting TM5 for clarity. Remdesivir and selected side chains of the receptor are shown as sticks and colored gray and green, respectively. Black dashed lines indicate hydrogen bonds. b, c Effects of the indicated mutations on remdesivir (b)- or UT2 (c)- mediated receptor activation. The y-axis (⊿pEC50 value) represents the relative activation potency of each mutant receptor compared to the WT receptor. ⊿pEC50 value = pEC50 mutant − pEC50 WT. The ⊿pEC50 cutoff value was set to −1, as indicated by dashed lines. EC50 values were determined by the TGFα-shedding assay. **p < 0.01, ***p < 0.001 vs. WT by one-way ANOVA followed by Dunnett's multiple comparison test. Data are shown as means ± SEM (n ≥ 3).

To validate the docking model, we mutated these UTS2R residues (T3047.42×41, N2977.35×34, and M1343.36) to other amino acids and tested whether these mutations affect remdesivir-mediated UTS2R activation (Fig. 2b, c, Supplementary Fig. 2a, b). In line with the in silico simulation, mutation of the three residues almost completely abolished the activation potency of remdesivir to UTS2R, indicating that remdesivir needs to contact UTS2R at multiple residues to achieve stable binding. Both nucleoside and McGuigan moiety are essential for receptor binding, with N2977.35×34 interacting with the McGuigan moiety and M1343.36 and T3047.42×41 interacting with the nucleoside moiety. This structural observation explains why mutating a single amino acid strikingly reduces the remdesivir's activity towards UTS2R. This explanation is further supported by the result that UTS2R does not respond to any remdesivir metabolites, from which the McGuigan moiety is metabolically removed. In contrast, these mutations only partially reduced the UT2-mediated UTS2R activation. Meanwhile, contrary to these three residues, we found that D1303.32 has an inhibitory effect on remdesivir–UTS2R binding. The D1303.32N mutation abolished UT2-mediated UTS2R activation but significantly upregulated remdesivir-mediated receptor activation (Fig. 2b, c, Supplementary Fig. 2a, b). According to the docking model, the D1303.32 residue is localized near the nucleobase moiety of remdesivir (Fig. 2a). We speculate that the negative charge of D1303.32 residue may induce an electron repulsion effect32, leading to destabilization of the interaction between D1303.32 residue and the nucleobase. These results demonstrate that remdesivir-mediated UTS2R activation is mediated by specific binding between UTS2R and its McGuigan prodrug moiety and nucleoside base, which differs from that of UTS2R and the endogenous ligand.

To determine whether remdesivir-mediated UTS2R activation induces intracellular signaling transduction, we stimulated UTS2R-expressing HEK293 cells with remdesivir and examined the phosphorylation status of extracellular signal-regulated kinase (ERK)1/2. Application of remdesivir for up to 72 h evoked long-lasting and dose-dependent phosphorylation of ERK1/2 (Fig. 3a, Supplementary Fig. 3a). Importantly, remdesivir-mediated ERK phosphorylation was abolished by the UTS2R antagonist (Fig. 3a, Supplementary Fig. 3a). The response of ERK1/2 induced by remdesivir was similar to that induced by UT2 (Supplementary Fig. 3b).

a Serum-starved HEK293 cells overexpressing UTS2R were stimulated with the indicated concentrations of remdesivir for 5 min with or without urantide, a UTS2R antagonist, and the lysates subjected to western blotting analysis. ERK1 and ERK2 activation ratios (pERK1/ERK1 and pERK2/ERK2) were calculated with data normalized to the vehicle. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001 by Tukey's multiple comparisons test. Data are represented as means ± SEM (n = 3). b Left, Temporal correlation between the field potential duration and the QT interval on the surface ECG. Right, Schematic of the multielectrode array (MEA) platform. c Representative field potential waveform in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) treated with 1 µM remdesivir in the presence or absence of urantide for 72 h. d Effect of remdesivir and urantide on field potential prolongation in hiPSC-CMs. *p < 0.05, **p < 0.01 by two-way ANOVA followed by Šídák's multiple comparisons tests. Data are represented as means ± SEM (n = 3). e Perforated patch-clamp to record spontaneous action potentials of hiPS-CMs in current-clamp models. The hiPSC-CMs were treated with 10 µM remdesivir in the presence or absence of 50 µM urantide for 72 h. *p < 0.05 and **p < 0.01 by Tukey's multiple comparisons test. Data are represented as means ± SEM (n ≥ 3).

UT2 and its receptor UTS2R are widely expressed in tissues, with relatively high expression in cardiovascular systems33 (Supplementary Fig. 3c, d). Prompted by the cardiotoxic potential of remdesivir10,11,12,13, we assessed the impact of remdesivir on cardiomyocyte functions. We examined the effect of remdesivir on the field potential (FP) using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs)34, in which the expression level of UTS2R is comparable to that of the human heart (Supplementary Fig. 3e). The FP duration (FPD) correlates closely with the QT interval on an electrocardiogram (ECG)35 (Fig. 3b left). Notably, the prolongation of the QT interval has been linked with the occurrence of severe and potentially fatal arrhythmias, and QT prolongation is the leading cause of drug-induced cardiovascular toxicity36. For the assessment of FPD, the use of a multielectrode array (MEA) platform is well-accepted for its ability to monitor the electrophysiology of cardiomyocytes at the cell population level (Fig. 3b right)35. An MEA analysis revealed that hiPSC-CMs treated with remdesivir showed a prolonged FPD of 1.32 ± 1.38% at 24 h, 5.60 ± 1.60% at 48 h, and 15.57 ± 3.49% at 72 h, respectively. Notably, the prolongation was significantly suppressed by the UTS2R antagonist (Fig. 3c, d). Indeed, the application of UTS2R antagonist with remdesivir showed almost no FPD delay at 24 h (−1.71 ± 1.60%) and 48 h (−0.61 ± 0.11%), while the reversal effect was still significant but became partial at 72 h (Fig. 3d).

In addition to FPD, we evaluated QT prolongation by perforated patch–clamp and recorded spontaneous action potentials of hiPS-CMs. We measured the APD90 (action potential duration at 90% repolarization). Remdesivir significantly prolonged the APD90 in hiPS-CMs, and this prolongation was markedly attenuated by the administration of urantide (Fig. 3e). Thus, these results elucidate a previously unknown mechanism of the reported proarrhythmic risks of remdesivir11,12,13, which is at least partially dependent on UTS2R. Next, we assessed the effects of remdesivir on cardiac contractility. To this end, we used neonatal rat cardiomyocytes (NRCMs) and evaluated the contraction force. Under a constant pacing protocol, NRCMs with chronic remdesivir treatment showed reduced contractility, which was significantly attenuated by the UTS2R antagonist (Fig. 4a).

a, b Left, Representative waveform of contractility under pacing in neonatal rat cardiomyocytes (NRCMs) by 1 µM remdesivir with or without a urantide or b pertussis toxin (PTX), a Gαi/o inhibitor, and YM-254890, a Gαq/11 inhibitor, at 48 h. Right, the effect of remdesivir with a urantide or b PTX and YM-254890 on contractility under pacing in NRCMs. **p < 0.01, ***p < 0.001, ****p < 0.0001 by Tukey's multiple comparisons test. Data are represented as means ± SEM (n = 3). c, d Representative western blot for phosphorylation of c ERK1/2 and d protein kinase B (AKT). Serum-starved HEK293 cells overexpressing UTS2R were stimulated with the indicated concentrations of remdesivir for 48 h. For Gi/o protein inhibition, cells were incubated with PTX for at least 18 h at 150 ng/mL. ERK1, ERK2, and AKT activation ratios were calculated with data normalized to the vehicle. *p < 0.05, **p < 0.01, ***p < 0.001 by Tukey's multiple comparisons test. Data are shown as means ± SEM (n = 3). e, f The effects of remdesivir on the contractility of adult mouse cardiomyocytes. Isolated adult mouse cardiomyocytes were treated with remdesivir (10 µM) for 30 min, and the percentile change in cell area (e) and the ratio of cardiomyocyte shortening (f) was measured during electrical pacing. ****p < 0.001 by unpaired t-test. Data are shown as means ± SEM (a total of 54 cells from 3 mice).

Heterotrimeric G proteins, including Gαs, Gαi/o, Gαq/11, and Gα12/13, are the downstream effectors of GPCRs. Among those, the Gαi/o family has been implicated in myocardial contractility and heart rate via the modulation of ion channels37. Since UTS2R is coupled to Gαi/o and Gαq29, we sought to determine which Gα protein is involved in the remdesivir-mediated decrease in myocardial contraction. Gαi/o inhibitor pertussis toxin (PTX), but not Gαq/11 inhibitor YM-254890, completely blocked the effect of remdesivir and restored the peak contractions of NRCMs (Fig. 4b). A previous study demonstrated how activating Gαi/o in NRCMs leads to liberation of Gβγ, which intracellularly activates PI3K, leading to the activation of AKT and ERK1/238. In line with this observation, Gαi/o inhibition reduced the remdesivir-induced phosphorylation of ERK1/2 and AKT (Fig. 4c, d). These results suggest that remdesivir reduces cardiomyocyte contractility through the Gαi/o-dependent AKT/ERK signal transduction pathway. Collectively, our results indicate that remdesivir itself can function as an exogenous ligand of UTS2R. Moreover, we identified the proarrhythmic and negative inotropic potential of remdesivir, both of which are UTS2R dependent.

Because neonatal cardiomyocytes are not terminally differentiated, we further investigated the effect of remdesivir on mature cardiomyocytes isolated from adult mouse hearts (Supplementary Fig. 3f). Contraction of cardiomyocytes was recorded under pacing conditions, and the degree of contraction was evaluated by morphological analysis. Similar to neonatal cardiomyocytes, remdesivir significantly impaired the contraction of adult cardiomyocytes (Fig. 4e, f).

A previous study suggested that remdesivir-related cardiotoxicity can be caused by mitochondrial dysfunction39 since the active form of remdesivir shows an inhibitory effect toward mitochondrial RNA polymerase (mtRNAP) at a high dose40. However, treatment with remdesivir at 10 μM, which is equivalent to the maximum plasma concentration following remdesivir administration in humans27, did not affect the steady-state levels of mitochondrial respiratory complex proteins (Supplementary Fig. 3g, h).

To understand the impact of genetic variance on the susceptibility to remdesivir–UTS2R signaling in humans, we extract single-nucleotide variant (SNV) information from a genome database that includes the 14KJPN Genome Reference Panel, which is a large-scale population-based genomic database constructed from the DNA sequence of 14,000 Japanese individuals41, followed by functional assay on receptor activation. A total of 2178 variants are reported in the UTS2R locus, of which 139 variants are missense variants (Supplementary Data 2). A considerable number of missense SNVs listed in the 14KJPN are also reported in the gnomAD database (https://gnomad.broadinstitute.org), which contains SNVs from broader ethnicities.

Accordingly, we generated 110 missense mutants corresponding to the human SNVs in the UTS2R gene. We excluded 29 missense mutations in the 5′-region of UTS2R as the confidence score for the structural prediction of the N-terminus was relatively low. Among the 110 missense SNVs, 44 SNVs showed a decrease in the sensitivity to remdesivir compared to the WT receptor (⊿pEC50 <−0.3, which corresponds to over a twofold EC50 increase compared to WT receptor; Fig. 5a upper panel). Meanwhile, 47 SNVs displayed a decreased sensitivity to UT2 compared to WT receptor, of which 18 SNVs overlapped with those that showed a decreased sensitivity to remdesivir (Fig. 5a lower panel, Supplementary Fig. 4a). Notably, we found four missense SNVs (G681.49C, D1303.32G, V15934.54M, and A249ICL3G) that can increase the receptor sensitivity toward remdesivir compared to the WT receptor (⊿pEC50 > 0.3, which corresponds to a < 0.5-fold EC50 decrease compared to the WT receptor; Fig. 5a–c). Furthermore, among these four gain-of-function remdesivir-sensitive UTS2R SNVs, the G681.48C and D1303.32G mutants conversely exhibited a decrease in the sensitivity toward UT2, while the V15934.54M and A249ICL3 mutants showed a moderate or insignificant increase in UT2 sensitivity (⊿pEC50 < 0.3; Fig. 5c, Supplementary Fig. 4a). Collectively, the results suggest that individuals with a G681.49C, D1303.32G, V15934.54M, or A249ICL3G mutation in the UTS2R gene are sensitive to remdesivir, possibly making them more susceptible to UTS2R-mediated cardiotoxicity, although the allele frequencies of the gain-of-function variants are low (Supplementary Fig. 4b).

a Effects of 110 missense SNVs from the UTS2R gene on (upper panel) remdesivir-mediated and (lower panel) UT2-mediated receptor activation, measured by the TGFα-shedding assay. The y-axis (⊿pEC50 value) represents the relative activation potency of each mutant receptor compared to the WT receptor. The ⊿pEC50 cutoff value was set to −1, as indicated by dashed lines. Light blue bands represent the range of −0.3 < ⊿pEC50 < 0.3, which corresponds to the range of an 0.5-fold to 2-fold (less than 2-fold) change in the EC50 of the mutant receptor compared to the WT receptor. All experiments were performed in triplicate, and the data are expressed as the means. b Snake plot diagram of UTS2R showing locations of mutations and their effects on the potency of remdesivir (modified from www.gpcrdb.org). c Effects of the selected mutations on (left) remdesivir-mediated and (right) UT2-mediated receptor activation. EC50 values are determined by the TGFα-shedding assay. "ns" p > 0.05, and ****p < 0.0001 vs. WT by one-way ANOVA followed by Dunnett's multiple comparisons tests. Data are shown as means ± SEM (n ≥ 3).

Following the completion of the Adaptive COVID-19 Treatment Trial 1 (ACTT-1)7, which demonstrated remdesivir's superiority over placebo for improving time to recovery in hospitalized COVID-19 patients, remdesivir has been one of the most commonly prescribed medications for patients hospitalized for COVID-19 infection. Importantly, the ACTT-1 investigators reported that 0.2% of patients receiving remdesivir, but not patients receiving placebo, showed arrhythmias (other than atrial fibrillation, supraventricular tachycardia, ventricular tachycardia, and ventricular fibrillation), although these cardiovascular effects were not considered as adverse effects. However, the percentage of patients experiencing arrhythmias could be underestimated since early clinical trials are insufficiently powered to detect uncommon adverse events42. Indeed, a large retrospective pharmacovigilance cohort study, which used medical records from more than 130 countries and 20 million patients, reported that cardiac arrest, bradycardia, and hypotension are associated with remdesivir use10. The elevation of plasma remdesivir concentration is associated with an increase in FP duration with decreased Na+ peak amplitudes and spontaneous beating rates, which might potentially induce prolonged QT interval and torsade de point43,44. Importantly, the cardiac adverse effects appear to be directedly caused by remdesivir, since they were reported to resolve within 24–48 h of discontinuing remdesivir45,46. To date, the precise mechanism underlying the remdesivir's cardiac side effects has remained unclear, with no means of predicting the populations susceptible to its cardiac side effects and no specific treatment for the cardiotoxicity. Therefore, there is a pressing need to understand how remdesivir induces cardiac dysfunction.

Using an unbiased and large-scale GPCR screening strategy, we identified that remdesivir, but not molnupiravir and favipiravir, can selectively activate UTS2R. The interaction between remdesivir and UTS2R was further supported by the competitive binding assay using the biotin–UT2 peptide. Interestingly, UTS2R is highly expressed in heart tissue, including cardiomyocytes. Activation of UTS2R by Urotensin-II (UT2), the endogenous ligand of UTS2R, has been implicated in cardiac dysfunction. For example, the plasma level of UT2 and expression level of UTS2R in cardiomyocytes are elevated in patients with end-stage congestive heart failure47. In line with these previous studies, we found that activation of UTS2R by remdesivir at a concentration of 1 μM induced electrical abnormalities and contraction force impairment in cultured cardiomyocytes, both of which resemble the reported cardiac side effects in humans. Furthermore, these adverse effects were effectively blocked by antagonizing UTS2R or inhibiting its downstream signaling. Clinically, remdesivir is administrated intravenously at a dose of 200 mg once, followed by 100 mg daily for a total of 5–10 days in adults and children ≥40 kg. The estimated peak plasma concentration of remdesivir is 9.03 μM in healthy adults and can be higher in patients with renal or hepatic impairment because of its renal and biliary excretion48. Our results suggest that the clinical dosage of remdesivir is sufficient to activate UTS2R, and patients with renal or hepatic impairment may be at high risk for adverse events mediated through the remdesivir–UTS2R axis. Notably, our GPCR assay clearly showed that remdesivir has no effect on adenosine receptors, of which the activation can also impact cardiac functions. Thus, the activation of UTS2R is likely responsible for the cardiac side effects.

An important finding of this study is that SNVs in the coding region of UTS2R have large impacts on its response to remdesivir. This result is in line with previous studies that show SNVs in GPCR genes are associated with the pathophysiology of various diseases, including CV diseases, and affect therapeutic outcomes49. While 40% of SNVs in UTS2R (44/110) showed at least a twofold reduction in potency toward remdesivir compared to the WT receptor, and 56% (62/110) showed no remarkable change (within the range of an 0.5-fold to 2-fold change), we identified four gain-of-function SNVs that showed more than two-fold increase in response to remdesivir. Notably, D1303.32G was a variant at the same amino acid for which we identified a gain-of-function mutation (D1303.32N) using in silico modeling, indicating the importance of D1303.32 residue in remdesivir recognition.

Activation of GPCR often induces recruitment of β-arrestin to the receptor, followed by receptor internalization to the intracellular compartment. This internalization terminates GPCR activation and promotes secondary signaling pathways50. Interestingly, unlike the endogenous ligand UT2, which effectively induces β-arrestin recruitment (Fig. 1c)51, remdesivir-mediated UTS2R activation did not induce β-arrestin recruitment. Thus, our results suggest that remdesivir is a G-protein-biased ligand (Fig. 1c, Supplementary Fig. 1b). Such biased activation of UTS2R by remdesivir may have an important impact on cardiac function in a way that allows prolonged activation of UTS2R without β-arrestin-mediated shutdown and thereby an exaggeration of downstream signaling and enhanced cardiotoxic effects. However, it is also conceivable that the modest potency of remdesivir could have hindered the detection of β-arrestin response. Further study using a new methodology for sensitive detection of the β-arrestin response is needed to evaluate the bias effect of remdesivir.

Upon ligand binding, GPCRs, including UTS2R, initiate conformational change that induces the activation of heterotrimeric G proteins and dissociation of Gα and Gβγ subunit complexes. Gα proteins include Gαs, Gαi/o, Gαq/11, and Gα12/13 proteins, which are responsible for downstream signaling transduction. Members of the Gαi/o family are widely distributed, including in the cardiac system, where they are highly expressed and act to regulate myocardial contractility and heart rate via modulation of ion channels37. For example, resting heart rate is controlled by cholinergic signals mediated through muscarinic M2 Gαi-coupled receptors. These effects occur by inhibition of adenylyl cyclase (AC) and by Gβγ-inhibition of a potassium channel in the sinoatrial node. The transduction mechanism of UTS2R is the coupling and activation of Gαi/o as well as Gαq29, consistent with the UT2–UTS2R axis being implicated in CV regulation through complex signaling pathways, both physiologically and pathologically. Our results using NRCMs clearly showed that the remdesivir-mediated decrease in myocardial contraction is dependent on Gαi/o, but not Gαq/11 (Fig. 4b). Since NRCMs resemble the phenotype of atrial myocytes52, the decrease in myocardial contractility under constant pacing by remdesivir indicates impairment of calcium handling53, which is consistent with the slowing of the heart rate, a distinctive cardiovascular side effect of remdesivir10. Additionally, we showed that remdesivir significantly impaired the contractility of adult cardiomyocytes (Fig. 4e, f).

Despite the discovery that remdesivir can activate UTS2R and cause cardiotoxicity, the lack of clinical evidence is a major limitation of this study. Yet, as the allele frequencies of the gain-of-function variants are low, (Supplementary Fig. 4b) and usage of remdesivir is expected to decline due to the recent prevalence of the Omicron COVID-19 variant with its milder symptoms, a large-scale clinical study on the association between remdesivir sensitivity and genomic variance is challenging to execute. Additionally, we have not determined the precise molecular mechanisms that induce the UTS2R-dependent proarrhythmic risk of remdesivir. Possible explanations are the impaired regulation of gene expression or trafficking of ERG potassium channels, which are essential for electrical activity in the heart54. Alternatively, the chronic and cumulative cardiotoxic effects of the remdesivir–UTS2R axis are consistent with a downstream impact on translation and transcription. Although we do not exclude the possibility that remdesivir might affect mitochondrial metabolism through UTS2R signal transduction (Supplementary Fig. 3f), the current data suggest that the remdesivir–UTS2R axis is the major off-target of remdesivir.

In conclusion, to our knowledge, this is the first report showing that remdesivir is a selective agonist of UTS2R and that remdesivir-mediated UTS2R activation underlies drug-mediated cardiotoxicity. Furthermore, we discovered that specific SNVs in UTS2R can increase the sensitivity to remdesivir. Thus, this study provides mechanistic insights into remdesivir-mediated cardiac side effects and the therapeutic opportunity to prevent aversive events in the future.

The overall goal of this study was to explore the molecular mechanisms of remdesivir-related cardiotoxicity in anti-COVID-19 therapy. We first performed the GPCR screening using major anti-COVID-19 drugs as ligands (n = 3 per GPCR). We validated the results of GPCR screening by concentration-response analysis and determined the EC50 and Emax values of remdesivir against UTS2R (n = 16). We also performed docking simulation using AlphaFold and determined the key amino acid residues that were essential for remdesivir–UTS2R binding, as demonstrated by a mutagenesis study (n ≥ 3). We next examined the cardiotoxic potential of remdesivir. Using hiPS-derived cardiomyocytes, we found that the proarrhythmic risk of remdesivir is largely mediated by UTS2R (n = 3). Moreover, using NRCMs, we showed that the decrease in contraction force by remdesivir is mediated by the UTS2R-Gαi/o axis (n = 3). Finally, the effects of SNVs in the UTS2R gene on remdesivir-mediated receptor activation were evaluated using a GPCR assay (n = 3 per SNV).

Remdesivir (GS-5734), GS-441524, and sofosbuvir (GS-7977) were purchased from Selleck Chemicals. Favipiravir was purchased from Tokyo Chemical Industry. Fibronectin was purchased from Sigma-Aldrich. All other reagents were of analytical grade and obtained from commercial sources.

To search for SNVs in UTS2R, we used a public database of 14,000 healthy Japanese individuals, called 14KJPN, released by Tohoku University's Tohoku Medical Megabank Organization (ToMMo, https://jmorp.megabank.tohoku.ac.jp/)41, which covers variant frequency data.

To measure the activation of the GPCR, a transforming growth factor-α (TGFα) shedding assay was performed26. Briefly, HEK293A cells were seeded in a 12-well culture plate. A total of 348 pCAG plasmids encoding an untagged GPCR, including the UTS2R construct, were prepared (Supplementary Data 1). The UTS2R mutants, including 110 missense mutants corresponding to the human SNVs in the UTS2R gene, were generated by introducing single-point mutations using a KOD-Plus-Mutagenesis kit. Cells were transfected using a polyethylenimine (PEI) reagent (2.5 µl of 1 mg/ml per well hereafter; Polysciences) with these UTS2R plasmids (100 ng per well hereafter), together with the plasmids encoding alkaline phosphatase (AP)-tagged TGFα (AP-TGFα; 250 ng). UTS2R activation assays were performed with endogenous G protein unless otherwise specified. Chimeric Gα subunit proteins (mixture of Gαq/s, Gαq/i1, Gαq/i3, Gαq/o, Gαq/z, Gαq/12, Gαq/13, and Gα16; each 10 ng) were used for initial screening of 348 GPCRs (Fig. 1a). After a 24-h culture, the transfected cells were harvested and collected by centrifugation. Cells were suspended in Hank's Balanced Salt Solution (HBSS) containing 5 mM HEPES (pH 7.4) and seeded in a 96-well plate. After a 30-min incubation, the test compound was added to the cells. After 1-h incubation, the conditioned media were transferred to an empty 96-well plate. The AP reaction solution (a mixture of 10 mM p-nitrophenylphosphate (p-NPP), 120 mM Tris-HCl (pH 9.5), 40 mM NaCl, and 10 mM MgCl2) was added to plates containing cells and conditioned media. The absorbance at a wavelength of 405 nm was measured using a microplate reader (Molecular Devices) before and after 1-h incubation of the plates at room temperature. Compound-induced AP-TGFα release was calculated by subtracting spontaneous AP-TGFα release signal from compound-induced AP-TGFα release signal, and the percentages were fitted to a four-parameter sigmoidal concentration–response curve using Prism 9 software (GraphPad Prism), and the EC50 and Emax values were obtained.

NanoBiT enzyme complementation-based β-arrestin recruitment assay was performed55. Briefly, HEK293A cells were seeded in a 6-well culture plate and transfected using a PEI reagent (5 µl of 1 mg/ml per well hereafter) with a mixture of plasmids consisting of 500 ng ssHA-FLAG-UTS2R-SmBiT construct (N-terminal hemagglutinin signal sequence followed by a FLAG epitope tag, plus C-terminal SmBiT with a 15-amino acid flexible linker), 100 ng LgBiT-β-arrestin construct (N-terminal LgBiT with a 15-amino acid flexible linker) and 400 ng empty plasmid. After 24-h culture, the transfected cells were harvested, suspended in 2 ml of assay buffer (Hanks’ Balanced Salt Solution containing 5 mM HEPES (pH 7.4) and 0.01% bovine serum albumin), and seeded in a 96-well plate at a volume of 80 µl per well. Coelenterazine (20 µl of 50 µM diluted in the assay buffer) was added to the cell plates, followed by incubation at room temperature for 2 h in the dark. After measurement of baseline luminescence using a Spectra Max L Microplate Reader (Molecular Devices), remdesivir or UT2 (20 µl of 6× concentration) was added. Luminescence was measured every 20 s after compound addition. The average luminescent signal over 5–10 min was normalized to an initial value. The fold-change values were further normalized to that of the vehicle-treated condition and were fitted to a four-parameter sigmoidal concentration-response curve to obtain EC50 values using Prism 9 software.

To affinity pulldown of UTS2R, a UT2 peptide with an N-terminal biotin tag and a GSSG spacer was synthesized (Peptide Institute, INC)28. HEK293 cells were transfected with the plasmid encoding FLAG-UTS2R and were serum-starved for 6 h before experimentation. Cells were washed twice with ice-cold PBS. Hypotonic solution (50 mM HEPES, pH 7.5 with protease inhibitor cocktails) was added to the cell plate, and the cells were harvested by scraping. Cells were homogenized using a dounce homogenizer (IKA EUROSTAR 20 digital) at 2,000 rpm for 20 strokes at 4 °C, and further homogenized using a tissue rupture (Qiagen) for 5 s at 4 °C. The homogenate was centrifuged at 800 x g for 10 min at 4 °C to remove unbroken cells and the nucleus. The supernatant was further centrifuged at 30,000×g for 30 min at 4° (CP80NX, himac). The resulting pellet was used as crude membrane fraction and was suspended in hyperosmotic solution (50 mM HEPES, pH 7.5, 250 mM NaCl, 2 mM EDTA with protease inhibitor).

Afterward, the crude membrane suspension was incubated in the presence or absence of remdesivir for 1-h, followed by incubation with biotin-UT2 overnight. The biotin-UT2-bound UTS2R complex was solubilized using a hyperosmotic buffer containing detergent (N-Dodecyl-β-D-maltopyranoside; DDM, 5× CMC) and pulled down using streptavidin magnetic beads. The suitable beads for the experiment were selected from 5 types (Supplementary Fig. 1d, Dynabeads Streptavidin Trial Kit, Thermo Fisher, and Streptavidin magnetic beads, NEB) and denatured with SDS buffer. The lysates were then subjected to SDS-PAGE followed by a western blot using an anti-FLAG antibody to detect FLAG-UTS2R (Sigma M2 monoclonal antibody, 1:2000).

A human UTS2R structure was prepared from the AlphaFold Protein Structure Database56,57 (AlphaFold Monomer v2.0) by trimming the lid-like N-terminus region (1–42) with a low pLDDT score. Hydrogen atoms were added to the trimmed receptor with the program Reduce58. Remdesivir was docked in the orthosteric pocket of the receptor by AutoDockFR59 with 50 genetic algorithm evolutions and a maximum of 2,500,000 evaluations. Docking poses were clustered at the 2 Å cutoff, and the top five clusters were inspected based on the criterion of hydrogen bonding between the receptor and remdesivir. Structures were visualized and analyzed with CueMol and PyMOL (Schrödinger, Inc.). Files used in this study were deposited to Zenodo.

For blots of pERK, ERK, pAKT, AKT, and β-actin, UTS2R-transfected and serum-starved HEK293A cells were stimulated with 0.1 or 1 μM remdesivir or 10 or 100 nM UT2. Before stimulation with these compounds, cells were incubated with urantide (Peptide Institute) for 30 min at 100 nM (for UTS2R inhibition) or with 150 ng/mL pertussis toxin (PTX; Wako) overnight (for Gαi/o inhibition). For blots of mitochondrial proteins (Total OXPHOS, MT-COI, TFAM, and VDAC), 1 or 10 μM remdesivir or 1 or 10 μM GS-441524 was added to UTS2R-transfected HEK293A cells for 48 h. Whole-cell lysates were prepared in RIPA Lysis Buffer (Thermo Fischer Scientific) containing protease inhibitors (Thermo Fischer Scientific) and phosphatase inhibitors (Nacalai Tesque). Cell lysates were sonicated and clarified by centrifugation. Protein concentrations were measured using the BCA Protein Assay Kit (Pierce) and were adjusted to 1 mg/mL, and 10 μg of proteins were subjected to western blotting for each experiment. Samples were separated using SDS polyacrylamide gels and transferred onto polyvinylidene difluoride (PVDF) membranes. The blots were blocked with 5% non-fat milk in TBS-T (Tris-buffered saline with 0.1% Tween-20) for 1 h, then probed using various primary antibodies. The primary antibodies and working dilutions used were as follows: pERK (AB_331646), 1:2000; ERK (AB_330744), 1:2000; pAKT (AB_329825), 1:1000; AKT (AB_329827), 1:1000; βactin (AB_10697039), 1:10000; total OXPHOS (AB_2756818), 1:10000; MTCO1 (AB_2084810), 1:2000; TFAM (AB_10841294), 1:1000; VDAC (AB_2272627), 1:1000. Target antigens were incubated with appropriate HRP-conjugated secondary antibodies (Cell Signaling) and were visualized by an ECL substrate (GE Healthcare). Phosphorylated protein levels were normalized against the total levels of the target protein (Figs. 3a and 4c, d, Supplementary Fig. 3a, b). For measuring mitochondrial proteins (Supplementary Fig. 3g), OXPHOS subunits, VDAC, MTCO1, and TFAM were detected by specific antibodies on different PVDF membranes, respectively. Protein levels of OXPHOS subunits, MTCO1, and TFAM were normalized to the level of VDAC.

To assess the human UTS2R (hUTS2R) mRNA expression in normal adult human tissues, commercially available cDNAs originating from various human tissues were purchased from TaKaRa (detailed sample information is given in Supplementary Data 3), and 2 ng of cDNA were subjected to quantitative PCR (qPCR) analysis (Rotor-Gene Q, Qiagen), as recommended by the manufacturers. To assess mouse UTS2R (mUTS2R) mRNA expression in normal adult mouse tissues, C57BL/6J male mice were purchased from Japan Clea, and RNA was isolated from multiple tissues. RNA samples were reverse-transcribed into cDNA and were subjected to qPCR analysis. The primers used for the relative quantification are shown in Supplementary Data 4. The expression levels of hUTS2R and mUTS2R were normalized to housekeeping genes of G3PDH and 18 s rRNA, respectively.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) were purchased from Fujifilm Cellular Dynamics, Inc. (CDI; iCell cardiomyocytes 2.0). To prepare the probe (MED-PG515A, Alpha MED Sciences), the recording areas were coated with fibronectin and dissolved in Dulbecco's phosphate-buffered saline to create a 50 µg/mL solution. The recording area was covered with 1–2 µL fibronectin solution, and the probes were incubated at 37 °C for at least 1 h. Cells were thawed and suspended in iCell Cardiomyocytes Plating Medium (CDI) and plated onto the probes at a density of 3.5×104 cells in a 2 µL plating medium. The cells were incubated at 37 °C in 5% CO2 for 3–4 h before filling each probe with 1 mL iCell cardiomyocytes maintenance medium (CDI), which was used as the culture medium. The medium was changed every 2–3 days, with the cells cultured in the probes for 5–14 days to obtain a sheet of cardiomyocytes with spontaneous and synchronous electrical activity.

FP recording was carried out35,60. Before measuring FPs, iCell cardiomyocyte sheets were equilibrated for at least 4 h in a fresh culture medium using a 5% CO2 incubator at 37 °C. After equilibration, the probes were transferred to the multi-electrode array (MEA) system (MED64, Alpha med Scientific) and incubated in a humidified 5% CO2 atmosphere. The stability and constancy of the waveforms were confirmed by monitoring the signals for at least 30 min. Once the FPs reached a stable state, they were recorded for 10 min at the baseline and then 24, 48, and 72 h after drug treatment. The stock solutions of drugs were prepared in dimethyl sulfoxide (DMSO) or distilled water at a 1000-fold target concentration. The stock solution was diluted in a culture medium to a final concentration of 0.1% DMSO.

The FP duration (FPD) was defined as the duration from the first to the second peak in FP35. The FPD was corrected for the beat rate (inter-spike interval (ISI)) with Fridericia's formula (Eq. 1), which was the primary method of correction in this study:

The ISI and FPD values were averaged from the last 30 beats or 30 beats at time points exhibiting a stable ISI and FPD.

The perforated patch-clamp technique was used to record spontaneous action potentials in current-clamp models of hiPS-CMs at 25 °C, using an EPC-10 patch–clamp amplifier (HEKA, Lambrecht, Germany)61. The hiPS-CMs typically generate spontaneous automaticity in normal Tyrode solution. Spontaneous action potentials were recorded using a fire-polished patch pipette (resistance, 2.0–4.0 MΩ) filled with a pipette solution containing 70 mM potassium aspartate, 50 mM KCl, 10 mM KH2PO4, 1 mM MgSO4, 3 mM adenosine triphosphate (ATP) (disodium salt; Sigma Chemical Company, St Louis), 5 mM EGTA, 5 mM HEPES, and 0.1 mM Li2-guanosine triphosphate (GTP) (Sigma Chemical Company) (pH adjusted to 7.2 with KOH). Amphotericin B (Wako Pure Chemical Industries) was added to the pipette solution at a concentration of 100 mg/mL. The glass coverslip (3 × 5 mm2) containing hiPS-CMs plated in 35-mm culture dishes was placed into a recording chamber with normal Tyrode solution."

Neonatal rat cardiomyocytes (NRCMs) were isolated from Sprague-Dawley rat (Japan SLC Inc.) pups on postnatal days 1–362. The ventricles of the pups were dissected and minced on ice after terminal anesthesia (1.5% sevoflurane inhalation) and euthanasia by cervical dislocation. The minced tissue was pre-digested in 0.05% trypsin-EDTA (Gibco) overnight at 4 °C and then digested in 1 mg/mL collagenase type2 (Worthington) in PBS for 30 min at 37 °C while shaking the flask (125 rpm). The dissociated cells were filtered through a cell strainer (70 μm, Falcon) and centrifuged for 2 min at 180×g. After removing the supernatant, cell pellets were resuspended in Dulbecco's modified Eagle's medium 12 (DMEM) supplemented with 10% FBS and 1% penicillin and streptomycin plated in a 10-cm culture dish. Cells were incubated at 37 °C in a humidified atmosphere (5% CO2, 95% air) for 90 min. Floating NRCMs were collected and plated to matrigel-coated culture dishes. After 24 h, the culture medium was changed to serum-free DMEM and incubated for more than 2 days before experiments. All animal experiments were reviewed and approved by the ethics committees at the National Institutes of Natural Sciences or the Animal Care and Use Committee at Kyushu University. All reported procedures conform to the NIH Guide for the Care and Use of Laboratory Animals.

NRCMs (4 × 105 cells/well) were seeded on 3.5-cm glass-bottomed dishes in serum-free DMEM. Cells were treated with remdesivir (1 μM) for 48 h at 37 °C and 5%CO2. Urantide (100 nM) was added at the same time. Pertussis toxin (PTX; 150 ng/mL) was treated 18 h before adding remdesivir, and YM-254890 (1 μM) was treated 30 minutes before remdesivir addition. Microscopic images of NRCMs were recorded for 13 s under pacing conditions. The pacing of NRCM was stimulated at 4 V every second, with a frequency of 10 ms, using an electrical stimulator (NIHON KOHDEN, SEN-3301). Imaging data were acquired at six frames per second using a BZ-X800 microscope (Keyence) and analyzed using Fiji software63.

C57BL/6J mice for the isolation of adult mouse cardiomyocytes were purchased from Japan SLC Inc. (Shizuoka, Japan). Isolation of adult mouse ventricular cardiomyocytes was performed64. Male mice were anesthetized by inhalation with isoflurane. The heart was quickly excised, and the aorta was clamped with a small vascular clamp. The heart was antegradely perfused with isolation buffer (less than 3 mL) containing 130 mM NaCl, 5.4 mM KCl, 0.5 mM MgCl2, 0.33 mM NaH2PO4, 25 mM HEPES, 22 mM glucose, and 50 μU/mL bovine insulin (Sigma) (pH 7.4, adjusted with NaOH) supplemented with 0.4 mM EGTA, followed by the perfusion of enzyme solution (10 mL) (isolation buffer containing 1 mg/mL collagenase type 2 (Worthington), 0.06 mg/mL trypsin (Sigma), 0.06 mg/mL protease (Sigma), and 0.3 mM CaCl2). Thereafter, the LV chamber was removed and cut into small pieces in an enzyme solution containing 0.2% bovine serum albumin (BSA) and 0.7 mM CaCl2. The tissue-cell suspension was filtered through a cell strainer (100 μm, Falcon) and centrifuged for 3 min at 50×g. After removing the supernatant, the cell pellet was resuspended in isolation buffer supplemented with 0.2% BSA and 1.2 mM CaCl2 and incubated for 7 min at 37 °C. After centrifugation (50×g, 3 min), the cells were resuspended in Tyrode's solution A (containing 140 mM NaCl, 5.4 mM KCl, 1.8 mM CaCl2, 0.5 mM MgCl2, 0.33 mM NaH2PO4, 5.0 mM HEPES, and 5.5 mM glucose (pH 7.4)) supplemented with 0.2% BSA. Cells were seeded onto Matrigel (Corning)-coated 35 mm glass base dishes (IWAKI) and incubated at 37 °C for 30 min. Cells were used within 1–6 h after isolation.

Isolated adult mouse cardiomyocytes were treated with remdesivir (10 μM) for 30 min at 37 °C. Microscopic images of cardiomyocytes were recorded for 13 seconds under pacing conditions. The pacing of cardiomyocytes was stimulated at 30 V every second with a frequency of 10 msec using an electrical stimulator (NIHON KOHDEN, SEN-3301). Imaging data were analyzed using Fiji software.

Statistical analyses were performed with Prism Software (GraphPad), and methods are described in the figure legends. Symbols are mean values, and error bars denote the standard error of the mean (SEM). Representative results from at least three independent experiments are shown for every figure unless stated otherwise in the figure legends.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

All data are available in the main text or the supplementary materials. Source data underlying figures are provided in Supplementary Data 1. Uncropped blots are presented in Supplementary Fig. 5. Molecular docking simulation files were deposited to Zenodo.

Seley-Radtke, K. L. & Yates, M. K. The evolution of nucleoside analogue antivirals: a review for chemists and non-chemists. Part 1: early structural modifications to the nucleoside scaffold. Antiviral Res. 154, 66–86 (2018).

Article CAS PubMed PubMed Central Google Scholar

Chien, M. et al. Nucleotide analogues as inhibitors of SARS-CoV-2 polymerase, a key drug target for COVID-19. J. Proteome Res. 19, 4690–4697 (2020).

Article CAS PubMed PubMed Central Google Scholar

Alanazi, A. S., James, E. & Mehellou, Y. The ProTide prodrug technology: where next? ACS Med. Chem. Lett. 10, 2–5 (2019).

Article CAS PubMed Google Scholar

Sheahan, T. P. et al. Broad-spectrum antiviral GS-5734 inhibits both epidemic and zoonotic coronaviruses. Sci. Transl. Med. 9, eaal3653 (2017).

Article PubMed PubMed Central Google Scholar

Siegel, D. et al. Discovery and synthesis of a phosphoramidate prodrug of a Pyrrolo[2,1-f][triazin-4-amino] adenine C-nucleoside (GS-5734) for the treatment of Ebola and emerging viruses. J. Med. Chem. 60, 1648–1661 (2017).

Article CAS PubMed PubMed Central Google Scholar

Warren, T. K. et al. Therapeutic efficacy of the small molecule GS-5734 against Ebola virus in rhesus monkeys. Nature 531, 381–385 (2016).

Article CAS PubMed PubMed Central Google Scholar

Beigel, J. H. et al. Remdesivir for the treatment of Covid-19—final report. N. Engl. J. Med. 383, 1813–1826 (2020).

Article CAS PubMed Google Scholar

Gottlieb, R. L. et al. Early remdesivir to prevent progression to severe Covid-19 in outpatients. N. Engl. J. Med. 386, 305–315 (2022).

Article CAS PubMed Google Scholar

Coronavirus Disease 2019 (COVID-19) Treatment Guidelines. (2022).

Jung, S. Y. et al. Cardiovascular events and safety outcomes associated with remdesivir using a World Health Organization international pharmacovigilance database. Clin. Transl. Sci. 15, 501–513 (2022).

Article CAS PubMed Google Scholar

Haghjoo, M. et al. Effect of COVID-19 medications on corrected QT interval and induction of torsade de pointes: Results of a multicenter national survey. Int. J. Clin. Pract. 75, e14182 (2021).

Article CAS PubMed Google Scholar

Liu, D. et al. Adverse cardiovascular effects of anti-COVID-19 drugs. Front. Pharmacol. 12, 699949 (2021).

Article CAS PubMed PubMed Central Google Scholar

Gupta, A. K., Parker, B. M., Priyadarshi, V. & Parker, J. Cardiac adverse events with remdesivir in COVID-19 infection. Cureus 12, e11132 (2020).

PubMed PubMed Central Google Scholar

Hu, W. J. et al. Pharmacokinetics and tissue distribution of remdesivir and its metabolites nucleotide monophosphate, nucleotide triphosphate, and nucleoside in mice. Acta Pharmacol. Sin. 42, 1195–1200 (2021).

Article CAS PubMed Google Scholar

Gordon, C. J., Tchesnokov, E. P., Schinazi, R. F. & Götte, M. Molnupiravir promotes SARS-CoV-2 mutagenesis via the RNA template. J. Biol. Chem. 297, 100770 (2021).

Article CAS PubMed PubMed Central Google Scholar

Jayk Bernal, A. et al. Molnupiravir for oral treatment of Covid-19 in nonhospitalized patients. N. Engl. J. Med. 386, 509–520 (2022).

Article PubMed Google Scholar

Furuta, Y. et al. T-705 (favipiravir) and related compounds: Novel broad-spectrum inhibitors of RNA viral infections. Antiviral Res. 82, 95–102 (2009).

Article CAS PubMed PubMed Central Google Scholar

Padhi, A. K., Dandapat, J., Saudagar, P., Uversky, V. N. & Tripathi, T. Interface-based design of the favipiravir-binding site in SARS-CoV-2 RNA-dependent RNA polymerase reveals mutations conferring resistance to chain termination. FEBS Lett. 595, 2366–2382 (2021).

Article CAS PubMed PubMed Central Google Scholar

Wen, W. et al. Efficacy and safety of three new oral antiviral treatment (molnupiravir, fluvoxamine and Paxlovid) for COVID-19: a meta-analysis. Ann. Med. 54, 516–523 (2022).

Article CAS PubMed PubMed Central Google Scholar

Hung, D. T. et al. The efficacy and adverse effects of favipiravir on patients with COVID-19: a systematic review and meta-analysis of published clinical trials and observational studies. Int. J. Infect. Dis. 120, 217–227 (2022).

Article CAS PubMed PubMed Central Google Scholar

Jacobson, K. A. & Gao, Z. G. Adenosine receptors as therapeutic targets. Nat Rev Drug Discov 5, 247–264 (2006).

Article CAS PubMed PubMed Central Google Scholar

Ogawa, A. et al. N(6)-methyladenosine (m(6)A) is an endogenous A3 adenosine receptor ligand. Mol. Cell 81, 659–674.e657 (2021).

Article CAS PubMed Google Scholar

Rengo, G., Lymperopoulos, A. & Koch, W. J. Future g protein-coupled receptor targets for treatment of heart failure. Curr. Treat. Options Cardiovasc. Med. 11, 328–338 (2009).

Article PubMed Google Scholar

Sato, M. & Ishikawa, Y. Accessory proteins for heterotrimeric G-protein: implication in the cardiovascular system. Pathophysiology 17, 89–99 (2010).

Article CAS PubMed Google Scholar

Wacker, D., Stevens, R. C. & Roth, B. L. How ligands illuminate GPCR molecular pharmacology. Cell 170, 414–427 (2017).

Article CAS PubMed PubMed Central Google Scholar

Inoue, A. et al. TGFα shedding assay: an accurate and versatile method for detecting GPCR activation. Nat. Methods 9, 1021–1029 (2012).

Article CAS PubMed Google Scholar

Summary-compassionate-use-remdesivir-gilead_en.pdf. (2020).

Du, A. T. et al. Ligand-supported purification of the urotensin-II receptor. Mol. Pharmacol. 78, 639–647 (2010).

Article CAS PubMed Google Scholar

Castel, H. et al. The G protein-coupled receptor UT of the neuropeptide urotensin II displays structural and functional chemokine features. Front. Endocrinol. 8, 76 (2017).

Article Google Scholar

Ballesteros, J. A. & Weinstein, H. Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors. Methods Neurosci. 25, 366–428 (1995).

Article CAS Google Scholar

Biswal, H. S., Gloaguen, E., Loquais, Y., Tardivel, B. & Mons, M. Strength of NH···S hydrogen bonds in methionine residues revealed by gas-phase IR/UV spectroscopy. J. Phys. Chem. Lett. 3, 755–759 (2012).

Article CAS PubMed Google Scholar

Pokkuluri, P. R. et al. Factors contributing to decreased protein stability when aspartic acid residues are in beta-sheet regions. Protein Sci. 11, 1687–1694 (2002).

Article CAS PubMed PubMed Central Google Scholar

Pereira-Castro, J., Brás-Silva, C. & Fontes-Sousa, A. P. Novel insights into the role of urotensin II in cardiovascular disease. Drug Discov. Today 24, 2170–2180 (2019).

Article CAS PubMed Google Scholar

Mummery, C. L. et al. Differentiation of human embryonic stem cells and induced pluripotent stem cells to cardiomyocytes: a methods overview. Circ. Res. 111, 344–358 (2012).

Article CAS PubMed PubMed Central Google Scholar

Yanagida, S., Satsuka, A., Hayashi, S., Ono, A. & Kanda, Y. Comprehensive cardiotoxicity assessment of COVID-19 treatments using human-induced pluripotent stem cell-derived cardiomyocytes. Toxicol. Sci. 183, 227–239 (2021).

Article CAS PubMed Google Scholar

Cubeddu, L. X. QT prolongation and fatal arrhythmias: a review of clinical implications and effects of drugs. Am. J. Ther. 10, 452–457 (2003).

Article PubMed Google Scholar

Wettschureck, N. & Offermanns, S. Mammalian G proteins and their cell type specific functions. Physiol. Rev. 85, 1159–1204 (2005).

Article CAS PubMed Google Scholar

Nishida, M. et al. G alpha(i) and G alpha(o) are target proteins of reactive oxygen species. Nature 408, 492–495 (2000).

Article CAS PubMed Google Scholar

Kwok, M. et al. Remdesivir induces persistent mitochondrial and structural damage in human induced pluripotent stem cell derived cardiomyocytes. Cardiovasc. Res. 118, 2652–2664 (2022).

Tchesnokov, E. P., Feng, J. Y., Porter, D. P. & Götte, M. Mechanism of inhibition of Ebola virus RNA-dependent RNA polymerase by remdesivir. Viruses 11, 326 (2019).

Article CAS PubMed PubMed Central Google Scholar

Tadaka, S. et al. jMorp updates in 2020: large enhancement of multi-omics data resources on the general Japanese population. Nucleic Acids Res. 49, D536–d544 (2021).

Article CAS PubMed Google Scholar

Lasser, K. E. et al. Timing of new black box warnings and withdrawals for prescription medications. J. Am. Med. Assoc. 287, 2215–2220 (2002).

Article Google Scholar

Michaud, V. et al. Risk assessment of drug-induced long QT syndrome for some COVID-19 repurposed drugs. Clin. Transl. Sci. 14, 20–28 (2021).

Article CAS PubMed Google Scholar

Nabati, M. & Parsaee, H. Potential cardiotoxic effects of remdesivir on cardiovascular system: a literature review. Cardiovasc. Toxicol. 22, 268–272 (2022).

Article CAS PubMed Google Scholar

Day, L. B., Abdel-Qadir, H. & Fralick, M. Bradycardia associated with remdesivir therapy for COVID-19 in a 59-year-old man. Can. Med. Assoc. J. 193, E612–e615 (2021).

Article CAS Google Scholar

Sanchez-Codez, M. I., Rodriguez-Gonzalez, M. & Gutierrez-Rosa, I. Severe sinus bradycardia associated with Remdesivir in a child with severe SARS-CoV-2 infection. Eur. J. Pediatr. 180, 1627 (2021).

Article CAS PubMed PubMed Central Google Scholar

Douglas, S. A., Tayara, L., Ohlstein, E. H., Halawa, N. & Giaid, A. Congestive heart failure and expression of myocardial urotensin II. Lancet 359, 1990–1997 (2002).

Article CAS PubMed Google Scholar

Adamsick, M. L. et al. Remdesivir in patients with acute or chronic kidney disease and COVID-19. J. Am. Soc. Nephrol. 31, 1384–1386 (2020).

Article CAS PubMed PubMed Central Google Scholar

Insel, P. A., Tang, C. M., Hahntow, I. & Michel, M. C. Impact of GPCRs in clinical medicine: monogenic diseases, genetic variants and drug targets. Biochim. Biophys. Acta 1768, 994–1005 (2007).

Article CAS PubMed Google Scholar

Lefkowitz, R. J. & Shenoy, S. K. Transduction of receptor signals by beta-arrestins. Science 308, 512–517 (2005).

Article CAS PubMed Google Scholar

Esposito, G. et al. EGFR trans-activation by urotensin II receptor is mediated by β-arrestin recruitment and confers cardioprotection in pressure overload-induced cardiac hypertrophy. Basic Res. Cardiol. 106, 577–589 (2011).

Article CAS PubMed Google Scholar

Korhonen, T., Hänninen, S. L. & Tavi, P. Model of excitation-contraction coupling of rat neonatal ventricular myocytes. Biophys. J. 96, 1189–1209 (2009).

Article CAS PubMed PubMed Central Google Scholar

Sala, L. et al. MUSCLEMOTION: a versatile open software tool to quantify cardiomyocyte and cardiac muscle contraction in vitro and in vivo. Circ Res. 122, e5–e16 (2018).

Article CAS PubMed PubMed Central Google Scholar

Sanguinetti, M. C. & Tristani-Firouzi, M. hERG potassium channels and cardiac arrhythmia. Nature 440, 463–469 (2006).

Article CAS PubMed Google Scholar

Hisano, Y. et al. Lysolipid receptor cross-talk regulates lymphatic endothelial junctions in lymph nodes. J. Exp. Med. 216, 1582–1598 (2019).

Article CAS PubMed PubMed Central Google Scholar

Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589 (2021).

Article CAS PubMed PubMed Central Google Scholar

Varadi, M. et al. AlphaFold Protein Structure Database: massively expanding the structural coverage of protein-sequence space with high-accuracy models. Nucleic Acids Res. 50, D439–d444 (2022).

Article CAS PubMed Google Scholar

Word, J. M., Lovell, S. C., Richardson, J. S. & Richardson, D. C. Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation. J. Mol. Biol. 285, 1735–1747 (1999).

Article CAS PubMed Google Scholar

Ravindranath, P. A., Forli, S., Goodsell, D. S., Olson, A. J. & Sanner, M. F. AutoDockFR: advances in protein-ligand docking with explicitly specified binding site flexibility. PLoS Comput. Biol. 11, e1004586 (2015).

Article PubMed PubMed Central Google Scholar

Ando, H. et al. A new paradigm for drug-induced torsadogenic risk assessment using human iPS cell-derived cardiomyocytes. J. Pharmacol. Toxicol. Methods 84, 111–127 (2017).

Article CAS PubMed Google Scholar

Ding, W. G. et al. Identification of verapamil binding sites within human Kv1.5 channel using mutagenesis and docking simulation. Cell Physiol. Biochem. 52, 302–314 (2019).

Article CAS PubMed Google Scholar

Kitajima, N. et al. TRPC3 positively regulates reactive oxygen species driving maladaptive cardiac remodeling. Sci. Rep. 6, 37001 (2016).

Article CAS PubMed PubMed Central Google Scholar

Matsuda, Y., Takahashi, K., Kamioka, H. & Naruse, K. Human gingival fibroblast feeder cells promote maturation of induced pluripotent stem cells into cardiomyocytes. Biochem. Biophys. Res. Commun. 503, 1798–1804 (2018).

Article CAS PubMed Google Scholar

Omatsu-Kanbe, M., Yoshioka, K., Fukunaga, R., Sagawa, H. & Matsuura, H. A simple antegrade perfusion method for isolating viable single cardiomyocytes from neonatal to aged mice. Physiol. Rep. 6, e13688 (2018).

Article PubMed PubMed Central Google Scholar

Download references

We are grateful for suggestions and technical support from all members of the Department of Modomics Biology & Medicine, Tohoku University. We thank Y. Takahata and H. Miyamoto for technical assistance and Natalie D. DeWitt for critical reading and language editing. This work was supported by funding provided by JSPS KAKENHI (Grant nos. JP20K18371, JP22H04628, and JP22H02813 to A.O., JP21H05269 and JP22H02772 to M.N., JP21H04791 and JP21H051130 to A.I., JP21H02659 and JP21H05265 to F.Y.W., and JP21H02634 to Y. Kanda. This work was also supported by the AMED (Grant nos. JP21mk0101189 to Y. Kanda, JP22zf0127007 to A.I., 21he2302008j0002 to F.Y.W., BINDS JP22ama121031 to M.N. and JP20am0101095 to A.I., LEAP JP20gm0010004 to A.I.), FOREST from JST (Grant nos. JPMJFR220K to A.O., JPMJFR215T to A.I., JPMJFR205Y to F.Y.W.), Moonshot R&D JPMJMS2023 from JST to A.I., CREST JPMJCR2024 (20348438) from J.S.T. to M.N., ERATO JPMJER2002 from J.S.T. to F.Y.W., Daiichi Sankyo Foundation of Life Science to A.I., Takeda Science Foundation to A.O. and A.I., Uehara Memorial Foundation to A.O. and A.I., Inamori Foundation to A.O., the Mitsubishi Foundation to A.O., Terumo Life Science Foundation to A.O., and the Naito Foundation to A.O.

These authors contributed equally: Akiko Ogawa, Seiya Ohira.

Department of Modomics Biology and Medicine, Institute of Development, Aging and Cancer (IDAC), Tohoku University, Sendai, Miyagi, 980-8575, Japan

Akiko Ogawa, Seiya Ohira & Fan-Yan Wei

Graduate School of Medicine, Tohoku University, Sendai, Miyagi, 980-8575, Japan

Seiya Ohira

Department of Physiology, Graduate School of Pharmaceutical Sciences, Kyushu University, Fukuoka, 812-8582, Japan

Yuri Kato, Xinya Mi, Yukina Ishii & Motohiro Nishida

Laboratory of Molecular & Cellular Biochemistry, Graduate School of Pharmaceutical Sciences, Tohoku University, 6-3, Aoba, Aramaki, Aoba-ku, Sendai, Miyagi, 980-8578, Japan

Tatsuya Ikuta & Asuka Inoue

Division of Pharmacology, National Institute of Health Sciences, Kanagawa, 210-9501, Japan

Shota Yanagida & Yasunari Kanda

Division of Pharmaceutical Sciences, Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan

Shota Yanagida

National Institute for Physiological Sciences and Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki, 444-8787, Japan

Motohiro Nishida

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

The authors confirm their contribution to the paper as follows: study conceptualization and design: A.O., M.N., A.I., and F.Y.W.; data collection: A.O., S.O., Y. Kato, T.I., S.Y., X.M., Y.I., and A.I.; analysis and interpretation of results: A.O., Y. Kato, T.I., S.Y., X.M., Y. Kanda, M.N., A.I., and F.Y.W. draft paper preparation: A.O., Y. Kato, T.I., S.Y., X.M., A.I., F.Y.W. project supervision: Y. Kanda, M.N., A.I., F.Y.W. All authors reviewed the results and approved the final version of the paper.

Correspondence to Motohiro Nishida, Asuka Inoue or Fan-Yan Wei.

The authors declare no competing interests.

Communications Biology thanks Dipayan Chaudhuri and the other anonymous reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: George Inglis.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Ogawa, A., Ohira, S., Kato, Y. et al. Activation of the urotensin-II receptor by remdesivir induces cardiomyocyte dysfunction. Commun Biol 6, 511 (2023). https://doi.org/10.1038/s42003-023-04888-x

Download citation

Received: 10 November 2022

Accepted: 28 April 2023

Published: 12 May 2023

DOI: https://doi.org/10.1038/s42003-023-04888-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.