Prerpints.org logo
Preprint
Article

This version is not peer-reviewed.

Repurposing Drugs for Senotherapeutic Effect: Potential Senomorphic Effects of Female Synthetic Hormones

A peer-reviewed article of this preprint also exists.

Submitted:

26 January 2024

Posted:

29 January 2024

You are already at the latest version

Abstract
Repurposing previously approved drugs may fast track the route to clinic for potential senotherapeutics and improves the inefficiency of the clinical drug development pipeline. We carried out a repurposing screen of 240 clinically approved molecules in human primary dermal fibroblasts for effects on CDKN2A expression. Molecules demonstrating effects on CDKN2A expression underwent secondary screening for senescence-associated beta galactosidase (SAB) activity, based on effect size, direction and/or molecule identity. Selected molecules then underwent a more detailed assessment of senescence phenotypes including proliferation, apoptosis, DNA damage, senescence-associated secretory phenotype (SASP) expression and regulators of alternative splicing. A selection of the molecules demonstrating effects on senescence were then used in a new bioinformatic structure-function screen to identify common structural motifs. 90 molecules displayed altered CDKN2A expression at one or other dose, of which 15 also displayed effects on SAB positivity in primary human dermal fibroblasts. Of these, three were associated with increased SAB activity, and 11 with reduced activity. The female synthetic sex hormones; diethylstilboestrol, ethynyl estradiol and levonorgestrel; were all associated with a reduction in aspects of the senescence phenotype in male cells, with no effects visible in female cells. Finally, we identified that the 30 compounds that decreased CDKN2A activity the most had a common substructure linked to this function. Our results suggest that several drugs licenced for other indications may warrant exploration as future senotherapies, but that different donors and potentially different sexes may respond differently to senotherapeutic compounds. This underlines the importance of consideration of donor-related characteristics when designing drug screening platforms.
Keywords: 
Senescence
;  
structure-function screen
;  
synthetic hormone
;  
sex differences
;  
sex-specific
;  
senomorphic
Subject: 
Biology and Life Sciences  -   Biochemistry and Molecular Biology

Introduction

Senescence is a hallmark of ageing, and an emerging therapeutic target [1,2]. Senescence may appear as part of natural development, but during ageing, it is induced by replicative exhaustion, or by cellular stressors such as DNA damage, oncogenes and other forms of cellular stress [3–7]. Despite the original definition that senescence is irreversible, recent research indicates that the senescence phenotypes can be reversed by some classes of drugs [8,9]. Senotherapeutics (compounds that target senescence) include those that attenuate the deleterious characteristics of senescent cells (senomorphics) and drugs that cause preferential lysis of senescent cells (senolytics) [8–10]. Clearance of senescent cells significantly extends lifespan, improves mobility and fur condition in mouse models of progeria, and improves multiple aspects of functionality in aged wild-type mice [11–15]. Senolysis has also been seen to confer additional health benefits in humans; combinations of the senolytic drugs dasatinib and quercetin compounds are currently in trials for diabetic kidney disease and idiopathic pulmonary fibrosis (IPF) [16,17].
Attenuating the senescent state using senomorphic approaches may also be useful. It is possible to uncouple features of senescence such as reversal of senescence-associated beta galactosidase (SAB) staining from other aspects such as proliferation; such effects are often dose-dependent [18]. The ideal senotherapeutic candidate would be able to reverse senescence and attenuate the senescence-associated secretory phenotype (SASP) (a senomorphic effect), but would not necessarily elicit re-entry to cell cycle, since rejuvenated cells may still carry a mutation load. Conversely, any compounds that are identified as increasing senescence might represent potential oncodrugs. Forcing cancerous cells to enter a senescent state might provide a better tolerated oncotherapeutic approach and provide an opportunity to selectively target the resulting cells with senolytic drugs subsequently.
It is likely that some known and licensed drugs have some senomorphic or senolytic capacity. The drug development pipeline is inefficient, with only 15.3% of drugs in phase 1 clinical trials in the US advancing to gain FDA-approval [19]. Repurposing drugs which are already approved for clinical use represents a tactic which avoids the problems with the leaky pipeline of drug development. For example, trametinib, a MEK inhibitor currently used as a cancer treatment, exhibits a biphasic dose response, affecting different aspects of senescence depending on dose [20]. Panels of small molecules for drug repurposing studies can be procured and customised commercially, giving plenty of opportunity to adapt drug repurposing screens for different indications.
Bioinformatic approaches can also be used to complement wet laboratory screening. Structure-function associations may be of particular interest in the context of a screen for senescence. If a certain structure is associated with a senomorphic or senolytic function, then this provides an opportunity to identify potentially useful compounds from public drug databases by screening them for the structure. This strategy could offer the discovery of novel drugs in a quicker way than traditional pharmaceutical discovery processes. Similarly, any structural association with specific senescence related functions may provide mechanistic insight into the cellular processes at hand.
We aimed to screen a range of compounds for effects on aspects of the senescence phenotype using in vitro screens in primary human dermal fibroblasts and bioinformatic structure-function analysis. We identified several existing clinically approved molecules as having capacity to attenuate aspects of the senescence phenotype in a sex-specific manner. Finally, we have worked up a structure-function screening pipeline and identified a molecular substructure that is associated with alterations in CDKN2A expression (a biomarker of senescence that encodes p14 and p16) or SAB positivity. Our work indicates that repurposing studies augmented by bioinformatic or machine learning approaches may prove a rich vein of research for the identification of new classes of senotherapeutic molecules, but donor characteristics such as sex and individual genetics can influence senescence outcome and should be accounted for in study design.

Materials & Methods

Drug Panel, Screen Design and Preparation

A selection of 240 compounds were chosen from the MedChemExpress FDA-Approved Drug Library Plus panel of 2278 compounds (MedChemTronica, Stockholm, Sweden). We selected drugs that target known senescence/cell fate pathways, that cover a variety of other cellular functions (including apoptosis and autophagy), and a variety of commonly prescribed or household medicines. Drug identities and targets are given in Supplementary Table S1 and Supplementary Figure S1 (summary graphs were produced using Microsoft Office). Most compounds are FDA-approved with the remainder approved by the EMA or other countries. Compounds were all supplied at 10 mM concentrations in DMSO. Concordant with similar screens [18,21], compounds were diluted to 1 µM or 10 µM in Gibco™ ultrapure RNase/DNase free water suitable for preparation of cell culture media and laboratory reagents (A1287301, Gibco™, Billings, USA). Vehicle DMSO controls were prepared in the same manner. On the basis of our first pass results, an additional synthetic female hormone (Levonorgestrel; also known as D-Norgestrel) which targets the progesterone receptor was also added to secondary screen. A flow chart (produced using Microsoft Office) illustrating our screening approach is given in Figure 1.

Cells used in this study

Normal human dermal fibroblast (nHDF) cells from one male and one female donor were commercially sourced with full ethical permission granted at source (Promocell, Heidelberg, catalogue number C-12302, lot numbers 445Z026.3 (male) and 467Z026.3 (female). Both donors were Caucasian. The male donor was 36 years old at the time of donation, and the female donor was 28 years old. The cells were taken from the abdomen of the male donor, whereas the female donor’s cells were taken from the breast. Cells were grown in DMEM 1 g/l glucose + phenol red (31885023, Gibco™, Billings, USA), 10% human serum (H3667, Sigma Aldrich, St. Louis, USA) and 1% 10,000 U/ml penicillin - 10,000 µg/ml streptomycin (15140122, Gibco™, Billings, USA). Cells were grown in antibiotic-free media for 48 to 72 hours before seeding, and all treatments were performed without the presence of antibiotics in the medium.

Primary screen

Tissue culture and drug treatment conditions

Male nHDF cells had average cumulated population doublings (cPDL) of 38.91 (range of 34.41 - 40.13 cPDL) at the time of seeding for the primary screen. For this, cells were seeded out in 96-well plates at 6000 cells per well (a density of 1880 cells/cm2) and treated 72 hours after seeding. Dose and incubation times were informed by previous work from our research group and the literature as discussed in the introduction [18,21],22]. On the day of treatment, the medium was removed and replaced with 135 µl of fresh medium and 15 µl of the appropriate stock solution of each drug or control. The drug or control was applied for 24 hours before two washes in DPBS (catalogue number 14190136, Gibco™, Billings, USA) and performing the RNA extraction.

Quantification of CDKN2A expression

