GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study.

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

In silico design and cell-based evaluation of two dual anti breast cancer compounds targeting Bcl-2 and GPER.

According to WHO statistics, breast cancer (BC) disease represents about 2.3 million diagnosed and 685,000 deaths globally. Regarding histological classification of BC, the Estrogen (ER) and Progesterone (PR) receptors negative-expression cancer, named Triple-Negative BC (TNBC), represents the most aggressive type of this disease, making it a challenge for drug discovery. In this context, our research group, applying a well-established Virtual Screening (VS) protocol, in addition to docking and molecular dynamics simulations studies, yielded two ligands identified as 6 and 37 which were chemically synthesized and evaluated on MCF-7 and MDA-MB-231 cancer cell lines. Strikingly, 37 assayed on MDA-MB-231 (a TNBC cell model) depicted an outstanding value of 18.66 µM much lower than 65.67 µM yielded by Gossypol Bcl-2 inhibitor whose main disadvantage is to produce multiple toxic effects. Highlighted above, enforce the premise of the computational tools to find new therapeutic options against the most aggressive forms of breast cancer, as the results herein showed.

Targeting Several Biologically Reported Targets of Glioblastoma Multiforme by Assaying 2D and 3D Cultured Cells.

Glioblastoma multiforme (GBM) is account for 70% of all primary malignancies of the central nervous system. The median survival of human patients after treatment is around 15 months. There are several biological targets which have been reported that can be pursued using ligands with varied structures to treat this disease. In our group, we have developed several ligands that target a wide range of proteins involved in anticancer effects, such as histone deacetylase (HDACs), G protein-coupled estrogen receptor 1 (GPER), estrogen receptor-beta (ERß) and NADPH oxidase (NOX), that were screened on bidimensional (2D) and tridimensional (3D) GBM stem cells like (GSC). Our results show that some HDAC inhibitors show antiproliferative properties at 21-32 µM. These results suggest that in this 3D culture, HDACs could be the most relevant targets that are modulated to induce the antiproliferative effects that require in the future further experimental studies.

Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells.

The implementation of chemo- and bioinformatics tools is a crucial step in the design of structure-based drugs, enabling the identification of more specific and effective molecules against cancer without side effects. In this study, three new compounds were designed and synthesized with suitable absorption, distribution, metabolism, excretion and toxicity (ADME-tox) properties and high affinity for the G protein-coupled estrogen receptor (GPER) binding site by in silico methods, which correlated with the growth inhibitory activity tested in a cluster of cancer cell lines. Docking and molecular dynamics (MD) simulations accompanied by a molecular mechanics/generalized Born surface area (MMGBSA) approach yielded the binding modes and energetic features of the proposed compounds on GPER. These in silico studies showed that the compounds reached the GPER binding site, establishing interactions with a phenylalanine cluster (F206, F208 and F278) required for GPER molecular recognition of its agonist and antagonist ligands. Finally, a 3-(4,5-dimethylthiazol-2-yl)2,5-diphenyltetrazolium bromide (MTT) assay showed growth inhibitory activity of compounds 4, 5 and 7 in three different cancer cell lines-MIA Paca-2, RCC4-VA and Hep G2-at micromolar concentrations. These new molecules with specific chemical modifications of the GPER pharmacophore open up the possibility of generating new compounds capable of reaching the GPER binding site with potential growth inhibitory activities against nonconventional GPER cell models.

Structural and energetic basis for novel epicatechin derivatives acting as GPER agonists through the MMGBSA method.

(-)-Epicatechin (Epi) has been demonstrated to activate pathways involved in GPER-stimulated nitric oxide (NO) production via endothelial NO synthase, known as the eNOS/NO pathway. Previous studies combining synthesis of four Epi derivatives demonstrated that Epi and Epi-prop, Epi-4-prop and Epi-5-prop were able to bind GPER by acting as GPER agonists, whereas docking studies allowed observation of structural details of the binding of these derivatives at the GPER binding site. However, due to the nature of past studies, the theoretical methods employed did not allow observation of structural and energetic details linked to ligand binding at the GPER binding site. In this contribution, we explore the structural and energetic changes coupling the binding of Epi and its four derivatives to GPER. To this end, MD simulations on the microsecond scale (1 µs) with an MMGBSA approach were used for each GPER-ligand complex. Energetic analysis demonstrated that incorporation of several aliphatic chains to Epi contributed to increasing the affinity towards the GPER binding site, thus helping to explain the experimental evidence. Structural analysis demonstrated that Epi, Epi-4-prop and Epi-5-prop share more similar interactions at GPER binding sites with similar conformational behavior than with Epi-prop and Epi-Ms. However, Epi-prop had additional residues that could explain its different but related biological effects.

