Three-dimensional model of the human aromatase enzyme and density functional parameterization of the iron-containing protoporphyrin IX for a molecular dynamics study of heme-cysteinato cytochromes

Developing new anticancer agents: Design, synthesis, biological evaluation and in silico study of several functionalized pyrimidine-5-carbonitriles as small molecules modulators targeting breast cancer

Committed to our growing effort addressed toward the development of potent anti-breast cancer candidates, new 4-hydrazinylpyrimidine-5-carbonitriles featuring a morpholinyl or piperidinyl moiety at the position-2 and derivatized with various functionalities at the hydrazinyl group were designed through structure optimization, and their antiproliferative potency against two human breast cancer (BC) cell lines, relative to the reference drug 5-FU, was evaluated. Compounds showing remarkable cytotoxic activity versus the hormone dependent MCF-7 cell line (IC50 = 1.62 ± 0.06 µM- 9.88 ± 0.38 µM) and the non-hormone dependent MDA-MB-231 cell line (IC50 = 3.26 ± 0.14 µM-12.93 ± 0.55 µM) were further tested by multiple assays for clarification of their potential activity. Promising derivatives revealing low damage to healthy cells were subject to enzymatic inhibitory assessment against ARO and EGFR and their activities compared to letrozole and erlotinib respectively. Compounds 3c, 6a as well as compounds 4c, 4d proved to be good inhibitors of the ARO and EGFR enzymes respectively. Active compounds were also evaluated for their underlying mode of action by further investigation for CDK, Hsp90, PI3K inhibition and compared to normal MCF-10A cells and assessed for their enhancement of the caspase 9 levels. Additionally, cell cycle analysis and apoptotic induction were performed. They demonstrated remarkable activities in the previous assays and emanated as leads as anti-breast cancer candidates. Eventually, molecular docking analysis revealed that hit compounds 3c, 4c, 4d, and 6a could bind favorably to the proposed in silico models of various protein-ligand interactions. Therefore, our promising top candidates, by demonstrating appreciable anti-breast cancer activities, present valuable prospects for optimization, potency enhancement and future application.

Pharmacophore-based virtual screening, molecular docking, and molecular dynamics investigation for the identification of novel, marine aromatase inhibitors

Breast cancer remains a leading cause of mortality among women worldwide. Our current research focuses on identifying effective therapeutic agents by targeting the human aromatase enzyme. Aromatase inhibitors (AIs) have been effective in treating postmenopausal breast cancer but face challenges such as drug resistance and long-term side effects like cognitive decline and osteoporosis. Natural products, especially from marine organisms, are emerging as potential sources for new drug candidates due to their structural diversity and pharmacological properties. This study aims to discover marine natural products capable of inhibiting human aromatase by combining ligand-based and structure-based pharmacophore models for virtual screening against the Comprehensive Marine Natural Products Database. From the initial virtual screening of more than 31,000 compounds, 1,385 marine natural products were identified as possible candidates. Following initial molecular docking analysis, only four compounds managed to pass the criteria this research has introduced to confirm strong binding affinity to aromatase. All four compounds yielded acceptable binding affinities, with CMPND 27987 having the highest -10.1 kcal/mol. All four hits were subjected to molecular dynamics, and CMPND 27987 was further confirmed to be the most stable at the protein's active site, with an MM-GBSA free binding energy of -27.75 kcal/mol. Our in silico studies indicate that CMPND 27987 interacts effectively within the binding site of the human aromatase, maintaining high affinity and stability. Based on these findings, we propose that CMPND 27987 could hold significant potential for further lead optimization and drug development.

Design, synthesis and molecular docking study of new pyrimidine-based hydrazones with selective anti-proliferative activity against MCF-7 and MDA-MB-231 human breast cancer cell lines

Efforts were directed on the design, synthesis and evaluation of the anticancer activity of some pyrimidine-based hydrazones against two breast cancer cell lines, MCF-7 and MDA-MB-231. Preliminary screening results revealed that some candidates scrutinized for their antiproliferative activities exhibited IC50 values of 0.87 μM-12.91 μM in MCF-7 and 1.75 μM-9.46 μM in MDA-MB-231 cells, indicating almost equal activities on both cell lines and better growth inhibition activities than those of the positive control 5-fluorouracil (5-FU) which displayed IC50 values of 17.02 μM and 11.73 μM respectively. Selectivity of the significantly active compounds was estimated against MCF-10A normal breast cells when compounds 7c, 8b, 9a and 10b exhibited superior activity for cancerous cells than for normal cells when compound 10b presented the best selectivity Index (SI) with respect to both MCF-7 and MDA-MB-231 cancer cells in comparison to the reference drug 5-FU. Mechanisms of their actions were explored by inspecting activation of caspase-9, annexin V staining and cell cycle analysis. It was noticed that compounds 7c, 8b, 8c 9a-c and 10b produced an increase in caspase-9 levels in MCF-7 treated cells with 10b inducing the highest elevation (27.13 ± 0.54 ng/mL) attaining 8.26-fold when compared to control MCF-7 which was higher than that of staurosporine (19.011 ± 0.40 ng/mL). The same compounds boosted caspase-9 levels in MDA-MB-231 treated cells when an increase in caspase-9 concentration reaching 20.40 ± 0.46 ng/mL (4.11-fold increase) was observed for compound 9a. We also investigated the role of these compounds for their increasing apoptosis ability against the 2 cell lines. Compounds 7c, 8b and 10b tested on MCF-7 cells displayed pre-G1 apoptosis and arrested cell cycle in particular at the S and G1 phases. Further clarification of their effects was made by modulating their related activities as inhibitors of ARO and EGFR enzymes when 8c and 9b showed 52.4% and 58.9% inhibition activity relative to letrozole respectively and 9b and 10b showed 36% and 39% inhibition activity of erlotinib. Also, the inhibition activity was verified by docking into the chosen enzymes.