RNA was extracted from treated cells using the PureLink™ Pro 96 RNA Purification Kit (catalogue number 12173–011A, Fisher Scientific, Pittsburgh, USA) according to manufacturer’s instructions and eluted in a volume of 45 µl of RNase-free water. RNA quality and quantity was sampled using the Thermo Scientific™ Nanodrop 8000 Spectrophotometer (ThermoFisher, Waltham, Massachusetts, USA). Two test plates had compromised RNA quality and were excluded from the analysis. The maximum RNA volume possible based on reaction volume constraints was reverse-transcribed using the High-Capacity cDNA Reverse Transcription Kit (catalogue number 4368813, ThermoFisher, Waltham, Massachusetts, USA) according to the manufacturer’s instructions. Reverse transcription was carried out on the Applied Biosystems™ Veriti™ 96-Well Fast Thermal Cycler platform. Cycling conditions were: 25°C for 10 minutes, 37°C for 120 minutes, 85°C for 5 minutes and a 4°C hold step. Real-time quantitative PCR (RT-qPCR) reactions were carried out on the Quantstudio 12K platform (Applied Biosystems™, Birchwood, UK) as 5 µl reactions on 384-well plates. Cycling conditions were: 50°C for 2 minutes, 95°C for 10 minutes, followed by 50 cycles of 95°C for 15 seconds and 60°C for 1 minute. Each reaction contained 1 µl of cDNA product, 900 nM each primer and 250 nM probe and TaqMan™ Universal Mastermix II. CDKN2A expression was assessed using the Hs00923894_m1 TaqMan™ Gene Expression Assay (FAM) (catalogue 4331182, ThermoFisher, Waltham City, USA). Endogenous housekeeper control genes were PGK1 (assay ID HS99999906_m1), PPIA (assay ID Hs04194521_s1) and UBC (assay ID Hs05002522_g1), empirically determined to represent the most stable baseline accordingly to the RefFinder webtool [23]. All were procured from ThermoFisher (Waltham, City, USA). Assays were run in two biological and three technical replicates for each compound/control at both treatment doses. Relative gene expression levels were calculated using the comparative CT technique relative to the geometric mean expression level of the three housekeeping genes [24]. Levels were normalised to the average of the vehicle control on each plate and were expressed as natural log to aid against skew of data. The mean ± three standard deviations was used to provide upper and lower bounds for prioritisation of compounds for follow up.

Secondary screen

To assess induction of senescence, experiments were carried out using early passage male cells (cPDL = 32.69), whereas work to assess potential reduction in senescent cell load was carried out using later passage cells (cPDL = 40.77 - 43.2). Later passage cells were assessed at the point that they had slowed to half their original division speed. For assessment of SAB activity, cells were seeded in 12-well plates at an average seeding density of 6226 cells/cm2. Cells were grown for 24 to 48 hours before treatment as in the primary screen and each compound was applied for 24 hours prior to assessment of SAB activity, which was carried out using the Senescence Cells Histochemical Staining Kit (Merck, Gillingham, UK), according to manufacturer’s instructions. Five images per biological replicate were imaged at 10 × magnification using a Zeiss AxioCam ERC55 PrimoVert microscope and later counted manually using ImageJ 1.47v software (US National Institute of Health, Bethesda, Maryland, USA) [25]. Differences in SAB staining between test compounds and controls were assessed by one-way ANOVA with an uncorrected Fisher’s LSD post hoc test and graphed using GraphPad Prism version 9.4.1 for Windows (GraphPad Software, San Diego, California USA, www.graphpad.com).

In depth characterisation of female synthetic hormone compounds

Tissue culture and dosing regime