Synthesis and In Vitro Evaluation of Tetrahydroquinoline Derivatives as Antiproliferative Compounds of Breast Cancer via Targeting the GPER.

BACKGROUND: Some reports have demonstrated the role of the G Protein-coupled Estrogen Receptor (GPER) in growth and proliferation of breast cancer cells. OBJECTIVE: In an effort to develop new therapeutic strategies against breast cancer, we employed an in silico study to explore the binding modes of tetrahydroquinoline 2 and 4 to be compared with the reported ligands G1 and G1PABA. METHODS: This study aimed to design and filter ligands by in silico studies determining their Lipinski's rule, toxicity and binding properties with GPER to achieve experimental assays as anti-proliferative compounds of breast cancer cell lines. RESULTS: In silico studies suggest as promissory two tetrahydroquinoline 2 and 4 which contain a carboxyl group instead of the acetyl group (as is needed for G1 synthesis), which add low (2) and high hindrance (4) chemical moieties to explore the polar, hydrophobic and hindrance effects. Docking and molecular dynamics simulations of the target compounds were performed with GPER to explore their binding mode and free energy values. In addition, the target small molecules were synthesized and assayed in vitro using breast cancer cells (MCF-7 and MDA-MB-231). Experimental assays showed that compound 2 decreased cell proliferation, showing IC50 values of 50µM and 25µM after 72h of treatment of MCF-7 and MDA-MB-231 cell lines, respectively. Importantly, compound 2 showed a similar inhibitory effect on proliferation as G1 compound in MDA-MB-231 cells, suggesting that both ligands reach the GPER-binding site in a similar way, as was demonstrated through in silico studies. CONCLUSION: A concentration-dependent inhibition of cell proliferation occurred with compound 2 in the two cell lines regardless of GPER.

Structural insight into the binding mechanism of ATP to EGFR and L858R, and T790M and L858R/T790 mutants.

The L858R mutation in EGFR is particularly responsive to small tyrosine kinase inhibitors (TKIs) such as gefitinib and erlotinib. This efficacy decreases due to drug resistance conferred by a second mutation, T790M, which subsequently produces a double mutant, L858R/T790M. Although this resistance was initially attributed to steric blocking by the T790M mutation, experimental studies have demonstrated that differences in the binding affinities of TKIs to T790M and L858R/T790M mutants are more a result of the increased sensitivity of these mutants to ATP than to a decrease in the affinity to TKIs. Regrettably, detailed information at the atomic level on the origins of the increased binding affinity of mutants for ATP is lacking. In this study, we have combined structural data and molecular dynamics simulations with the MMGBSA approach to determine how the L858R, T790M and L858R/T790 mutations impact the binding mechanism of ATP with respect to wild-type EGFR. Structural and energetic analyses provided novel information that helps to explain the increased affinity of ATP to T790M and L858R/T790 mutants with respect to L858R and wild-type systems. In addition, it was observed that dimerization of the wild-type and mutant systems exerts dissimilar effects on the ATP binding affinity characteristic of negative cooperativity. Communicated by Ramaswamy H. Sarma.

Selection of a GPER1 Ligand via Ligand-based Virtual Screening Coupled to Molecular Dynamics Simulations and Its Anti-proliferative Effects on Breast Cancer Cells.