Developing multitarget coumarin based anti-breast cancer agents: synthesis and molecular modeling study

A new series of 7-substituted coumarin scaffolds containing a methyl ester moiety at the C4-position were synthesized and tested for their in vitro anti-proliferative activity against MCF-7 and MDA-MB-231 breast cancer cell lines using Doxorubicin (DOX) as reference. Compounds 2 and 8 showed noticeable selectivity against MCF-7 with IC50 = 6.0 and 5.8 µM, respectively compared to DOX with IC50 = 5.6 µM. Compounds 10, 12, and 14 exhibited considerable selectivity against Estrogen Negative cells with IC50 = 2.3, 3.5, and 1.9 µM, respectively) compared to DOX with (IC50 = 7.3 µM). The most promising compounds were tested as epidermal growth factor receptor and aromatase (ARO) enzymes inhibitors using erlotinib and exemestane (EXM) as standards, respectively. Results proved that compound 8 elicited the highest inhibitory activity (94.73% of the potency of EXM), while compounds 10 and 12 displayed 97.67% and 81.92% of the potency of Erlotinib, respectively. Further investigation showed that the promising candidates 8, 10, and 12 caused cell cycle arrest at G0-G1 and S phases and induced apoptosis. The mechanistic pathway was confirmed by elevating caspases-9 and Bax/Bcl-2 ratio. A set of in silico methods was also performed including docking, bioavailability ADMET screening and QSAR study.

Benzoxazole-appended piperidine derivatives as novel anticancer candidates against breast cancer

Novel series of benzoxazole-appended piperidine derivatives were planned, synthesized and screened against two breast cancer cell lines. Considerable antiproliferative activity was observed for screened compounds (IC₅₀ = 33.32 ± 0.2 µM to 7.31 ± 0.43 µM and 1.66 ± 0.08 µM to 12.10 ± 0.57 µM) against MCF-7 and MDA-MB-231 cell lines respectively being more potent than doxorubicin (IC₅₀ = 8.20 ± 0.39 µM and 13.34 ± 0.63 µM respectively). Active compounds were submitted for enzyme inhibition assays when 4d and 7h demonstrated potent EGFR inhibition (0.08 ± 0.002 µM and 0.09 ± 0.002 µM respectively) compared to erlotinib (0.11 ± 0.003 µM). However, no one compound displayed effective ARO inhibition activity as tested compounds were less active than letrozole. Apoptosis inducing ability results implied that apoptosis was provoked by significant stimulation of caspase-9 protein levels (4.25-7.04-fold) upon treatment of MCF-7 cells with 4a, 7h, 9, 12e and 12f. Alternatively, MDA-MB-231 cells treated with 4d, 7a, 12b and 12c considerably increased caspase-9 levels (2.32-4.06-fold). Cell cycle arrest and annexin-V/Propidium iodide assays further confirmed apoptosis when tested compounds arrested cell cycle at various phases and demonstrated high annexin V binding affinity. Docking outcomes proved valuable binding affinities for compounds 4d and 7h to EGFR enzyme while compounds 4a and 12e, upon docking into the active site of ARO, failed to interact with heme, suggesting their inabilities to act as AIs. Therefore, these benzoxazoles can act as promising candidates exhibiting EGFR inhibition and apoptosis-promoting properties.

Cytotoxic Potential of Alternaria tenuissima AUMC14342 Mycoendophyte Extract: A Study Combined with LC-MS/MS Metabolic Profiling and Molecular Docking Simulation

Breast, cervical, and ovarian cancers are among the most serious cancers and the main causes of mortality in females worldwide, necessitating urgent efforts to find newer sources of safe anticancer drugs. The present study aimed to evaluate the anticancer potency of mycoendophytic Alternaria tenuissima AUMC14342 ethyl acetate extract on HeLa (cervical cancer), SKOV-3 (ovarian cancer), and MCF-7 (breast adenocarcinoma) cell lines. The extract showed potent effect on MCF-7 cells with an IC50 value of 55.53 μg/mL. Cell cycle distribution analysis of treated MCF-7 cells revealed a cell cycle arrest at the S phase with a significant increase in the cell population (25.53%). When compared to control cells, no significant signs of necrotic or apoptotic cell death were observed. LC-MS/MS analysis of Alternaria tenuissima extract afforded the identification of 20 secondary metabolites, including 7-dehydrobrefeldin A, which exhibited the highest interaction score (-8.0156 kcal/mol) in molecular docking analysis against human aromatase. Regarding ADME pharmacokinetics and drug-likeness properties, 7-dehydrobrefeldin A, 4'-epialtenuene, and atransfusarin had good GIT absorption and water solubility without any violation of drug-likeness rules. These findings support the anticancer activity of bioactive metabolites derived from endophytic fungi and provide drug scaffolds and substitute sources for the future development of safe chemotherapy.

