Research Article |
Corresponding author: Vivekanda A. Chatpalliwar ( vchatpalliwar@yahoo.co.in ) Corresponding author: Mithun Rudrapal ( rsmrpal@gmail.com ) Corresponding author: Johra Khan ( j.khan@mu.edu.sa ) Academic editor: Alexander Zlatkov
© 2023 Bhagyashri J. Warude, Sandip N. Wagh, Vivekanda A. Chatpalliwar, Merve Yildirim, Ismail Celik, Mithun Rudrapal, Johra Khan, Sampath Chinnam, Aniket A. Garud, Vishnu S. Neharkar.
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
Warude BJ, Wagh SN, Chatpalliwar VA, Yildirim M, Celik I, Rudrapal M, Khan J, Chinnam S, Garud AA, Neharkar VS (2023) Design, docking, MD simulation and in-silco ADMET prediction studies of novel indole-based benzamides targeting estrogen receptor alfa positive for effective breast cancer therapy. Pharmacia 70(2): 307-316. https://doi.org/10.3897/pharmacia.70.e100356
|
Breast cancer is one of the most common malignancies in women, afflicting millions of lives each year. Our current study suggests that the development of the most promising 7-substituted -1-(4-(piperidine-1-yl methoxy)benzyl)-1H-indole-3-carboxamide derivatives results in potent anticancer agents through in-silico investigations. The molecular docking was performed against estrogen receptor alpha (ER-α) positive (PDB ID: 3UUD) of breast cancer cells to anticipate the binding modes of the designed compounds and the likely mode of action. The interactions between the ligands and amino acid residues were thoroughly elucidated. The stability of the docked protein-ligand complexes was further confirmed by 100 ns molecular simulations methods. From in-silico studies, indole-based benzamides exhibited satisfactory physicochemical, drug-likeness and toxicity properties. To conclude, the most promising substituted benzamide analogs on the indole ring could serve as a possible modulator against ER-α positive breast cancer.
breast cancer, estrogen receptor alpha, indole scaffold, benzamide, bazedoxifene, docking, molecular dynamics, swiss admet
Breast cancer (BC) affects women globally at any age after puberty with increasing incidence in the future. Human breast cancer is the second largest cause of death in women. In 2020, there was 2.3 million women diagnosed with BC, with 685 000 deaths worldwide. As of the end of 2020, there were 7.8 million active instances of BC in women over the previous five years. Around 50% of BC develops in women due to BC risk factors other than gender, such as being female or being over the age of 40. Obesity, radiation exposure, excessive alcohol and tobacco consumption, reproductive disorders, and a family history of BC are all risk factors for BC (
The work presented here is based on the structure-based drug design (
Bazedoxifene (BSD, 1H-indo-5-ol, 1-[[4-[2(hexahydro-1H-azepin-1-yl)ethoxy]methyl] acetic acid, 2-(-4-hydroxyphenlyl)-3-methyl) is an indole derivative and third-generation SERM, which acts as an estrogen receptor antagonist in breast cancer (
In designing, the scaffold of indole-based analogues involves substitution at the 7th position with benzamide, 3rd position with amide group and 3՛ position with different alkyl group with ester and alkyl halide functional group and 1st position substituted with 1-((4ethylphenoxy)methyl)piperidine and 1-((4ethylphenoxy)ethyl)piperidine is mentioned in Table
The molecular docking software, AutoDock Vina (Virtual screening tool) was employed for the docking study, and the Biovia Discovery Studio visualizer was used to study the 2D and 3D interactions of the ligand-receptor complex after docking. Molecular docking analyses were performed via the CB-Dock server (http://clab.labshare.cn/cb-dock/php/) (
The molecular dynamics (MD) simulation study was performed with Gromacs v2021.2 (
The ADMET data of ligands were evaluated on pkCSM (http://biosig.unimelb.edu.au/pkcsm/prediction/) and SwissADME (http://www.swissadme.ch/index.php/) software (
Molecular docking is a computational tool used to monitor the formation of a stable protein-ligand complex between the active site of the protein and ligand molecule (
AutoDock Vina molecular docking interaction energies, parameters, and contact residues of the compounds B73aIII, B73aV, and cocrystal ligand estradiol (EST) with ER (PDB ID: 3UUD).
Compounds | Vina score | Cavity volume (Å3) | Center (x, y, z) | Docking size (x, y, z) | Contact residues | |
---|---|---|---|---|---|---|
H Bonds | Hydrophobic | |||||
B73aIII | -8.7 | 1383 | 16, 5, 7 | 27, 27, 27 | Trp393 (2.87 Å), Arg394 (3.09 Å) | Glu323, Pro324, Pro325, Ile326, Glu353, His356, Met357, Ile386, Leu387, Gly390, Leu391, Glu397, His398, Leu403, Phe445, Lys449 |
B73aV | -8.0 | 322 | -3, 14, -5 | 29, 29, 29 | Trp393 (2.97 Å), Arg394 (3.36 Å) | Leu320, Glu323, Pro324, Pro325, Ile326, Glu353, His356, Met357, Trp360, Ile386, Leu387, Gly390, Leu391, Trp393, Arg394, Gly442, Glu443, Phe445, Lys449 |
EST | -10.9 | 1383 | 16, 5, 7 | 27, 20, 20 | Arg394 (3.19 Å), His524 (2.95Å) | Met343, Leu346, Thr347, Leu349, Ala350, Glu353, Leu384, Leu387, Met388, Gly390, Leu391, Phe404, Met421, Ile424, Leu428, Gly521, Leu525 |
Fig.
Amino acids located in the active site of the ER interact with the ligand in different types of bonds. It has been shown that an amino acid can form different bonds with different parts of the ligand. Conventional hydrogen bond, van der Waals, carbon-hydrogen bond, π-anion, π-donor hydrogen bond, π-sigma, alkyl, π-alkyl bonds were observed. On the contrary, it was observed that an unfavorable hydrophobic bond was formed between the Ile356 amino acid of receptor molecule and the ligand B73aIII. In Fig.
EST, B73aIII, and B73aV molecules were performed in MD simulations for examination of protein-ligand complex stability (
In ADMET studies, the molecular weights of ligands, LogP values, the characteristics of drug-likeness, the bioavailability scores, absorption, distribution, metabolism, excretion, and toxicity properties were evaluated (
Descriptive properties of ligands examined in pkCSM and SwissADME software and estimated values of drug-likeness.
Descriptor | Value | ||
---|---|---|---|
B73aIII | B73aV | EST | |
Molecular weight | 531.056 | 554.647 | 272.388 |
LogP | 5.6018 | 4.6921 | 3.2651 |
Rotatable bonds | 9 | 10 | 0 |
Acceptors | 5 | 7 | 2 |
Donors | 2 | 2 | 1 |
Surface area | 226.679 | 238.381 | 120.496 |
Drug-likeness | B73aIII | B73aV | EST |
Lipinski | Yes; 1 violation: MW>500 | Yes; 1 violation: MW>500 | Yes; 0 violation |
Ghose | No; 2 violations: MW>480, MR>130 | No; 3 violations: MW>480, MR>130, #atoms>70 | Yes |
Veber | Yes | No; 1 violation: Rotors>10 | Yes |
Egan | Yes | Yes | Yes |
Muegge | Yes | Yes | Yes |
Bioavailability score | 0.55 | 0.55 | 0.55 |
Table
Property | Model name | Predicted Value | Unit | ||
---|---|---|---|---|---|
B73aIII | B73aV | EST | |||
Absorption | Water solubility | -5.332 | -5.065 | -4.356 | Numeric (log mol/L) |
Absorption | Caco-2 permeability | 0.987 | 0.71 | 1.513 | Numeric (log Papp in 10-6 cm/s) |
Absorption | Intestinal absorption (human) | 89.409 | 87.114 | 97.263 | Numeric (% Absorbed) |
Absorption | Skin permeability | -2.785 | -2.864 | -3 | Numeric (log Kp) |
Absorption | P-glycoprotein substrate | Yes | Yes | No | Categorical (Yes/No) |
Absorption | P-glycoprotein I inhibitor | Yes | Yes | Yes | Categorical (Yes/No) |
Absorption | P-glycoprotein II inhibitor | Yes | Yes | No | Categorical (Yes/No) |
Distribution | VDss (human) | 0.108 | -0.17 | 0.403 | Numeric (log L/kg) |
Distribution | Fraction unbound (human) | 0.01 | 0.059 | 0.185 | Numeric (Fu) |
Distribution | BBB permeability | -0.826 | -1.024 | 0.128 | Numeric (log BB) |
Distribution | CNS permeability | -1.993 | -2.413 | -2.232 | Numeric (log PS) |
Metabolism | CYP2D6 substrate | No | No | No | Categorical (Yes/No) |
Metabolism | CYP3A4 substrate | Yes | Yes | Yes | Categorical (Yes/No) |
Metabolism | CYP1A2 inhibitor | No | No | No | Categorical (Yes/No) |
Metabolism | CYP2C19 inhibitor | No | No | Yes | Categorical (Yes/No) |
Metabolism | CYP2C9 inhibitor | Yes | No | No | Categorical (Yes/No) |
Metabolism | CYP2D6 inhibitor | No | No | No | Categorical (Yes/No) |
Metabolism | CYP3A4 inhibitor | Yes | Yes | No | Categorical (Yes/No) |
Excretion | Total Clearance | 1.172 | 0.991 | 1.025 | Numeric (log ml/min/kg) |
Excretion | Renal OCT2 substrate | No | No | Yes | Categorical (Yes/No) |
Toxicity | AMES Toxicity | No | No | No | Categorical (Yes/No) |
Toxicity | Max. Tolerable dose (human) | 0.467 | 0.313 | -0.723 | Numeric (log mg/kg/day) |
Toxicity | hERG I inhibitor | No | No | No | Categorical (Yes/No) |
Toxicity | hERG II inhibitor | Yes | Yes | Yes | Categorical (Yes/No) |
Toxicity | Oral Rat Acute Toxicity (LD50) | 2.466 | 2.542 | 1.779 | Numeric (mol/kg) |
Toxicity | Oral Rat Chronic Toxicity (LOAEL) | 2.214 | 2.323 | 1.893 | Numeric (log mg/kg_bw/day) |
Toxicity | Hepatotoxicity | Yes | Yes | No | Categorical (Yes/No) |
Toxicity | Skin sensitivity | No | No | No | Categorical (Yes/No) |
Toxicity | T. Pyriformis Toxicity | 0.31 | 0.316 | 1.154 | Numeric (log ug/L) |
Toxicity | Minnow toxicity | -0.697 | -0.165 | 0.583 | Numeric (log mM) |
In Fig.
(A1) Chemical structure of the B73aIII (A2) Bioavailability radar diagram of the B73aIII, (A3) ‘BOILED-Egg’ images of the B73aIII in SwissADME software, (B1)Chemical structure of the ligand B73aV, (B2) Bioavailability radar diagram of the B73aV, (B3) ‘BOILED-Egg’ images of the B73aV in SwissADME software, (C1) Chemical structure of the cocrystal ligand EST, (C2) Bioavailability radar diagram of the EST, (C3) ‘BOILED-Egg’ images of the EST in SwissADME software. (Boiled-egg images indicate BBB (Blood-Brain barrier) permeability, HIA (gastrointestinal absorption), PGP+ and PGP- (P glycoprotein)).
In this study, novel estrogen receptor (ER) inhibitors based upon the indole-based benzamide scaffold were designed and developed by molecular modeling methods. The 7-substituted-1-(4-(piperidin-1-ylmethoxy)benzyl)-1H-indole-3-carboxamide derivatives would have the potential to modulate ER-α in breast cancer. Molecular docking, molecular dynamics and ADMET studies revealed the binding modes, protein-ligand stability and in-silico pharmacokinetic/toxicities of indole-based benzamide derivatives. Finally, it is concluded that the newly designed indole-based benzamides can be used as selective estrogen receptor modulators (SERMs) against ER-α positive breast cancer.
The authors thank the Deanship of Scientific Research at Majmaah University, Saudi Arabia for supporting this work under the project number (R-2023-158). The molecular dynamics numerical calculations reported in this paper were partially performed at TUBITAK ULAKBIM in TURKEY, High Performance and Grid Computing Center (TRUBA resources).