Based on the results of our primary and secondary screen, we selected three female synthetic hormones for follow up due to evidence of effects on senescence. Cells for this work had an average cPDL of 39.46 at the time of seeding for male cells and 33.68 for female cells and were assessed as being late passage at the point that they had slowed to half their original division speed. Cells were seeded at approximately 7,200 cells/cm2 in a 12-well plate for the SAB assay, at ~ 6,000 cells/cm2 in a 12-well plate on 13 mm coverslips for immunocytochemical staining, at ~ 7,000 cells/cm2 in a 24-well plate on 13 mm coverslips for the TUNEL assay experiments and at ~ 14,000 cells/cm2 in a 6-well plate for RNA extractions. Cells were treated with either a DMSO vehicle control (J66650.AD, Thermo Scientific Alfa Aesar), or a 10 µM dose of diethylstilboestrol, ethynyl estradiol or levonorgestrel (Catalogue numbers HY-14598, HY-B0216 or HY-B0257 respectively (MedChemExpress, Stockholm, Sweden). Fresh medium was added to the plates before the addition of the treatment stock. Cells were treated for 24 hours before staining or harvesting.

Quantification of senescent cell load using SAB staining

Cultures were stained for SAB activity using the Senescence Cells Histochemical Staining kit (CS0030, Merck) according to manufacturer’s instructions. After 24 hours of staining, cells were imaged at 10 × magnification using a Zeiss AxioCam ERC55 PrimoVert. Five images per biological replicate were captured and later counted manually using ImageJ 1.47v software (US National Institute of Health, Bethesda, Maryland, USA) [25]. Differences in SAB staining between test compounds and controls were assessed by one-way ANOVA with an uncorrected Fisher’s LSD post hoc test and graphed using GraphPad Prism version 9.4.1 for Windows (GraphPad Software, San Diego, California USA, www.graphpad.com).

Quantification of cellular proliferation and DNA damage repair using immunocytochemical staining for Ki67 and γH2AX

Following two washes in DPBS (14190136, Gibco™), the cells were fixed with 4% paraformaldehyde and stored in DPBS. Prior to staining, the cells were washed again in DPBS and blocked with ADST [antibody diluent solution - triton: DPBS, 0.1M L-Lysine (303341000, Thermo Scientific™), 1% w/v Human Serum Albumin Fraction V (12668-10GM, Sigma-Aldrich), Triton X-100 (A16046.AP, Thermo Scientific Alfa Aesar)] and 5% human serum (H3667, Sigma Aldrich) for 30 minutes. Antibodies were commercially derived from Abcam: Rb anti-Ki67 (ab15580, ab16667), Ms anti-γH2AX (ab26350), Alexa Fluor ® 555 Goat pAb to Rb (ab150078, ab150086) and Alexa Fluor ® 488 Goat pAb to Ms (ab150117). Primary antibodies were applied overnight at 2.5 µg/ml (suspended in ADST with 2% human serum). Secondary antibodies were applied at 5 µg/ml and 4′,6-diamidino-2-phenylindole (DAPI, D1306, Invitrogen™) at 1 µg/ml (suspended in ADST with 2% human serum) were applied for 1 hour. Then the coverslips were mounted with Dako mounting medium (S302380-2, Agilent). Five representative images per coverslip were captured at 10 × magnification using a Leica DM4 B Upright Microscope and cells were manually scored positive or negative for each parameter manually using the Leica Application Suite X 2019 3.7.1.21655v software (Leica Microsystems, Wetzlar, Germany). Differences in cell kinetic parameters between treated and control cells were assessed by one-way ANOVA with an uncorrected Fisher’s LSD post hoc test and graphed using GraphPad Prism version 9.4.1 for Windows (GraphPad Software, San Diego, California USA, www.graphpad.com).

Quantification of apoptosis using TUNEL assay

Cells were washed in DPBS (14190136, Gibco™), before the cells were fixed with 4% paraformaldehyde, washed again and stored in DPBS. The Click-iT® TUNEL Alexa Fluor® Imaging Assay (C10245, ThermoFisher) was performed according to the manufacturer’s instructions using additional DPBS, bovine serum albumin (BSA) fraction V fatty acid-free (10775835001, Roche), and Triton X-100 (A16046.AP, Thermo Scientific Alfa Aesar). In the same manner as for the other immunofluorescently-stained cells, the Leica DM4 B Upright Microscope at 10 × magnification was used to capture five images per coverslip. The cells in the images were later counted manually using Leica Application Suite X 2019 3.7.1.21655v software (Leica Microsystems, Wetzlar, Germany). Differences in TUNEL staining between test compounds and controls were assessed by one-way ANOVA with an uncorrected Fisher’s LSD post hoc test and graphed using GraphPad Prism version 9.4.1 for Windows (GraphPad Software, San Diego, California USA, www.graphpad.com).

Quantitative RT-qPCR assessment of gene expression

RNA was extracted from cells using TRI Reagent Solution (AM9738, Invitrogen™) according to the manufacturer’s instructions, with the addition of 10 mM MgCl2 (AM9530G, Invitrogen™) before phase separation (to aid in RNA recovery [26]) and 1.2 µl of 15 mg/ml GlycoBlue™ Coprecipitant (AM9516, Invitrogen™) prior to washing (to aid in pellet visualisation). RNA was resuspended in 20 µl 1 × TE buffer, pH 8.0 (BP2473-500, Fisher Bioreagents) and assessed for concentration and quality using the Thermo Scientific™ Nanodrop 8000 Spectrophotometer (Thermo Fisher Scientific, Waltham, Massachusetts, USA). RNA was reverse-transcribed at 10 ng/µl in a 20 µl reaction, using the High-Capacity cDNA Reverse Transcription Kit (4368813, Applied Biosystems™) according to the manufacturer’s instructions. Reverse transcription was performed on an Applied Biosystems™ Veriti™ 96-Well Fast Thermal Cycler with the following cycling conditions: 25°C for 10 minutes, 37°C for 120 minutes, 85°C for 5 minutes and a 4°C hold step. 12.5 ng of cDNA was pre-amplified according to the manufacturer’s instructions using TaqMan™ PreAmp Master Mix (4384266, Applied Biosystems™) and pooled TaqMan™ Gene Expression Assays (FAM) (4331182, TaqMan®). Transcripts encoding factors associated with apoptosis, the senescence-associated secretory phenotype (SASP) and regulators of alternative splicing were assessed by RT-qPCR. A table of genes assessed is provided in Table 1. Using the Applied Biosystems™ Veriti™ 96-Well Fast Thermal Cycler, the cycling conditions were: 95°C for 10 minutes, 14 cycles of [95°C for 15 seconds, 60°C for 4 minutes], 99°C for 10 minutes, and a 4°C hold step. The pre-amplified cDNA products were diluted by a factor of 10 in 1 × TE buffer, pH 8.0 (BP2473-500, Fisher Bioreagents). RT-qPCR was performed in three biological and three technical replicates on the Quantstudio 12K platform (Applied Biosystems™) as 5 µl reactions on 384-well plates. The cycling conditions were: 50°C for 2 minutes, 95°C for 10 minutes, followed by 50 cycles of 95°C for 15 seconds and 60°C for 1 minute. 1 µl of diluted, pre-amplified cDNA product was used per reaction with 0.25 µl of Taqman™ Gene Expression Assay (equating to 900 nM primer and 250 nM probe). Gene expression was calculated using the comparative CT technique [24] relative to the geometric mean of five housekeeping genes (GUSB, IDH3B, PGK1, PPIA and UBC) empirically selected for stability as described above [23], and normalised to expression levels in the respective cell type’s vehicle-treated control. Results were assessed for statistical significance using a one-way ANOVA with an uncorrected Fisher’s LSD post hoc test and graphed using GraphPad Prism version 9.4.1 for Windows (GraphPad Software, San Diego, California USA, www.graphpad.com).

Bioinformatic assessment of structure-function relationships

Structural information on each compound tested was obtained from the supplier (MedChemTronica, Stockholm). The SMILES (simplified molecular input line entry system) data was transformed into SDF (structure data file) data for analysis with ChemmineR and fmcsR packages in Rstudio software version 4.1.0 [27–30]. Tanimoto coefficients, measures of structural similarity [27], were computed for each pair of compounds and used to construct a matrix of intragroup comparisons of structural similarity. Computing resource limitations imposed a maximum of 30 compounds per test group. The average Tanimoto coefficient across the matrix of functionally related compounds was compared against the average Tanimoto coefficient for a control group of non-functionally related compounds using an unpaired t test in GraphPad Prism version 9.4.1 for Windows (GraphPad Software, San Diego, California USA, www.graphpad.com). A dendrogram was constructed for groups of interest to illustrate the structural similarity between compounds. If a structure-function relationship was suggested, the exact maximum common substructure was computed for the two least similar compounds (as identified in the dendrogram) to identify the maximum common substructure across the whole test group.

Methodological validation

Making a statistical comparison between similarity matrices represents a novel use for the matrix outputs of ChemmineR. A significant difference between intragroup average Tanimoto coefficients indicates that compounds in the test group are more structurally similar than the control group. When the test group contains only compounds with a particular function of interest, a significant difference may suggest a structure-function relationship. The maximum common substructure of the group may therefore suggest (or be incorporated within) a substructure which is associated with the function of interest. To validate this approach, a group of compounds that share a known functionally related substructure was compared against a control group. Validation compounds were selected from the MedChemExpress FDA-Approved Drug Library Plus panel of 2278 compounds (MedChemTronica, Stockholm, Sweden).
Given that our in vitro screens had already highlighted some oestrogenic compounds and that the provided drug library information identified compounds that target the oestrogen receptor, we decided to validate the approach by trying to identify a known structure-function relationship using compounds that target the oestrogen receptor. These compounds are known to share substructures which are linked with their function of targeting the oestrogen receptor. The validation test group consisted of 30 compounds versus a control group of 30 functionally unrelated compounds. The number of oestrogenic compounds in the validation test group was varied to assess the sensitivity of the method: using 30, 10, 4, 3 and 2 oestrogenic compounds in a group of other non-oestrogenic compounds totalling 30 for comparison against the control group of 30 functionally unrelated compounds.

Structure-function analysis of in vitro screen results

The first two test groups were the compounds that had either increased or decreased CDKN2A gene expression the most (averaged across both doses). The third test group was a selection of compounds that had decreased SAB activity in the screen. Control compounds acting as a non-functionally associated control group were selected based on the compounds with the least effect on CDKN2A expression and were matched to the number of compounds in each test group. 78 individual compounds were used in the study. Test groups of compounds are given in Supplementary Figure S2, Supplementary Figure S3 and Supplementary Figure S4.

Results

Primary and secondary screens

We identified 90 molecules that altered CDKN2A gene expression, with 20 increasing senescent cell load and 70 decreasing senescent cell load by more than mean ± three standard deviations of the control treatments (Table 2). 32 compounds were selected for secondary screening based on effect size, widespread usage or due having different effect directionality between doses from the primary screen. Of these, 11 compounds elicited a reduction in SAB positivity, and three caused an increase in SAB positivity (Figure 2; Table 3). Compounds causing a statistically significant decrease in SAB activity included the non-steroidal anti-inflammatory drug (NSAID) aspirin, the cancer drugs cabozantinib, and carmofur, the antihistamine chlorpheniramine (maleate), the 11β-hydroxylase inhibitor metyrapone, the antipsychotic penfluridol, the ammonia lowering drug sodium-4-phenylbutyrate and the synthetic female sex hormones diethylstilboestrol, ethynyl estradiol and levonorgestrel. Most of these effects were evident at 10 µM concentration, though aspirin and penfluridol had effects at a lower concentration of 1 µM. Although the cancer drug, sunitinib, caused a significant decrease in SAB activity at 10 µM in both early and late passage cells, the drug caused mass cell death rather than acting as a senotherapeutic. Compounds demonstrating induction of senescence included the anticancer agents doxorubicin, homoharringtonine and imatinib.

Potential donor characteristic-specific differences in cellular senescence kinetics in response to treatment with female synthetic sex hormones.

Female synthetic sex hormones were prominent across both screens, so we examined the effects of these compounds in more detail in senescent male and female primary dermal fibroblasts. We identified that all three synthetic female hormones caused a decrease in SAB activity in male fibroblasts (a 30%, 32% and 51% decrease in stained cells for diethylstilboestrol, ethynyl estradiol and levonorgestrel; p = 0.0122, 0.0083 and 0.0002 respectively). Notably, these effects were not evident in the female cells (Figure 3A; Table 4).
Effects on proliferation were minimal, with only diethylstilboestrol demonstrating a 45% decrease in proliferation (p = 0.0289) in male cells (Figure 3B; Table 4). Levels of γH2AX, (indicating DNA damage repair) were very low in all cell types and treatments, reflected in very low levels of cell death in the culture as measured by TUNEL assay; an average of 2.1% of cells had evidence of double strand breaks with no significant difference noted between any of the experimental groups (Figure 3C; Table 4). We also noted some donor-specific changes in apoptotic markers (Figure 3; Table 5; Table 6). We observed an 82% and a 91% increase in BCL2 expression following treatment with ethynyl estradiol or levonorgestrel in female cells (p = 0.0351 and 0.0214). BCL2 expression was unchanged in male cells. Conversely, CASP3 was increased by 51% in response to levonorgestrel (p = 0.0397) in male cells but was unchanged in female cells.
Furthermore, the treatments only affected SASP factor expression in the female cells, appearing to be mildly pro-inflammatory (Figure 4; Table 5; Table 6). Diethylstilboestrol caused an 89% increase in IL6 expression (p = 0.0057) and increased IL8 expression by 74% (p = 0.0062). Ethynyl estradiol caused a 54% increase in expression of CXCL1 (p = 0.0174). Levonorgestrel caused a very large effect in CXCL10 expression (a 14-fold increase, p = 0.0162) but it is important to note that gene expression of CXCL10 in the controls was very low. Levonorgestrel also elicited a 40% and an 42% increase in CXCL1 and IL12A expression respectively (p = 0.0033 and 0.0148) No other SASP markers were altered in either the male or the female cells.
Splicing factor dysregulation is known to be a driver of senescence, and targeted restoration of splicing factor expression yields senomorphic effects [31]. We noted differences in the expression of splicing factor genes between the male and the female primary dermal fibroblasts in response to synthetic female sex hormones (Figure 5; Table 5; Table 6). Diethylstilboestrol, ethynyl estradiol and levonorgestrel induced a 55%, 50% and 62% increase in HNRNPK expression in male cells (p = 0.0029, 0.0060 and 0.0012 respectively), whereas diethylstilboestrol caused a 73% decrease in female cells (p = 0.0003). In male cells, ethynyl estradiol induced a 24% increase in SRSF6 expression (p = 0.0449) and levonorgestrel induced a 30% increase in TRA2B expression (p = 0.0361), but no effects on expression of either gene were observed in female cells. In female cells, ethynyl estradiol and levonorgestrel induced a 78% and 72% increase in the gene expression of the spliceosomal component, NOVA1 (p = 0.0065 and 0.0109 respectively), and induced a 44% and 42% increase in PNISR expression (p = 0.0372 and 0.0452 respectively), whereas male cells were unaffected.

A common substructure was identified for compounds that decreased CDKN2A

A structure-function analysis was used to identify any substructure associated with compounds that were grouped by their functionality from the screens for CDKN2A gene expression and/or SAB activity. Prior to this analysis, methodological validation of the bioinformatic statistical approach indicated that the method is not very sensitive. Only the methodological validation test group consisting of 30 oestrogen receptor-targeting compounds (mean Tanimoto coefficient ± standard error of the mean (SEM); 0.3634 ± 0.007869, n = 900) versus 30 control compounds (0.1886 ± 0.005936, n = 900) was significantly different (p < 0.0001). Variations of the 30-compound validation test group with fewer oestrogenic compounds and more functionally unrelated compounds were not significant when compared against the control group: 10 oestrogenic compounds (0.1877 ± 0.006082, n = 900, p = 0.9167), five oestrogenic compounds (0.1820 ± 0.006138, n = 900, p = 0.4418), three oestrogenic compounds (0.1885 ± 0.006182, n = 900, p = 0.9985), two oestrogenic compounds (0.1861 ± 0.006148, n = 900, p = 0.7690), and a group with zero oestrogenic compounds (0.1851 ± 0.006323, n = 900, p = 0.6895).
78 individual compounds (and the associated data on their effects on CDKN2A expression and/or SAB activity) were used to provide input for a structure-function analysis (Supplementary Figure 2, Supplementary Figure 3 and Supplementary Figure 4). The first test group consisted of the 30 compounds that decreased CDKN2A expression the most (averaged across both doses). Structures for this group are shown in Supplementary Figure 2. This test group (0.2107 ± 0.005580, n = 900) was significantly structurally different from the control group (0.1755 ± 0.006050, n = 900) with a p < 0.0001. A dendrogram was constructed to visualise the structural similarity of the compounds in the group in Figure 6A. The maximum common substructure between the two least similar compounds in the second test group is shown in red in Figure 6B.
The number of compounds that increased SAB activity (n = 3) was too small to be suitable for this type of analysis, however the number of compounds that decreased SAB activity was appropriate (n = 11). Sunitinib was omitted from the analysis as it caused mass cell death rather than acting to reduce senescence. Test group two therefore consisted of a group of ten compounds that decreased SAB activity (Supplementary Figure 3), and their Tanimoto coefficients were compared against the coefficients of a control group of ten compounds that did not have an effect on CDKN2A expression. The average Tanimoto coefficient was not significantly different in the test group (0.2928 ± 0.02544, n = 100) compared to the control group (0.2524 ± 0.02900, n = 100, p = 0.2964).
The third test group (Supplementary Figure 3) comprised the eight compounds that increased CDKN2A expression (above the mean ± 3 SDs criterion) when averaged across both doses. The average Tanimoto coefficient of this group (0.2817 ± 0.03492: mean ± SEM, n = 64) had no significant difference when compared against the average Tanimoto coefficient of eight compounds that had no effect on CDKN2A expression (0.2482 ± 0.03798, n = 64, p = 0.5169).

Discussion

We carried out a drug repurposing screen on 240 FDA-approved molecules for effects on cellular senescence phenotypes. We identified 90 compounds that have effects on CDKN2A expression in human primary dermal fibroblasts, 11 of which bring about a reduction in senescent cell load and 3 of which increase senescent cell load as measured by SAB activity. Three of the compounds that reduce senescent cell load are the synthetic female sex hormones diethylstilboestrol, ethynyl estradiol and levonorgestrel, which exert senotherapeutic effects in male dermal fibroblasts, but not in female cells, where their effects are mildly inflammatory. Finally, we have identified a chemical substructure associated with reduced CDKN2A expression and/or SAB activity. Our findings are important for future research into drugs to target the molecular basis of ageing, as they indicate that some senotherapeutic effects may be specific to certain donor characteristics e.g. sex, which has major implications for therapeutic screening cascades and eventual population level treatment options.
Several of the compounds that we identified as having effects on senescence phenotypes are frequently prescribed or are common household drugs, some of which have also been previously linked with pathways associated with ageing. Aspirin, for example has been shown to extend lifespan in mice [32], metyrapone is an 11β-hydroxylase inhibitor known to activate autophagy [33,34] and penfluridol, a potent antipsychotic medication, has been shown to increase lifespan in Drosophila melanogaster [35]. Several known senotherapeutic compounds (dasatinib (hydrochloride), metformin (hydrochloride), resveratrol and trametinib) were not amongst the largest effects on CDKN2A, suggesting that smaller effects could also be worth examining in similar screens in the future. In the primary screen, we observed that some drug classes had interesting effects on senescence, suggesting that more research is needed into these effects as it may be that certain drugs are more (or less) suitable for use in older patients due to their effects on senescence. Unsurprisingly, we detected effects on senescence kinetics for drugs used in the treatment of cancer, but perhaps less predictably, we also detected effects for antidepressant drugs, anticonvulsant drugs, and female synthetic sex hormones.
Our study has identified a maximum common structural motif that was present even in molecules with very little other structural similarity. This compares well with work in the literature from Olascoaga-Del Angel et al., where several chemotypes associated with senomorphic or senolytic properties were identified [36]. The maximum common substructure that we identified was also common across 11 of the 13 structures in their larger-scale analysis [36]. This finding is strengthened when we consider that the new approach used for the identification of structural similarities was not very sensitive, as noted during the methodology validation.
Diethylstilboestrol, ethynyl estradiol and levonorgestrel were all associated with a decrease in CDKN2A expression in male cells. These drugs are commonly used in hormone replacement therapy or contraceptives [37]. Female hormones are associated with protective benefits in ageing [38–40], and there is some evidence of sex differences in senescence-associated phenotypes [41–51]. It is clear that being biologically female offers protective benefits against ageing [45,52], and the two main female hormones, oestrogen and progesterone, are known to be involved in many ageing and senescence-related pathways [46,53–56]. The typical nuclear receptors for these two hormones, the oestrogen receptors (ERα and ERβ) and progesterone receptors (PR-A and PR-B) are involved in the same pathways [55,57]. There is comparatively little information about the senotherapeutic properties of synthetic female sex hormones in humans [47,58,59]; most research has been carried out in mouse models treated with synthetic oestrogens [56,60–62].
We found differences in SAB positivity, expression of splicing factors and expression of mRNAs encoding SASP proteins between male and female cells in response to female sex hormones. Sex differences in drug responses are not uncommon, and a sexual dimorphism has been reported in mice in response to senotherapeutics [63,64]. Recently, the NIA Interventions Testing Program in mice has revealed sex differences in effects on longevity in response to 17-α-estradiol and aspirin [65,66]. Anthropometric parameters such as bodyweight, fat distribution and differences in pharmacokinetics and pharmacodynamics means that women are more sensitive to some drugs, have altered clearance kinetics and may experience more drug interactions [67]. In humans, oestrogen and progesterone are endogenous to both sexes, but differ in their circulating levels [68,69]. Unlike progesterone, there are many forms of oestrogen: estrone (E1), estradiol (E2), estriol (E3) and other minor oestrogens, but the major oestrogen is E2. This has two isoforms: 17α-estradiol and the more potent and biologically-most relevant 17β-estradiol [57,70]. Oestrogens are discussed more often than progesterones in relation to senescence, but in this study levonorgestrel, a progesterone, had a larger effect on senescence than the oestrogens. Diethylstilboestrol decreased proliferation in male cells, which is at odds with oestrogen’s often growth-inducing effects, e.g. during the female pubertal growth spurt [68]. At the present time, it is not clear whether the observed sex differences arise from differences in bioavailability, or from an undescribed non-canonical role of the hormones over and above canonical oestrogen/progesterone signalling, particularly given the senomorphic effect occurs with treatment of either a synthetic oestrogen or a progesterone. The classical signalling pathways for both oestrogen and progesterone feature the hormone and its respective nuclear receptor(s) acting as ligand-activated transcription factors. The complex binds to hormone responsive elements (HREs) in the genome to control gene expression. There are many HREs across the genome, for example there are over 70,000 oestrogen-responsive-elements identified [71]. Both hormones can act via other pathways, including membrane bound GPCRs. Activation of their respective GPCRs can activate cell fate pathways such as Ras/Raf/MEK/ERK and PI3K/Akt, as well as cross-signalling with classical hormonal signalling pathways [72,73]. Both pathways have previously been implicated in senescence in human cells [20] and in lifespan in invertebrate models [74]. Differing expression, activity and/or sensitivity of receptors between the sexes might also be factoring into the senotherapeutic effect observed in this study. Another consideration is that the female fibroblasts used in our study were donated by a pre-menopausal woman: it is possible that cells from women who are undergoing or have gone through the menopause may have differing responses to synthetic female hormones, or indeed they may have a similar effect compared to the effect seen in the male cells.
Translating the findings of repurposing screens into the clinic needs careful consideration. When considering these compounds in vivo, dosage is also a factor. Many compounds associated with senomorphic effect display biphasic dose responses, which may arise from the autoregulatory relationships between the affected genes and pathways [20]. It is therefore possible that repeated exposure and/or higher/lower dosage may have different effects in a systemic setting. It is also possible that the effects may be tissue specific. Repurposing drugs identified to have senotherapeutic effect may also not be clinically feasible as severe side effects may alter the risk-benefit relationship for milder age-related diseases. The three female synthetic hormones identified in this study do not currently offer a potential clinical application as a mainstream senotherapeutic drug as the effect is not observed in females who routinely take the medicines, and males taking the hormones may have feminising side-effects.
In conclusion, our work demonstrates the utility of repurposing screens, combined with bioinformatic structure-function analyses to identify chemical structures that may be suitable for eventual senotherapeutic benefit. Our study suggests that the sexual dimorphisms in senomorphic/geroprotective effects in animal models may also exist in human cells. We identify several compounds of interest for future senotherapeutic research in the screen including the three female synthetic hormones. We use a new approach to also identify a chemical substructure associated with a decrease in senescence. Our work also highlights the need for patient characteristics such as biological sex to be taken into consideration even in early in vitro pre-clinical work; high throughput screening cascades are often carried out using a single clone of a well characterised transformed cell line, and other senotherapeutic compounds may be sex-specific. This statement could equally be applied to other individual anthropometric or genetic characteristics. Biological sex in in vitro experiments can cause dimorphic effects and this should be considered more regularly when designing experiments, particularly in the process of investigating senotherapeutic compounds. The easiest cell type may not always be the best candidate for such screens. However, provided studies are designed appropriately to factor in donor characteristics such as sex, repurposing remains a potent mechanism for identifying new jobs for old drugs.

Ethical and conflict of interest statements

Ethical statement:

This work was performed with commercially sourced cell lines from donors with progeroid syndromes, and as such, ethical requirements were fulfilled at source.

Author Contributions

LRB designed the experiment, performed laboratory work, data analysis and prepared the manuscript. RF performed laboratory work and helped prepare the manuscript. LWH provided overall oversight and managed the project.

Funding

The authors gratefully acknowledge funding from Animal Free Research UK. No reagents containing animal-derived biomaterials were purchased with this grant.

Acknowledgements

This research was supported by the National Institute for Health and Care Research (NIHR) Exeter Biomedical Research Centre (BRC). The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care.

Conflicts of Interest

LWH is an inventor on patent PCT/GB2019/052125, and is a founder, chief scientific officer, and research and development lead for SENISCA, Ltd. LRB and RF have no financial interest to declare in relation to the content of this manuscript.

References

  1. M. Y. Terzi, M. Izmirli, and B. Gogebakan, “The cell fate: senescence or quiescence,”. Mol. Biol. Rep. 2016, 43, 1213–1220. [CrossRef]
  2. J. Tigges et al., “The hallmarks of fibroblast ageing,”. Mech. Ageing Dev. 2014, 138, 26–44. [CrossRef] [PubMed]
  3. D. Muñoz-Espín et al., “Programmed cell senescence during mammalian embryonic development,”. Cell 2013, 155, 1104. [CrossRef] [PubMed]
  4. M. Storer et al., “Senescence Is a Developmental Mechanism that Contributes to Embryonic Growth and Patterning,”. Cell 2013, 155, 1119–1130. [CrossRef] [PubMed]
  5. L. Hayflick and P. S. Moorhead, “The serial cultivation of human diploid cell strains,”. Exp. Cell Res. 1961, 25, 585–621. [CrossRef] [PubMed]
  6. Toussaint et al., From the Hayflick mosaic to the mosaics of ageing. Role of stress-induced premature senescence in human ageing,”. Int. J. Biochem. Cell Biol. 2002, 34, 1415–1429. [CrossRef]
  7. M. Serrano, A. W. Lin, M. E. McCurrach, and D. Beach, “Oncogenic ras Provokes Premature Cell Senescence Associated with Accumulation of p53 and p16 INK4a,”. Cell 1997, 88, 593–602. [CrossRef] [PubMed]
  8. L. J. Niedernhofer and P. D. Robbins, “Senotherapeutics for healthy ageing,”. Nat. Publ. Gr. 2018. Available online: www.nature.com/nrd (accessed on 1 January 2020).
  9. B. P. Lee and L. W. Harries, “Senotherapeutic drugs: A new avenue for skincare?,”. Plast. Reconstr. Surg. 2021, 148, 21S–26S. [CrossRef] [PubMed]
  10. J. L. Kirkland, T. Tchkonia, Y. Zhu, L. J. Niedernhofer, and P. D. Robbins, “The Clinical Potential of Senolytic Drugs,”. J. Am. Geriatr. Soc. 2017, 65, 2297–2301. [CrossRef]
  11. D. J. Baker et al., “Clearance of p16 Ink4a-positive senescent cells delays ageing-associated disorders,”. Nature 2011, 479, 232–236. [CrossRef]
  12. D. J. Baker et al., “Naturally occurring p16Ink4a-positive cells shorten healthy lifespan,”. Nature 2016, 530, 184–189. [CrossRef]
  13. Y. Gu et al., “The tyrosine kinase inhibitor Dasatinib reduces cardiac steatosis and fibrosis in obese, type 2 diabetic mice,”. Cardiovasc. Diabetol. 2023, 22, 22. [CrossRef]
  14. Krzystyniak et al., “Combination of dasatinib and quercetin improves cognitive abilities in aged male Wistar rats, alleviates inflammation and changes hippocampal synaptic plasticity and histone H3 methylation profile,”. Aging (Albany. NY). 2022, 14, 572–595. [CrossRef]
  15. M. Dungan et al., “Deletion of SA β-Gal+ cells using senolytics improves muscle regeneration in old mice,”. Aging Cell 2022, 21. [CrossRef]
  16. L. J. Hickson et al., “Senolytics decrease senescent cells in humans: Preliminary report from a clinical trial of Dasatinib plus Quercetin in individuals with diabetic kidney disease,”. EBioMedicine 2019, 47, 446–456. [CrossRef]
  17. J. N. Justice et al., “Senolytics in idiopathic pulmonary fibrosis: Results from a first-in-human, open-label, pilot study,”. EBioMedicine 2019, 40, 554–563. [CrossRef]
  18. E. Latorre, R. Torregrossa, M. E. Wood, M. Whiteman, and L. W. Harries, “Mitochondria-targeted hydrogen sulfide attenuates endothelial senescence by selective induction of splicing factors HNRNPD and SRSF2,”. Aging (Albany. NY). 2018, 10, 1666–1681. [CrossRef]
  19. M. Hay, D. W. Thomas, J. L. Craighead, C. Economides, and J. Rosenthal, “Clinical development success rates for investigational drugs,”. Nat. Biotechnol. 2014, 32, 40–51. [CrossRef]
  20. E. Latorre, E. L. Ostler, R. G. A. Faragher, and L. W. Harries, “FOXO1 and ETV6 genes may represent novel regulators of splicing factor expression in cellular senescence,”. FASEB J. 2019, 33, 1086–1097. [CrossRef]
  21. E. Latorre et al., “Small molecule modulation of splicing factor expression is associated with rescue from cellular senescence,”. BMC Cell Biol. 2017, 18. [CrossRef]
  22. V. Barra, R. F. Chiavetta, S. Titoli, I. M. Provenzano, P. S. Carollo, and A. Di Leonardo, “Specific Irreversible Cell-Cycle Arrest and Depletion of Cancer Cells Obtained by Combining Curcumin and the Flavonoids Quercetin and Fisetin,”. Genes (Basel). 2022, 13, 1125. [CrossRef]
  23. F. Xie, P. Xiao, D. Chen, L. Xu, and B. Zhang, “miRDeepFinder: A miRNA analysis tool for deep sequencing of plant small RNAs,”. Plant Mol. Biol. 2012, 80, 75–84. [CrossRef] [PubMed]
  24. M. W. Pfaffl, “A new mathematical model for relative quantification in real-time RT–PCR,”. Nucleic Acids Res. 2001, 29, 2002–2007. [CrossRef]
  25. C. A. Schneider, W. S. Rasband, and K. W. Eliceiri, “NIH Image to ImageJ: 25 years of image analysis,”. Nat. Methods 2012, 9, 671–675. [CrossRef] [PubMed]
  26. Y. K. Kim, J. Yeo, B. Kim, M. Ha, and V. N. Kim, “Short Structured RNAs with Low GC Content Are Selectively Lost during Extraction from a Small Number of Cells,”. Mol. Cell 2012, 46, 893–895. [CrossRef] [PubMed]
  27. Y. Cao, A. Charisi, L. C. Cheng, T. Jiang, and T. Girke, “ChemmineR: a compound mining framework for R,”. Bioinformatics 2008, 24, 1733–1734. [CrossRef]
  28. T. W. H. Backman, Y. Cao, and T. Girke, “ChemMine tools: An online service for analyzing and clustering small molecules,”. Nucleic Acids Res. 2011, 39, 486–491. [CrossRef]
  29. Y. Wang, T. W. H. Backman, K. Horan, and T. Girke, “FmcsR: Mismatch tolerant maximum common substructure searching in R,”. Bioinformatics 2013, 29, 2792–2794. [CrossRef]
  30. RStudio Team, “RStudio: Integrated Development Environment for R.” RStudio, PBC, Boston, Massachusetts, 2020. Available online: http://www.rstudio (accessed on 20 October 2020).
  31. C. Holly et al., “Changes in splicing factor expression are associated with advancing age in man,”. Mech. Ageing Dev. 2013, 134, 356–366. [CrossRef]
  32. R. Strong et al., “Nordihydroguaiaretic acid and aspirin increase lifespan of genetically heterogeneous male mice,”. Aging Cell 2008, 7, 641–650. [CrossRef]
  33. H. Jahn, M. Schick, F. Kiefer, M. Kellner, A. Yassouridis, and K. Wiedemann, “Metyrapone as Additive Treatment in Major Depression: A Double-blind and Placebo-Controlled Trial,”. Arch. Gen. Psychiatry 2004, 61, 1235–1244. [CrossRef]
  34. M. Holczer, M. Márton, A. Kurucz, G. Bánhegyi, and O. Kapuy, “A Comprehensive Systems Biological Study of Autophagy-Apoptosis Crosstalk during Endoplasmic Reticulum Stress,”. Biomed Res. Int. 2015, 2015. [CrossRef]
  35. H. L. Kang, S. Benzer, and K. T. Min, “Life extension in Drosophila by feeding a drug,”. Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 838–843. [CrossRef] [PubMed]
  36. K. S. Olascoaga-Del Angel, H. Gutierrez, M. Königsberg, J. Pérez-Villanueva, and N. E. López-Diazguerrero, “Exploring the fuzzy border between senolytics and senomorphics with chemoinformatics and systems pharmacology,”. Biogerontology 2022, 23, 453–471. [CrossRef] [PubMed]
  37. FSRH, “FSRH Clinical Guideline: Combined Hormonal Contraception,”. 2020. Available online: https://www.fsrh (accessed on 21 October 2020).
  38. C. J. Pike, “Sex and the development of Alzheimer’s disease,”. J. Neurosci. Res. 2017, 95, 671–680. [CrossRef] [PubMed]
  39. B. Tramunt et al., “Sex differences in metabolic regulation and diabetes susceptibility,”. Diabetologia 2020, 63, 453–461. [CrossRef] [PubMed]
  40. V. Regitz-Zagrosek and G. Kararigas, “Mechanistic pathways of sex differences in cardiovascular disease,”. Physiol. Rev. 2017, 97, 1–37. [CrossRef] [PubMed]
  41. Ruggierii, S. Anticoli, A. D’ambrosio, L. Giordani, and M. Mora, “The influence of sex and gender on immunity, infection and vaccination,”. Ann. Ist. Super. Sanita 2016, 52, 198–204. [CrossRef]
  42. Borrás, J. Sastre, D. García-Sala, A. Lloret, F. V. Pallardó, and J. Viña, “Mitochondria from females exhibit higher antioxidant gene expression and lower oxidative damage than males,”. Free Radic. Biol. Med. 2003, 34, 546–552. [CrossRef] [PubMed]
  43. S. Groß, U.-D. Mmel, M. Klintschar, and F. Bartel, “Germline Genetics of the p53 Pathway Affect Longevity in a Gender Specific Manner,”. Curr. Aging Sci. 2014, 7. [CrossRef]
  44. Y. Han, S. A. Wennersten, J. M. Wright, R. W. Ludwig, E. Lau, and M. P. Y. Lam, “Proteogenomics reveals sex-biased aging genes and coordinated splicing in cardiac aging,”. Am. J. Physiol. Hear. Circ. Physiol. 2022, 323, H538–H558. [CrossRef]
  45. S. Hägg and J. Jylhävä, “Sex differences in biological aging with a focus on human studies,”. Elife 2021, 10, 1–27. [CrossRef]
  46. J. B. Rubin et al., “Sex differences in cancer mechanisms,”. Biol. Sex Differ. 2020, 11, 1. [CrossRef]
  47. M. Ng and L. N. Hazrati, “Evidence of sex differences in cellular senescence,”. Neurobiol. Aging, vol. 2022, 120, 88–104. [CrossRef]
  48. M. J. Yousefzadeh et al., “Tissue specificity of senescent cell accumulation during physiologic and accelerated aging of mice,”. Aging Cell 2020, 19, e13094. Available online: https://pubmed.ncbi.nlm.nih.gov/31981461/ (accessed on 25 October 2020).
  49. M. Waskar et al., “Drosophila melanogaster p53 has developmental stage-specific and sex-specific effects on adult life span indicative of sexual antagonistic pleiotropy,”. Aging (Albany. NY). 2009, 1, 903–936. [CrossRef]
  50. M. Rall-Scharpf, T. W. P. Friedl, S. Biechonski, M. Denkinger, M. Milyavsky, and L. Wiesmüller, “Sex-specific differences in DNA double-strand break repair of cycling human lymphocytes during aging,”. Aging. 2021, 13, 21066–21089. [CrossRef] [PubMed]
  51. R. Trzeciak et al., “Age, sex, and race influence single-strand break repair capacity in a human population NIH Public Access,”. Radic Biol Med 2008, 45, 1631–1641. [CrossRef] [PubMed]
  52. T. L. Rochelle, D. K. Y. Yeung, M. H. Bond, and L. M. W. Li, “Predictors of the gender gap in life expectancy across 54 nations,”. Psychol. Heal. Med. 2015, 20, 129–138. [CrossRef]
  53. Y. Sasaki, Y. Ikeda, T. Miyauchi, Y. Uchikado, Y. Akasaki, and M. Ohishi, “Estrogen-sirt1 axis plays a pivotal role in protecting arteries against menopause-induced senescence and atherosclerosis,”. J. Atheroscler. Thromb. 2020, 27, 47–59. [CrossRef] [PubMed]
  54. M. Vernier and V. Giguère, “Aging, senescence and mitochondria: The PGC-1/ERR axis,”. J. Mol. Endocrinol. 2021, 66, R1–R14. [CrossRef] [PubMed]
  55. C. H. Diep, N. J. Charles, C. B. Gilks, S. E. Kalloger, P. A. Argenta, and C. A. Lange, “Progesterone receptors induce FOXO1-dependent senescence in ovarian cancer cells,”. Cell Cycle 2013, 12, 1433–1449. [CrossRef]
  56. D. E. Harrison et al., “17-a-estradiol late in life extends lifespan in aging UM-HET3 male mice; nicotinamide riboside and three other drugs do not affect lifespan in either sex,”. Aging Cell 2021, 00, e13328. [CrossRef]
  57. H. N. Wilkinson and M. J. Hardman, “The Role of Estrogen in Cutaneous Ageing and Repair,”. Maturitas 2017, 103, 60–64. [CrossRef] [PubMed]
  58. S. J. Russell and C. R. Kahn, “Endocrine regulation of ageing,”. Nat. Rev. Mol. Cell Biol. 2007, 8, 681–691. [CrossRef] [PubMed]
  59. W. van den Beld, J. M. Kaufman, M. C. Zillikens, S. W. J. Lamberts, J. M. Egan, and A. J. van der Lely, “The physiology of endocrine systems with ageing,”. Lancet Diabetes Endocrinol. 2018, 6, 647–658. [CrossRef]
  60. X. Dai et al., “Estradiol-induced senescence of hypothalamic astrocytes contributes to aging-related reproductive function declines in female mice,”. Aging. 2020, 12, 6089–6108. [CrossRef]
  61. M. Garratt, K. A. Lagerborg, Y. M. Tsai, A. Galecki, M. Jain, and R. A. Miller, “Male lifespan extension with 17-α estradiol is linked to a sex-specific metabolomic response modulated by gonadal hormones in mice,”. Aging Cell 2018, 17, e12786. [CrossRef]
  62. C. H. Song, N. Kim, R. H. Nam, S. I. Choi, H. N. Lee, and Y. J. Surh, “17β-Estradiol supplementation changes gut microbiota diversity in intact and colorectal cancer-induced ICR male mice,”. Sci. Rep. 2020, 10, 12283. [CrossRef]
  63. Y. Fang et al., “Sexual dimorphic metabolic and cognitive responses of C57BL/6 mice to Fisetin or Dasatinib and quercetin cocktail oral treatment,”. GeroScience 2023. [CrossRef]
  64. L. Saemann, P. Naujoks, L. Hartrumpf, S. Pohl, A. Simm, and G. Szabó, “Sex-Specific Protection of Endothelial Function after Vascular Ischemia/Reperfusion Injury by the Senomorphic Agent Ruxolitinib,”. Int. J. Mol. Sci. 2023, 24, 11727. [CrossRef] [PubMed]
  65. National Institute on Aging, “Intervention Testing Program: Supported Interventions.,”. 2023. Available online: https://www.nia.nih.gov/research/dab/in- terventions-testing-program-itp/supported-interventions (accessed on 25 January 2024).
  66. N. Jiang and J. F. Nelson, “Sex Differences in Mouse Longevity and Responses to Geroprotective Drugs: Implications for Human Intervention,”. Public Policy Aging Rep. 2023, 33, 120–124. [CrossRef] [PubMed]
  67. P. Soldin and D. R. Mattison, “Sex Differences in Pharmacokinetics and Pharmacodynamics,”. Clin. Pharmacokinet. 2009, 48, 143. [CrossRef] [PubMed]
  68. Institute of Medicine (US) Committee on Understanding the Biology of Sex and Gender, T. M. Wizemann, and M.-L. Pardue, “Exploring the Biological Contributions to Human Health: Does Sex Matter?: Sex Begins in the Womb,” Washington (DC): National Academies Press (US). 2001. Available online: https://www.ncbi.nlm.nih.gov/books/NBK222286/.
  69. M. Oettel and A. K. Mukhopadhyay, “Progesterone: the forgotten hormone in men?,”. Aging Male 2004, 7, 236–257. [CrossRef] [PubMed]
  70. S. Hiller-Sturmhöfel and A. Bartke, “The endocrine system: an overview,”. Alcohol Heal. Res World 1998, 22, 153–164, [Online] Available: https://wwwncbinlmnihgov/pubmed/15706790.
  71. N. Fuentes and P. Silveyra, “Estrogen receptor signaling mechanisms,”. Adv. Protein Chem. Struct. Biol. 2019, 116, 135. [CrossRef]
  72. P. Ranganathan, N. Nadig, and S. Nambiar, “Non-canonical Estrogen Signaling in Endocrine Resistance,”. Front. Endocrinol. 2019, 10, 708. [CrossRef]
  73. P. Thomas, “Membrane Progesterone Receptors (mPRs, PAQRs): Review of Structural and Signaling Characteristics,”. Cells 2022, 11. [CrossRef]
  74. C. Slack, N. Alic, A. Foley, M. Cabecinha, M. P. Hoddinott, and L. Partridge Correspondence, “The Ras-Erk-ETS-Signaling Pathway Is a Drug Target for Longevity,”. Cell 2015, 162, 72–83. [CrossRef]
Preprints 97420 g001
Figure 1. Flowchart describing experimental design.
Figure 1. Flowchart describing experimental design.
Preprints 97420 g001
Preprints 97420 g002
Figure 2. Assessment of effects on senescent cell load using senescence-associated beta galactosidase (SAB) activity. A. Effect of treatments on SAB activity in cells at a late passage with higher levels of SAB activity. B. Effect of treatments on SAB activity in cells at an early passage with low levels of SAB activity. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests. (ns) = not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Figure 2. Assessment of effects on senescent cell load using senescence-associated beta galactosidase (SAB) activity. A. Effect of treatments on SAB activity in cells at a late passage with higher levels of SAB activity. B. Effect of treatments on SAB activity in cells at an early passage with low levels of SAB activity. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests. (ns) = not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Preprints 97420 g002
Preprints 97420 g003
Figure 3. Senescence kinetics for senescent male and female primary dermal fibroblasts. Percentage of cells stained for A. senescence-associated beta galactosidase (SAB), B. Ki67, a marker of proliferation and C. Terminal deoxynucleotidyl transferase dUTP nick end labelling (TUNEL), a marker of DNA damage, in female (F) and male (M) dermal fibroblast cells treated with synthetic female hormones at 10 μM or a DMSO-only control. Gene expression of markers for apoptosis, D. BCL2 and E. CASP3, in female (F) and male (M) dermal fibroblast cells. n = 3 for all groups. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests is reported: (ns) not significant, * p < 0.05, ** p< 0.01, *** p < 0.001 and **** p < 0.0001.
Figure 3. Senescence kinetics for senescent male and female primary dermal fibroblasts. Percentage of cells stained for A. senescence-associated beta galactosidase (SAB), B. Ki67, a marker of proliferation and C. Terminal deoxynucleotidyl transferase dUTP nick end labelling (TUNEL), a marker of DNA damage, in female (F) and male (M) dermal fibroblast cells treated with synthetic female hormones at 10 μM or a DMSO-only control. Gene expression of markers for apoptosis, D. BCL2 and E. CASP3, in female (F) and male (M) dermal fibroblast cells. n = 3 for all groups. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests is reported: (ns) not significant, * p < 0.05, ** p< 0.01, *** p < 0.001 and **** p < 0.0001.
Preprints 97420 g003
Preprints 97420 g004
Figure 4. Gene expression of genes encoding senescence-associated secretory phenotype (SASP) factors in female (F) and male (M) dermal fibroblast cells treated with synthetic female hormones at 10 μM or a DMSO-only control. Graph demonstrating the effect of synthetic female sex hormones on A. IL6, B. IL8, C. IL12A, D. CXCL1 and E. CXCL10 expression. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests is reported: (ns) not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Figure 4. Gene expression of genes encoding senescence-associated secretory phenotype (SASP) factors in female (F) and male (M) dermal fibroblast cells treated with synthetic female hormones at 10 μM or a DMSO-only control. Graph demonstrating the effect of synthetic female sex hormones on A. IL6, B. IL8, C. IL12A, D. CXCL1 and E. CXCL10 expression. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests is reported: (ns) not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Preprints 97420 g004
Preprints 97420 g005
Figure 5. Splicing factor expression following treatment with synthetic female sex hormones. Graph demonstrating the effect of synthetic female sex hormones on A. HNRNPK, B. NOVA1, C. PNISR, D. SRSF6 and E. TRA2B expression. n = 3 for all groups. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests is reported: (ns) not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Figure 5. Splicing factor expression following treatment with synthetic female sex hormones. Graph demonstrating the effect of synthetic female sex hormones on A. HNRNPK, B. NOVA1, C. PNISR, D. SRSF6 and E. TRA2B expression. n = 3 for all groups. Error bars show standard error of the mean (SEM), and statistical significance of p values computed using one-way ANOVA with uncorrected Fisher’s LSD post hoc tests is reported: (ns) not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Preprints 97420 g005
Preprints 97420 g006
Figure 6. Figure 6. Structure function analysis of compounds that decreased CDKN2A expression. A. Dendrogram constructed using the Tanimoto coefficient to show structural similarity of compounds tested that decreased CDKN2A gene expression. B. Maximum common substructure of the two least structurally similar compounds that decreased CDKN2A gene expression.
Figure 6. Figure 6. Structure function analysis of compounds that decreased CDKN2A expression. A. Dendrogram constructed using the Tanimoto coefficient to show structural similarity of compounds tested that decreased CDKN2A gene expression. B. Maximum common substructure of the two least structurally similar compounds that decreased CDKN2A gene expression.
Preprints 97420 g006
Table 1. Gene name and TaqManTM Gene Expression assay IDs used for characterisation experiments.
Table 1. Gene name and TaqManTM Gene Expression assay IDs used for characterisation experiments.
Preprints 97420 i001Preprints 97420 i002
Table 2. Fold change in CDKN2A (arbitrary units, relative to control) by compound and dose in the initial senescence screen. All effects listed here were more than three standard deviations above or below the mean of control treatments.
Table 2. Fold change in CDKN2A (arbitrary units, relative to control) by compound and dose in the initial senescence screen. All effects listed here were more than three standard deviations above or below the mean of control treatments.
Drug Name Dose (µM) Fold change in CDKN2A
Tucidinostat 10 2.048
Doxifluridine 10 1.559
Doxorubicin (hydrochloride) 10 1.498
Bromhexine (hydrochloride) 10 1.167
Homoharringtonine 10 1.160
Chlorambucil 10 1.133
Aspirin 10 1.072
Amoxapine 10 1.034
Doxorubicin (hydrochloride) 1 0.969
Imatinib 10 0.948
Montelukast (sodium) 10 0.888
Atorvastatin (hemicalcium salt) 10 0.822
Ribociclib 10 0.820
Baricitinib (phosphate) 10 0.820
Irinotecan (hydrochloride) 10 0.804
Levoleucovorin (calcium) 10 0.798
Epirubicin (hydrochloride) 10 0.790
Cobimetinib 10 0.773
Homoharringtonine 1 0.765
Decitabine 10 0.744
Sunitinib 10 0.722
Temozolomide 10 0.700
Silibinin 10 -0.686
Diacerein 10 -0.694
Vinorelbine (ditartrate) 1 -0.713
Alpelisib 10 -0.717
Ethamsylate 10 -0.734
Diethylstilboestrol 1 -0.753
Altretamine 10 -0.782
Panobinostat 1 -0.791
Sertraline (hydrochloride) 1 -0.805
Deferoxamine (mesylate) 10 -0.822
Balsalazide 1 -0.852
Pexidartinib 1 -0.890
Bexarotene 10 -0.894
Clofarabine 10 -0.897
Caffeic acid 10 -0.903
Pazopanib (hydrochloride) 10 -0.909
Aspirin 1 -0.916
Dexamethasone 1 -0.917
Pazopanib 10 -0.921
Rucaparib (phosphate) 10 -0.984
Glasdegib 1 -1.005
Aceglutamide 10 -1.020
Trimethoprim 10 -1.021
Crizotinib (hydrochloride) 10 -1.051
Acalabrutinib 1 -1.069
Zidovudine 10 -1.080
Citalopram (hydrobromide) 10 -1.094
Topotecan (hydrochloride) 10 -1.111
Rucaparib (phosphate) 1 -1.126
Alpelisib 1 -1.153
Sertraline (hydrochloride) 10 -1.154
Erlotinib 1 -1.157
Triclabendazole 10 -1.168
Nefopam (hydrochloride) 10 -1.174
Altretamine 1 -1.184
Bortezomib 1 -1.212
Nefopam (hydrochloride) 1 -1.217
Penfluridol 10 -1.230
Clioquinol 10 -1.241
Ethynyl estradiol 1 -1.259
Panobinostat 10 -1.260
Clofibrate 1 -1.272
Mizoribine 10 -1.291
Belinostat 10 -1.330
Valpromide 10 -1.351
Bosutinib 1 -1.354
Berberine (chloride hydrate) 10 -1.367
Nelarabine 1 -1.403
Acalabrutinib 10 -1.405
Tofacitinib (citrate) 10 -1.412
Erdosteine 1 -1.470
Bortezomib 10 -1.475
Bosutinib 10 -1.478
Osalmid 1 -1.493
Topotecan (hydrochloride) 1 -1.515
Bezafibrate 10 -1.523
Orotic acid 10 -1.532
Methylthiouracil 1 -1.551
Chlorpheniramine (maleate) 10 -1.559
Nitisinone 1 -1.561
Teniposide 10 -1.577
Sulfasalazine 10 -1.584
Pemetrexed (disodium hemipenta hydrate) 1 -1.702
Nifuroxazide 10 -1.705
Osalmid 10 -1.716
Nicotinamide 1 -1.717
Erlotinib 10 -1.741
Bendazol 1 -1.820
Bexarotene 1 -1.835
5-Azacytidine 1 -1.837
Nelarabine 10 -1.893
Clofarabine 1 -1.905
Niraparib 10 -1.927
Mycophenolic acid 10 -1.963
5-Azacytidine 10 -2.022
Chlorzoxazone 1 -2.045
Metyrapone 1 -2.066
Dimethyl fumarate 10 -2.099
Dexamethasone 10 -2.209
Dimethyl fumarate 1 -2.227
Chromocarb 10 -2.277
Penfluridol 1 -2.460
Bendazol 10 -2.486
Methylthiouracil 10 -2.527
Ethynyl estradiol 10 -2.684
Abemaciclib (methanesulfonate) 10 -2.768
Conivaptan (hydrochloride) 10 -2.908
Sunitinib 1 -2.926
Diethylstilbestrol 10 -3.068
Dronedarone 1 -4.099
Sodium 4-phenylbutyrate 10 -4.861
Cabozantinib 10 -7.875
Metyrapone 10 -8.532
Abemaciclib (methanesulfonate) 1 -11.417
Cabozantinib 1 -11.571
Carmofur 10 -11.805
Balsalazide 10 -11.887
Chlorzoxazone 10 -12.035
Table 3. Results from a screen for senescence-associated beta galactosidase (SAB) activity. The mean percentages of cells stained for SAB were compared against the appropriate experimental control for each batch of the screen. Assays 1-5 were performed on later passage fibroblasts to investigate potential reductions in senescence. Assays 6-7 were performed on earlier passage fibroblasts to investigate potential increases in senescence. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported: (ns) not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Table 3. Results from a screen for senescence-associated beta galactosidase (SAB) activity. The mean percentages of cells stained for SAB were compared against the appropriate experimental control for each batch of the screen. Assays 1-5 were performed on later passage fibroblasts to investigate potential reductions in senescence. Assays 6-7 were performed on earlier passage fibroblasts to investigate potential increases in senescence. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported: (ns) not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.
Treatment Mean SEM P Significance
Assay 1 Control 10µM 44.17 7.506 - -
5-Azacytidine 10µM 36.99 6.191 0.2496 ns
Caffeic Acid 10µM 31.67 3.34 0.0553 ns
Chlorpheniramine (maleate) 10µM 29.33 3.805 0.0264 *
Diethylstilboestrol 10µM 30.98 1.507 0.0445 *
Ethynyl estradiol 10µM 30.1 1.317 0.0337 *
Levonorgestrel 10µM 21.54 0.8541 0.002 **
Assay 2 Control 10µM 40.05 9.082 - -
Amoxapine 10µM 28.92 9.597 0.4353 ns
Bendazol 10µM 23.67 6.348 0.2568 ns
Citalopram (hydrobromide) 10µM 33.56 11.83 0.6466 ns
Methylthiouracil 10µM 33.69 12.09 0.6531 ns
Sertraline (hydrochloride) 10µM 26.81 10.88 0.3556 ns
Valpromide 10µM 23.84 7.242 0.2615 ns
Assay 3 Control 10µM 27.52 3.686 - -
Balsalazide 10µM 22.99 2.48 0.3251 ns
Carmofur 10µM 16.7 2.985 0.0288 *
Chlorzoxazone 10µM 19.91 2.438 0.109 ns
Conivaptan (hydrochloride) 10µM 22.33 2.988 0.2627 ns
Metyrapone 10µM 15.36 0.7593 0.0161 *
Sodium-4-Phenylbutyrate 10µM 16.16 4.995 0.0228 *
Assay 4 Control 1µM 19.19 4.546 - -
Abemaciclib (methanesulfonate) 1µM 12.09 1.888 0.1025 ns
Cabozantinib 1µM 8.13 0.2987 0.0159 *
Dronedarone 1µM 14.12 2.768 0.234 ns
Nicotinamide 1µM 12.11 1.795 0.1034 ns
Penfluridol 1µM 10.48 1.586 0.0495 *
Assay 4 Control 10µM 16.11 4.794 - -
Dexamethasone 10µM 10.26 2.553 0.1728 ns
Assay 5 Control 1µM 39.1 8.275 - -
Aspirin 1µM 32.61 2.587 0.3269 ns
Sunitinib 1µM 40.42 1.816 0.8398 ns
Assay 5 Control 10µM 36.06 3.345 - -
Aspirin 10µM 21.39 4.997 0.0396 *
Sunitinib 10µM 2.563 2.563 0.0002 ***
Assay 6 Control 1µM 3.163 0.5069 - -
Aspirin 1µM 4.287 0.5053 0.4963 ns
Sunitinib 1µM 3.493 1.016 0.2333 ns
Assay 6 Control 10µM 4.223 0.6868 - -
Aspirin 10µM 6.227 1.665 0.8404 ns
Sunitinib 10µM 0 0 0.0199 *
Imatinib 10µM 18.67 2.064 <0.0001 ****
Assay 7 Control 10µM 4.07 0.8632 - -
Bromhexine (hydrochloride) 10µM 5.65 1.818 0.3679 ns
Doxifluridine 10µM 6.04 1.637 0.2654 ns
Doxorubicin (hydrochloride) 10µM 15.84 0.1804 <0.0001 ****
Ethynyl estradiol 10µM 5.33 0.8184 0.4704 ns
Homoharringtonine 10µM 13.11 1.609 0.0001 ***
Tucidinostat 10µM 3.13 0.27 0.5886 ns
Table 4. Statistics of the percentage of cells stained for biomarkers in female (F) and male (M) dermal fibroblast cells treated with female synthetic hormones at 10 μM or a DMSO-only control. Biomarkers for senescence (senescence-associated beta galactosidase (SAB), proliferation (Ki67) and DNA damage (γH2AX and TUNEL)) are assessed. Although some cells stained for it, γH2AX staining was negligible across all experimental groups. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported. Significant p values > 0.05 are emboldened. n = 3 for all groups.
Table 4. Statistics of the percentage of cells stained for biomarkers in female (F) and male (M) dermal fibroblast cells treated with female synthetic hormones at 10 μM or a DMSO-only control. Biomarkers for senescence (senescence-associated beta galactosidase (SAB), proliferation (Ki67) and DNA damage (γH2AX and TUNEL)) are assessed. Although some cells stained for it, γH2AX staining was negligible across all experimental groups. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported. Significant p values > 0.05 are emboldened. n = 3 for all groups.
Preprints 97420 i003
Table 5. Gene expression data in female (F) dermal fibroblast cells treated with female synthetic hormones at 10 μM or a DMSO-only control. Genes relating to apoptosis, senescence, senescence-associated secretory phenotype (SASP) factors, splicing factors and spliceosomal components are assessed. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported. Significant p values > 0.05 are emboldened. n = 3 for all groups.
Table 5. Gene expression data in female (F) dermal fibroblast cells treated with female synthetic hormones at 10 μM or a DMSO-only control. Genes relating to apoptosis, senescence, senescence-associated secretory phenotype (SASP) factors, splicing factors and spliceosomal components are assessed. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported. Significant p values > 0.05 are emboldened. n = 3 for all groups.
Preprints 97420 i004Preprints 97420 i005
Table 6. Gene expression data in male (M) dermal fibroblast cells treated with female synthetic hormones at 10 μM or a DMSO-only control. Genes relating to apoptosis, senescence, senescence-associated secretory phenotype (SASP) factors, splicing factors and spliceosomal components are assessed. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported. Significant p values > 0.05 are emboldened. n = 3 for all groups.
Table 6. Gene expression data in male (M) dermal fibroblast cells treated with female synthetic hormones at 10 μM or a DMSO-only control. Genes relating to apoptosis, senescence, senescence-associated secretory phenotype (SASP) factors, splicing factors and spliceosomal components are assessed. The mean ± standard error of the mean (SEM) and p values from one-way ANOVAs with Fisher’s post hoc test are reported. Significant p values > 0.05 are emboldened. n = 3 for all groups.
Preprints 97420 i006Preprints 97420 i007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
bsky
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated