- Open Access
Large-scale in silico identification of drugs exerting sex-specific effects in the heart
Journal of Translational Medicine volume 16, Article number: 236 (2018)
Major differences exist between men and women in both physiology and pathophysiology. Dissecting the underlying processes and contributing mechanisms of sex differences in health and disease represents a crucial step towards precision medicine. Considering the significant differences between men and women in the response to pharmacotherapies, our aim was to develop an in silico model able to predict sex-specific drug responses in a large-scale.
For this purpose, we focused on cardiovascular effects because of their high morbidity and mortality. Our model predicted several drugs (including acebutolol and tacrine) with significant differences in the heart between men and women. To validate the sex-specific drug responses identified by our model, acebutolol was selected to lower blood pressure in spontaneous hypertensive rats (SHR), tacrine was used to assess cardiac injury in mice and metformin as control for a non-sex-specific response.
As our model predicted, acebutolol exhibited a stronger decrease in heart rate and blood pressure in female than male SHRs. Tacrine lowered heart rate in male but not in female mice, induced higher plasma cTNI level and increased cardiac superoxide (DHE staining) generation in female than male mice, indicating stronger cardiac toxicity in female than male mice. To validate our model in humans, we employed two Chinese cohorts, which showed that among patients taking a beta-receptor blocker (metoprolol), women reached significantly lower diastolic blood pressure than men.
We conclude that our in silico model could be translated into clinical practice to predict sex-specific drug responses, thereby contributing towards a more appropriate medical care for both men and women.
There are significant differences between men and women in physiology, in the development and outcome of disease, as well as in the response to pharmacotherapies [1,2,3]. However, the role and effects of biological sex have been poorly investigated on the basis of the assumption that men and women respond similarly to disease and therapeutic interventions [3, 4]. Dissecting the differences between male and female individuals in health and disease is being recognized as critically important for improving the efficiency, efficacy, accuracy and precision of medical care [3, 5].
Among others, genetic mechanisms are expected to contribute to sex differences in physiology, pathophysiology and drug responses [3, 6]. In the post-genomic era, the fast development of high-throughput techniques has enabled the identification of a large number of genes that are under sex-specific regulation in several tissues . For example, using transcriptomic approaches, projects, such as the Genotype-Tissue Expression (GTEx), have identified hundreds of genes that are differentially expressed between the sexes in various tissues [8, 9]. These studies have made important contributions towards the understanding of the biology of sex differences, but methods to infer sex-specific medicine based on functional genomic datasets are still limited.
Public databases [10, 11] have accumulated enormous profusion of functional genomic datasets, which may facilitate the research and development of novel and innovative in silico models for the prediction of sex-specific drug responses. On the basis of this, we propose here an in silico model, SexBiasedDrug, which infers drugs that exert different effects between male and female individuals. For this purpose, we focused on cardiovascular effects because of their high morbidity and mortality, based on acquired big functional genomic data responses to many drugs. SexBiasedDrug first calculates the male-biased genes (MGs) and female-biased genes (FGs) based on gene expression profiles from the public GEO database  and then determines the Connectivity Map  drugs that show significant associations to MGs and FGs. These resulting drugs are the ones exerting sex-specific effects. Finally, we confirmed the accuracy of our model using experimental animals treated with acebutolol, tacrine and metformin and show its clinical relevance employing two Chinese cohorts.
Generation of the in silico model
We used the normalized human cardiac gene expression profiles (GEO Accession Number: GSE57345-GPL11532 ) from the GEO database . In the original study, Liu et al., studied the changes in cardiac gene expression between failing and non-failing hearts. We previously presented a method to identify sex-biased genes from public gene expression datasets . Here, we used this approach to identify sex-biased genes in the healthy heart. For this purpose, we only extracted the 135 gene expression profiles from non-failing (normal control) heart samples and classified them into two groups (male group and female group) based on the sex annotation. There was no difference in the age between the male and female groups (data not shown, Wilcoxon test). Next, we compared the expression level between the two groups for each gene using multiple t tests. The male/female fold change was calculated. The P values were corrected by FDR. We considered genes with FDR ≤ 0.05 to be sex-biased.
We used gene expression signatures following treatment with more than 1000 drugs from the Connectivity Map , which is a database with thousands of gene expression profiles for various drug treatments and was originally used for drug-repurposing. Here, we used the data in the Connectivity Map for predicting sex-biased responsive drugs. For a given drug treatment, genes with fold change ≥ 2.0 and ≤ − 2.0 between treatment and control were considered upregulated and downregulated, respectively. Next, we evaluated the association of sex-biased genes and a given drug by the Fisher’s exact test, which was performed based on two kinds of classifications, i.e. sex-biased gene and drug-regulated gene. Additional file 1 lists the top 20 drugs that show the most significant sex-specific effects in the heart. We developed a computer script using R and Java, SexBiasedDrug, for this purpose. Finally, it should be noted that a bigger/smaller fold change does not always mean more/less significant. Therefore, methods to identify better gene signatures could be proposed in the future.
All animal procedures complied with the Animal Management Rule of the Ministry of Health, People’s Republic of China (Document No. 55, 2001) and the Care and Use of Laboratory Animals published by the US National Institutes of Health (NIH Publication No. 85-23, updated 2011). The care and use of laboratory animals was approved by the Laboratory Animal Ethics Committee of Peking University. Male and female Sprague–Dawley (SD) rats (200–220 g), Spontaneously Hypertensive Rats (SHR) (9 weeks old), C57BL/6J mice (20–25 g) were purchased from the Animal Center of Peking University Health Science Center. Animals were kept under temperature 23 ± 2 °C, humidity range of 40–70% and 12-h light/dark cycles. Food and water were given ad libitum. Acebutolol hydrochloride, isoproterenol, tacrine, metformin, and dihydroethidium (DHE) were obtained from Sigma-Aldrich. Other chemicals and reagents were of analytical grade.
For the animal experiments with acebutolol, eight male and eight female SHRs were treated with acebutolol (50 mg/kg/day, gavage), control groups were given gavage of equal volume saline. The changes of blood pressure were recorded via a noninvasive tail-cuff apparatus (Softron, Japan). Blood pressure was calculated as the average value of three successive measurements. Animals were anesthetized with isoflurane and then pentobarbital sodium (50 mg/kg) was given by intraperitoneal injection. Heart rate changes were recorded by electrocardiogram (ECG) transduction sensor connecting with Powerlab 4S (Adinstrument, Australia).
For tacrine and metformin treatments, 10 male and 10 female C57BL/6J mice were administered tacrine (5 mg/kg/day, gavage) for 4 weeks; or treated with metformin (30 mg/kg/day, gavage). Control mice were given equal volume of saline. Blood samples were collected from the angular artery. Serum cardiac troponin I (cTNI) was measured by an ELISA kit (R&D System Inc.) according to the manufacturer’s protocol. Mice were sacrificed and the heart was quickly removed. Frozen cardiac sections (6 μm) were stained with dihydroethidium (DHE, 5 μM) at room temperature for 60 min and images were acquired with a fluorescence microscope (90i, Nikon, Japan).
Consecutive patients were recruited at two Chinese sites. One diabetes cohort was obtained from Luhe Hospital at Beijing city with 35 patients receiving metoprolol (9 men and 26 women). The other cohort was from the Jidong community for health management at Tangshan city with 153 patients (86 men and 67 women) receiving metoprolol. The basal clinical characteristics are shown in Tables 1 and 2, respectively. The investigation conforms with the principles outlined in the Declaration of Helsinki and patient written informed consent was obtained.
All data are expressed as mean ± SD. Differences among groups were analyzed by two-way ANOVA (multiple groups) or two-tailed t test (two groups). The significance level was set at 0.05.
An in silico model for inferring drugs that exert sex-specific effects
We have developed an in silico pipeline, SexBiasedDrug, which predicts drugs exerting different effects between male and female individuals. For this purpose, we focused on public datasets (GEO: GSE57345-GPL11532)  of non-failing and failing human hearts (Fig. 1a). We analyzed 135 gene expression profiles from non-failing cardiac samples for sex-specific gene expression (72 male vs. 63 female samples). There was no significant age difference between the groups (data not shown). We identified 134 sex-biased genes (Additional file 2); 79 male-biased genes (MGs, which are higher in males than in females) and 55 female-biased genes (FGs, which are higher in females than in males). We next determined the associations between the sex-biased genes (MGs and FGs) and the deregulated genes (upregulated genes, UGs, and downregulated genes, DGs) induced by drugs in the Connectivity Map database one by one using Fisher’s exact test (Fig. 1a). Additional file 1 shows the top 20 drugs that are predicted to exert significantly different effects between men and women (Fig. 1b).
Experimental validation of the in silico model
To confirm the accuracy of our in silico pipeline, SexBiasedDrug, we selected two positive predictions for sex-specific effects, i.e. acebutolol and tacrine, and one negative prediction for sex-specific effects, i.e. metformin. Acebutolol is a selective β1 receptor blocker used mainly for the treatment of hypertension and arrhythmias . As our in silico model predicted, following a 4-week treatment of spontaneously hypertensive rats (SHR) with acebutolol, the heart rate was significantly lower in female than male SHRs (Fig. 2a). Similarly, systolic blood pressure (SBP) and diastolic blood pressure (DBP) were significantly lower in female than male SHR treated with acebutolol (Fig. 2b). These data together confirmed our in silico prediction that acebutolol exerts sex-specific effects.
Tacrine is a cholinesterase inhibitor, which was aimed for the treatment of Alzheimer’s disease but has been discontinued from the US market . Given that tacrine could exert cardiac-toxic effects leading to bradycardia, we investigated heart rate and cardiac toxicity after tacrine treatment, using metformin as negative control for the prediction of sex-specific effects. In agreement with our in silico predictions, after 3 weeks of treatment of healthy mice, tacrine significantly decreased heart rate in male mice but not in female mice (Fig. 2c). In contrast, metformin did not affect the heart rate in both male and female mice (Fig. 2c). We further found that plasma cardiac troponin I (cTNI) levels (Fig. 2d) and cardiomyocyte superoxide generation, as assessed by dihydroethidium (DHE) staining (Fig. 2e), were induced to a significantly higher extent in female than male mice (Additional file 3: Figure S1). Metformin also slightly increased cardiomyocyte superoxide level, but there was no significant difference between male and female animals (Fig. 2e, Additional file 3: Figure S1). The changes in cTNI levels and superoxide generation indicate that tacrine exerted more severe toxicity in female than male mice. These data together further support the accuracy of our in silico model for the prediction of sex-specific drug effects.
Clinical validation of the in silico model
To demonstrate the clinical relevance of our model, we performed a retrospective study comparing the blood pressure changes between male and female patients receiving metoprolol in two Chinese communities. In China, acebutolol had been little used because of its side-effects and it was replaced with metoprolol. The chemical structure of metoprolol is close to that of acebutolol in chemical structure (Additional file 4: Figure S2), and metoprolol has a closest similar decreasing HR effect on heart rate with acebutolol . Here, we performed a retrospective study to compare the blood pressure changes taken metoprolol between male and female from two Chinese northern communities. In a diabetes cohort from Beijing city, with a total of 35 patients who receiving metoprolol were recruited (basal clinical data see Table 1), female patients exhibited lower DBP than male patients (Table 1, P = 0.001) although with larger BMI (P = 0.035). Our data showed that metoprolol treatment had more anti-hypertension effect in obese females than lean males, these differences showed the stronger evidence for sex-response of metoprolol to support our in silico model. In the second cohort with 153 patients (86 males vs. 67 female patients) who were from the Jidong community at Tangshan city receiving metoprolol, the DBP decrease by the drug was more prominent in the female than male patients (Table 2, P = 0.0004). Together, these data demonstrate the accuracy and clinical relevance of our prediction model.
Using accumulating big datasets in public databases, we have developed a computational model for the prediction of drugs exerting sex-specific effects. Dissecting the differences between the sexes in health and disease is crucial for the improvement of the efficiency, efficacy, accuracy and precision of medical care. However, this has not received much attention in both medical care and clinical research. Thousands of drugs have been approved , but systematic evaluation for sex-specific effects has been lacking. Computational approaches and frameworks for these purposes are still not readily available. The model presented here addresses this major gap. The accuracy of this model was confirmed by animal experimentation with two positive predictions for sex-specific effects, i.e. acebutolol and tacrine, and one negative prediction for sex-specific effects, i.e. metformin. Importantly, the clinical relevance of our model was demonstrated in studies with male and female patients receiving metoprolol, which showed that a significantly lower DBP was achieved in female than male patients.
In fact, in our attempt to validate the model in humans, we tried to employ publicly available clinical data to support the predictions. Sutandar et al.  reported the acebutolol plus hydrochlorothiazide effects for the treatment of essential hypertension. Since detailed data were provided in the study for 11 male and 11 female patients with acebutolol treatment, we performed a meta-analysis for any differences between men and women, but we did not find significant effects in SBP (Additional file 5: Figure S3a), DBP (Additional file 5: Figure S3b) or HR (Additional file 5: Figure S3c). Analysis of absolute heart rate values revealed significant effects of acebutolol in both sexes. However, this study did not exclude the interference of combination with other drugs, such as diuretics, and the sample size was rather to small significant sex differences. To this extent, in multicenter clinical trials with acebutolol alone for the treatment of essential hypertension, SBP and DBP were lowered effectively at rest but not exercise [17,18,19,20]. Another multicenter clinical trial reported significant effects of acebutolol and confirming its anti-hypertensive actions . However, none of these studies compared the responses to acebutolol between men and women. For our clinical validation, we chose metoprolol, due to its similar structure to acebutolol and that its use has been limited in China because of its side-effects. In first cohort, obese female patients got the more lowering blood pressure effect than lean male patients (overweight and obesity is a risk factor for cardiovascular diseases and diabetes), the great differences supported the sex-response of metoprolol. The second cohort’s study also showed the same sex-response of metoprolol in age, gender, BMI and biochemical parameters matched patients. Our data show for the first time that the response to this drug differs significantly between men and women, which is also a strong clinical confirmation for our model.
Our study was limited by scarce sex-specific genomic and transcriptomic data, especially, data on sex-specific responses to first-line drugs. This lack of data limited our prediction model in assessing sex-specific responses to further drugs. Another limitation is that gene expression profile based approaches do not address direct protein targets of the interested drugs. The expression profile-derived “gene signature” represent one type of “molecular phenotype” but not the drug targets. Although this limitation existed, this class of approaches has been successfully applied in a number of studies in recent years. In addition, genetic or genomic variations in drug target DNA regions could lead to changing in drug effects. Therefore, integrating GWAS data and eQTL data into this study could be useful to improve the prediction accuracy.
In summary, we developed an in silico model for the prediction of drugs exerting sex-specific effects. On the basis of this, possible responses or therapeutic effects of drugs to male and female patients using our prediction model may easily be evaluated, to avoid side-effects or adjust treatments accordingly. We put forward that our model may contribute to more appropriate and personalized pharmacological intervention approaches.
Gaignebet L, Kararigas G. En route to precision medicine through the integration of biological sex into pharmacogenomics. Clin Sci. 2017;131:329–42.
Pagidipati NJ, Peterson ED. Acute coronary syndromes in women and men. Nat Rev Cardiol. 2016;13:471–80.
Regitz-Zagrosek V, Kararigas G. Mechanistic pathways of sex differences in cardiovascular disease. Physiol Rev. 2017;97:1–37.
Morselli E, Frank AP, Santos RS, Fatima LA, Palmer BF, Clegg DJ. Sex and gender: critical variables in pre-clinical and clinical medical research. Cell Metab. 2016;24:203–9.
Kararigas G, Seeland U, Barcena de Arellano ML, Dworatzek E, Regitz-Zagrosek V. Why the study of the effects of biological sex is important. Commentary. Ann Ist Super Sanita. 2016;52:149–50.
Gilks WP, Abbott JK, Morrow EH. Sex differences in disease genetics: evidence, evolution, and detection. Trends Genet. 2014;30:453–63.
Ellegren H, Parsch J. The evolution of sex-biased genes and sex-biased gene expression. Nat Rev Genet. 2007;8:689–98.
Mele M, Ferreira PG, Reverter F, DeLuca DS, Monlong J, Sammeth M, Young TR, Goldmann JM, Pervouchine DD, Sullivan TJ, et al. Human genomics. The human transcriptome across tissues and individuals. Science. 2015;348:660–5.
Trabzuni D, Ramasamy A, Imran S, Walker R, Smith C, Weale ME, Hardy J, Ryten M, North American Brain Expression C. Widespread sex differences in gene expression and splicing in the adult human brain. Nat Commun. 2013;4:2771.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991–5.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, et al. The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313:1929–35.
Liu Y, Morley M, Brandimarto J, Hannenhalli S, Hu Y, Ashley EA, Tang WH, Moravec CS, Margulies KB, Cappola TP, et al. RNA-seq identifies novel myocardial gene expression signatures of heart failure. Genomics. 2015;105:83–9.
Guo S, Zhou Y, Zeng P, Xu G, Wang G, Cui Q. Identification and analysis of the human sex-biased genes. Brief Bioinform. 2016;19:188–98.
Law V, Knox C, Djoumbou Y, Jewison T, Guo AC, Liu Y, Maciejewski A, Arndt D, Wilson M, Neveu V, et al. DrugBank 4.0: shedding new light on drug metabolism. Nucleic Acids Res. 2014;42:D1091–7.
Newman RJ. Comparison of propranolol, metoprolol, and acebutolol on insulin-induced hypoglycaemia. Br Med J. 1976;2:447–9.
Sutandar H, Sugeng I, Utama H, Suryaatmaja M, Sukaman. The efficacy and acceptability of the combination of acebutolol and hydrochlorothiazide in the treatment of essential hypertension. Curr Med Res Opin. 1984;9:323–8.
Ashton WL. An open, multicentre study of acebutolol, given as a single daily dose, in the management of hypertension. Curr Med Res Opin. 1977;5:279–83.
Baker PG, Goulton J. A multicentre study of twice daily acebutolol (‘Sectral’) in the treatment of hypertension in general practice. Curr Med Res Opin. 1979;6:50–61.
DiBianco R, Singh SN, Shah PM, Newton GC, Miller RR, Nahormek P, Costello RB, Laddu AR, Gottdiener JS, Fletcher RD. Comparison of the antianginal efficacy of acebutolol and propranolol. A multicenter, randomized, double-blind placebo-controlled study. Circulation. 1982;65:1119–28.
DiBianco R, Singh S, Singh JB, Katz RJ, Bortz R, Gottdiener JS, Spodick DH, Laddu AR, Fletcher RD. Effects of acebutolol on chronic stable angina pectoris. A placebo-controlled, double-blind, randomized crossover study. Circulation. 1980;62:1179–87.
Ashton WL. An open, multicentre study of acebutolol in hypertension. Curr Med Res Opin. 1976;4:442–54.
QC and BG designed and supervised this project. CH implemented the predictive model. CC performed the animal experiments. KL, BH, YZ, and YF analyzed clinical data. QC, BG, GK, JY, and GX wrote the paper. All authors read and approved the final manuscript.
We thank members in Qinghua Cui’s lab and Bing Geng’s lab for valuable discussions.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The care and use of laboratory animals was approved by the Laboratory Animal Ethics Committee of Peking University.
This work was supported by the Special Project on Precision Medicine under the National Key R&D Program (2016YFC0903003), the National Natural Science Foundation of China (81670462; 81422006) and the Recovery Medical Science Foundation.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.