BACKGROUND: Recent reports have demonstrated the role of the G Protein-Coupled Estrogen Receptor 1 (GPER1) on the proliferation of breast cancer. The coupling of GPER1 to estrogen triggers cellular signaling pathways related to cell proliferation. OBJECTIVE: Develop new therapeutic strategies against breast cancer. METHOD: We performed in silico studies to explore the binding mechanism of a set of G15 /G1 analogue compounds. We included a carboxyl group instead of the acetyl group from G1 to form amides with several moieties to increase affinity on GPER1. The designed ligands were submitted to ligand-based and structure-based virtual screening to get insights into the binding mechanism of the best designed compound and phenol red on GPER1. RESULTS: According to the in silico studies, the best molecule was named G1-PABA ((3aS,4R,9bR)-4-(6- bromobenzo[d][1,3]dioxol-5-yl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-carboxylic acid). It was synthesized and assayed in vitro in breast cancer (MCF-7 and MDA-MB-231) and normal (MCF-10A) cell lines. Experimental studies showed that the target compound was able to decrease cell proliferation, IC50 values of 15.93 µM, 52.92 µM and 32.45 µM in the MCF-7, MDA-MB-231 and MCF-10A cell lines, respectively, after 72 h of treatment. The compound showed better IC50 values without phenol red, suggesting that phenol red interfere with the G1-PABA action at GPER1, as observed through in silico studies, which is present in MCF-7 cells according to PCR studies and explains the cell proliferation effects. CONCLUSION: Concentration-dependent inhibition of cell proliferation occurred with G1-PABA in the assayed cell lines and could be due to its action on GPER1.

Aromatic Regions Govern the Recognition of NADPH Oxidase Inhibitors as Diapocynin and its Analogues.

Oxidative stress is related to the pathogenesis and progress of several human diseases. NADPH oxidase (NOX), and mainly the NOX2 isoform, produces superoxide anions (O2•- ). To date, it is known that NOX2 can be inhibited by preventing the assembly of its subunits, p47phox and p22phox. In this work, we analyzed the binding to NOX2 of the apocynin dimer, diapocynin (C1), a known NOX2 inhibitor, and of 18 designed compounds (C2-C19) which have chemical relationships to C1, by in silico methods employing a p47phox structure from the Protein Data Bank (PDB code: 1WLP). C1 and six of the designed compounds were recognized in the region where p22phox binds to p47phox and makes π-π interactions principally with W193, W263, and Y279, which form an aromatic-rich region. C8 was chosen as the best compound according to the in silico studies and was synthesized and evaluated in vitro. C8 was able to prevent the production of reactive oxygen species (ROS) similar to C1. In conclusion, targeting the aromatic region of p47phox through π-interactions is important for inhibiting NOX activity.

Antihyperglycemic Activity of the Leaves from Annona cherimola Miller and Rutin on Alloxan-induced Diabetic Rats.

BACKGROUND: Annona cherimola, known as "chirimoya" has been reported in Mexican traditional medicine for the treatment of diabetes. OBJECTIVE: The aims of the present study were to validate and assess the traditional use of A. cherimola as an antidiabetic agent. MATERIALS AND METHODS: The ethanol extract from A. cherimola (300 mg/kg, EEAc), subsequent fractions (100 mg/kg), and rutin (30 mg/kg) were studied on alloxan-induced type 2 diabetic (AITD) and normoglycemic rats. In addition, oral glucose tolerance test (OGTT) and oral sucrose tolerance test (OSTT) were performed in normoglycemic rats. Molecular docking technique was used to conduct the computational study. RESULTS: Bioassay-guided fractionation of EEAc afforded as major antihyperglycemic compound, rutin. EEAc attenuated postprandial hyperglycemia in acute test using AITD rats (331.5 mg/dL) carrying the glycemic levels to 149.2 mg/dL. Rutin after 2 h, attenuated postprandial hyperglycemia in an acute assay using AITD rats such as EEAc, with maximum effect (150.0 mg/dL) being seen at 4 h. The antihyperglycemic activities of EEAc and rutin were comparable with acarbose (151.3 mg/dL). In the subchronic assay on AITD rats, the EEAc and rutin showed a reduction of the blood glucose levels since the 1st week of treatment, reaching levels similar to normoglycemic state (116.9 mg/kg) that stayed constant for the rest of the assay. OGTT and OSTT showed that EEAc and rutin significantly lowered blood glucose levels in normoglycemic rats at 2 h after a glucose or sucrose load such as acarbose. Computational molecular docking showed that rutin interacted with four amino acids residues in the enzyme α-glucosidase. CONCLUSION: The results suggest that rutin an α-glucosidase inhibitor was responsible in part of the antihyperglycemic activity of A. cherimola. Its in vivo antihyperglycemic activity is in good agreement with the traditional use of A. cherimola for the treatment of diabetes. SUMMARY: The ethanol extract from Annona cherimola (300 mg/kg, EEAc), subsequent fractions (100 mg/kg) and rutin (30 mg/kg) were studied on alloxan-induced type 2 diabetic (AITD) and normoglycemic rats. The results suggest that rutin; an α-glucosidase inhibitor was responsible in part of the antihyperglycemic activity of A. cherimola. Its in vivo antihyperglycemic activity is in good agreement with the traditional use of A. cherimola for the treatment of diabetes. Abbreviations Used: EEAc: The ethanol extract from Annona cherimola, AITD: Alloxan-induced type 2 diabetic rats, OGTT: Oral glucose tolerance test, OSTT: Oral sucrose tolerance test, DM: Diabetes mellitus.

Antiprotozoal Constituents from Annona cherimola Miller, a Plant Used in Mexican Traditional Medicine for the Treatment of Diarrhea and Dysentery.

BACKGROUND: Annona cherimola Miller (Annonaceae) is a medicinal plant frequently recommended in Mexican traditional medicine for the treatment of gastrointestinal disorders such as diarrhea and dysentery. OBJECTIVE: This work was undertaken to obtain information that support the traditional use of A. cherimola, on pharmacological basis using in vitro and computational experiments. MATERIAL AND METHODS: Bioassay-guided fractionation of the ethanol extract of the leaves of A. cherimola afforded five phenolic compounds: caffeic acid, quercetin, kaempferol, nicotinflorin, and rutin. RESULTS: The in vitro antiprotozoal assay showed that kaempferol was the most potent antiamoebic and antigiardial compound with IC50 values of 7.9 µg/mL for Entamoeba histolytica and 8.7 µg/mL for Giardia lamblia. Computational molecular docking study showed that kaempferol interacted in a region different than metronidazole in the enzyme pyruvate: ferredoxin oxidoreductase (PFOR). CONCLUSION: Considering that PFOR is a target of metronidazole; kaempferol may be a lead compound for the development of novel antiprotozoal agent. Also, these findings give support to the use of A. cherimola in the traditional medicine from México for the treatment of diarrhea and dysentery. SUMMARY: Bioassay-guided fractionation of the ethanol extract of the leaves of Annona cherimola afforded five phenolic compounds: caffeic acid, quercetin, kaempferol, nicotinflorin and rutin. The in vitro antiprotozoal assay showed that kaempferol was the most potent antiamoebic and antigiardial compound with IC50 values of 7.9 µg/mL for Entamoeba histolytica and 8.7 µg/mL for Giardia lamblia. Computational molecular docking study showed that kaempferol interacted in a region different that metronidazole in the enzyme pyruvate: ferredoxin oxidoreductase. Abbreviations used: PFOR:Pyruvate:ferredoxin oxidoreductase, G: lamblia: Giardia lamblia, E: histolytica: Entamoeba histolytica.

Design, in silico studies, synthesis and in vitro evaluation of oseltamivir derivatives as inhibitors of neuraminidase from influenza A virus H1N1.

Since the neuraminidase (NA) enzyme of the influenza A virus plays a key role in the process of release of new viral particles from a host cell, it is often a target for new drug design. The emergence of NA mutations, such as H275Y, has led to great resistance against neuraminidase inhibitors, including oseltamivir and zanamivir. Hence, we herein designed a set of derivatives by modifying the amine and/or carboxylic groups of oseltamivir. After being screened for their physicochemical (Lipinski's rule) and toxicological properties, the remaining compounds were submitted to molecular and theoretical studies. The docking simulations provided insights into NA recognition patterns, demonstrating that oseltamivir modified at the carboxylic moiety and coupled with anilines had higher affinity and a better binding pose for NA than the derivatives modified at the amine group. Based on these theoretical studies, the new oseltamivir derivatives may have higher affinity to mutant variants and possibly to other viral subtypes. Accordingly, two compounds were selected for synthesis, which together with their respective intermediates were evaluated for their cytotoxicity and antiviral activities. Their biological activity was then tested in cells infected with the A/Puerto Rico/916/34 (H1N1) influenza virus, and virus yield reduction assays were performed. Additionally, by measuring neuraminidase activity with the neuraminidase assay kit it was found that the compounds produced inhibitory activity on this enzyme. Finally, the infected cells were analysed with atomic force microscopy (AFM), observing morphological changes strongly suggesting that these compounds interfered with cellular release of viral particles.

Docking and QSAR Studies of Aryl-valproic Acid Derivatives to Identify Antiproliferative Agents Targeting the HDAC8.

BACKGROUND: Histone deacetylase 8 (HDAC8) is a plausible target for the development of novel anticancer drugs using a metal-chelating group and hydrophobic moieties as pharmacophores. It is known that valproic acid (administered as its salt, sodium valproate; VPANa+) is an HDAC8 inhibitor characterized by its hydrophobic chains. Nevertheless, VPA is hepatotoxic and VPA analogues might be explored for less hepatotoxic antiproliferative compounds. METHOD: In this work, docking and QSAR studies of 500 aryl-VPA derivatives as possible HDAC8 inhibitors were performed in order to explore and select potential anti-proliferative compounds. Docking results identified π-π, hydrogen bonds as the most important noncovalent interactions between HDAC8 (PDB: 3F07) and the ligands tested, whereas Belm4 was the best QSAR descriptor and classified as a 2D-BCUT descriptor. RESULT: Based on theoretical studies, compound DAVP042 was synthesized and evaluated in vitro for its antiproliferative activities on several cancer cell lines (A549-lung, MCF-7-breast, HCT116-colon and U937- lymphoid tissue) in comparison to VPA, as well as for its inhibitory activity on HDAC8 using in vitro models. DAVP042 demonstrated to have antiproliferative activity on all cancer cell lines employed, not only suggesting that this compound should be further studied, but also demonstrating that the methodology herein employed is appropriated to identify new therapeutic candidates.

N-(2-hydroxyphenyl)-2-propylpentanamide, a valproic acid aryl derivative designed in silico with improved anti-proliferative activity in HeLa, rhabdomyosarcoma and breast cancer cells.

Epigenetic alterations are associated with cancer and their targeting is a promising approach for treatment of this disease. Among current epigenetic drugs, histone deacetylase (HDAC) inhibitors induce changes in gene expression that can lead to cell death in tumors. Valproic acid (VPA) is a HDAC inhibitor that has antitumor activity at mM range. However, it is known that VPA is a hepatotoxic drug. Therefore, the aim of this study was to design a set of VPA derivatives adding the arylamine core of the suberoylanilide hydroxamic acid (SAHA) with different substituents at its carboxyl group. These derivatives were submitted to docking simulations to select the most promising compound. The compound 2 (N-(2-hydroxyphenyl)-2-propylpentanamide) was the best candidate to be synthesized and evaluated in vitro as an anti-cancer agent against HeLa, rhabdomyosarcoma and breast cancer cell lines. Compound 2 showed a better IC50 (µM range) than VPA (mM range) on these cancer cells. And also, 2 was particularly effective on triple negative breast cancer cells. In conclusion, 2 is an example of drugs designed in silico that show biological properties against human cancer difficult to treat as triple negative breast cancer.

Understanding the molecular basis of agonist/antagonist mechanism of GPER1/GPR30 through structural and energetic analyses.

The G-protein coupled receptors (GPCRs) represent the largest superfamily of membrane proteins in charge to pass the cell signaling after binding with their cognate ligands to the cell interior. In breast cancer, a GPCR named GPER1 plays a key role in the process of growth and the proliferation of cancer cells. In a previous study, theoretical methods were applied to construct a model of GPER1, which later was submitted to molecular dynamics (MD) simulations to perform a docking calculation. Based on this preceding work, it is known that GPER1 is sensitive to structural differences in its binding site. However, due to the nature of that past study, conformational changes linked to the ligand binding were not observed. Therefore, in this study, in order to explore the conformational changes coupled to the agonist/antagonist binding, MD simulations of about 0.25µs were performed for the free and bound states, summarizing 0.75µs of MD simulation in total. For the bound states, one agonist (G-1) and antagonist (G-15) were chosen since is widely known that these two molecules cause an impact on GPER1 mobility. Based on the conformational ensemble generated through MD simulations, we found that despite G-1 and G-15 being stabilized by similar map of residues, the structural differences between both ligands impact the hydrogen bond pattern not only at the GPER1 binding site but also along the seven-helix bundle, causing significant differences in the conformational mobility along the extracellular and cytoplasmic domain, and to a lesser degree in the curvatures of helix 2, helix 3 and helix 7 between the free and bound states, which is in agreement with reported literature, and might be linked to microscopic characteristics of the activated-inactivated transition. Furthermore, binding free energy calculations using the MM/GBSA method for the bound states, followed by an alanine scanning analysis allowed us to identify some important residues for the complex stabilization.

The effects of (-)-epicatechin on endothelial cells involve the G protein-coupled estrogen receptor (GPER).

We have provided evidence that the stimulatory effects of (-)-epicatechin ((-)-EPI) on endothelial cell nitric oxide (NO) production may involve the participation of a cell-surface receptor. Thus far, such entity(ies) has not been fully elucidated. The G protein-coupled estrogen receptor (GPER) is a cell-surface receptor that has been linked to protective effects on the cardiovascular system and activation of intracellular signaling pathways (including NO production) similar to those reported with (-)-EPI. In bovine coronary artery endothelial cells (BCAEC) by the use of confocal imaging, we evidence the presence of GPER at the cell-surface and on F-actin filaments. Using in silico studies we document the favorable binding mode between (-)-EPI and GPER. Such binding is comparable to that of the GPER agonist, G1. By the use of selective blockers, we demonstrate that the activation of ERK 1/2 and CaMKII by (-)-EPI is dependent on the GPER/c-SRC/EGFR axis mimicking those effects noted with G1. We also evidence by the use of siRNA the role that GPER has on mediating ERK1/2 activation by (-)-EPI. GPER appears to be coupled to a non Gαi/o or Gαs, protein subtype. To extrapolate our findings to an ex vivo model, we employed phenylephrine pre-contracted aortic rings evidencing that (-)-EPI can mediate vasodilation through GPER activation. In conclusion, we provide evidence that suggests the GPER as a potential mediator of (-)-EPI effects and highlights the important role that GPER may have on cardiovascular system protection.

Three amino acid derivatives of valproic acid: design, synthesis, theoretical and experimental evaluation as anticancer agents.

Valproic acid (VPA) is extensively used as an anticonvulsive agent and as a treatment for other neurological disorders. It has been shown that VPA exerts an anti-proliferative effect on several types of cancer cells by inhibiting the activity of histone deacetylases (HDACs), which are involved in replication and differentiation processes. However, VPA has some disadvantages, among which are poor water solubility and hepatotoxicity. Therefore, the aim of the present study was to design and synthesize three derivatives of VPA to improve its physicochemical properties and anti-proliferative effects. For this purpose, the amino acids aspartic acid, glutamic acid and proline were added to the molecular structure of VPA. Docking and molecular dynamics simulations were used to determine the mode of recognition of these three derivatives by different conformations of HDAC8. This receptor was used as the specific target because of its high affinity for this type of substrate. The results demonstrate that, compared to VPA, the test compounds bind to different sites on the enzyme and that hydrogen bonds and hydrophobic interactions play key roles in this difference. The IC50 values of the VPA derivatives, experimentally determined using HeLa cells, were in the mM range. This result indicates that the derivatives have greater antiproliferative effects than the parent compound. Hence, these results suggest that these amino acid derivatives may represent a good alternative for anticancer treatment.