Repurposing statins as a potential ligand for estrogen receptor alpha via molecular docking

 

Khandelwal Alisha, Sharma Tripti

School of Pharmaceutical Sciences, Siksha ‘O’ Anusandhan Deemed to be University,

Bhubaneswar - 751003, Odisha, India.

*Corresponding Author E-mail: triptisharma@soa.ac.in

 

ABSTRACT:

Computational drug repurposing is the strategy for drug development which remarkably reduces the cost and development time. Research suggests that breast cancer development in women have been associated with cholesterol and its transporters. Cholesterol lowering drugs can be repurposed as potential therapeutic agents to prevent high cholesterol in estrogen receptor positive- breast cancer. The objective of this study was to carryout in-silico molecular docking of HMG-CoA reductase inhibitors (statins) with estrogen α receptor (3ERT) to repurpose the statins as breast cancer inhibitors.  Molecular docking studies were performed to explore the mechanism of interactions between the statins and human estrogen α receptor. Docking results revealed that statins bind to the hydrophobic pocket of the estrogen α receptor with high binding affinity. The docking scores were compared with the standard drug 4- hydroxy tamoxifen. The study helped to compare the interactions amongst different statins with the receptor and the energy values produced were ranging from -8.5 to -5.5 kcal/mol.  Molinspiration web servers was used to calculate the physiochemical properties and ADMET of the statins. Simvastatin showed better interaction amongst the docked statins with best protein ligand interactions, it was found to exhibit higher docking score of -8.5 kcal/mol. Therefore, we conclude that statins can be employed as an alternative drug for treatment of breast cancer.

 

KEYWORDS: Molecular docking, HMG-CoA reductase inhibitors, Estrogen α receptor, Breast cancer.

 

 


INTRODUCTION:

Breast cancer is the second most prevalent malignancy in women around the world. According to 2019 breast cancer statistics report more than 2,68,600 new cases of breast cancer were diagnosed in women1. Steroid hormone estrogen or 17β-estradiol is responsible for development of mammary gland and is one of the major risk factor for breast cancer. Estrogen receptor alpha (ERα) is target for ER-positive breast cancer and about 70% of breast cancer expresses ERα receptor2.

 

Statins reduces blood cholesterol by competitive inhibition of HMG-CoA reductase, the rate limiting step in the mevalonate pathway3. Apart from lipid lowering effects, statins also produce vasoprotective and cardioprotective effects4.

 

 

Statins are reported to reduce tumor cell growth and proliferation by inducing cell cycle arrest, inducing apoptosis and by negative impact on the tumor vasculature through suppression of angiogenesis5. Many clinical studies have demonstrated influence of statins on liver cancer6, colorectal7, lungs8,9, prostate10,11,12 and breast cancer13,14,15.

 

Drug repurposing is a drug discovery strategy for identification of a novel clinical use for an existing, approved or investigational drug. Since most of the drug interacts with multiple targets, these interactions could result in toxicities or side effects. The resulting side effects can be a new medical indication of the drug and used for treating other diseases16. The objective of this study was to evaluate the binding interactions of statins against ERα receptor through molecular docking.

 

MATERIAL AND METHODS:

Retrieval of target protein:

