- 1Key Laboratory of Medical Electrophysiology, Ministry of Education and Medical Electrophysiological Key Laboratory of Sichuan Province, (Collaborative Innovation Center for Prevention of Cardiovascular Diseases), Institute of Cardiovascular Research, Southwest Medical University, Luzhou, China
- 2Biological Physics Group, Department of Physics and Astronomy, The University of Manchester, Manchester, United Kingdom
- 3Beijing Institute of Artificial Intelligence, Beijing, China
Remdesivir (RDV) is the first drug approved by the FDA for clinical treatment of hospitalized patients infected with COVID-19 because it has been shown to have good antiviral activity against a variety of viruses, including Arenaviridae and Coronaviridae viral families. However, it has been reported that its clinical treatment leads to the symptoms of sick sinus syndrome such as sinus bradycardia, conduction block, and sinus arrest, but the electrophysiological mechanism of its specific cardiac adverse events is still unclear. We report complementary, experimental, studies of its electrophysiological effects. In wireless cardiac telemetry experiments in vivo and electrocardiographic studies in ex vivo cardiac preparations, RDV significantly caused sinus bradycardia, sinus atrial block, and prolongation of the QT interval in guinea pigs. Dose-dependent effects of RDV on the electrical activities of sinoatrial node (SA node) preparations of guinea pigs were characterised by multielectrode, optical RH237 voltage mapping. These revealed reversibly reduced sinoatrial conduction time (SACT), increased AP durations (APDs), and decreased the pacemaking rate of the SA node. Patch-clamp experiments showed that RDV significantly inhibited the If current of HCN4 channels, resulting in a significant decrease in the spontaneous firing rate of SA node cells, which may underlie the development of sick sinus node syndrome. In addition, RDV significantly inhibits IKr currents in hERG channels, leading to prolongation of the QT interval and playing a role in bradycardia. Therefore, these findings provide insights into the understanding the bradycardia effect of RDV, which may be used as basic theoretical guidance for the intervention of its adverse events, and prompt safety investigations of RDV’s cardiac safety in the future.
1 Introduction
Remdesivir (RDV) is an adenosine triphosphate analogue first synthetically developed as a candidate for the treatment of hepatitis C virus and first described in the literature as a potential treatment for Ebola in 2016 (Warren et al., 2016). The mechanism of action of RDV suggests a broad antiviral activity spectrum. To date, in vitro studies have demonstrated its efficacy against multiple viral families including Filoviridae, Paramyxoviridae, Arenaviridae, Flaviviridae, Pneumoviridae, and Coronaviridae (Malin et al., 2020). Initial observations of RDV’s efficacy against the Coronaviridae family in 2017 sparked significant interest in its potential as a treatment for COVID-19 (Sheahan et al., 2017). However, concerns have been raised about its clinical safety in treating COVID-19. Nucleosides have long been associated with cardiovascular side effects, altering cardiomyocyte action potentials and potentially leading to adverse cardiac effects such as conduction block, atrial fibrillation, and ventricular fibrillation (Ahmad et al., 2015; Antonioli et al., 2023). Additionally, in vitro experiments have demonstrated that RDV can cause mitochondrial damage in cardiomyocytes (Bjork and Wallace, 2021; Kwok et al., 2022).
As COVID-19 escalates into an unparalleled global public health crisis, there is an urgent demand for efficacious drugs to treat clinically severe cases (Jeong et al., 2023). Research indicates that RDV can target the non-specific chain terminator of RdRp (RNA-dependent RNA polymerase) in SARS-CoV-2, as well as its related strains SARS-CoV and MERS-CoV, thereby exhibiting antiviral effects (De Wit et al., 2020). And studies have demonstrated clinical benefits of RDV treatment in rhesus monkeys infected with SARS-CoV-2 (Williamson et al., 2020). Consequently, RDV became the first drug to receive emergency authorization from the U.S. FDA for the treatment of hospitalized patients with COVID-19 under these conditions. Among 2,186 U.S. adult patients with laboratory-confirmed SARS-CoV-2 infection, RDV showed potential benefits for COVID-19 treatment (Rivera et al., 2020), and it also demonstrated a reduction in recovery time for hospitalized COVID-19 patients with lower respiratory tract infections (Williamson et al., 2020). However, with the clinical use of RDV, a significant number of cardiac adverse reactions have been reported. In a prospective longitudinal study involving COVID-19 patients, those receiving RDV treatment had an overall incidence of bradycardia (heart rate < 60 bpm) of 27%, with 19% of them experiencing severe bradycardia (heart rate < 50 bpm) (Hajimoradi et al., 2023). Al-Jammali et al. (2022) conducted a comprehensive review of published literature and a meta-analysis study focusing on RDV’s impact on SARS-CoV-2 patients, particularly those who experienced arrhythmias post RDV administration. Their findings revealed an elevated risk of sinus bradycardia associated with RDV treatment compared to control groups. The meta-analysis indicated a bradycardia incidence of 34.07% in RDV-treated patients, notably higher than the control group’s 18.13% incidence. Furthermore, a joint study by the FDA and the Academy of Medical Toxicology reported that among 760 patients, 220 developed sinus bradycardia post RDV use, with 181 of these cardiac adverse reactions directly attributed to RDV. The observed cardiac adverse events included sinus bradycardia, prolonged QRS/QT interval, sinus arrest, and cardiac arrest (Devgun et al., 2023). In summary, RDV’s association with sick sinus syndrome is confirmed, yet the precise electrophysiological mechanism triggering adverse cardiac reactions remains unclear.
Currently, research on RDV primarily focuses on clinical observations, with limited studies on its mechanisms. In this study, we first utilized wireless telemetry electrocardiography and ex vivo cardiac electrophysiology to examine RDV’s impact on in vivo and ex vivo guinea pig cardiac electrophysiology, revealing consistent results at both animal and ion cardiac levels with known clinical observations of the drug. Subsequently, we studied RDV’s effects on sinoatrial node (SA node) tissue pacing, action potential conduction, and APD duration under sinus rhythm-driven conditions, which have proven increasingly valuable in elucidating mechanisms of cardiac adverse reactions and possess significance in drug safety assessments, providing insights not attainable solely from in vitro single-cell single-channel studies (Strauss et al., 2019). Finally, we investigated RDV’s influence on acutely isolated SA node cells and the pivotal pacemaker ion channel HCN4 channel, while also examining RDV’s impact on hERG channel IKr currents in HEK293 cells, elucidating the drug’s actions at cellular and ion channel levels. Therefore, through the investigation of RDV’s impact on guinea pig cardiac electrophysiology, we can elucidate the clinical mechanisms behind RDV-induced sick sinus node syndrome and its implications for drug safety.
2 Methods
The experiments included the following: (1) Telemetric electrocardiography in healthy guinea pigs to investigate RDV’s impact on guinea pig cardiac electrophysiology in vivo. (2) Ex vivo electrocardiography studies on RDV’s effects on guinea pig ex vivo cardiac electrophysiology under Langendorff perfusion, further complementing the in vivo electrocardiography results. (3) Optical mapping and electrocardiographic techniques to study RDV’s effects on pacemaking signals, action potential conduction and duration in acutely isolated SA node tissues. (4) Patch clamp techniques to investigate RDV’s effects on acute isolated guinea pig SA node cells and the key ion channels involved in pacemaking. In experiments, we ensured the consistency of variables such as animal weight, procedural methods, experimental conditions, and analytical techniques. By minimizing variability across different experimental protocols, we enhanced the reliability of our results. Finally, we analysed findings from various levels to evaluate the electrophysiological mechanisms by which the drug of RDV induces sick sinus syndrome.
2.1 Ethics and study approval
All experimental procedures were approved by the Institutional Animal Care and Use Committee of Southwest Medical University and complied with the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health. Healthy female guinea pigs weighing 250–350 g were provided by Chongqing Tengxin Huafu Experimental Animal Sales Co., Ltd. They were housed at the Animal Center of Southwest Medical University with ad libitum access to food and water. Animals were euthanized humanely by intraperitoneal injection of 3% pentobarbital sodium (50 mg/kg).
2.2 Wireless cardiac telemetry experiments
The guinea pigs were placed in an anesthesia chamber and anesthetized with 3%–4% isoflurane. Once adequately anesthetized, the guinea pigs were transferred to the surgical table in a prone position and maintained under anesthesia using a face mask. The placement of the wireless sensor on the animal’s back, at the surgical incision site, was determined based on the guinea pig’s body size and electrode length. After local disinfection, the skin was incised with scissors, and two subcutaneous tunnels were created using a blunt needle from the pre-placed sensor site to the lead position. The electrode lead was then threaded through the tunnels, and the electrode ring was sutured onto the subcutaneous fascia at the lead position. Following successful electrode placement, wound closure and disinfection were completed, and the guinea pigs were housed in individual cages for recovery. They were used for electrocardiogram experiments 1 week after surgery. A backpack was installed to facilitate real-time electrocardiogram recording, capturing baseline ECGs for 24 h before and after intraperitoneal administration of RDV. Throughout the experiment, the animals had ad libitum access to food and water. In the in vivo electrocardiogram experiment, following clinical drug absorption patterns and treatment regimens, intraperitoneal drug administration was used to study drug-induced cardiac adverse reactions. The dosing regimen for this animal experiment spanned 5 days, with recordings taken for 24 h before drug administration, 24 h after the first dose, 24 h after the fifth dose, and 24 h after a recovery period of 3 days post-treatment cessation.
2.3 Langendorff-perfused isolated hearts
Guinea pigs (female; 250–350 g) were humanely killed by intraperitoneal injection with pentobarbital sodium (50 mg/kg). Hearts rapidly excised after thoracotomy were mounted onto a Langendorff perfusion system and perfused with a modified Tyrode’s solution (119 NaCl, 25 NaHCO3, 4 KCl, 1.2 KH2PO4, 1 MgCl2, 1.8 CaCl2·2H2O, and 10 D-glucose (in mM) equilibrated with 5% CO2 and 95% O2) with the flow rate of 8 mL/min at 37°C. Hearts were perfused and monitored for stability for 20 min before experimental procedures commenced.
2.4 Measurement of heart rate in Langendorff perfusion-maintained ex vivo hearts
ECG was recorded via two silver electrodes attached to the ventricular apex and to the aortic cannula. An equilibration period of approximately 20 min was allowed to ensure stable ECG recordings. The hearts were then successively perfused with RDV at concentrations of 0.1, 0.3, 1, 3 and 10 μM for 15 min with intervening washout periods of 15 min. It should be noted that the effect of RDV on heart rate was almost reversed during the washout period.
2.5 Optical mapping of sinoatrial node region
After the Langendorff-perfused hearts reached steady state, contraction artifacts were mimimized using blebbistatin (10 μM). RH237 (1 μg/mL) were perfused to enable membrane potential at 37°C. The hearts were dissected to prepare preparations containing left and right atria in 37°C bath. ECG electrodes were placed on the right and left auricles to record the pacemaking rate of the SA node tissue. An equilibration period of approximately 10 min was allowed to ensure the stable pacemaking recordings. RH237 (1 μg/mL) were perfused to enable membrane potential again. Two 530 nm LEDs were used to illuminate the heart after their emissions were bandpass filtered (wavelengths 530 ± 20 nm) to minimize stray excitation light reaching the dyes. The fluorescence light was passed through a 550 nm long-pass filter and then a dichroic mirror with a cutoff of 638 nm. Fluorescence light with wavelengths above 638 nm was passed through a 700 nm long-pass filter and then imaged by the camera for recording voltage signals. The cameras of the optical mapping system, LED lights, and ECG recording were simultaneously driven using OMapRecord 4.0 software. For the analysis of optical mapping signals and generation of isochronal maps, data were semiautomatically processed using Electromap software. Furthermore, the interval between the initial pacing center point of the sinoatrial node and crista terminalis, during which the action potential travels, was designated as the sinus atrial conduction time (SACT).
2.6 SA node cell preparation
Single SA node cells were isolated using an enzymatic dispersion procedure similar to that described previously (Guo et al., 1997; Kojima et al., 2012). Briefly, guinea-pigs were deeply anaesthetized with sodium pentobarbital (i.p., 120 mg kg−1), and then the chest cavity was opened under artificial respiration. The ascending aorta was cannulated in situ and the heart was then excised and retrogradely perfused via the aortic cannula on a Langendorff apparatus at 37°C, initially for 4 min with normal Tyrode solution containing (in mM) 140 NaCl, 5.4 KCl, 1.8 CaCl2, 0.5 MgCl2, 0.33 NaH2PO4, 5.5 glucose and 5 HEPES (pH adjusted to 7.4 with NaOH), and then for 4 min with nominally Ca2+-free Tyrode solution (made by simply omitting CaCl2 from the normal Tyrode solution). This was followed by 8–12 min of perfusion with nominally Ca2+-free Tyrode solution containing 0.4 mg mL−1 collagenase (Wako Pure Chemical Industries, Osaka, Japan). All these solutions were oxygenated with 100% O2. The digested heart was then removed from the Langendorff apparatus, and the SA node region, bordered by the crista terminalis, the intra-atrial septum, and superior and inferior vena cava, was dissected out and cut perpendicular to the crista terminalis into small strips measuring about 0.5 mm in width. These SA node tissue strips were further digested at 37°C for 20 min with nominally Ca2+-free Tyrode solution containing 1.0 mg mL−1 collagenase and 0.1 mg mL−1 elastase. Finally, the enzymatically digested SA node strips were gently agitated in a high-K+, low-Cl− Kraftbrühe (KB) solution, containing (in mM) 70 K-glutamate, 30 KCl, 10 KH2PO4, 1 MgCl2, 20 taurine, 0.3 EGTA, 10 glucose and 10 HEPES (pH adjusted to 7.2 with KOH) (Isenberg and Klockner, 1982). The isolated cells thus obtained were then stored at 4°C in the KB solution for experimental use within 8 h.
2.7 Whole-cell patch-clamp recordings
Perforated and conventional (ruptured) whole-cell patch-clamp techniques were used to record the spontaneous action potentials, If currents and IKr currents in the current- and voltage-clamp modes respectively (Hamill et al., 1981; Horn and Marty, 1988). An EPC-10 patch-clamp amplifier was used for these recordings. The fire-polished patch electrodes had resistance of 3.0–6.0 MΩ when filled with the pipette solution. In experiments, the cells adhered to the Round Coverslip were transferred to the recording chamber (volume of 0.5 mL) mounted on the stage of an inverted microscope. The chamber was perfused with standard Tyrode’s solution at a temperature of 35°C–37°C. The cells included cultured HEK293 cells stably transfected with hERG channels and acutely isolated SA node cells. The SA node cells dissociated from the entire SA node region comprised morphologically heterogeneous cells, such as spindle- and spider-shaped cells. Because spindle-shaped cells closely resemble the primary pacemaker cells in the SA node (Belardinelli et al., 1988), such cells displaying regular contraction were selected for the present experiments to assess the effect of RDV on SA node electrical activity.
Spontaneous action potentials were recorded in normal Tyrode solution with a pipette solution containing (in mM) 70 potassium aspartate, 50 KCl, 10 KH2PO4, 1 MgSO4 and 5 HEPES (pH adjusted to 7.2 with KOH). Amphotericin B (150 μg·mL−) was added to the pipette solution just before use and measurements were started 5–10 min after giga-seal formation. If was recorded in normal Tyrode solution with a K+-rich pipette solution containing (in mM) 70 potassium aspartate, 50 KCl, 10 KH2PO4, 1 MgSO4, 5 ATP, 0.1 GTP, 5 EGTA, 1.2 CaCl2 and 5 HEPES (pH adjusted to 7.2 with KOH). If was activated by hyperpolarizing voltage-clamp steps applied from a holding potential of −40 mV to test potentials in 10 mV decrements from −40 to −130 mV for a duration of 2.5 s, followed by a return to −40 mV. IKr was recorded in normal Tyrode solution with a pipette solution containing (in mM): KOH 110, KCl 40, K2ATP 5, MgCl2 5, EGTA 5, HEPES 10 (pH was adjusted to 7.2 by KOH). The hERG current was assessed using the patch-clamp technique with the membrane potential clamped at −80 mV. Depolarizing voltage steps were applied in 10 mV increments from −60 mV to +60 mV for a duration of 3 s, followed by a return to −50 mV.
2.8 Statistical analysis
In vivo electrocardiographic data were analyzed using LabChart 7.0 (ADInstruments). Ex vivo cardiac electrograms were analyzed using Spike 2 software (CED). Data from optogenetic experiments were analyzed using ElectroMap, with action potential durations (APDs) measured as the time required for 50% and 80% repolarization. Raw current data from patch clamp experiments were analyzed using Clampfit 10.2 (Molecular Devices, United States). All data are presented as mean ± standard deviation (SD). Comparisons between the two groups were performed using the unpaired Student’s t-test. One-way ANOVA, followed by post hoc tests, were performed to compare multiple groups using GraphPad Prism8 software (GraphPad Software, San Diego, CA, United States). Statistical significance was set at p < 0.05. In the figures, the designations for p values are *p < 0.05, **p < 0.01 and ***p < 0.001 respectively.
3 Results
3.1 Wireless cardiac telemetry experiments in guinea pigs
Firstly, we investigated the impact of RDV on in vivo guinea pig cardiac electrophysiology using wireless telemetry electrocardiography. Referencing clinical RDV treatment concentrations and regimens for COVID-19 patients, we administered RDV via intraperitoneal injection to guinea pigs (n = 7) for 5 days, with an initial dose of 23.6 mg/kg on day 1 and 11.8 mg/kg on days 2–5. We recorded the guinea pigs’ normal cardiac electrophysiology before administration, as well as on the first and fifth days of administration, and the third day after the last dose. In guinea pigs, the normal heart rate in control conditions ranges from 200 to 360 beats per minute (bpm). According to Haq et al. (2022), young healthy guinea pigs exhibit a heart rate not lower than 230 bpm. Therefore, heart rates below 200 bpm can be classified as bradycardia.
Figure 1A shows that compared to baseline, RDV administration significantly reduced heart rates in guinea pigs on the first and fifth days, with heart rates returning to normal after cessation of RDV. Figure 1B illustrates that RDV led to severe sinus bradycardia in guinea pigs, with the lowest heart rate dropping to 60 BPM. Additionally, we conducted statistical analysis on guinea pig heart rates 24 h before and after administration. As shown in Figure 1C, as compared to controls, intraperitoneal RDV significantly decreased sinus heart rates on the first and fifth days. Additionally, two guinea pigs experienced sudden cardiac arrest and death following the conclusion of the administration period. Therefore, the in vivo cardiac electrophysiology results in guinea pigs suggest that RDV significantly inhibits sinus node pacing function, leading to a marked decrease in sinus heart rate. It may also cause sinus node organic lesions, preventing recovery to the optimal state before administration and potentially resulting in sudden cardiac arrest after discontinuation.
Figure 1. Effect of RDV on electrocardiography in guinea pigs. (A) Representative electrocardiogram of guinea pigs before and after RDV administration and recovery. (B) Representative segment showing the severe decrease in heart rate of guinea pigs induced by RDV. (C) Heart rate statistics of guinea pigs before and after RDV administration (n = 7). (D) Representative electrocardiogram of guinea pigs before and after RDV administration. (E) Statistical chart of PR interval. (F) Statistical chart of QRS. (G) Statistical chart of RR interval. (H) Statistical chart of P duration. (I) Statistical chart of QT interval. (J) Statistical chart of QTc (n = 6). (K) RDV causes second-degree type I and second-degree type II cardiac conduction block. (L) RDV induces junctional escape rhythm.
Subsequently, we analyzed electrocardiographic parameters (Figure 1D), revealing that RDV had no significant effect on the PR interval (Figure 1E) or QRS duration (Figure 1F). However, RDV notably prolonged the RR interval (Figure 1G), P duration after the fifth day of administration (Figure 1H), QT interval (Figure 1I), and QTc (Figure 1J) after the fifth day of administration. Prolongation of the RR interval, P duration, and QT interval suggests that RDV not only affects sinus node pacing function but also influences sinus atrial conduction, with concomitant prolonged ventricular repolarization. Among these effects, RDV-induced suppression to SA node pacing function appears to be the primary factor contributing to sinus bradycardia.
Finally, we observed that RDV not only induced severe sinus bradycardia with arrhythmia but also elicited characteristic electrocardiographic responses representative of sick sinus syndrome. These included RDV-induced manifestations of marked second-degree type I and type II sinus atrioventricular (AV) block in guinea pigs (Figure 1K), further indicating severe sinus AV conduction block. Moreover, RDV also resulted in junctional escape rhythms in guinea pigs (Figure 1L), where P waves appeared after QRS complexes, underscoring RDV’s profound impact on SA node function, leading to sinus arrest and initiation of escape pacemaker activity. This pacemaker first activated the ventricles, then retrogradely activated the atria, while maintaining a bradycardic state throughout, suggesting RDV’s suppressive effect on cardiac pacing. In addition, some guinea pigs experienced profound bradycardia, culminating in sudden cardiac arrest.
In summary, these findings demonstrate that RDV disrupts SA node pacing function, leading to sinus atrioventricular conduction block, QT prolongation, resulting in severe bradycardia and arrhythmia, ultimately culminating in sick sinus syndrome.
3.2 Electrocardiographic studies in ex vivo cardiac preparations
Next, we excluded the influence of other factors such as neural and hormonal effects and investigated the impact of RDV on ex vivo guinea pig cardiac electrophysiology using Langendorff perfusion.
Administering various concentrations of RDV (0.1 μM, 0.3 μM, 1 μM, 3 μM, 10 μM) to ex vivo guinea pig hearts via Langendorff perfusion, we monitored electrocardiographic changes pre- and post-administration. Our results revealed that compared to the control, both 3 μM and 10 μM RDV significantly slowed ex vivo guinea pig heart rates, with recovery observed upon washout (Figures 2A,B). Conversely, lower concentrations of 0.1 μM, 0.3 μM, and 1 μM RDV showed no significant impact on ex vivo cardiac rhythm (Supplementary Table S1). Noticeably, there exist significant discrepancies in heart rates between ex vivo and in vivo guinea pig hearts, attributed to the inability of ex vivo experimental conditions to fully replicate the physiological state of the heart in situ, including the regulatory influences such as neural and hormonal inputs. Consequently, ex vivo cardiac conditions cannot precisely replicate the exact physiological status of the in vivo heart, thus producing inherent baseline differences between these two conditions (Wang et al., 2021).
Figure 2. Effect of RDV on electrocardiography in Langendorff-perfused guinea pig isolated hearts. (A) Representative electrocardiograms of guinea pig isolated hearts before and after administration of 3 μM RDV, 10 μM RDV and after wash. (B) A statistical graph depicting the changes in heart rate of guinea pig isolated hearts before and after the administration of 3 μM and 10 μM RDV (n = 5). (C) Statistical chart of PR interval. (D) Statistical chart of QRS. (E) Statistical chart of RR interval. (F) Statistical chart of P duration. (G) Statistical chart of QTc (n = 5).
Subsequently, we conducted further analysis of ex vivo guinea pig cardiac electrophysiological parameters. Our results indicated that different concentrations of RDV had no significant impact on PR interval (Figure 2C) or QRS duration (Figure 2D). However, 3 μM and 10 μM RDV significantly prolonged the RR interval (Figure 2E), with only 10 μM RDV causing significant prolongation of P duration (Figure 2F), and QTc (Figure 2G).
The consistent findings between ex vivo guinea pig cardiac electrophysiology and in vivo guinea pig cardiac electrophysiology indicate that RDV primarily inhibits SA node pacing function, leading to sinus atrioventricular conduction block, resulting in sick sinus syndrome, accompanied by a certain degree of QTc prolongation.
3.3 Optical mapping studies in ex vivo SA node preparations
Next, we employed optical mapping techniques to investigate the effects of RDV on acutely dissected guinea pig SA node tissue. The dissected tissue consisted of cardiac tissue from the left and right atria and atrial appendages, with the SA node located as a translucent structure between the superior and inferior vena cava in the right atrium and the crista terminalis (Figure 3A). Recording electrodes were placed on both sides of the left and right atrial appendages, and the optical mapping data acquisition area covered the entire left and right atria (Figure 3B).
Figure 3. Effect of RDV on pacemaking rate and AP conduction in SA node preparations of guinea pig. (A) Schematic representation of the SA node (SAN) preparations after dissection. RAA, Right Atrial Appendage; LAA, Left Atrial Appendage; SVC, Superior Vena Cava; IVC, Inferior Vena Cava; RV, right ventricle; LV, Left Ventricle; AVJ, Artioventricular Junction. (B) Schematic of the optical mapping area and ECG recording configuration. (C) Representative isochronal conduction maps of SA node tissue before and after administration of 3 μM and 10 μM RDV. (D) Schematic of conduction time from the SA node site to the crista terminalis (SACT), CT, Crista terminalis. (E) Statistical graph of SACT after the administration of 3 μM (n = 5) and 10 μM (n = 6) of RDV. (F) A statistical graph depicting the changes in pacemaking rate of sinus node before and after the administration of 3 μM RDV (n = 5). (G) A statistical graph depicting the changes in pacemaking rate of sinus node before and after the administration of 10 μM RDV (n = 6).
Figure 3C illustrates voltage mapping measurements of action potential initiation and conduction under sinus rhythm conditions. The results of sinoatrial conduction time (Figures 3D,E) indicated sinoatrial conduction block under the influence of 3 μM and 10 μM RDV, while SA node pacing rates (Figures 3F,G) showed a significant decrease, along with visible conduction heterogeneity and slowing (Figure 3C). Finally, optical action potential maps (Figures 4A–F) demonstrated that 3 μM and 10 μM RDV prolonged the duration of action potentials at 50% and 80% repolarization (APD50 and APD80) in SA node tissue.
Figure 4. Effect of RDV on APD50 and APD80 in SA node tissue. (A) Effect of 3 μM RDV on APD50 and APD80 in SA node tissue. (B) Effect of 10 μM RDV on APD50 and APD80 in SA node tissue. (C) APD50 statistics after administration of 3 μM RDV (n = 5). (D) APD80 statistics after administration of 3 μM RDV (n = 5). (E) APD50 statistics after administration of 10 μM RDV (n = 6). (F) APD80 statistics after administration of 10 μM RDV (n = 6).
3.4 Patch-clamp studies of spontaneous activity and HCN4 channel characteristics in SA node cells
Next, we acutely isolated guinea pig SA node cells and utilized patch-clamp techniques to investigate the effects of RDV on spontaneous action potentials and the key pacemaker ion channel HCN4 channel If current. The results revealed that 3 μM RDV significantly decreased the pacing rate of SA node cells (Figures 5A,B). Additionally, the original current traces of the HCN4 channel (Figure 5C) and current traces at −130 mV (Figure 5D) showed that 3 μM RDV inhibited the HCN4 channel, leading to a significant reduction in If current. Moreover, the inhibition rate of If current by 3 μM RDV exceeded 40% (Figures 5E,F). Therefore, our study has confirmed that RDV significantly suppresses the crucial pacemaker current If current in SA node cells, contributing to a decrease in the pacing rate of sinoatrial node cells.
Figure 5. Effect of RDV on the pacemaking rate and HCN4 channel characteristics of SA node cells. (A) Representative plot of the effect of administration of 3 μM RDV on the Pacemaking rate of SA node cells. (B) Statistical plot of the effect of 3 μM RDV on SA node cell beating rate (n = 5). (C) Representative plots of If current changes before and after administration of 3 μM RDV. (D) Representative plots of the effect of 3 μM RDV on If current at a hyperpolarizing voltage of − 130 mV. (E) I-V curves of If current before and after 3 μM RDV administration (n = 5). (F) Suppression of current at hyperpolarizing voltage (n = 5).
3.5 Patch-clamp studies of hERG channel characteristics in HEK293 cells
Finally, due to the observed QT prolongation in both clinical and our in vivo and in vitro electrocardiographic data, we further investigated the effects of RDV on hERG channel IKr current using patch-clamp techniques in HEK293 cells. The results demonstrated that both 3 μM and 10 μM RDV significantly inhibited the hERG channel current compared to the control group (Figures 6A,B). Moreover, the I-V curve analysis indicated that the inhibition of the hERG channel by RDV was voltage-dependent (Figure 6C). Additionally, as illustrated in Figure 6D, the inhibition of the hERG channel by RDV displayed concentration-dependence, with average inhibition rates of 32.7% and 34.9% for IKr current at 3 μM and 10 μM RDV, respectively. These findings suggest that RDV can lead to inhibition of the hERG channel, and the reduction in IKr current due to hERG channel inhibition can result in QT interval prolongation.
Figure 6. Effect of RDV on IKr current of hERG channels in HEK293. (A) Representative plots of IKr current changes before and after administration of 3 μM and 10 μM RDV. (B) Effect of administration of 3 μM RDV and 10 μM RDV on IKr current at 60 mV voltage. (C) The I-V curve graph of IKr current after the action of 3 μM (n = 5) and 10 μM RDV (n = 7). (D) Inhibition of hERG channels by different concentrations of RDV.
4 Discussion
This study investigated the effects of RDV on guinea pig hearts and isolated guinea pig heart electrophysiology under Langendorff perfusion to assess whether the obtained results align with clinical observations. Furthermore, we explored the impact of RDV on SA node tissue pacing, action potential signal propagation, pacing of SA node cells, and the HCN4 ion channel in acute isolated SA node cells, thereby further studying the electrophysiological mechanism of RDV-induced sick sinus syndrome at the tissue, cellular, and channel levels. Additionally, we examined the effects of RDV on hERG channel IKr current in HEK293 cells to confirm its association with clinical QT interval prolongation. This study represents the first comprehensive investigation into the cardiac electrophysiological effects of RDV, incorporating animal electrocardiography, ex vivo cardiac electrophysiology, SA node tissue pacing, action potential conduction, SA node cell pacing, and ion channel alterations. Finally, we correlate the observed cardiac electrophysiological responses induced by RDV with its potential to cause sick sinus syndrome clinically. The diagnosis of sick sinus syndrome, the latter of which is based on drug-induced bradycardia, conduction block, QT interval prolongation, and cardiac arrest, among other electrocardiographic abnormalities observed in in vivo electrocardiography, ex vivo cardiac electrophysiology, and sinoatrial node tissue conduction, involves symptoms such as sinus bradycardia, sinus arrest or pause, sinoatrial exit block, chronotropic incompetence, or bradycardia-tachycardia syndrome (Brignole et al., 2013; Kusumoto et al., 2019). The investigation into the electrophysiological mechanism of RDV-induced sick sinus syndrome stems from numerous clinical case reports of adverse cardiac events, including bradycardia, conduction blocks, QT interval prolongation, and cardiac arrest (Nabati and Parsaee, 2022; Nie et al., 2022; Wu et al., 2022; Ai et al., 2023; Devgun et al., 2023; Kanagala et al., 2023), highlighting the critical need for mechanistic research to elucidate the cardiac adverse effects of RDV in COVID-19 patients.
Our study yielded the following findings: Telemetry electrocardiography in guinea pigs revealed that RDV induced sick sinus syndrome symptoms, including sinus bradycardia, sinoatrial conduction block, junctional escape rhythm, and sinus arrest, as well as long QT syndrome. Prolongation of P duration also indicated RDV-induced sinoatrial conduction block, consistent with clinical observations. Additionally, electrocardiographic parameters showed RDV-induced QT interval prolongation, aligning with clinical findings. Moreover, ex vivo electrocardiography results similarly demonstrated that μM RDV caused bradycardia, sinoatrial conduction block, and QT interval prolongation, with heart rate recovery upon washout, consistent with our in vivo and clinical findings, thus providing further evidence of the adverse cardiac effects of RDV. Furthermore, we employed optical mapping techniques for the first time to study the effects of RDV on action potential conduction and duration at the sinoatrial node tissue level. Tissue-level results indicated that RDV inhibited sinoatrial node pacing and prolonged action potential conduction time and duration, further elucidating the reasons for RDV-induced sick sinus syndrome at the ex vivo cardiac and animal levels. Lastly, we investigated for the first time at the acute isolated sinoatrial node cell level the effects of RDV on native cell pacing and the key pacing ion channel HCN4 channel, revealing significant inhibition of sinoatrial node cell pacing by RDV and significant inhibition of the key pacing ion channel HCN4 channel’s If current. Additionally, RDV significantly inhibited hERG channel IKr current. Thus, further ion channel and cellular evidence confirmed that RDV causes a decrease in sinoatrial node cell pacing by inhibiting the key pacing current If, leading to decreased sinoatrial node tissue pacing accompanied by sinoatrial conduction tissue and action potential duration prolongation, ultimately resulting in symptoms of sick sinus syndrome such as sinus bradycardia, sinoatrial conduction block, junctional escape rhythm, and sinus arrest. Furthermore, RDV significantly inhibits hERG channel IKr current, leading to QT prolongation, presenting another cause of bradycardia.
Previous studies on RDV have primarily focused on clinical cases and a few in vitro investigations. One pediatric case with severe acute COVID-19 showed sinus bradycardia during RDV treatment (Al-Jammali et al., 2022). In a prospective longitudinal study of COVID-19 patients, the overall incidence of bradycardia (heart rate < 60 bpm) was 27% among those treated with RDV, with 19% of them experiencing severe bradycardia (heart rate < 50 bpm) (Hajimoradi et al., 2023). Furthermore, a study by the US FDA and the Medical Toxicology Academy confirmed RDV-induced adverse cardiac reactions including sinus bradycardia, QRS/QT interval prolongation, sinus arrest, and cardiac arrest (Devgun et al., 2023). These previous findings underscore the cardiac adverse effects of RDV in clinical settings. Our in vivo and ex vivo electrocardiography results further demonstrate that RDV can induce sick sinus syndrome and prolong QT intervals.
Previous studies by Wolfes et al. (2024) and Pilote et al. (2022) have also demonstrated that RDV can cause a decrease in heart rate in isolated hearts, accompanied by QT interval prolongation and shortened ventricular APD90. Note that their research focused on the ventricular region of the heart, specifically on QT interval prolongation. However, based on our in vivo electrocardiographic and ex vivo heart findings in guinea pigs, QT interval prolongation is not the primary cause of RDV-induced cardiac adverse effects. Instead, dysfunction of the sinoatrial node (SAN) is the major factor. Therefore, our study primarily investigates the conduction and APD of excitation in the SA node region, revealing that RDV leads to reduced SAN pacing, sinoatrial conduction block, and APD prolongation. Moreover, this study represents the first to utilize optical mapping techniques to investigate RDV’s effects on action potential conduction and duration at the sinoatrial node tissue level, providing valuable insights into the mechanisms of cardiac adverse reactions and enhancing drug safety assessment by offering experimental insights not achievable solely through in vitro single-cell single-channel studies (Guo et al., 2009). These results further elucidate the pathogenesis of RDV-induced sick sinus syndrome clinically. Furthermore, our investigation at the acute isolated sinoatrial node cell level and ion channel level demonstrates convincing evidence of RDV’s specific effects on sinoatrial node cells and channels.
The study by Szendrey et al. (2021) demonstrates that in stably transfected HEK293 cells, RDV not only lacks an effect on the hERG current but also does not block the KCNQ1 + KCNE1 current. However, until recently, Al-Moubarak et al. (2021) demonstrated, contrary to Szendrey et al., that RDV indeed blocks the hERG current in HEK293 transfected cells. This result is consistent with our patch clamp results, showing significant inhibition of hERG current by RDV, which may account for as a major cause the QT interval prolongation and a contributing factor to bradycardia. However, based on the in vivo telemetry electrocardiography results and ex vivo electrocardiography results, pathological sick sinus syndrome appears to be the main cause of bradycardia. Therefore, close attention should be paid to the cardiac side effects of RDV in clinical use, and appropriate ion channel modulation can be employed to prevent the occurrence of severe sick sinus syndrome.
4.1 Limitations
This study primarily investigated the effects of RDV on sinoatrial node cells at the cellular level, focusing on pacemaking and the If current of the key pacemaking channel HCN4. However, the action potential generation of sinoatrial node cells involves not only the phase 4 primary hyperpolarization-activated current If, but also currents from Nav1.5 channels, T-type and L-type Ca2+ channels, as well as the repolarization phase IKs current (Donald and Lakatta, 2023; Fan et al., 2023). Therefore, effects of RDV on other ion channel currents warrants further investigations. Additionally, unlike typical COVID-19 patients treated with RDV, this study used healthy guinea pigs, which cannot fully reflect the complexity of RDV use in critically ill COVID-19 patients. Although this could be addressed with an in vivo COVID-19 guinea pig model, it is currently challenging due to safety considerations and the absence of relevant guinea pig experimental models. Nevertheless, the findings from this study contribute to understanding the cardiac safety profile of RDV when used alone, providing foundational mechanistic insights that could guide future clinical use of RDV and aid in mitigating cardiac adverse reactions associated with its clinical application.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by the Institutional Animal Care and Use Committee of Southwest Medical University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
SL: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing. LY: Formal Analysis, Investigation, Methodology, Writing–review and editing. YX: Data curation, Formal Analysis, Investigation, Methodology, Writing–review and editing. HZ: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (No. 61803318).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2024.1436727/full#supplementary-material
References
Ahmad T., Yin P., Saffitz J., Pockros P. J., Lalezari J., Shiffman M., et al. (2015). Cardiac dysfunction associated with a nucleotide polymerase inhibitor for treatment of hepatitis C. Hepatology 62 (2), 409–416. doi:10.1002/hep.27488
Ai M. Y., Chang W. L., Yang C. J. (2023). Remdesivir-induced bradycardia and mortality in SARS-CoV-2 infection, potential risk factors assessment: a systematic review and meta-analysis. J. Clin. Med. 12 (24), 7518. doi:10.3390/jcm12247518
Al-Jammali S., Al-Zakhari R., Sheets N., Mahtani A., Stefanishina V., Isber N. (2022). Bradyarrhythmia after remdesivir administration in SARS-CoV-2: a review of literature and meta-analysis of observational studies in epidemiology. Cardiol. Res. 13 (3), 135–143. doi:10.14740/cr1377
Al-Moubarak E., Sharifi M., Hancox J. C. (2021). In silico exploration of interactions between potential COVID-19 antiviral treatments and the pore of the hERG potassium channel-A drug antitarget. Front. Cardiovasc. Med. 8, 645172. doi:10.3389/fcvm.2021.645172
Antonioli L., Fornai M., Pellegrini C., Pacher P., Haskó G. (2023). Adenosine signaling as target in cardiovascular pharmacology. Curr. Opin. Pharmacol. 71, 102393. doi:10.1016/j.coph.2023.102393
Belardinelli L., Giles W. R., West A. (1988). Ionic mechanisms of adenosine actions in pacemaker cells from rabbit heart. J. Physiology 405, 615–633. doi:10.1113/jphysiol.1988.sp017352
Bjork J. A., Wallace K. B. (2021). Remdesivir; molecular and functional measures of mitochondrial safety. Toxicol. Appl. Pharmacol. 433, 115783. doi:10.1016/j.taap.2021.115783
Brignole M., Auricchio A., Baron-Esquivias G., Bordachar P., Boriani G., Breithardt O. A., et al. (2013). 2013 ESC Guidelines on cardiac pacing and cardiac resynchronization therapy: the Task Force on cardiac pacing and resynchronization therapy of the European Society of Cardiology (ESC). Developed in collaboration with the European Heart Rhythm Association (EHRA). Eur. Heart J. 34 (29), 2281–2329. doi:10.1093/eurheartj/eht150
Devgun J. M., Zhang R., Brent J., Wax P., Burkhart K., Meyn A., et al. (2023). Identification of bradycardia following remdesivir administration through the US food and drug administration American college of medical Toxicology COVID-19 toxic pharmacovigilance Project. JAMA Netw. open 6 (2), e2255815. doi:10.1001/jamanetworkopen.2022.55815
De Wit E., Feldmann F., Cronin J., Jordan R., Okumura A., Thomas T., et al. (2020). Prophylactic and therapeutic remdesivir (GS-5734) treatment in the rhesus macaque model of MERS-CoV infection. Proc. Natl. Acad. Sci. U. S. A. 117 (12), 6771–6776. doi:10.1073/pnas.1922083117
Donald L., Lakatta E. G. (2023). What makes the sinoatrial node tick? A question not for the faint of heart. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 378 (1879), 20220180. doi:10.1098/rstb.2022.0180
Fan W., Sun X., Yang C., Wan J., Luo H., Liao B. (2023). Pacemaker activity and ion channels in the sinoatrial node cells: MicroRNAs and arrhythmia. Prog. Biophysics Mol. Biol. 177, 151–167. doi:10.1016/j.pbiomolbio.2022.11.005
Guo J., Mitsuiye T., Noma A. (1997). The sustained inward current in sino-atrial node cells of Guinea-pig heart. Pflügers Arch. 433 (4), 390–396. doi:10.1007/s004240050293
Guo L., Dong Z., Guthrie H. (2009). Validation of a Guinea pig Langendorff heart model for assessing potential cardiovascular liability of drug candidates. J. Pharmacol. Toxicol. Methods 60 (2), 130–151. doi:10.1016/j.vascn.2009.07.002
Hajimoradi M., Sharif Kashani B., Dastan F., Aghdasi S., Abedini A., Naghashzadeh F., et al. (2023). Remdesivir associated sinus bradycardia in patients with COVID-19: a prospective longitudinal study. Front. Pharmacol. 13, 1107198. doi:10.3389/fphar.2022.1107198
Hamill O. P., Marty A., Neher E., Sakmann B., Sigworth F. J. (1981). Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflugers Arch. 391 (2), 85–100. doi:10.1007/BF00656997
Haq K. T., Cooper B. L., Berk F., Roberts A., Swift L. M., Posnack N. G. (2022). Demographic and methodological heterogeneity in electrocardiogram signals from Guinea pigs. Front. Physiology 13, 925042. doi:10.3389/fphys.2022.925042
Horn R., Marty A. (1988). Muscarinic activation of ionic currents measured by a new whole-cell recording method. J. General Physiol. 92 (2), 145–159. doi:10.1085/jgp.92.2.145
Isenberg G., Klockner U. (1982). Calcium tolerant ventricular myocytes prepared by preincubation in a “KB medium”. Pflugers Arch. 395 (1), 6–18. doi:10.1007/BF00584963
Jeong Y. J., Wi Y. M., Park H., Lee J. E., Kim S. H., Lee K. S. (2023). Current and emerging knowledge in COVID-19. Radiology 306 (2), e222462. doi:10.1148/radiol.222462
Kanagala S. G., Dholiya H., Jhajj P., Patel M. A., Gupta V., Gupta S., et al. (2023). Remdesivir-induced bradycardia. South. Med. J. 116 (3), 317–320. doi:10.14423/SMJ.0000000000001519
Kojima A., Kitagawa H., Omatsu-Kanbe M., Matsuura H., Nosaka S. (2012). Inhibitory effects of sevoflurane on pacemaking activity of sinoatrial node cells in Guinea-pig heart. Br. J. Pharmacol. 166 (7), 2117–2135. doi:10.1111/j.1476-5381.2012.01914.x
Kusumoto F. M., Schoenfeld M. H., Barrett C., Edgerton J. R., Ellenbogen K. A., Gold M. R., et al. (2019). 2018 ACC/AHA/HRS guideline on the evaluation and management of patients with bradycardia and cardiac conduction delay: executive summary: a report of the American college of cardiology/American heart association task force on clinical practice guidelines, and the heart rhythm society. J. Am. Coll. Cardiol. 74 (7), 932–987. doi:10.1016/j.jacc.2018.10.043
Kwok M., Lee C., Li H. S., Deng R., Tsoi C., Ding Q., et al. (2022). Remdesivir induces persistent mitochondrial and structural damage in human induced pluripotent stem cell derived cardiomyocytes. Cardiovasc. Res. 118, 2652–2664. doi:10.1093/cvr/cvab311
Malin J. J., Suárez I., Priesner V., Fätkenheuer G., Rybniker J. (2020). Remdesivir against COVID-19 and other viral diseases. Clin. Microbiol. Rev. 34 (1), e001622–20. doi:10.1128/CMR.00162-20
Nabati M., Parsaee H. (2022). Potential cardiotoxic effects of remdesivir on cardiovascular system: a literature review. Cardiovasc. Toxicol. 22 (3), 268–272. doi:10.1007/s12012-021-09703-9
Nie Z., Sun T., Zhao F. (2022). Safety and efficacy of antiviral drugs for the treatment of COVID-19: a systematic review. Infect. Drug Resist. 15, 4457–4466. doi:10.2147/IDR.S362946
Pilote S., Simard C., Drolet B. (2022). Remdesivir (VEKLURY) for treating COVID-19: Guinea pig ex vivo and in vivo cardiac electrophysiological effects. J. Cardiovasc. Pharmacol. 80 (4), 616–622. doi:10.1097/FJC.0000000000001321
Rivera D. R., Peters S., Panagiotou O. A., Shah D. P., Kuderer N. M., Hsu C. Y., et al. (2020). Utilization of COVID-19 treatments and clinical outcomes among patients with cancer: a COVID-19 and cancer consortium (CCC19) cohort study. Cancer Discov. 10 (10), 1514–1527. doi:10.1158/2159-8290.CD-20-0941
Sheahan T. P., Sims A. C., Graham R. L., Menachery V. D., Gralinski L. E., Case J. B., et al. (2017). Broad-spectrum antiviral GS-5734 inhibits both epidemic and zoonotic coronaviruses. Sci. Transl. Med. 9 (396), eaal3653. doi:10.1126/scitranslmed.aal3653
Strauss D. G., Gintant G., Li Z., Wu W., Blinova K., Vicente J., et al. (2019). Comprehensive in vitro proarrhythmia assay (CiPA) update from a cardiac safety research consortium/health and environmental sciences institute/FDA meeting. Ther. Innovation Regul. Sci. 53 (4), 519–525. doi:10.1177/2168479018795117
Szendrey M., Guo J., Li W., Yang T., Zhang S. (2021). COVID-19 drugs chloroquine and hydroxychloroquine, but not azithromycin and remdesivir, block hERG potassium channels. J. Pharmacol. Exp. Ther. 377, 265–272. doi:10.1124/jpet.120.000484
Wang G., Lu C. J., Trafford A. W., Tian X., Flores H. M., Maj P., et al. (2021). Electrophysiological and proarrhythmic effects of hydroxychloroquine challenge in Guinea-pig hearts. ACS Pharmacol. Transl. Sci. 4 (5), 1639–1653. doi:10.1021/acsptsci.1c00166
Warren T. K., Jordan R., Lo M. K., Ray A. S., Mackman R. L., Soloveva V., et al. (2016). Therapeutic efficacy of the small molecule GS-5734 against Ebola virus in rhesus monkeys. Nature 531 (7594), 381–385. doi:10.1038/nature17180
Williamson B. N., Feldmann F., Schwarz B., Meade-White K., Porter D. P., Schulz J., et al. (2020). Clinical benefit of remdesivir in rhesus macaques infected with SARS-CoV-2. Nature 585 (7824), 273–276. doi:10.1038/s41586-020-2423-5
Wolfes J., Kirchner L., Doldi F., Wegner F., Rath B., Eckardt L., et al. (2024). Electrophysiological profile of different antiviral therapies in a rabbit whole-heart model. Cardiovasc. Toxicol. 24, 656–666. doi:10.1007/s12012-024-09872-3
Keywords: Remdesivir, cardiac electrophysiology, sick sinus syndrome, patch clamp, COVID-19
Citation: Li S, Yue L, Xie Y and Zhang H (2024) Electrophysiological and sick sinus syndrome effects of Remdesivir challenge in guinea-pig hearts. Front. Physiol. 15:1436727. doi: 10.3389/fphys.2024.1436727
Received: 22 May 2024; Accepted: 04 July 2024;
Published: 13 August 2024.
Edited by:
Drew Nassal, The Ohio State University, United StatesReviewed by:
Thomas Hund, The Ohio State University, United StatesSamantha Salvage, University of Cambridge, United Kingdom
Copyright © 2024 Li, Yue, Xie and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Henggui Zhang, henggui.zhang@manchester.ac.uk