New functionalized 6-thienylpyrimidine-5-carbonitriles as antiproliferative agents against human breast cancer cells

6-Thienylpyrimidine-5-carbonitrile derivatives were synthesized and screened for their in vitro antiproliferative activities against two human breast cancer cell lines in comparison to 5-fluorouracil as a reference. Compounds 2, 3a-c, and 6b evolved as the most active congeners against both cell lines, while others showed selectivity for only one cell line. Compound 2 exerted its effect through inhibition of the epidermal growth factor receptor (EGFR), while 6b showed less aromatase inhibitory activity than letrozole. The rest of the tested compounds did not show significant inhibition, and it can be assumed that they exert their antiproliferative activity through different target mechanisms. In addition, caspase-9 protein activation assays, cell cycle analysis using flow cytometry, and annexin V-fluorescein isothiocyanate-propidium iodide (FITC/PI) dual staining assays were performed for the most active compounds. All the tested compounds were found to be potent pyrimidine derivatives able to initiate apoptosis in MCF-7 and MDA-MB-231 cells.

New thieno[3,2-d]pyrimidine-based derivatives: Design, synthesis and biological evaluation as antiproliferative agents, EGFR and ARO inhibitors inducing apoptosis in breast cancer cells

An array of newly synthesized thieno[3,2-d]pyrimidine-based derivatives and thienotriazolopyrimidines hybridized with some pharmacophoric anticancer fragments were designed, synthesized and assessed for their in vitro antiproliferative activity against MCF-7 and MDA-MB-231 breast cancer cell lines using erlotinib and pictilisib as reference standards in the MTT assay. In general, many compounds were endowed with considerable antiproliferative activity (IC₅₀ = 0.43-1.31 µM). Some of the tested compounds, namely 3c, 5b, 5c, 9d, 10, 11b and 13 displayed remarkable antiproliferative activity against both cell lines. Meanwhile, compounds 2c-e, 3b, 4a, 5a, 9c and 15b showed noticeable selectivity against MCF-7 cells while compounds 2b, 3a, 4b, 6a-c, 7, 8, 9b and 12 exhibited considerable selectivity against MDA-MB-231 cells. Further mechanistic evidences for their anticancer activities were provided by screening the most potent compounds against MCF-7 and/or MDA-MB-231 cells for EGFR and ARO inhibitory activities using erlotinib and letrozole as reference standards respectively. Results proved that, in general, tested compounds were better EGFRIs than ARIs. In addition, significant overexpression in caspase-9 level in treated MCF-7 breast cell line samples was observed for all tested compounds with the 4-fluorophenylhydrazone derivative 2d exhibiting the highest activation. In treated MDA-MB-231 breast cell line samples, 11b was found to highly induce caspase-9 level thereby inducing apoptosis. Cell cycle analysis and Annexin V-FITC/PI assay were also assessed for active compounds where results indicated that all tested compounds induced preG1 apoptosis and cell cycle arrest at G2/M phase. Compound 9d, as an inhibitor of ARO, was observed to downregulate the downstream signaling proteins HSP27 and p-ERK in MCF-7 cells. Furthermore, compound 11b downregulated EGFR expression as well as the downstream signaling protein p-AKT. Docking experiments on EGFR and ARO enzymes supported their in vitro results. Thus, the thienotriazolopyrimidines 11b and 12 showing good EGFR inhibition and the thieno[3,2-d]-pyrimidine derivatives 3b and 9d, eliciting the best ARO inhibition activity, can be considered as new candidates as anti-breast cancer agents that necessitate further development.

Reconsidering Aromatase for Breast Cancer Treatment: New Roles for an Old Target

The current therapeutic approach for the treatment of hormone dependent breast cancer includes interference with estrogen receptors via either selective modulators or estrogens deprivation, by preventing their biosynthesis with aromatase inhibitors. Severe side effects and acquired resistance are drawbacks of both drug classes, and the efforts to overcome these issues still allow for research in this field to be animated. This review reports on recent findings that have opened new avenues for reconsidering the role of aromatase enzymes (and estrogen receptors) leading to the possibility of looking at well-known targets in a new perspective.

Aromatase Inhibitors for the Treatment of Breast Cancer: A Journey from the Scratch

Background:

Estrogens are essential for the growth of breast cancer in the case of premenopausal as well as in postmenopausal women. However, most of the breast cancer incidences are reported in postmenopausal women and the concurrent risk surges with an increase in age. Since the enzyme aromatase catalyses essential steps in estrogen biosynthesis, Aromatase Inhibitors (AIs) are effective targeted therapy in patients with Estrogen Receptor positive (ER+) breast cancer. AIs are more effective than Selective Estrogen Receptor Modulators (SERMs) because they block both the genomic and nongenomic activities of ER. Till date, first, second and third-generation AIs have been approved by the FDA. The third-generation AIs, viz. Letrozole, Anastrozole, Exemestane, are currently used in the standard treatment for postmenopausal breast cancer.

Methods:

Data were collected from Medline, PubMed, Google Scholar, Science Direct through searching of keywords: ‘aromatase’, ‘aromatase inhibitors’, ‘breast cancer’, ‘steroidal aromatase inhibitors’, ‘non-steroidal inhibitors’ and ‘generations of aromatase inhibitors’.

Results:

In the current scenario of breast cancer chemotherapy, AIs are the most widely used agents which reveal optimum efficacy along with the least side effects. Keeping in view the prominence of AIs in breast cancer therapy, this review covered the detailed description of aromatase including its role in the biosynthesis of estrogen, biochemistry, gene expression, 3D-structure, and information of reported AIs along with their role in breast cancer treatment.

Conclusion:

AIs are the mainstream solution of the ER+ breast cancer treatment regimen with the continuous improvement of human understanding of the importance of a healthy life of women suffering from breast cancer.

History of Estrogen: Its Purification, Structure, Synthesis, Biologic Actions, and Clinical Implications

This mini-review summarizes key points from the Clark Sawin Memorial Lecture on the History of Estrogen delivered at Endo 2018 and focuses on the rationales and motivation leading to various discoveries and their clinical applications. During the classical period of antiquity, incisive clinical observations uncovered important findings; however, extensive anatomical dissections to solidify proof were generally lacking. Initiation of the experimental approach followed later, influenced by Claude Bernard's treatise "An Introduction to the Study of Experimental Medicine." With this approach, investigators began to explore the function of the ovaries and their "internal secretions" and, after intensive investigations for several years, purified various estrogens. Clinical therapies for hot flashes, osteoporosis, and dysmenorrhea were quickly developed and, later, methods of hormonal contraception. Sophisticated biochemical methods revealed the mechanisms of estrogen synthesis through the enzyme aromatase and, after discovery of the estrogen receptors, their specific biologic actions. Molecular techniques facilitated understanding of the specific transcriptional and translational events requiring estrogen. This body of knowledge led to methods to prevent and treat hormone-dependent neoplasms as well as a variety of other estrogen-related conditions. More recently, the role of estrogen in men was uncovered by prismatic examples of estrogen deficiency in male patients and by knockout of the estrogen receptor and aromatase in animals. As studies became more extensive, the effects of estrogen on nearly every organ were described. We conclude that the history of estrogen illustrates the role of intellectual reasoning, motivation, and serendipity in advancing knowledge about this important sex steroid.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

Coumarin: A Natural, Privileged and Versatile Scaffold for Bioactive Compounds

Many naturally occurring substances, traditionally used in popular medicines around the world, contain the coumarin moiety. Coumarin represents a privileged scaffold for medicinal chemists, because of its peculiar physicochemical features, and the versatile and easy synthetic transformation into a large variety of functionalized coumarins. As a consequence, a huge number of coumarin derivatives have been designed, synthesized, and tested to address many pharmacological targets in a selective way, e.g., selective enzyme inhibitors, and more recently, a number of selected targets (multitarget ligands) involved in multifactorial diseases, such as Alzheimer's and Parkinson's diseases. In this review an overview of the most recent synthetic pathways leading to mono- and polyfunctionalized coumarins will be presented, along with the main biological pathways of their biosynthesis and metabolic transformations. The many existing and recent reviews in the field prompted us to make some drastic selections, and therefore, the review is focused on monoamine oxidase, cholinesterase, and aromatase inhibitors, and on multitarget coumarins acting on selected targets of neurodegenerative diseases.

Comprehensive and Automated Linear Interaction Energy Based Binding-Affinity Prediction for Multifarious Cytochrome P450 Aromatase Inhibitors

Cytochrome P450 aromatase (CYP19A1) plays a key role in the development of estrogen dependent breast cancer, and aromatase inhibitors have been at the front line of treatment for the past three decades. The development of potent, selective and safer inhibitors is ongoing with in silico screening methods playing a more prominent role in the search for promising lead compounds in bioactivity-relevant chemical space. Here we present a set of comprehensive binding affinity prediction models for CYP19A1 using our automated Linear Interaction Energy (LIE) based workflow on a set of 132 putative and structurally diverse aromatase inhibitors obtained from a typical industrial screening study. We extended the workflow with machine learning methods to automatically cluster training and test compounds in order to maximize the number of explained compounds in one or more predictive LIE models. The method uses protein-ligand interaction profiles obtained from Molecular Dynamics (MD) trajectories to help model search and define the applicability domain of the resolved models. Our method was successful in accounting for 86% of the data set in 3 robust models that show high correlation between calculated and observed values for ligand-binding free energies (RMSE < 2.5 kJ mol-1), with good cross-validation statistics.

Quantum mechanical/molecular mechanical and docking study of the novel analogues based on hybridization of common pharmacophores as potential anti-breast cancer agents

In an attempt to identify some new potential leads as anti-breast cancer agents, novel hybrid compounds were designed by molecular hybridization approach. These derivatives were structurally derived from hybrid benzofuran-imidazole and quinazolinone derivatives, which had shown good cytotoxicity against the breast cancer cell line (MCF-7). Since aromatase enzyme (CYP19) is highly expressed in the MCF-7 cell line, the binding of these novel hybrid compounds to aromatase was investigated using the docking method. In this study, due to the positive charge on the imidazole ring of the designed ligands and also, the presence of heme iron in the active site of the enzyme, it was decided to optimize the ligand inside the protein to obtain more realistic atomic charges for it. Quantum mechanical/molecular mechanical (QM/MM) method was used to obtain more accurate atomic charges of ligand for docking calculations by considering the polarization effects of CYP19 on ligands. It was observed that the refitted charge improved the binding energy of the docked compounds. Also, the results showed that these novel hybrid compounds were adopted properly within the aromatase binding site, thereby suggesting that they could be potential inhibitors of aromatase. The main binding modes in these complexes were through hydrophobic and H bond interactions showing agreement with the basic physicochemical features of known anti aromatase compounds. Finally, the complex structures obtained from the docking study were used for single point QM/MM calculations to obtain more accurate electronic interaction energy, considering the electronic polarization of the ligand by its protein environment.

Ligand- and Structure-Based Drug Design of Non-Steroidal Aromatase Inhibitors (NSAIs) in Breast Cancer

Aromatase is a multienzyme complex overexpressed in breast cancer and responsible for estrogen production. It is the potential target for designing anti-breast cancer drugs. Ligand and Structure-Based Drug Designing approaches (LBDD and SBDD) are involved in development of active and more specific Nonsteroidal Aromatase Inhibitors (NSAIs). Different LBDD and SBDD approaches are presented here to understand their utility in designing novel NSAIs. It is observed that molecules should possess a five or six membered heterocyclic nitrogen containing ring to coordinate with heme portion of aromatase for inhibition. Moreover, one or two hydrogen bond acceptor features, hydrophobicity, and steric factors may play crucial roles for anti-aromatase activity. Electrostatic, van der Waals, and p-p interactions are other important factors that determine binding affinity of inhibitors. HQSAR, LDA-QSAR, GQSAR, CoMFA, and CoMSIA approaches, pharmacophore mapping followed by virtual screening, docking, and dynamic simulation may be effective approaches for designing new potent anti-aromatase molecules.

Characterization of differences in substrate specificity among CYP1A1, CYP1A2 and CYP1B1: an integrated approach employing molecular docking and molecular dynamics simulations

Recent trends in new drug discovery of anticancer drugs have made oncologists more aware of the fact that the new drug discovery must target the developing mechanism of tumorigenesis to improve the therapeutic efficacy of antineoplastic drugs. The drugs designed are expected to have high affinity towards the novel targets selectively. Current research highlights overexpression of CYP450s, particularly cytochrome P450 1A1 (CYP1A1), in tumour cells, representing a novel target for anticancer therapy. However, the CYP1 family is identified as posing significant problems in selectivity of anticancer molecules towards CYP1A1. Three members have been identified in the human CYP1 family: CYP1A1, CYP1A2 and CYP1B1. Although sequences of the three isoform have high sequence identity, they have distinct substrate specificities. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics, design novel antitumour compounds that could be specifically metabolized by only CYP1A1 to mediate their antitumour activity and elucidate the reasons for differences in substrate specificity profile among the three proteins. In the present study, we employed a combination of computational methodologies: molecular docking and molecular dynamics simulations. We utilized eight substrates for elucidating the difference in substrate specificity of the three isoforms. Lastly, we conclude that the substrate specificity of a particular substrate depends upon the type of the active site residues, the dynamic motions in the protein structure upon ligand binding and the physico-chemical characteristics of a particular ligand. Copyright © 2016 John Wiley & Sons, Ltd.

Recent Progress in the Discovery of Next Generation Inhibitors of Aromatase from the Structure-Function Perspective

Human aromatase catalyzes the synthesis of estrogen from androgen with high substrate specificity. For the past 40 years, aromatase has been a target of intense inhibitor discovery research for the prevention and treatment of estrogen-dependent breast cancer. The so-called third generation aromatase inhibitors (AIs) letrozole, anastrozole, and the steroidal exemestane were approved in the U.S. in the late 1990s for estrogen-dependent postmenopausal breast cancer. Efforts to develop better AIs with higher selectivity and lower side effects were handicapped by the lack of an experimental structure of this unique P450. The year 2009 marked the publication of the crystal structure of aromatase purified from human placenta, revealing an androgen-specific active site. The structure has reinvigorated research activities on this fascinating enzyme and served as the catalyst for next generation AI discovery research. Here, we present an account of recent developments in the AI field from the perspective of the enzyme's structure-function relationships.

Celebrating 75 years of oestradiol

Oestrogens exert important effects on the reproductive as well as many other organ systems in both men and women. The history of the discovery of oestrogens, the mechanisms of their synthesis, and their therapeutic applications are very important components of the fabric of endocrinology. These aspects provide the rationale for highlighting several key components of this story. Two investigators, Edward Doisy and Alfred Butenandt, purified and crystalized oestrone nearly simultaneously in 1929, and Doisy later discovered oestriol and oestradiol. Butenandt won the Nobel Prize for this work and Doisy's had to await his purification of vitamin K. Early investigators quickly recognized that oestrogens must be synthesized from androgens and later investigators called this process aromatization. The aromatase enzyme was then characterized, its mechanism determined, and its structure identified after successful crystallization. With the development of knock-out methodology, the precise effects of oestrogen in males and females were defined and clinical syndromes of deficiency and excess described. Their discovery ultimately led to the development of oral contraceptives, treatment of menopausal symptoms, therapies for breast cancer, and induction of fertility, among others. The history of the use of oestrogens for postmenopausal women to relieve symptoms has been characterized by cyclic periods of enthusiasm and concern. The individuals involved in these studies, the innovative thinking required, and the detailed understanding made possible by evolving biologic and molecular techniques provide many lessons for current endocrinologists.

Design of novel potential aromatase inhibitors via hybrid pharmacophore approach: docking improvement using the QM/MM method

QM/MM and docking methods were used for designing novel hybrid aromatase inhibitors incorporating benzofuran, imidazole and quinazolinone moieties.

Ligand- and Structure-Based Drug Design of Non-Steroidal Aromatase Inhibitors (NSAIs) in Breast Cancer

Aromatase is a multienzyme complex overexpressed in breast cancer and responsible for estrogen production. It is the potential target for designing anti-breast cancer drugs. Ligand and Structure-Based Drug Designing approaches (LBDD and SBDD) are involved in development of active and more specific Nonsteroidal Aromatase Inhibitors (NSAIs). Different LBDD and SBDD approaches are presented here to understand their utility in designing novel NSAIs. It is observed that molecules should possess a five or six membered heterocyclic nitrogen containing ring to coordinate with heme portion of aromatase for inhibition. Moreover, one or two hydrogen bond acceptor features, hydrophobicity, and steric factors may play crucial roles for anti-aromatase activity. Electrostatic, van der Waals, and p-p interactions are other important factors that determine binding affinity of inhibitors. HQSAR, LDA-QSAR, GQSAR, CoMFA, and CoMSIA approaches, pharmacophore mapping followed by virtual screening, docking, and dynamic simulation may be effective approaches for designing new potent anti-aromatase molecules.

Theoretical study on phthalocynine-Fe(II)-based fluorescent sensors for cyanide anion

Three fluorescent sensors bearing phthalocynine- Fe ( II ) moiety were designed specifically for detecting cyanide anions, and investigated by DFT/TDDFT method. Comparison of the geometrical and photophysical properties of these sensor molecules, equipped with H -, carbamoyl and phthalimino groups, provided a deep insight into the sensor–cyanide interactions. The binding energy calculation shows that all the three sensors have good selectivity to the cyanide anion. Especially, frontier molecular orbital analysis confirmed that there was a photoinduced electron transfer (PET) process in the sensor with phthalimino group upon the addition of cyanide anion. This process could cause the fluorescence change. As a result, the sensor with phthalimino group displayed several favorable sensing properties.

Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure–activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein–ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

Aromatase, estrone sulfatase, and 17β-hydroxysteroid dehydrogenase: structure-function studies and inhibitor development

Aromatase, estrone sulfatase, and 17β-hydroxysteroid dehydrogenase type 1 are involved in the key steps of 17β-estradiol biosynthesis. Structure-function studies of aromatase, estrone sulfatase and 17β-hydroxysteroid dehydrogenase type 1 are important to evaluate the molecular basis of the interaction between these enzymes and their inhibitors. Selective and potent inhibitors of the three enzymes have been developed as antiproliferative agents in hormone-dependent breast carcinoma. New treatment strategies for hormone-dependent breast cancer are discussed.

Molecular Docking of Aromatase Inhibitors

Aromatase is an enzyme that plays a critical role in the development of estrogen receptor positive breast cancer. As aromatase catalyzes the aromatization of androstenedione to estrone, a naturally occurring estrogen, it is a promising drug target for therapeutic management. The undesirable effects found in aromatase inhibitors (AIs) that are in clinical use necessitate the discovery of novel AIs with higher selectivity, less toxicity and improving potency. In this study, we elucidate the binding mode of all three generations of AI drugs to the crystal structure of aromatase by means of molecular docking. It was demonstrated that the docking protocol could reliably reproduce the interaction of aromatase with its substrate with an RMSD of 1.350 Å. The docking study revealed that polar (D309, T310, S478 and M374), aromatic (F134, F221 and W224) and non-polar (A306, A307, V370, L372 and L477) residues were important for interacting with the AIs. The insights gained from the study herein have great potential for the design of novel AIs.

Design, synthesis, and biological evaluation of imidazolyl derivatives of 4,7-disubstituted coumarins as aromatase inhibitors selective over 17-α-hydroxylase/C17-20 lyase

The design, synthesis, and biological evaluation of a series of new aromatase (AR, CYP19) inhibitors bearing an imidazole ring linked to a 7-substituted coumarin scaffold at position 4 (or 3) are reported. Many compounds exhibited an aromatase inhibitory potency in the nanomolar range along with a high selectivity over 17-α-hydroxylase/C17-20 lyase (CYP17). The most potent AR inhibitor was the 7-(3,4-difluorophenoxy)-4-imidazolylmethyl coumarin 24 endowed with an IC(50) = 47 nM. Docking simulations on a selected number of coumarin derivatives allowed the identification of the most important interactions driving the binding and clearly indicated the allowed and disallowed regions for appropriate structural modifications of coumarins and closely related heterocyclic molecular scaffolds.

Molecular docking and QSAR studies of aromatase inhibitor androstenedione derivatives

Objectives

Aromatase (CYP19) inhibitors have emerged as promising candidates for the treatment of estrogen-dependent breast cancer. In this study, a series of androstenedione derivatives with CYP19 inhibitory activity was subjected to a molecular docking study followed by quantitative structure–activity relationship (QSAR) analyses in search of ideal physicochemical characteristics of potential aromatase inhibitors.

Methods

The QSAR studies were carried out using both two-dimensional (topological, and structural) and three-dimesional (spatial) descriptors. We also used thermodynamic parameters along with 2D and 3D descriptors. Genetic function approximation (GFA) and genetic partial least squares (G/PLS) were used as chemometric tools for QSAR modelling.

Key findings

The docking study indicated that the important interacting amino acids in the active site were Met374, Arg115, Ile133, Ala306, Thr310, Asp309, Val370, Leu477 and Ser478. The 17-keto oxygen of the ligands is responsible for the formation of a hydrogen bond with Met374 and the remaining parts of the molecules are stabilized by the hydrophobic interactions with the non-polar amino acids. The C2 and C19 positions in the ligands are important for maintaining the appropriate orientation of the molecules in the active site. The results of docking experiments and QSAR studies supported each other.

Conclusions

The developed QSAR models indicated the importance of some Jurs parameters, structural parameters, topological branching index and E-state indices of different fragments. All the developed QSAR models were statistically significant according to the internal and external validation parameters.

Molecular docking and QSAR study on steroidal compounds as aromatase inhibitors

In order to develop more potent, selective and less toxic steroidal aromatase (AR) inhibitors, molecular docking, 2D and 3D hybrid quantitative structure-activity relationship (QSAR) study have been conducted using topological, molecular shape, spatial, structural and thermodynamic descriptors on 32 steroidal compounds. The molecular docking study shows that one or more hydrogen bonds with MET374 are one of the essential requirements for the optimum binding of ligands. The QSAR model obtained indicates that the aromatase inhibitory activity can be enhanced by increasing SIC, SC_3_C, Jurs_WNSA_1, Jurs_WPSA_1 and decreasing CDOCKER interaction energy (ECD), IAC_Total and Shadow_XZfrac. The predicted results shows that this model has a comparatively good predictive power which can be used in prediction of activity of new steroidal aromatase inhibitors.

Brain chemistry: how does P450 catalyze the O-demethylation reaction of 5-methoxytryptamine to yield serotonin?

Density functional theory has been applied to elucidate the mechanism of the O-demethylation reaction that generates serotonin from 5-methoxytryptamine (5-MT); a process that is efficiently catalyzed by P450 CYP2D6. Two substrates, the neutral 5-MT and the protonated 5-MTH(+), were used to probe the reactivity of CYP2D6 compound I. Notably, the H-abstraction process is found to be slightly more facile for 5-MT. However, our DFT augmented by docking results show that the amino acid Glu216 in the active site holds the NH(3)(+) tail of the 5-MTH(+) substrate in an upright conformation and thereby controls the regioselectivity of the bond activation. Thus, the substrate protonation serves an important function in maximizing the yield of serotonin. This finding is in accord with experimental conclusions that 5-MTH(+) serves as the substrate for the CYP2D6 enzyme. The study further shows that the H-abstraction follows two-state reactivity (TSR), whereas the rebound path may involve more states due to the appearance of both Fe(IV) and Fe(III) electromers during the reaction of 5-MTH(+).

Progress in aromatase research and identification of key future directions

The IX International Aromatase Conference focused upon key developments in research related to the aromatase enzyme that had occurred since the last meeting. A session took place at the conclusion of conference discussing key areas for future research and issues currently facing researchers in the field. While significant progress on understanding structural elements of the enzyme and regulatory mechanisms of both the gene and protein provides an excellent basis for development of improved aromatase inhibitors and exploration of the important problem of aromatase inhibitor resistance, significant challenges remain. Increasing the speed with which findings are translated into clinical practice and finding an appropriate balance between basic and translational research were identified as areas which require further attention before the next meeting in 2010.

History of aromatase: saga of an important biological mediator and therapeutic target

Aromatase is the enzyme that catalyzes the conversion of androgens to estrogens. Initial studies of its enzymatic activity and function took place in an environment focused on estrogen as a component of the birth control pill. At an early stage, investigators recognized that inhibition of this enzyme could have major practical applications for treatment of hormone-dependent breast cancer, alterations of ovarian and endometrial function, and treatment of benign disorders such as gynecomastia. Two general approaches ultimately led to the development of potent and selective aromatase inhibitors. One targeted the enzyme using analogs of natural steroidal substrates to work out the relationships between structure and function. The other approach initially sought to block adrenal function as a treatment for breast cancer but led to the serendipitous finding that a nonsteroidal P450 steroidogenesis inhibitor, aminoglutethimide, served as a potent but nonselective aromatase inhibitor. Proof of the therapeutic concept of aromatase inhibition involved a variety of studies with aminoglutethimide and the selective steroidal inhibitor, formestane. The requirement for even more potent and selective inhibitors led to intensive molecular studies to identify the structure of aromatase, to development of high-sensitivity estrogen assays, and to "mega" clinical trials of the third-generation aromatase inhibitors, letrozole, anastrozole, and exemestane, which are now in clinical use in breast cancer. During these studies, unexpected findings led investigators to appreciate the important role of estrogens in males as well as in females and in multiple organs, particularly the bone and brain. These studies identified the important regulatory properties of aromatase acting in an autocrine, paracrine, intracrine, neurocrine, and juxtacrine fashion and the organ-specific enhancers and promoters controlling its transcription. The saga of these studies of aromatase and the ultimate utilization of inhibitors as highly effective treatments of breast cancer and for use in reproductive disorders serves as the basis for this first Endocrine Reviews history manuscript.

Molecular characterization of aromatase

Aromatase is an estrogen synthetase. Estrogens are female sex hormones involved in the development and growth of breast tumors. Elucidation of the structure-function relationship of aromatase has been of significant interest since its inhibitors have shown great promise in fighting breast cancer. Aromatase belongs to the cytochrome P450 family and forms an electron-transfer complex with its partner, NADPH-cytochrome P450 reductase. Because of the membrane-bound character and heme-binding instability, no crystal structure of aromatase has been reported so far. Much remains to be investigated, including the 3-dimensional structure of aromatase, interaction between aromatase and reductase, catalytic mechanism of estrogen synthesis by aromatase, and the binding mechanism of aromatase inhibitors. This review presents current knowledge about structural and functional characteristics of aromatase to address unsolved mysteries about this enzyme.

Synthesis and biochemical studies of 17-substituted androst-3-enes and 3,4-epoxyandrostanes as aromatase inhibitors

A series of 5alpha-androst-3-enes and 3alpha,4alpha-epoxy-5alpha-androstanes were synthesized and tested for their abilities to inhibit aromatase in human placental microsomes. In these series the original C-17 carbonyl group was replaced by hydroxyl, acetyl and hydroxyimine groups. Inhibition kinetic analysis on the most potent steroid of these series revealed that it inhibits the enzyme in a competitive manner (IC(50)=6.5 microM). The achieved data pointed out the importance of the C-17 carbonyl group in the D-ring of the studied steroids as a structural feature required to reach maximum aromatase inhibitory activity. Further, at least one carbonyl group (C-3 or C-17) seems to be essential to effective aromatase inhibition.

Molecular mechanisms of aromatase inhibition by new A, D-ring modified steroids

A recent approach for treatment and prevention of estrogen-dependent breast cancer focuses on the inhibition of aromatase, the enzyme that catalyzes the final step of estrogen biosynthesis. Some synthetic steroids, such as formestane and exemestane, resembling the natural enzyme substrate androstenedione, revealed to be potent and useful aromatase inhibitors (AIs) and were approved for the treatment of estrogen-dependent breast cancer in postmenopausal women. Recently, we found that five newly synthesized steroids with chemical features in the A- and D-rings considered important for drug-receptor interaction efficiently inhibit aromatase derived from human placental microsomes. In this work, these steroids showed a similar pattern of anti-aromatase activity in several aromatase-expressing cell lines. 5alpha-androst-3-en-17-one and 3alpha,4alpha-epoxy-5alpha-androstan-17-one were revealed to be the most potent inhibitors. These compounds induced a time-dependent inhibition of aromatase, showing to be irreversible AIs. The specific interactions of these compounds with aromatase active sites were further demonstrated by site-directed mutagenesis studies and evaluated by computer-aided molecular modeling. Both compounds were able to suppress hormone-dependent proliferation of MCF-7aro cells in a dose-dependent manner. These findings are important for the elucidation of a structure-activity relationship on aromatase, which may help in the development of new AIs.