- 1Section on Fibrotic Disorders, National Institute on Alcohol Abuse and Alcoholism, National Institutes of Health, Rockville, MD, United States
- 2Laboratory of Physiologic Studies, National Institute on Alcohol Abuse and Alcoholism, National Institutes of Health, Rockville, MD, United States
- 3Section on Medicinal Chemistry, National Institute on Alcohol Abuse and Alcoholism, National Institutes of Health, Rockville, MD, United States
Scleroderma, or systemic sclerosis, is a multi-organ connective tissue disease resulting in fibrosis of the skin, heart, and lungs with no effective treatment. Endocannabinoids acting via cannabinoid-1 receptors (CB1R) and increased activity of inducible NO synthase (iNOS) promote tissue fibrosis including skin fibrosis, and joint targeting of these pathways may improve therapeutic efficacy. Recently, we showed that in mouse models of liver, lung and kidney fibrosis, treatment with a peripherally restricted hybrid CB1R/iNOS inhibitor (MRI-1867) yields greater anti-fibrotic efficacy than inhibiting either target alone. Here, we evaluated the therapeutic efficacy of MRI-1867 in bleomycin-induced skin fibrosis. Skin fibrosis was induced in C57BL/6J (B6) and Mdr1a/b-Bcrp triple knock-out (KO) mice by daily subcutaneous injections of bleomycin (2 IU/100 µL) for 28 days. Starting on day 15, mice were treated for 2 weeks with daily oral gavage of vehicle or MRI-1867. Skin levels of MRI-1867 and endocannabinoids were measured by mass spectrometry to assess target exposure and engagement by MRI-1867. Fibrosis was characterized histologically by dermal thickening and biochemically by hydroxyproline content. We also evaluated the potential increase of drug-efflux associated ABC transporters by bleomycin in skin fibrosis, which could affect target exposure to test compounds, as reported in bleomycin-induced lung fibrosis. Bleomycin-induced skin fibrosis was comparable in B6 and Mdr1a/b-Bcrp KO mice. However, the skin level of MRI-1867, an MDR1 substrate, was dramatically lower in B6 mice (0.023 µM) than in Mdr1a/b-Bcrp KO mice (8.8 µM) due to a bleomycin-induced increase in efflux activity of MDR1 in fibrotic skin. Furthermore, the endocannabinoids anandamide and 2-arachidonylglycerol were elevated 2-4-fold in the fibrotic vs. control skin in both mouse strains. MRI-1867 treatment attenuated bleomycin-induced established skin fibrosis and the associated increase in endocannabinoids in Mdr1a/b-Bcrp KO mice but not in B6 mice. We conclude that combined inhibition of CB1R and iNOS is an effective anti-fibrotic strategy for scleroderma. As bleomycin induces an artifact in testing antifibrotic drug candidates that are substrates of drug-efflux transporters, using Mdr1a/b-Bcrp KO mice for preclinical testing of such compounds avoids this pitfall.
Introduction
Scleroderma, or systemic sclerosis (SSc), is a connective tissue disease with multiple clinical manifestations, including autoimmunity, vascular dysfunction, and tissue fibrosis (1), and a prevalence in the United States of around 240 cases per 1 million adults (2). Scleroderma is a complex, heterogeneous disease with clinical forms ranging from limited skin involvement (limited cutaneous systemic sclerosis) to forms with diffuse skin sclerosis and severe and often progressive internal organ involvement (diffuse cutaneous systemic sclerosis) (3). Pulmonary fibrosis and interstitial lung diseases (ILD) occur in about 60% of patients, contributing to mortality (4), while dermal fibrosis causes significant morbidity in scleroderma (5, 6). In the absence of approved therapies, there is an unmet need for identifying new targets and treatment strategies. Due to the complex and multifactorial pathogenesis of scleroderma, targeting multiple signaling pathways may be essential for effective treatment (7).
Inducible nitric oxide synthase (iNOS) is an enzyme encoded by the Nos2 gene and is responsible for generating pro-inflammatory reactive nitrogen species (8). The relevance of iNOS as a target is based on evidence for overproduction of NO in the pathogenesis of SSc (9, 10). In SSc, the expression of iNOS in the endothelium, smooth muscle cells, fibroblasts, macrophages and many other cell types is robustly induced by inflammatory mediators and cytokines and its activity is increased at inflammatory sites (8). The iNOS-mediated formation of NO is increased in inflammatory cells such as macrophages or activated fibroblasts (10). Immunohistological studies of scleroderma skin show that disease progression involves iNOS upregulation (11). Previous studies also demonstrate that SSc lung macrophages express high levels of iNOS and produce a high quantity of ONOO- anions (11). In SSc patients, increased production of NO is suggested by the increased expression of iNOS in endothelial cells, fibroblasts and mononuclear cells infiltrating the fibrotic skin (12) as well as in alveolar macrophages (13). The role of NO synthases and especially iNOS is elegantly dissected by the work of Cotton et al., which proposes an active role of iNOS-induced NO production in endothelial cell damage and advances the concept of iNOS inhibition as a viable therapeutic strategy for SSc (14).
An additional target that is becoming increasingly relevant in the modulation of fibrotic responses is the endocannabinoid system. Endocannabinoids are lipid-signaling molecules that act through cannabinoid receptors CB1 and CB2. Endocannabinoids acting via CB1R promote fibrosis in multiple organs including skin (15), liver (16–18), kidney (19), and heart (20). Besides, CB1R activation is pro-inflammatory (21). Increased CB1R activity has been linked to different forms of pulmonary fibrosis such as radiation-induced pulmonary fibrosis (22), idiopathic pulmonary fibrosis (23) and Hermansky-Pudlak syndrome pulmonary fibrosis (HPSPF) (24). Conversely, CB1R antagonism prevents fibroblast activation and exerts potent antifibrotic effects in skin fibrosis (25). The role of CB1R as a pro-fibrotic receptor has also been confirmed in fatty acid amide hydrolase knock-out mice, in which elevated levels of anandamide induced skin fibrosis in a CB1R-dependent manner (26).
Bleomycin is widely used to induce fibrosis in rodent models of fibrotic disorders. We have earlier reported that bleomycin induces drug-efflux associated ABC transporters in the lung, which limits exposure of the fibrotic tissue to drug candidates that are substrates of such transporters. Here we show that this pitfall could be avoided by establishing bleomycin-induced skin fibrosis in -Bcrp(-/-) triple knock-out mice and using this model to reveal the antifibrotic therapeutic efficacy of the peripherally restricted hybrid CB1R/iNOS antagonist MRI-1867, a known substrate of drug efflux transporters.
Materials and Methods
Chemicals
S-MRI-1867, referred to as MRI-1867, was synthesized as described previously (16). Rimonabant was obtained from the National Institute of Drug Abuse Drug Supply Program (Research Triangle Park, NC, USA). Pharmaceutical grade bleomycin was obtained from Hospira. All other chemicals were from Sigma-Aldrich.
Animals
All animal procedures were conducted in accordance with the rules and regulations of the Institutional Animal Care and Use Committee of the National Institutes of Alcohol Abuse and Alcoholism (NIAAA). C57BL/6J mice were purchased from the Jackson Laboratory (Bar Harbor, ME, USA). -Bcrp(-/-) mice were purchased from Taconic (Rensselaer, NY, USA). Animals were housed individually under a 12-hour light/dark cycle and fed a standard diet, ad libitum (Teklad NIH-31; Envigo, Huntingdon, UK).
Bleomycin-Induced Skin Fibrosis
This study used male C57BL/6J and -Bcrp(-/-) (KO) mice, ranging from 16-20 weeks of age, with an average initial body weight of 28.5 g. Mice received daily 100µL subcutaneous injections dorsal to each scapula for 28 days of either vehicle (sterile 0.9% saline), or 2 IU of bleomycin (bleo) to induce skin fibrosis (Figure 1A).
Figure 1 Bleomycin-induced skin fibrosis minimizes skin exposure and target engagement by MRI-1867 in C57BL6/J mice. (A) Schematic presentation of bleomycin-induced skin fibrosis development and therapeutic treatment regime. (B) Dermal thickness quantification for histological images. (C) Representative histological images from skin biopsies stained with H&E Ep, Epidermis; D, Dermis; Hy, Hypodermis and adipose tissue; Pc, Panninculus carnosus. Double arrow indicates dermal length that used for quantification. (D) Levels of MRI-1867 in control and bleomycin-treated skin biopsy specimens following acute or 14-day chronic administration of MRI-1867 at 10 mg/kg, respectively. (E) Levels of endocannabinoids AEA and 2AG in skin biopsies. Data represent means ± SEM from 3-10 mice per group. Data were analyzed by 1-way ANOVA followed by Dunnett’s multiple comparisons test. Significant difference from corresponding values in controls (saline and vehicle treated group) (*P < 0.05).
Therapeutic Intervention
10 mg/kg MRI-1867 or its drug-free vehicle (5% DMSO, 5% Tween 80 in 0.9% saline) were administered by oral gavage either as a single dose in control mice or starting on day 15 and ending on day 28 of bleomycin treatment in bleomycin-treated mice. Mice were sacrificed 1 h following the single dose (control) or the last daily dose of MRI-1867 (bleo-treated mice).
Dermal Thickness Histological Assessment
Skin tissues were fixed in 10% neutralized formalin solution for 24 hours, embedded in paraffin, and sectioned (4μm) onto glass slides. Sections were counterstained with hematoxylin and eosin and imaged with an Olympus BX41 microscope. Images were taken at 4x magnification to cover 6 mm diameter skin biopsy section from each skin section. Skin thickness was quantified by measuring the length of a straight line between the limits of dermal tissue using ImageJ software. Data points reflect the means of 7-13 independent measurements from the same mouse.
Hydroxyproline Measurements by LC-MS/MS
Skin fibrosis was quantified by measuring hydroxyproline (Hyp) content of skin biopsies consisting of epidermis and dermis, using LC-MS/MS as described previously (23) with slight modifications. Briefly, 15-30 mg skin tissue was homogenized in 600 μL of ice-cold 0.1 N perchloric acid (PCA) then 200 µl of homogenate was aliquoted and prepared for endocannabinoid and MRI-1867 measurements as detailed below. One mL 12 N HCl was added to the remaining 400 µL skin homogenate and the homogenate was hydrolyzed at 100°C for 4 hours. Hydrolyzed samples were vortexed and centrifuged at 10,000 g for 10 minutes, and 5 μL hydrolysate was diluted 200-fold by the addition of 990 μL of 0.1 N PCA and 5 μL of L-Proline-13C5,15N as internal standard. Liquid chromatography tandem mass spectrometry (LC-MS/MS) analyses were conducted on an Agilent 6410 triple quadrupole mass spectrometer (Agilent Technologies) coupled to an Agilent 1200 LC system. 4-Hydroxyproline was separated using an Intrada Amino Acid column, 50 × 3 mm, 3 μm (Imtakt) at 40°C. Mobile phases consisted of acetonitrile/tetrahydrofuran/25 mM ammonium formate/formic acid = 9:75:16:0.3 (v/v/v/v) (phase A) and acetonitrile/100 mM ammonium formate = 20:80 (v/v) (phase B). Gradient elution (600 μL/min) was initiated and held at 0% B for 3 minutes, followed by a linear increase to 17% B by 6.5 minutes. This was followed by a step increase to 100% B, which was held until 10 minutes after the gradient had begun, and then by a linear decrease to 0% B by 11 minutes, which was held until 13 minutes after the gradient had begun. The mass spectrometer was set for electrospray ionization operated in positive ion mode. The source parameters were as follows: capillary voltage, 4,000 V; gas temperature, 330°C; and drying gas, 8 L/min. Nitrogen was used as the nebulizing gas. Collision-induced dissociation (CID) was conducted using nitrogen. Hydroxyproline level was analyzed by multiple reaction monitoring. L-Proline-13C5,15N (Sigma, cat#608114) was used as the internal standard. The molecular ion and fragments for hydroxyproline were measured as follows: m/z 132.1→86 and 132.1→68 (CID energy: 8 V and 20 V, respectively). Skin levels of hydroxyproline were determined against a standard curve, using trans-4-hydroxy-L-proline as standard (Sigma-Aldrich). Values are expressed as nmol/mg wet tissue.
Endocannabinoid Extraction and Analysis
Skin homogenate (200 µL) described in the hydroxyproline measurement section was used and transferred in 0.5 mL of ice-cold methanol/Tris buffer (50 mM, pH 8.0), 1:1, containing 7 ng of [2H4] arachidonoyl ethanolamide ([2H4] AEA) as internal standard. Homogenates were extracted three times with CHCl3: MeOH (2:1, vol/vol), dried under nitrogen flow, and reconstituted with MeOH after precipitating proteins with ice-cold acetone. LC-MS/MS analyses were conducted on an Agilent 6410 triple quadrupole mass spectrometer (Agilent Technologies) coupled to an Agilent 1200 LC system. Analytes were separated using a Zorbax SB-C18 rapid-resolution HT column. Gradient elution mobile phases consisted of 0.1% formic acid in H2O (phase A) and 0.1% formic acid in MeOH (phase B). Gradient elution (250 μL/min) was initiated and held at 10% B for 0.5 min, followed by a linear increase to 85% B at 1 min and maintained until 12.5 min, then increased linearly to 100% B at 13 min and maintained until 14.5 min. The mass spectrometer was set for electrospray ionization operated in positive ion mode. The source parameters were as follows: capillary voltage, 4,000 V; gas temperature, 350°C; drying gas, 10 L/min; nitrogen was used as the nebulizing gas. Collision-induced dissociation was performed using nitrogen. Levels of each compound were analyzed by multiple reaction monitoring. The molecular ion and fragment for each compound were measured as follows: m/z 348.3→62.1 for AEA, m/z 352.3→66.1 for [2H4] AEA, and m/z 379.3→91.1 for 2-arachidonoylglycerol (2-AG). Analytes were quantified using MassHunter Workstation LC/QQQ Acquisition and MassHunter Workstation Quantitative Analysis software (Agilent Technologies). Levels of AEA and 2-AG in the samples were measured against standard curves. Values are expressed as fmol/mg wet tissue.
Tissue Levels of MRI-1867
Two hundred μL of skin homogenate described above and 30 µl serum were used and extracted as described above for endocannabinoid extraction. MRI-1867 levels were determined by LC-MS/MS using an Agilent 6410 triple quadrupole mass spectrometer (Agilent Technologies) coupled to an Agilent 1200 LC system (Agilent Technologies). Levels of MRI-1867 were analyzed by multiple reaction monitoring. The molecular ion and fragments were measured as follows: m/z 548.1→145 and 548.1→257.1 for MRI-1867 (CID energy: 56 V and 24 V, respectively). The amounts of MRI-1867 in the samples were determined against a standard curve. Values are expressed as µM.
Immunohistochemistry
Immunohistochemistry was performed as previously described (27). Tissue sections were incubated overnight at 4°C in an optimized blocking solution: 5% blotting-grade milk (Bio-Rad, Hercules, CA) and 5% horse serum (Vector Laboratories, Burlingame, CA, USA) diluted in double distilled H2O. Rabbit monoclonal anti-P-Glycoprotein antibody (Abcam, Cat. ab170904) was diluted (1:1000 for all mouse tissues) in 5% milk and 5% horse serum as described above for overnight incubation at 4°C. The appropriate secondary antibody, anti-rabbit IgG made in horse (Cat. MP-7401, Vector Laboratories) was not diluted, and incubated at room temperature for 1 hour. Positive immunoreactivity was revealed via chromogenic detection with ImmPACT DAB Peroxidase (HRP) Substrate (SK-4105, Vector Laboratories) with incubating 3 minutes, then counterstained with hematoxylin (Gills Formula, Vector Laboratories) for five seconds and covered with a coverslip. Images were captured with a BX41 bright-field microscope (Olympus, Center Valley, PA, USA). Immunostaining intensity was quantified by using ImageJ software (NIH Public Domain) by a person blind to sample ID. Images were taken at 20x magnification from at least 5 randomly selected areas per skin specimen.
Statistical Analysis
Statistical analysis was performed using GraphPad Prism 8 (GraphPad Software Inc.). Normality test was performed by Prism 8 to determine whether samples show normal distribution. Then, one-way ANOVA followed by Dunnett’s multiple comparisons test was performed. P < 0.05 was considered significant. In multiple comparison post-hoc test, control groups were designated to address each statistical question as indicated in figure legends for statistical significance.
Results
Bleomycin Significantly Attenuates Skin Exposure of MRI-1867 in C57BL/6J Mice due to Increased Expression of P-Glycoprotein
One of the hallmarks of systemic sclerosis is skin thickening due to fibrosis. Bleomycin is commonly used as an exogenous inducer of fibrosis in murine models of skin and pulmonary fibrosis. In this study, skin fibrosis was induced by daily subcutaneous injection of bleomycin for 28 days as detailed in the methods and as depicted in Figure 1A. We recently demonstrated that bleomycin induces drug efflux transporters in murine lungs and that this mechanism compromises lung exposure to chemical compounds that happen to be substrates (27). To determine whether a similar mechanism is triggered by bleomycin in the skin, we first assessed the skin levels of the peripherally restricted hybrid CB1R/iNOS inhibitor MRI-1867, a substrate of P-gp (16), in a bleomycin-induced murine model of skin fibrosis.
Ten mg/kg is the maximally effective dose of MRI-1867 for peripheral CB1R antagonism (16), and this dose was previously shown to achieve dual-target inhibition of CB1R and iNOS in lung and kidney fibrosis (23, 28). MRI-1867 was administered by oral gavage either as a single dose in control mice or daily for 14 days in bleomycin-treated mice (Figure 1A). Bleomycin (2U/day for 28 days) induced skin fibrosis in C57BL/6J mice, as quantified by measuring dermal thickness (Figures 1B, C). Notably, the levels of MRI-1867 in fibrotic skin tissue were significantly lower (0.12 µM) than in healthy control skin (15 µM after single dosing) (Figure 1D). Skin levels of anandamide (AEA) and 2-arachidonoyl glycerol (2AG) were higher in fibrotic compared to normal skin (Figure 1E), suggesting an upregulated endocannabinoid system in the fibrotic tissue. However, chronic MRI-1867 administration did not reduce bleomycin-induced increase in endocannabinoids (Figure 1E), suggesting a lack of CB1R engagement by MRI-1867 in the fibrotic skin of C57BL/6J mice (Figure 1E), due to compromised skin exposure (Figure 1D). The dramatic loss in skin exposure to MRI-1867 might be attributed to bleomycin-induced over-expression and over-activity of P-gp, which was observed previously in bleomycin-induced pulmonary fibrosis (27). Indeed, P-gp protein expression was increased by bleomycin in skin biopsy specimens from C57BL/6J mice (Figures 2A, B).
Figure 2 Bleomycin induces drug efflux transporter P-gp in skin of C57BL6/J mice. (A) Representative histological images for P-gp immunohistochemistry staining in skin biopsies of control and bleomycin-induced fibrotic skin. A, adipocytes; BV, blood vessels; D, Dermis; PC, Panninculus carnosus. (B) Quantification of P-gp protein expression in skin biopsies. Data represent mean ± SEM from five subjects in each group. Data were analyzed by t-test for comparison of histological scoring. * (p <0.05) indicates significant difference from the control group.
Skin Exposure to MRI-1867 Was Recovered in Bleomycin-Induced Skin Fibrosis Using -Bcrp(-/-) Triple Knock-Out Mice
As bleomycin causes a ~100-fold reduction in skin exposure to MRI-1867 in C57BL/6J mice, this model is unsuitable for the preclinical testing of the antifibrotic potential of MRI-1867, a known substrate of drug efflux transporters. Instead, we decided to use Mdr1a/b-Bcrp triple knockout mice for this purpose as a way to bypass the artifact caused by increased activity of drug efflux transporters. First, we measured the levels of MRI-1867 in the fibrotic skin of bleomycin-exposed -Bcrp(-/-) triple KO mice after 14 days of chronic MRI-1867 treatment at 1, 3, 10, 30, 60 mg/kg doses (Figure 3). Levels of MRI-1867 dose-dependently increased in serum (Figure 3A). Importantly, skin exposure to MRI-1867 was much higher in -Bcrp(-/-) triple KO mice compared to wild-type mice (Figure 3B), such that the 10 mg/kg dose of MRI-1867 achieved a concentration of 8.8 µM in the fibrotic skin (Figure 3B) compared to 0.12 µM in wild-type C57BL6/J mice (Figure 1D). Maximum skin exposure in the triple KO mice was ~26 µM following chronic treatment with the 30 mg/kg dose (Figure 3B). However, skin exposure was not further increased with 60 mg/kg/day dosing, which might be due to an altered ADME/PK profile with the higher dose of MRI-1867. Therefore, the 30 mg/kg/day dose was selected to explore the maximum achievable efficacy of MRI-1867 in this skin fibrosis model and to establish the PK/PD relationship.
Figure 3 Dose-dependent systemic and skin exposure of MRI-1867 in -Bcrp(-/-) mice. (A) Serum and (B) Skin levels of MRI-1867 at 1 hour after the last dose of 14 days oral administration at 1, 3, 10, 30, and 60 mg/kg doses in bleomycin-treated -Bcrp(-/-) mice. Data represent means ± SEM from 4-10 mice per group.
MRI-1867 Significantly Attenuated Dermal Thickness and Skin Fibrosis in Bleomycin-Induced Skin Fibrosis in Bcrp-/- Knock-Out Mice
Daily subcutaneous bleomycin injections for 28 days significantly increased the levels of hydroxyproline (Figure 4A) and dermal thickness (Figures 4B, C), and endocannabinoids in the fibrotic skin (Figure 4D). Chronic daily oral administration of MRI-1867 for the last 14 days of the 28 day bleo treatment significantly attenuated bleomycin-induced hydroxyproline (Figure 4A), dermal thickness (Figures 4B, C), and endocannabinoid levels (Figure 4D) in the fibrotic skin of -Bcrp(-/-) mice, suggesting that endocannabinoid tone is reduced following treatment with MRI-1867 and that targeting CB1R is a putative target for fibrosis alleviation.
Figure 4 MRI-1867 (30 mg/kg) attenuates bleomycin-induced fibrosis and dermal thickness in -Bcrp(-/-) mice. (A) Hydroxyproline levels in biopsied skin. (B) Dermal thickness analyses from skin biopsies. (C) Representative histological images from skin biopsies stained with H&E. (D) AEA and 2AG levels in biopsied skin. Six mm skin biopsies are used for the analysis. Dermal thickness assessed from left scapular injection site biopsies, Hydroxyproline and endocannabinoid measurements performed using right scapular injection site biopsies from each mouse. Data represent means ± SEM from 3-13 mice per group. Data were analyzed by 1-way ANOVA followed by Dunnett’s multiple comparisons test. Significant difference from corresponding values in control (saline and vehicle treated group) (*P < 0.05) or from values in bleomycin and vehicle-treated group (#P < 0.05).
MRI-1867 Has Higher Antifibrotic Efficacy Than Rimonabant in Skin Fibrosis
We next compared the therapeutic efficacy of MRI-1867 and rimonabant at a 10 mg/kg dose, which was shown to provide equipotent CB1R antagonism (16, 23). In addition, we also tested MRI-1867 at 1 and 3 mg/kg doses to determine the minimum effective dose that provides anti-fibrotic efficacy in bleomycin-induced skin fibrosis. Ten mg/kg MRI-1867 significantly reduced bleomycin-induced dermal thickness (Figure 5A) and hydroxyproline content (Figure 5B). Furthermore, 10 mg/kg MRI-1867 significantly reduced dermal thickness compared to the rimonabant and vehicle (Figure 5A). However, rimonabant did not significantly reduce dermal thickness compared to the vehicle. At the dose of 1 mg/kg, MRI-1867 had no significant antifibrotic effect. While the 3 mg/kg dose significantly attenuated dermal thickness (Figure 5A), it did not significantly reduce hydroxyproline (Figure 5B). MRI-1867 dose-dependently attenuated bleomycin-induced increases in skin endocannabinoids. Both rimonabant and MRI-1867 at the 10 mg/kg dose significantly and comparably attenuated skin endocannabinoids, suggesting similar target engagement in fibrotic skin (Figure 5C).
Figure 5 Comparison of the anti-fibrotic effects of MRI-1867 and rimonabant in -Bcrp(-/-) mice. (A) Dermal thickness analysis, (B) Hydroxyproline levels, (C) AEA and 2AG levels from skin biopsies in control and bleomycin-injected Mdr1a/b(-/-) -Bcrp(-/-) mice. Rimonabant (10 mg/kg) or MRI-1867 (1, 3, 10 mg/kg) were orally administered daily for the last 14 days of bleomycin injections. Data represent means ± SEM from 4-10 mice per group. Data were analyzed by 1-way ANOVA followed by Dunnett’s multiple comparisons test. Significant difference from corresponding values in control (saline and vehicle treated group) (*P < 0.05), or from values in bleomycin and vehicle-treated group (#P < 0.05), or from values in bleomycin and rimonabant-treated group (+P < 0.05).
Discussion
We have discovered that using bleomycin to model skin fibrosis in C57BL/6J mice introduces an artifact related to upregulation of drug efflux transporters in skin tissue. This would confound studies aimed to test the pharmacokinetics and target engagement of drug candidates that happen to be substrates. This finding aligns with our previous work and the work of others in preclinical models of bleomycin-induced pulmonary fibrosis (27, 29), highlighting a limitation of bleomycin-induced fibrosis models. We also show that using -Bcrp(-/-) mice for preclinical testing of such compounds would avoid this pitfall. In lung specimens from IPF patients, efflux transporters such as P-gp and BCRP were expressed at the same level as in lung samples from appropriate controls (23), which indicates that the increased P-gp expression seen in the preclinical model is not part of the pathological process in human IPF. A recent study compared the frequencies of 3 single nucleotide polymorphisms (SNPs) in the ABCB1 gene, which encodes P-gp, and found no differences between patients with systemic sclerosis and their controls in a Polish population. Although a specific haplotype of these SNPs occurred significantly more frequently among patients than among their controls, there was no evidence presented for an association of this haplotype with altered gene or protein expression or transporter activity of ABCB1 (30).
Additionally, the present findings demonstrate that the dual-target inhibition of CB1R and iNOS by MRI-1867 is an effective anti-fibrotic strategy for scleroderma that warrants further study. This finding is in line with our previous studies showing that MRI-1867 can attenuate fibrosis in other organs as well, including the liver (16), kidney (28), and lungs (23, 24). This is consistent with mounting evidence that CB1R may be part of a core mechanism of fibrogenesis and that CB1R antagonism may have therapeutic potential in several fibrotic disorders, including chronic kidney (31, 32) and liver diseases (33) and cardiomyopathies (34, 35). Although MRI-1867 was more efficacious than rimonabant at equipotent doses for CB1R antagonism, we have not investigated the relative contribution of CB1R and iNOS inhibition, which may be subject to future studies. Recently, a structural analogue of MRI-1867 was identified as a β-arrestin-2 biased CB1R antagonist, whereas rimonabant was unbiased (36). Although we have not explored a potential signaling bias of CB1R activation in skin fibrosis development, the superior efficacy of MRI-1867 over rimonabant could not be attributed to functional selectivity since MRI-1867 does not display signaling bias in CB1R antagonism (unpublished information).
In addition to fibrosis, numerous studies have documented that an overactive endocannabinoid/CB1R system contributes to visceral obesity and its complications (37), including type-2 diabetes (21), and also play a role in the pathology of alcoholic liver disease (38) and viral hepatitis (39). Conversely, CB1R blockade has beneficial effects in preclinical models of these conditions as well as in overweight individuals with metabolic syndrome (40). However, brain-penetrant CB1R antagonists, such as rimonabant, cause psychiatric side effects due to the blockade of CB1R in the CNS, which had halted their therapeutic development. Non-brain-penetrant CB1R antagonists have recently been reported to retain the metabolic benefit of globally acting compounds without blocking CB1R in the CNS (21, 41–43). Thus, efforts to engage CB1Rs for mitigating fibrosis would require antagonists with limited brain exposure in order to avoid neuropsychiatric side effects, therefore peripheral dual-target CB1R antagonists might be an emerging therapeutic modality for metabolic and fibrotic disorders (44).
Previously it was shown that deletion of CB1R protected mice from bleomycin-induced skin fibrosis (25), which aligns with the current findings and supports the therapeutic potential of peripheral CB1R antagonism in skin fibrosis. Deletion of CB1R decreased the number of infiltrating T cells and macrophages in lesioned skin (25), suggesting critical roles of CB1R in leukocyte infiltration, inflammation, and fibroblast activation. Additionally, bone marrow transplantation from CB1R-/- mouse into CB1R+/+ mouse protected the CB1R+/+ mice from bleomycin-induced skin fibrosis development, which implicated CB1R expressing myeloid cell populations in inflammation-driven skin fibrosis development (25). Additionally, CB1R signaling in keratinocytes also regulates T-cell dependent inflammation in skin (45). It is important to note that CB1R is expressed in multiple cell types in skin, and its role in skin pathologies and inflammation can be context dependent (46). This suggests that activation status and potential paracrine regulation of endocannabinoid/CB1R system in different cells in the local pathologic microenvironment might be critical factors for the pro-inflammatory and profibrotic activity of CB1Rs in skin fibrosis. Indeed, this notion is supported by our finding that a significant loss of MRI-1867 exposure in the lesioned skin resulted in loss of its anti-fibrotic efficacy despite the high level of systemic exposure. Our finding demonstrated that CB1R inhibition is required at the site of action to result in anti-fibrotic efficacy. This could suggest that topical application of CB1R antagonists might be a therapeutic strategy in skin fibrosis. However, systemic administration should be the preferred therapeutic approach considering its potential therapeutic benefit in multi-organ involvement in systemic sclerosis (44).
A pro-inflammatory role of CB1R resulting in macrophage activation was established in the pancreas during diabetes and in the lung during pulmonary fibrosis (21, 23). Furthermore, interferon regulatory factor 5 (IRF5) was found to be an essential down-stream mediator of CB1R signaling in macrophages in diabetes (47) and transplantation of CB1R-/- bone marrow to pre-diabetic ZDF rats prevented β-cell loss and diabetic complications, supporting the pathogenic role CB1R-mediated IRF5 signaling. IRF5 is a master regulator of pro-inflammatory macrophages. Furthermore, IRF5 polymorphism increases the risk of systemic sclerosis whereas reduced expression of IRF5 increases survival (48, 49). Indeed, deletion of IRF5 protected mice from development of bleomycin-induced skin and pulmonary fibrosis (50), which makes IRF5 a potential therapeutic target in systemic sclerosis and scleroderma. Deletion of CB1R also attenuated bleomycin-induced increase in IRF5 in lungs and protected from pulmonary fibrosis (23). Thus, the intriguing possibility that CB1R-mediated IRF5 signaling may contribute to skin and pulmonary fibrosis development in systemic sclerosis and scleroderma, needs to be explored in future studies.
Interestingly, targeting CB1R may also be promising for symptom management. One of the most common symptoms of systemic sclerosis patients that affects quality of life is gastrointestinal dysmotility, which results in constipation (51). Previously, we found that MRI-1867 increases upper gastrointestinal motility in mice via peripheral CB1R inverse agonism (16, 52), which might compensate for constipation. In summary, the present findings introduce a polypharmacology approach to the treatment of skin fibrosis whereby simultaneous engagement of two therapeutic targets by a single molecule is harnessed for improved therapeutic efficacy. Clinical studies in scleroderma patients are warranted once MRI-1867 or related compounds become available for human studies.
Author Contributions
CZ and JP performed in vivo experiments, histology, immunohistochemistry, image acquisition and analysis, and participated in manuscript preparation. JA contributed immunohistochemistry. MRI synthesized and chemically analyzed MRI-1867 and contributed manuscript preparation and study concept. RC designed the study, planned experiments, performed mass spectrometry experiments, analyzed and interpreted data, and drafted the manuscript. GK reviewed data and finalized the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Intramural Research Programs of the NIAAA and by a Collaborative Research and Development Agreement (CRADA) with SCOPUS BIOPHARMA INC., NY, USA.
Conflict of Interest
This study received funding from SCOPUS BIOPHARMA INC., NY, USA. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication. RC, GK, and MRI are listed as coinventors on a US patent covering MRI-1867 and related compounds (patent no. US 9,765,031 B2).
The remaining 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.
Acknowledgments
The authors thank Ms. Judith Harvey-White for her technical assistance with liquid chromatography–tandem mass spectrometry experiments.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.744857/full#supplementary-material
References
1. Katsumoto TR, Whitfield ML, Connolly MK. The Pathogenesis of Systemic Sclerosis. Annu Rev Pathol (2011) 6:509–37. doi: 10.1146/annurev-pathol-011110-130312
2. Mayes MD. Scleroderma Epidemiology. Rheum Dis Clin North Am (2003) 29:239–54. doi: 10.1016/S0889-857X(03)00022-X
3. Moinzadeh P, Aberer E, Ahmadi-Simab K, Blank N, Distler JH, Fierlbeck G, et al. All Participating, Disease Progression in Systemic Sclerosis-Overlap Syndrome Is Significantly Different From Limited and Diffuse Cutaneous Systemic Sclerosis. Ann Rheum Dis (2015) 74:730–7. doi: 10.1136/annrheumdis-2013-204487
4. Steen VD, Medsger TA. Changes in Causes of Death in Systemic Sclerosis, 1972-2002. Ann Rheum Dis (2007) 66:940–4. doi: 10.1136/ard.2006.066068
5. de-Sa-Earp AP, do Nascimento AP, Carneiro SC, Porto LC, Monte-Alto-Costa A. Dermal Dendritic Cell Population and Blood Vessels Are Diminished in the Skin of Systemic Sclerosis Patients: Relationship With Fibrosis Degree and Disease Duration. Am J Dermatopathol (2013) 35:438–44. doi: 10.1097/DAD.0b013e3182712d1a
6. Akter T, Silver RM, Bogatkevich GS. Recent Advances in Understanding the Pathogenesis of Scleroderma-Interstitial Lung Disease. Curr Rheumatol Rep (2014) 16:411. doi: 10.1007/s11926-014-0411-1
7. Denton CP. Systemic Sclerosis: From Pathogenesis to Targeted Therapy. Clin Exp Rheumatol (2015) 33:S3–7.
8. Aktan F. iNOS-Mediated Nitric Oxide Production and Its Regulation. Life Sci (2004) 75:639–53. doi: 10.1016/j.lfs.2003.10.042
9. Dooley A, Bruckdorfer KR, Abraham DJ. Modulation of Fibrosis in Systemic Sclerosis by Nitric Oxide and Antioxidants. Cardiol Res Pract (2012) 2012:521958. doi: 10.1155/2012/521958
10. Sud A, Khullar M, Wanchu A, Bambery P. Increased Nitric Oxide Production in Patients With Systemic Sclerosis. Nitric Oxide (2000) 4:615–9. doi: 10.1006/niox.2000.0309
11. Failli P, Palmieri L, D'Alfonso C, Giovannelli L, Generini S, Rosso AD, et al. Effect of N-Acetyl-L-Cysteine on Peroxynitrite and Superoxide Anion Production of Lung Alveolar Macrophages in Systemic Sclerosis. Nitric Oxide (2002) 7:277–82. doi: 10.1016/S1089-8603(02)00120-9
12. Yamamoto T, Katayama I, Nishioka K. Nitric Oxide Production and Inducible Nitric Oxide Synthase Expression in Systemic Sclerosis. J Rheumatol (1998) 25:314–7.
13. Matucci Cerinic M, Kahaleh MB. Beauty and the Beast. The Nitric Oxide Paradox in Systemic Sclerosis. Rheumatol (Oxford) (2002) 41:843–7. doi: 10.1093/rheumatology/41.8.843
14. Cotton SA, Herrick AL, Jayson MI, Freemont AJ. Endothelial Expression of Nitric Oxide Synthases and Nitrotyrosine in Systemic Sclerosis Skin. J Pathol (1999) 189:273–8. doi: 10.1002/(SICI)1096-9896(199910)189:2<273::AID-PATH413>3.0.CO;2-4
15. Lazzerini PE, Natale M, Gianchecchi E, Capecchi PL, Montilli C, Zimbone S, et al. Adenosine A2A Receptor Activation Stimulates Collagen Production in Sclerodermic Dermal Fibroblasts Either Directly and Through a Cross-Talk With the Cannabinoid System. J Mol Med (Berl) (2012) 90:331–42. doi: 10.1007/s00109-011-0824-5
16. Cinar R, Iyer MR, Liu Z, Cao Z, Jourdan T, Erdelyi K, et al. Hybrid Inhibitor of Peripheral Cannabinoid 1 Receptors and Inducible Nitric Oxide Synthase Mitigates of Liver Fibrosis. JCI Insight (2016) 1(11):e87336. doi: 10.1172/jci.insight.87336
17. Teixeira-Clerc F, Julien B, Grenard P, Tran Van Nhieu J, Deveaux V, Li L, et al. CB1 Cannabinoid Receptor Antagonism: A New Strategy for the Treatment of Liver Fibrosis. Nat Med (2006) 12:671–6. doi: 10.1038/nm1421
18. Patsenker E, Stoll M, Millonig G, Agaimy A, Wissniowski T, Schneider V, et al. Cannabinoid Receptor Type I Modulates Alcohol-Induced Liver Fibrosis. Mol Med (2011) 17:1285–94. doi: 10.2119/molmed.2011.00149
19. Lin CL, Hsu YC, Lee PH, Lei CC, Wang JY, Huang YT, et al. Cannabinoid Receptor 1 Disturbance of PPARgamma2 Augments Hyperglycemia Induction of Mesangial Inflammation and Fibrosis in Renal Glomeruli. J Mol Med (Berl) (2014) 92:779–92. doi: 10.1007/s00109-014-1125-6
20. Slavic S, Lauer D, Sommerfeld M, Kemnitz UR, Grzesiak A, Trappiel M, et al. Cannabinoid Receptor 1 Inhibition Improves Cardiac Function and Remodelling After Myocardial Infarction and in Experimental Metabolic Syndrome. J Mol Med (Berl) (2013) 91:811–23. doi: 10.1007/s00109-013-1034-0
21. Jourdan T, Godlewski G, Cinar R, Bertola A, Szanda G, Liu J, et al. Activation of the Nlrp3 Inflammasome in Infiltrating Macrophages by Endocannabinoids Mediates Beta Cell Loss in Type 2 Diabetes. Nat Med (2013) 19:1132–40. doi: 10.1038/nm.3265
22. Bronova I, Smith B, Aydogan B, Weichselbaum RR, Vemuri K, Erdelyi K, et al. Protection From Radiation-Induced Pulmonary Fibrosis by Peripheral Targeting of Cannabinoid Receptor-1. Am J Respir Cell Mol Biol (2015) 53:555–62. doi: 10.1165/rcmb.2014-0331OC
23. Cinar R, Gochuico BR, Iyer MR, Jourdan T, Yokoyama T, Park JK, et al. Cannabinoid CB1 Receptor Overactivity Contributes to the Pathogenesis of Idiopathic Pulmonary Fibrosis. JCI Insight (2017) 2:e92281. doi: 10.1172/jci.insight.92281
24. Cinar R, Park JK, Zawatsky CN, Coffey NJ, Bodine SP, Abdalla J, et al. CB1R and iNOS Are Distinct Players Promoting Pulmonary Fibrosis in Hermansky–Pudlak Syndrome. Clin Trans Med (2021) 11:e471. doi: 10.1002/ctm2.471
25. Marquart S, Zerr P, Akhmetshina A, Palumbo K, Reich N, Tomcik M, et al. Inactivation of the Cannabinoid Receptor CB1 Prevents Leukocyte Infiltration and Experimental Fibrosis. Arthritis Rheum (2010) 62:3467–76. doi: 10.1002/art.27642
26. Palumbo-Zerr K, Horn A, Distler A, Zerr P, Dees C, Beyer C, et al. Inactivation of Fatty Acid Amide Hydrolase Exacerbates Experimental Fibrosis by Enhanced Endocannabinoid-Mediated Activation of CB1. Ann Rheum Dis (2012) 71:2051–4. doi: 10.1136/annrheumdis-2012-201823
27. Park JK, Coffey NJ, Bodine SP, Zawatsky CN, Jay L, Gahl WA, et al. Bleomycin Induces Drug Efflux in Lungs. A Pitfall for Pharmacological Studies of Pulmonary Fibrosis. Am J Respir Cell Mol Biol (2020) 62:178–90. doi: 10.1165/rcmb.2018-0147OC
28. Udi S, Hinden L, Ahmad M, Drori A, Iyer MR, Cinar R, et al. Dual Inhibition of Cannabinoid CB1 Receptor and Inducible NOS Attenuates Obesity-Induced Chronic Kidney Disease. Br J Pharmacol (2020) 177:110–27. doi: 10.1111/bph.14849
29. Togami K, Yamaguchi K, Tada H, Chono S. Assessment of Transporter-Mediated Efflux of Nintedanib Using In Vitro Cell Line Models of Idiopathic Pulmonary Fibrosis. Pharmazie (2020) 75:371–4. doi: 10.1691/ph.2020.0048
30. Baranska M, Rychlik-Sych M, Skretkowicz J, Dudarewicz M, Dziankowska-Bartkowiak B, Owczarek J, et al. Genotype and Haplotype Analysis of ABCB1 at 1236, 2677 and 3435 Among Systemic Sclerosis Patients. Autoimmunity (2017) 50:277–82. doi: 10.1080/08916934.2017.1329421
31. Lecru L, Desterke C, Grassin-Delyle S, Chatziantoniou C, Vandermeersch S, Devocelle A, et al. Cannabinoid Receptor 1 Is a Major Mediator of Renal Fibrosis. Kidney Int (2015) 88:72–84. doi: 10.1038/ki.2015.63
32. Udi S, Hinden L, Earley B, Drori A, Reuveni N, Hadar R, et al. Proximal Tubular Cannabinoid-1 Receptor Regulates Obesity-Induced CKD. J Am Soc Nephrol (2017) 28:3518–32. doi: 10.1681/ASN.2016101085
33. Teixeira-Clerc F, Julien B, Grenard P, Tran Van Nhieu J, Deveaux V, Li L, et al. CB1 Cannabinoid Receptor Antagonists: A Novel Approach for the Treatment of Liver Fibrosis. Med Sci (Paris) (2006) 22:683–5. doi: 10.1051/medsci/20062289683
34. Lin CY, Hsu YJ, Hsu SC, Chen Y, Lee HS, Lin SH, et al. CB1 Cannabinoid Receptor Antagonist Attenuates Left Ventricular Hypertrophy and Akt-Mediated Cardiac Fibrosis in Experimental Uremia. J Mol Cell Cardiol (2015) 85:249–61. doi: 10.1016/j.yjmcc.2015.06.010
35. Rajesh M, Batkai S, Kechrid M, Mukhopadhyay P, Lee WS, Horvath B, et al. Cannabinoid 1 Receptor Promotes Cardiac Dysfunction, Oxidative Stress, Inflammation, and Fibrosis in Diabetic Cardiomyopathy. Diabetes (2012) 61:716–27. doi: 10.2337/db11-0477
36. Liu Z, Iyer MR, Godlewski G, Jourdan T, Liu J, Coffey NJ, et al. Functional Selectivity of a Biased Cannabinoid-1 Receptor (CB1R) Antagonist. ACS Pharmacol Transl Sci (2021) 4:1175–87. doi: 10.1021/acsptsci.1c00048
37. Silvestri C, Di Marzo V. The Endocannabinoid System in Energy Homeostasis and the Etiopathology of Metabolic Disorders. Cell Metab (2013) 17:475–90. doi: 10.1016/j.cmet.2013.03.001
38. Jeong WI, Osei-Hyiaman D, Park O, Liu J, Batkai S, Mukhopadhyay P, et al. Paracrine Activation of Hepatic CB1 Receptors by Stellate Cell-Derived Endocannabinoids Mediates Alcoholic Fatty Liver. Cell Metab (2008) 7:227–35. doi: 10.1016/j.cmet.2007.12.007
39. Hezode C, Roudot-Thoraval F, Nguyen S, Grenard P, Julien B, Zafrani ES, et al. Daily Cannabis Smoking as a Risk Factor for Progression of Fibrosis in Chronic Hepatitis C. Hepatology (2005) 42:63–71. doi: 10.1002/hep.20733
40. Despres JP, Golay A, Sjostrom L. Effects of Rimonabant on Metabolic Risk Factors in Overweight Patients With Dyslipidemia. N Engl J Med (2005) 353:2121–34. doi: 10.1056/NEJMoa044537
41. Tam J, Cinar R, Liu J, Godlewski G, Wesley D, Jourdan T, et al. Peripheral Cannabinoid-1 Receptor Inverse Agonism Reduces Obesity by Reversing Leptin Resistance. Cell Metab (2012) 16:167–79. doi: 10.1016/j.cmet.2012.07.002
42. Cinar R, Godlewski G, Liu J, Tam J, Jourdan T, Mukhopadhyay B, et al. Hepatic Cannabinoid-1 Receptors Mediate Diet-Induced Insulin Resistance by Increasing De Novo Synthesis of Long-Chain Ceramides. Hepatology (2014) 59:143–53. doi: 10.1002/hep.26606
43. Tam J, Vemuri VK, Liu J, Batkai S, Mukhopadhyay B, Godlewski G, et al. Peripheral CB1 Cannabinoid Receptor Blockade Improves Cardiometabolic Risk in Mouse Models of Obesity. J Clin Invest (2010) 120:2953–66. doi: 10.1172/JCI42551
44. Cinar R, Iyer MR, Kunos G. The Therapeutic Potential of Second and Third Generation CB1R Antagonists. Pharmacol Ther (2020) 208:107477. doi: 10.1016/j.pharmthera.2020.107477
45. Gaffal E, Cron M, Glodde N, Bald T, Kuner R, Zimmer A, et al. Cannabinoid 1 Receptors in Keratinocytes Modulate Proinflammatory Chemokine Secretion and Attenuate Contact Allergic Inflammation. J Immunol (2013) 190:4929–36. doi: 10.4049/jimmunol.1201777
46. Toth KF, Adam D, Biro T, Olah A. Cannabinoid Signaling in the Skin: Therapeutic Potential of the "C(ut)annabinoid" System. Molecules (2019) 24:918. doi: 10.3390/molecules24050918
47. Jourdan T, Szanda G, Cinar R, Godlewski G, Holovac DJ, Park JK, et al. Developmental Role of Macrophage Cannabinoid-1 Receptor Signaling in Type 2 Diabetes. Diabetes (2017) 66:994–1007. doi: 10.2337/db16-1199
48. Sharif R, Mayes MD, Tan FK, Gorlova OY, Hummers LK, Shah AA, et al. IRF5 Polymorphism Predicts Prognosis in Patients With Systemic Sclerosis. Ann Rheum Dis (2012) 71:1197–202. doi: 10.1136/annrheumdis-2011-200901
49. Xu Y, Wang W, Tian Y, Liu J, Yang R. Polymorphisms in STAT4 and IRF5 Increase the Risk of Systemic Sclerosis: A Meta-Analysis. Int J Dermatol (2016) 55:408–16. doi: 10.1111/ijd.12839
50. Saigusa R, Asano Y, Taniguchi T, Yamashita T, Ichimura Y, Takahashi T, et al. Multifaceted Contribution of the TLR4-Activated IRF5 Transcription Factor in Systemic Sclerosis. Proc Natl Acad Sci USA (2015) 112:15136–41. doi: 10.1073/pnas.1520997112
51. Shreiner AB, Murray C, Denton C, Khanna D. Gastrointestinal Manifestations of Systemic Sclerosis. J Scleroderma Relat Disord (2016) 1:247–56. doi: 10.5301/jsrd.5000214
Keywords: endocannabinoids, skin fibrosis, bleomycin, ATP-binding cassette transporters, P-gp (P-glycoprotein), peripheral CB1R antagonist, cannabinoid (CB) receptor 1, polypharmacology
Citation: Zawatsky CN, Park JK, Abdalla J, Kunos G, Iyer MR and Cinar R (2021) Peripheral Hybrid CB1R and iNOS Antagonist MRI-1867 Displays Anti-Fibrotic Efficacy in Bleomycin-Induced Skin Fibrosis. Front. Endocrinol. 12:744857. doi: 10.3389/fendo.2021.744857
Received: 21 July 2021; Accepted: 09 September 2021;
Published: 28 September 2021.
Edited by:
Isabel Gonzalez Mariscal, Universidad de Málaga, SpainCopyright © 2021 Zawatsky, Park, Abdalla, Kunos, Iyer and Cinar. 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: Resat Cinar, resat.cinar@nih.gov
†These authors have contributed equally to this work