The X-ray crystal structure of human estrogen alpha receptor complex with 4-hydroxytamoxifen (4-HT) (PDBID:3ERT) at 1.9Ĺ resolution was retrieved from protein data bank (PDB) (https://www.rcsb.org/structure/3ert17. WHATIF server was used to repair the protein 3ERT. Hydrophobicity graphs were generated by Discovery Studio. Binding site analysis of 3ERT was performed using CASTp web server18.

 

Retrieval of ligands:

3D structures of eight statins Lovastatin (LS), Simvastatin (SS), Atorvastatin (AS), Pitavastatin (PTS), Cerivastatin (CS), Fluvastatin (FS), Rosuvastatin (RS), Pravastatin (PS) and standard drug 4-HT were retrieved from NCBI Pubchem compounds in SDF format and saved as mol file using Marvinsketch. (Table -1). ADMET properties of the selected compounds were predicted by molinspiration webserver.

 

Molecular docking:

Molecular docking of the ligands to the active binding site of the protein was performed using Autodock vina in PyRx virtual screening tool [19,20].  Protein pdb file (3ERT.pdb) was converted to pdbqt format file (3ERT.pdbqt) by removing water, adding hydrogen atom, kollman charges to the protein. The grid centre was positioned on the active site of the receptor 3ERT. The grid point for autogrid calculation was set to be X= 22.56, Y=-5.40, Z=21.93 with active site residue at the centre of the grid box. The Lamarckian Genetic Algorithm (LGA) was used to run docking simulations. Nine best poses are generated in the library after running vina wizard for selected 3ERT protein with each prepared ligand concurrently. The scoring, based on the docked energy of all the ligands were ranked version and the docking score obtained were then compared with 4-HT (Table 2). The ligands and protein were then submitted to swissdock for redocking. Protein ligand interactions were analyzed by Discovery Studio 4.5 and Chimera.

 

RESULTS AND DISCUSSIONS:

All the statins were found to strongly inhibit by occupying the active sites in the target protein 3ERT.


 

Figure 1: Protein ligand complex with target compound in binding site of Estrogen α receptor (3ERT) (a) Lovastatin (b) Simvastatin (c) Atrovastatin (d) Pitavastatin (e) Cerivastatin (f) Fluvastatin (g) Rosuvastatins (g) Pravastatin.


 

(e)

 

(d)

 
Table 1: 2D Structure of HMG-CoA reductase inhibitors (statins)

Sr. No

Compounds

Compound CID

2-D structure

1

Lovastatin

53232

 

2

Simvastatin

54454

 

3

Atorvastatin

60823

 

4

Pitavastatin

5282452

 

5

Cerivastatin

446156

 

6

Fluvastatin

446155

 

7

Rosuvastatin

446157

 

8

Pravastatin

54687

 

 

 

Table 2: List of docking interactions between the ligands 

Compound names

Docking score

(Pyrx)

(kcal mol-1)

Re docking

(swissdock)

(kcal mol-1)

Interacting residues

Lovastatin (LS)

-6.7

-8.49

Met 343, Leu 387, Met 388, Trp 383, Arg 394, Met 421, Leu 525.

Simvastatin (SS)

-8.5

-9.36

Leu 346, Met 388, Arg 394, Met 421, Ile 424, Leu 525.

Atrovastatin (AS)

-5.7

-7.89

Ile 326, Pro 406, Trp 393, Glu 322, Gly 390, Arg 394, Pro 324, Glu 442.

Pitavastatin (PTS)

-5.9

-7.71

Glu 323, Ile 326, Glu 353, Trp 393, Arg 394, Glu397, Pro 406, Phe 445, Lys 449.

Cerivastatin (CS)

-5.5

-7.32

Glu 323, Pro 324, Pro 325, Ile 326, Trp 393, Arg 394.

Fluvastatin (FS)

-6.9

-8.08

Met343, Ala 350, Glu353, Leu387, Leu 525

Rosuvastatins (RS)

-5.7

-7.26

  Glu 323, Glu 353, Arg 394, Trp 393, Lys 449, Pro, Ile 326, Glu 397

Pravastatin (PS)

-6.9

-8.47

 Leu 346, Met388, Arg394, Glu 353, Phe 404, Met 421, Leu 525

Tamoxifen

-9.0

-10.23

Ala45, Asp 46, Leu 49, Trp 78, Leu 220, Tyr 221, Lys 224, Cys 225, Val 228, Leu 231, Leu 234

 

Table 3: Drug likeness property (Lipinski’s rule of five) of statin

Compound names

Lipinski’s Properties

 

LV

MR

TPSA(Ĺ)

%ABS

HBD (≤10)

HBA

(≤10)

MW (g/Mol)

(≤ 500 g/mol)

logP

(≤5)

n

Lovastatin (LS)

1

5

404.55

4.34

7

0

113.92

72.84

83.87

Simvastatin (SS)

1

5

404.55

4.52

7

0

116.05

72.84

83.87

Atrovastatins (AS)

4

7

560.67

4.91

12

1

164.63

110.09

71.03

Pitavastatin (PTS)

3

5

421.47

3.91

8

0

117.82

90.65

77.72

Cerivastatin (CS)

3

6

459.56

4.45

11

0

127.82

99.88

74.54

Fluvastatin (FS)

2

4

409.50

4.44

8

0

119.66

62.46

87.45

Rosuvastatins (RS)

4

10

483.52

0.64

10

0

119.75

161.15

53.40

Pravastatin (PS)

4

7

424.53

2.15

11

0

114.04

124.29

66.12

Tamoxifen(4HT)

0

2

371.52

6.06

8

1

119.72

12.47

104.69

 


The 2-D sketch of the statins are drawn by using Marvin sketch (Table1)

 

The binding energy of the statins were in the range of -8.5 to 5.5kcal/mol and 4-HT taken as reference was -9 kcal/mol (Table 2).

 

Lipophillic statins exhibited absorption between 71.05 -87.45%, thus indicating better bioavailability as compared to hydrophilic statins. The no. of hydrogen bond donor and acceptor of all the prepared were ≤ 5 and 10 respectively. The MW of the ligand except AS were ≤ 500 and log P were ≤ 5 thus they are anticipated to have better absorption, distribution, permeation and metabolism. FS had higher % absorption 87.15 compared to other ligand. Drug likeliness properties of the statins are presented in Table3.

 

HBD: Hydrogen Bond Donor, HBA: Hydrogen Bond Acceptor, MW: molecular weight, n: No. of rotatable bonds, LV: Lipinski Violation, MR: Molar Refractivity, TPSA: Topological polar surface area, ABS: Absorption.

The binding pocket of 3ERT has key interacting amino acid residues Leu525, Gly521, Leu428, Met421, Gly420, Glu419, Phe404, Arg 394, Trp 383, Leu387, Glu353, Asp351, Ala350, Thr347, Leu346, Met343. Interaction of amino acids residues with the statin are described in Figure 2.


 


 


Figure 2: Investigation of binding mode of statins with Estrogen α receptor (3ERT) (a) Lovastatin (b) Simvastatin (c) Atrovastatin (d) Pitavastatin (e) Cerivastatin (f) Fluvastatin (g) Rosuvastatins (h) Pravastatin

 


Statins interacts with the amino acid residue of the protein by hydrogen bonding, vanderwaal interactions and π-alkyl interactions. Statins are reported to affect tumor cells directly by suppressing cancer proliferation, inducing apoptosis, antiangiogenic effects, anti-invasive and anti-metastaic effect. Lipophilic, high potency statins are found to be more effective at suppressing tumor cell growth21. This difference in antitumor activity of the statins may be accounted due to the difference in their physiochemical properties; lipophillic statins are more likely to cross the cellular membrane and hence more active. Statins are also found to enhance caspase 3 like activity and DNA fragmentation in MCF-7 cells leading to inhibition of proliferation of the cells22.

 

CONCLUSION:

Estrogen plays vital role in breast cancer development and estrogen receptor α is mainly responsible for breast cancer initiation and progression. There is a need for promising strategies for the design and synthesis of new therapeutic ligands which could selectively bind and inhibit estrogen dependent proliferation in estrogen α receptor. In this study we applied drug repurposing approach to find new ligands against estrogen α receptor, which can minimize the cost and risk associated in the development of new chemical entity. Lipophillic statins are found to be more active than hydrophilic statins.

 

According to the results obtained anti-cancerous drug 4-HT can be replaced by lipophilic statin SS due to least binding energy and higher absorption as compared to other statins. Lead modification of SS can be performed by bio-isosteric replacements and scaffold hopping. Further this needs proper investigation through experimental activity in an MCF-7 cell line to confirm the role of SS against breast cancer.

 

ACKNOWLEDGEMENT:

Authors are grateful to Dr. S.C. Si Dean SPS, for encouragements. Special thanks are due to honorable President, Prof. M.R. Nayak of the University for providing necessary facilities.

 

CONFLICT OF INTEREST:

The author(s) declares no conflicts of interest.

 

REFERENCE:

1.      American cancer society. About breast cancer. 2019. Available from: URL: https://www.cancer.org/cancer/breast-cancer.html

2.      Fishman J, Osborne MP, Telang NT. The role of estrogen in mammary carcinogenesis a. annals of the New York Academy of Sciences. 1995;768(1): 91-100.

3.      Endo A. A historical perspective on the discovery of statins. Proceedings of the Japan Academy, Series B. 2010;86(5): 484-93.

4.      Rosenson RS, Shott S, Lu L, Tangney CC. Hypertriglyceridemia and other factors associated with plasma viscosity. The American Journal of Medicine. 2001;110(6): 488-92.

5.      Beckwitt CH, Brufsky A, Oltvai ZN, Wells A. Statin drugs to reduce breast cancer recurrence and mortality. Breast Cancer Research. 2018 ; 20(1): 144.

6.      Shi M, Zheng H, Nie B, Gong W, Cui X. Statin use and risk of liver cancer: an update meta-analysis. BMJ Open. 2014; 4(9): e005399.

7.      Singh H, Mahmud SM, Turner D, Xue L, Demers AA, Bernstein CN. Long-term use of statins and risk of colorectal cancer: a population-based study. American Journal of Gastroenterology. 2009;104(12): 3015-23.

8.      Wang J, Li C, Tao H, Cheng Y, Han L, Li X, Hu Y. Statin use and risk of lung cancer: a meta-analysis of observational studies and randomized controlled trials. PloS one. 2013; 8(10).

9.      Tan M, Song X, Zhang G, Peng A, Li X, Li M, Liu Y, Wang C. Statins and the risk of lung cancer: a meta-analysis. PloS one. 2013; 8(2).

10.   Babcook MA, Joshi A, Montellano JA, Shankar E, Gupta S. Statin use in prostate cancer: an update. Nutrition and metabolic insights. 2016; 9: NMI-S38362.

11.   Murtola TJ, Tammela TL, Lahtela J, Auvinen A. Cholesterol-lowering drugs and prostate cancer risk: a population-based case-control study. Cancer Epidemiology and Prevention Biomarkers. 2007; 16(11): 2226-32.

12.   Yu O, Eberg M, Benayoun S, Aprikian A, Batist G, Suissa S, Azoulay L. Use of statins and the risk of death in patients with prostate cancer. Journal of Clinical Oncology. 2014;32(1): 5-11.

13.   Setoguchi S, Glynn RJ, Avorn J, Mogun H, Schneeweiss S. Statins and the risk of lung, breast, and colorectal cancer in the elderly. Circulation. 2007; 115(1):27.

14.   Ahern TP, Lash TL, Damkier P, Christiansen PM, Cronin-Fenton DP. Statins and breast cancer prognosis: evidence and opportunities. The Lancet Oncology. 2014; 15(10): e461-8.

15.   Undela K, Srikanth V, Bansal D. Statin use and risk of breast cancer: a meta-analysis of observational studies. Breast Cancer Research and Treatment. 2012;135(1): 261-9.

16.   Pushpakom S, Iorio F, Eyers PA, Escott KJ, Hopper S, Wells A, Doig A, Guilliams T, Latimer J, McNamee C, Norris A. Drug repurposing: progress, challenges and recommendations. Nature reviews Drug Discovery. 2019;18(1): 41-58.

17.   Williams CJ, Headd JJ, Moriarty NW, Prisant MG, Videau LL, Deis LN, Verma V, Keedy DA, Hintze BJ, Chen VB, Jain S. MolProbity: More and better reference data for improved all‐atom structure validation. Protein Science. 2018;27(1): 293-315.

18.   Tian W, Chen C, Lei X, Zhao J, Liang J. CASTp 3.0: computed atlas of surface topography of proteins. Nucleic Acids Research. 2018;46(W1): W363-7.

19.   Padmini R, Sitrarasi R, Razia M. Molecular docking studies of bioactive compounds from allium sativum against EML4-ALK receptor. Research Journal of Pharmacy and Technology. 2017; 10(11): 3741-7.

20.   Ramjith US, Muhammed S. Molecular docking study of novel Imidazo [2, 1-b]-1, 3, 4 thiadiazole derivatives. Research Journal of Pharmacy and Technology. 2013;6(6): 688-94

21.   Menter DG, Ramsauer VP, Harirforoosh S, Chakraborty K, Yang P, Hsi L, Newman RA, Krishnan K. Differential effects of pravastatin and simvastatin on the growth of tumor cells from different organ sites. PloS One. 2011; 6(12).

22.   Kotamraju S, Willams CL, Kalyanaraman B. Statin-induced breast cancer cell death: role of inducible nitric oxide and arginase-dependent pathways. Cancer Research. 2007;67(15): 7386-94.

 

 

 

 

 

 

 

Received on 12.04.2020           Modified on 08.06.2020

Accepted on 01.07.2020         © RJPT All right reserved

Research J. Pharm. and Tech. 2021; 14(7):3757-3762.

DOI: 10.52711/0974-360X.2021.00650