Anticancer mechanism of 7-α-hydroxyfrullanolide on microtubules and computational prediction of its target binding in triple-negative breast cancer cells

Background Triple-negative breast cancer (TNBC) responds poorly to the available drugs; thus, the mortality rate associated with TNBC remains high. 7-α-Hydroxyfrullanolide (7HF) possesses anticancer properties and arrests cells in the G2/M-phase via modulation of several proteins involved in the G2/M-phase transition, as well as the mitotic checkpoint in MDA-MB-468 (TNBC) cells. Microtubules (MTs) dynamically regulate cell division in the G2/M phase and are related to cancer cell stress response. However, antimitotic drug cytotoxicity to multiple cancer resistance developed in response to drugs are obstacles faced to date. Here, the activity and mechanism via which 7HF controls MTs dynamics was investigated in MDA-MB-468 cells. Methods 7HF uptake by MDA-MB-468 cells was assessed using spectrophotometry. The drug-like properties of 7HF were predicted using the Swiss-absorption, distribution, metabolism, and excretion (ADME) webtool. Then, the effect of 7HF treatment (6, 12, and 24 µM) on the dynamic arrangement of MTs was assessed for 1, 12, and 24 h using indirect immunofluorescence. Polymerization of α- and β-tubulin was assessed using different 7HF concentrations in a cell-free system for 1 h. Cell proliferation assay with bromodeoxyuridine plus propidium iodide staining and flow cytometry was performed at different 7HF concentrations and time points. The mechanism of action was assessed by detecting the expression of proteins, including Bub3, cyclin B1, p-Cdk1 (Tyr15), Rb, p-Rb (Ser780), Chk1, p-Chk1 (Ser345), Chk2, p-Chk2 (Ser516), and p-H2AX (Ser139), using western blotting. Molecular docking was used to predict the molecular interactions between 7HF and tubulins in MTs. Results We observed that 7HF was able to enter the MDA-MB-468 cells. The ADME webtool analysis predicted that it possesses the high passive permeation and gastrointestinal absorption properties of drugs. Various concentrations of 7HF disrupted the dynamic arrangement of spindle MTs by causing radial spindle array shrinkage and expansion of fibrous spindle density and radial array lengths in a time-dependent manner. 7HF reduced polymerization of α-, β-tubulin in dose-dependent manner. 7HF also triggered DNA damage response by inducing G2/M and G1 phase arrests in a concentration and time-dependent manner, which occurred due to the upregulation of Bub3, Chk1, p-Chk1 (Ser345), p-Cdk1 (Tyr15), and cyclin B1. According to molecular docking analysis, 7HF preferred to bind to β-tubulin over α-tubulin. The lactone, ketone, and hydroxyl groups of 7HF supported the 7HF-tubulin interactions. Hydrogen bonding with a hydrocarbon ring and salt bridge attractive forces were responsible for the binding versatility of 7HF. Conclusions This is the first study to investigate the molecular mechanism, MTs interacting sites, and the internalization and drug-like properties of 7HF in TNBC cells. The findings will be useful for developing 7HF-based treatment for patients with TNBC

Free Energy Simulations of a GTPase: GTP and GDP Binding to Archaeal Initiation Factor 2

International audienceArchaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, "ON" conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg(2+) ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg(2+):phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg(2+):phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2

Reorganization Energy for Internal Electron Transfer in Multicopper Oxidases.

We have calculated the reorganization energy for the intramolecular electron transfer between the reduced type 1 copper site and the peroxy intermediate of the trinuclear cluster in the multicopper oxidase CueO. The calculations are performed at the combined quantum mechanics and molecular mechanics (QM/MM) level, based on molecular dynamics simulations with tailored potentials for the two copper sites. We obtain a reorganization energy of 91-133 kJ/mol, depending on the theoretical treatment. The two Cu sites contribute by 12 and 22 kJ/mol to this energy, whereas the solvent contribution is 34 kJ/mol. The rest comes from the protein, involving small contributions from many residues. We have also estimated the energy difference between the two electron-transfer states and show that the reduction of the peroxy intermediate is exergonic by 43-87 kJ/mol, depending on the theoretical method. Both the solvent and the protein contribute to this energy difference, especially charged residues close to the two Cu sites. We compare these estimates with energies obtained from QM/MM optimizations and QM calculations in a vacuum and discuss differences between the results obtained at various levels of theory

Effects of rhodomyrtone on Gram-positive bacterial tubulin homologue FtsZ

Rhodomyrtone, a natural antimicrobial compound, displays potent activity against many Gram-positive pathogenic bacteria, comparable to last-defence antibiotics including vancomycin and daptomycin. Our previous studies pointed towards effects of rhodomyrtone on the bacterial membrane and cell wall. In addition, a recent molecular docking study suggested that the compound could competitively bind to the main bacterial cell division protein FtsZ. In this study, we applied a computational approach (in silico), in vitro, and in vivo experiments to investigate molecular interactions of rhodomyrtone with FtsZ. Using molecular simulation, FtsZ conformational changes were observed in both (S)- and (R)-rhodomyrtone binding states, compared with the three natural states of FtsZ (ligand-free, GDP-, and GTP-binding states). Calculations of free binding energy showed a higher affinity of FtsZ to (S)-rhodomyrtone (−35.92 ± 0.36 kcal mol−1) than the GDP substrate (−23.47 ± 0.25 kcal mol−1) while less affinity was observed in the case of (R)-rhodomyrtone (−18.11 ± 0.11 kcal mol−1). In vitro experiments further revealed that rhodomyrtone reduced FtsZ polymerization by 36% and inhibited GTPase activity by up to 45%. However, the compound had no effect on FtsZ localization in Bacillus subtilis at inhibitory concentrations and cells also did not elongate after treatment. Higher concentrations of rhodomyrtone did affect localization of FtsZ and also affected localization of its membrane anchor proteins FtsA and SepF, showing that the compound did not specifically inhibit FtsZ but rather impaired multiple divisome proteins. Furthermore, a number of cells adopted a bean-like shape suggesting that rhodomyrtone possibly possesses further targets involved in cell envelope synthesis and/or maintenance

How the Reorganization Free Energy Affects the Reduction Potential of Structurally Homologous Cytochromes

Differences in the reduction potential E-0 among structurally similar metalloproteins cannot always be fully explained on the basis of their 3-D structures. We investigate the molecular determinants to E-0 using the mixed quantum mechanics/molecular mechanics approach named perturbed matrix method (PMM); after comparison with experimental values to assess the reliability of our calculations, we investigate the relationship between the change in free energy upon reduction Delta A(0) and the reorganization energy. We find that the reduction potential of cytochromes can be regarded as the result of the sum of two terms, one being mostly dependent on the energy fluctuations within a limited range around the mean transition energy and the second being mostly dependent linearly on the difference Delta lambda = lambda(red) - lambda(ox) of the reorganization free energies for the ox -> red (lambda(red)) and for the red -> ox (lambda(ox)) relaxations