A DFT/PCM Study on the Affinity of Salinomycin to Bind Monovalent Metal Cations

Dudev, Todor; Cheshmedzhieva, Diana; Dorkov, Peter; Pantcheva, Ivayla

doi:10.3390/molecules27020532

Open AccessArticle

A DFT/PCM Study on the Affinity of Salinomycin to Bind Monovalent Metal Cations

¹

Laboratory of Computational Chemistry and Spectroscopy, Faculty of Chemistry and Pharmacy, “St. Kl. Ohridski” University of Sofia, 1164 Sofia, Bulgaria

²

Research & Development Department, Biovet Ltd., 4550 Peshtera, Bulgaria

³

Laboratory of Biocoordination and Bioanalytical Chemistry, Faculty of Chemistry and Pharmacy, “St. Kl. Ohridski” University of Sofia, 1164 Sofia, Bulgaria

^*

Authors to whom correspondence should be addressed.

Molecules 2022, 27(2), 532; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27020532

Submission received: 23 December 2021 / Revised: 8 January 2022 / Accepted: 11 January 2022 / Published: 14 January 2022

(This article belongs to the Special Issue Metal-Organic Complexes: Applications in Chemistry and Materials Science)

Download

Browse Figures

Review Reports Versions Notes

Abstract

:

The affinity of the polyether ionophore salinomycin to bind IA/IB metal ions was accessed using the Gibbs free energy of the competition reaction between SalNa (taken as a reference) and its rival ions: [M⁺-solution] + [SalNa] → [SalM] + [Na⁺-solution] (M = Li, K, Rb, Cs, Cu, Ag, Au). The DFT/PCM computations revealed that the ionic radius, charge density and accepting ability of the competing metal cations, as well as the dielectric properties of the solvent, have an influence upon the selectivity of salinomycin. The optimized structures of the monovalent metal complexes demonstrate the flexibility of the ionophore, allowing the coordination of one or two water ligands in SalM-W₁ and SalM-W₂, respectively. The metal cations are responsible for the inner coordination sphere geometry, with coordination numbers spread between 2 (Au⁺), 4 (Li⁺ and Cu⁺), 5/6 (Na⁺, K⁺, Ag⁺), 6/7 (Rb⁺) and 7/8 (Cs⁺). The metals’ affinity to salinomycin in low-polarity media follows the order of Li⁺ > Cu⁺ > Na⁺ > K⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺, whereas some derangement takes place in high-dielectric environment: Li⁺ ≥ Na⁺ > K⁺ > Cu⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺.

Keywords:

salinomycin; IA/IB metal ions; monovalent metal complex; DFT/PCM

Graphical Abstract

1. Introduction

Salinomycin is a polyketide ionophore extracted from Streptomyces albus [1]. It exerts an anticoccidial effect against all Eimeria species known as causative agents of coccidiosis in poultry as well as in the other farm animals [2,3,4]. It is also effective in treating bacterial, parasitic and some viral infections [5,6,7]. The biological activity of salinomycin is primarily based on its ability to form lipophilic complexes with the monovalent alkaline ions, making the membranes of target cells permeable to the ions mentioned. In this way, the intracellular cation concentration in the cells is disturbed, leading to a cascade of energy-consuming processes and ultimate cell death.

Recently salinomycin has re-attracted attention in the biomedical field due to its selective inhibition of cancer stem cells in a variety of types of cancer. The exact mode of action of the antibiotic is not yet known, but there are multiple pathways by which it may inhibit tumor growth, including interference with the Akt signaling pathway, Wnt/β-catenin, Hedgehog, and Notch pathways of cancer progression [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22].

Salinomycinic acid (SalH, Scheme 1) is a polyether polyalcohol monocarboxylic acid and its structure determines the ionophoric action of the antibiotic. The polyether chain contains a number of externally oriented alkyl substituents thus defining the overall lipophilic character of the drug. The weak interaction (H-bond formation) between the carboxylic and hydroxyl groups placed on the opposite ends of the molecule causes folding of the structure with the formation of a hydrophilic cavity due to the internally placed polyether oxygen atoms. This core is able to accommodate water molecules (SalH) or monovalent alkali ions (M⁺) preceded by the initial deprotonation of the carboxylic function (SalM). Compared to the polyether ionophorous analog monensin [23,24,25,26,27,28,29,30], only a single structural report on the salinomycin sodium complex (Figure 1) is available [31]. Moreover, although referred to as a potassium ionophore [32], in contrast to the sodium ionophore monensin, no systematic study on salinomycin’s selectivity towards the monovalent metal cations has been found in the literature.

Due to the very limited number of crystallographic structural data on metal complexes of salinomycin and lack of substantial knowledge on its affinity to bind metal cations, here, we highlight the factors which may influence the ability of the antibiotic to form coordination compounds with monovalent metal ions. For the first time, we employed the density functional theory (DFT) calculations combined with polarizable continuum model (PCM) computations to evaluate the free energies of complex formation between the salinomycinate anion and groups IA/IB metal ions. The quantum chemical study brought out the origin of the metal cation (radius, charge, accepting ability) and the dielectric properties of the environment as parameters that affect the coordination properties of the ionophore.

2. Results and Discussion

The metal(I)-loaded constructs of salinomycin were modeled according to the only available crystal structures of sodium salinomycinate (SalNa) [31]. Two types of complex species have been observed experimentally differing by the number of water molecules coordinated to Na⁺ ions in the binding site: a complex with one crystal water (SalNa-W₁) and its counterpart comprising two water molecules orbiting the metal cation (SalNa-W₂). Accordingly, the atomic coordinates of the respective experimental constructs were taken and used in subsequent geometry optimizations. Then, the Na⁺ cation (used as a reference) in the optimized structures was replaced by M⁺ cations and the respective SalM-W_1/2 structures were subjected to succeeding quantum-chemical calculations. The optimized structures of the entire series of SalNa/M–W₁/W₂ complexes are given in Figure 2, Figure 3 and Figure 4.

The affinity of salinomycin towards monovalent metal cations was evaluated modeling the competition Equation (1), where Na⁺, bound to the host ligand, is replaced by its competitor, M⁺ (M = Li, K, Rb, Cs, Cu, Ag, and Au) in an environment with a different dielectric constant ε:

[M⁺-solution] + [SalNa] → [SalM] + [Na⁺-solution].

(1)

[SalNa/M] and [M⁺/Na⁺-solution] represent the metal ion, coordinated to salinomycinate anion, and the corresponding ligand-free metal ion in solution, respectively. The outcome of the rivalry between the monovalent cations from groups IA and IB and their fellow congener Na⁺ for the ionophore was assessed by computing the Gibbs free energy change, ΔG^ε, of Equation (1).

2.1. Salinomycin Complexes with Alkali Metal Ions

Optimized structures of salinomycin complexes with Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺, containing one (W₁) and two (W₂) water molecules in the binding site, are shown in Figure 2 and Figure 3. The host salinomycin donates oxygen atoms from the ether, hydroxyl, carbonyl, and carboxylate groups in forming the metal complexes. The analysis of the structures presented reveals that the host ionophore is flexible and able to adapt to the specific physico-chemical properties of the guest metal ion species. The smallest member of the alkali group, Li⁺, adopts its preferred tetrahedral coordination sphere in both W₁ and W₂ series with a mean Li-O bond distance of 1.984/1.962 Å, respectively. The constructs with the heavier metals and one crystal water form penta- (SalNa-W₁), six- (SalK-W₁, SalRb-W₁), and eight-coordinated (SalCs-W₁) complexes with increased metal-oxygen bond distances going down the group: 2.407, 2.852, 2.985, and 3.211 Å for Na⁺, K⁺, Rb⁺, and Cs⁺, respectively. The main trend—increasing the coordination number (CN) and elongating the metal-oxygen bond distances with increasing the atomic number of the metal—is similar in the two-water SalM series.

However, as compared to the W₁ complexes, the metal coordination numbers differ in SalNa-W₂ (CN = 6), SalRb-W₂ (CN = 7) and SalCs-W₂ (CN = 7). The respective metal-oxygen bond lengths in the W₂ series are 2.511 (SalNa-W₂), 2.814 (SalK-W₂), 3.055 (SalRb-W₂), and 3.264 (SalCs-W₂) Å. In both the W₁ and W₂ series, the M-O bond elongation is compatible with the increased metal cation radius (cation bulkiness) with increasing the atomic number: 0.59 Å for tetra-coordinated lithium, 1.00/1.02 Å for penta-/hexa-coordinated sodium, 1.38 Å for hexa-coordinated potassium, 1.52/1.56 Å for hexa-/hepta-coordinated rubidium and 1.74 Å for octa-coordinated cesium [33] (Table 1). Water molecules stay firmly coordinated to the metal cation in all the complexes.

The thermodynamic outcome of the competition between Na⁺ (serving as a reference cation throughout the paper) and other metal species for binding salinomycin is encoded in the Gibbs free energies of metal substitution, given in Table 2. As seen, the trends of changes in both the W₁ and W₂ series are quite similar with numerical values in most cases very close to each other.

The smallest alkali cation, Li⁺, is the only one that can outcompete Na⁺ (in low dielectric medium for the W₁ series and throughout the entire dielectric region in the W₂ series), evidenced by the negative values of the respective ΔG^ε in Table 2. Note that the competitiveness of Li⁺ is much higher in the gas phase (very negative ΔG¹), where the electronic factors dominate the exchange reaction. However, the high desolvation penalty of the incoming Li⁺ (69, 95, 120 and 123 kcal/mol in ε = 2, 4, 32, and 78, respectively), which is higher than the solvation energy gain upon Na⁺ release (−54, −75, −95, and −98 kcal/mol [35]), attenuates the resulting ΔG^ε in condensed media. Furthermore, heavier alkali metals (K⁺, Rb⁺, and Cs⁺) are weaker competitors of Na⁺ and yield positive free energies of metal exchange which increase down the group. Thus, the affinity of the alkali metal cations for the host ionophore can be arranged in the following order: Li⁺ > Na⁺ > K⁺ > Rb⁺ > Cs⁺ (non-polar solvents) and Li⁺ ≥ Na⁺ > K⁺ > Rb⁺ > Cs⁺ (polar solvents). These results are fully in line with the respective thermodynamic characteristics of these metal species summarized in Table 1, Table 2, Table 3 and Table 4. The “winner” of the group, Li⁺, is characterized with the highest charge density and Lewis acidity, and the best charge accepting abilities. These characteristics gradually weaken for the rest of the cations, following the same affinity order as given above.

2.2. Salinomycin Complexes with Coinage Metal Ions

Competition between Na⁺ and Group IB monovalent cations (Cu⁺, Ag⁺, and Au⁺) was also investigated. The optimized structures of SalCu-W₁/W₂, SalAg-W₁/W₂, and SalAu-W₁/W₂ are shown in Figure 4.

Copper cation forms tetrahedral complexes in both W₁ and W₂ series with average Cu-O bond distances of 2.220/2.211 Å, respectively. Due to its high charge density and Lewis acidity, as well as strong charge accepting ability (Table 1, Table 2, Table 3 and Table 4), Cu⁺ ions outcompetes Na⁺ in both the gas phase and nonpolar solvents such as cyclohexane (ε ≈ 2). In more polar solvents, however, the trend changes: because of solvation effects (see above), Na⁺ complex becomes dominant over its Cu⁺ counterpart in higher-polarity solvents (positive ΔG⁴, ΔG³², and ΔG⁷⁸).

Expectedly, the metal coordination number in silver complexes increases (to 5 in SalAg-W₁ and 6 in SalAg-W₂) with concomitant increase in the respective ΔGs which stay on a positive ground throughout the entire dielectric range (Table 2). In contrast, Au⁺ coordination number sharply decreases to 2 (linear configuration) due to the strong relativistic effects [36]. Data collected in Table 2 imply that, although Au⁺ is more competitive than Ag⁺ (lower ΔGs for the former than the latter), it cannot outcompete Na⁺ (all positive ΔGs). Adding the results evaluated for the coinage cations to those of their alkali counterparts, the order of metals’ affinity to salinomycin in low-polarity solvents becomes: Li⁺ > Cu⁺ > Na⁺ > K⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺, whereas some reordering takes place in high-dielectric media: Li⁺ ≥ Na⁺ > K⁺ > Cu⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺.

2.3. Salinomycin vs. Monensin

Recently we studied the process of complexation of monensin—another representative of ionophore antibiotics, analogous to that of salinomycin—with the same series of monovalent metal cations, where the trends of changes in metal cation affinity to the host ligand were established [35]. Although generally similar, the two ionophores exhibit some variances in their complexation behavior, differing quite significantly in their flexibility. Monensin appears to be much more rigid than salinomycin forcing the incoming metal cations (regardless of their coordination preferences) to adopt a six-coordinated geometry with octahedrally arranged oxygen-containing groups. On the other hand, in salinomycin monovalent metal complexes, it is the cation that dictates the coordination geometry. As seen (above) the coordination numbers of different metals are scattered between 2 (Au⁺), 4 (Li⁺ and Cu⁺), 5/6 (Na⁺, K⁺, Ag⁺), 6/7 (Rb⁺) and 7/8 (Cs⁺). Along the same vein, salinomycin is flexible enough allowing the metal ions to accept and incorporate as a part of their coordination sphere water molecule(s). Differences in coordination pattern reflect on the relative affinity of metal cations to the two ionophores. Although the sequence is the same for non-polar solvents (Li⁺ > Cu⁺ > Na⁺ > K⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺), it is different in high dielectric media: Na⁺ > Li⁺ > K⁺ > Cu⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺ (monensin) and Li⁺ ≥ Na⁺ > K⁺ > Cu⁺ > Au⁺ > Ag⁺ > Rb⁺ > Cs⁺ (salinomycin). The calculations revealed that the change in solvent polarity may favor the binding of certain metal ions and this finding can be explored under laboratory conditions. The growth of crystals suitable for X-ray diffraction studies will bring new insights into the structural chemistry of salinomycin in the solid state.

3. Methods

3.1. Gibbs Free Energy Change for the Reaction Modeled

The cation exchange free energy for Equation (1) in a medium characterized by an effective dielectric constant ε can be calculated as a sum of electronic (ΔG¹) and solvation (ΔΔG^ε_solv) effects (Equation (2)):

ΔG^ε = ΔG¹ + ΔΔG^ε_solv

(2)

where ΔG¹ is the gas-phase free energy for Equation (1) (see below in next subsection), whereas ΔΔG^ε_solv is the difference in the solvation free energies ΔG^ε_solv of the products and reactants in the same process (Equation (3)).

ΔΔG^ε_solv = ΔG^ε_solv ([SalM]) + ΔG^ε_solv ([Na⁺-solution])
− ΔG^ε_solv ([SalNa]) − ΔG^ε_solv ([M⁺-solution]).

(3)

Equation (1) was modeled in different dielectric environments ranging from non-polar solvents such as cyclohexane (ε ≈ 2) and diethyl ether (ε ≈ 4) to highly polar counterparts such as methanol (ε ≈ 32) and water (ε ≈ 78). The positive Gibbs free energy change for Equation (1) suggests a Na⁺-selective ligand while the negative value implies an M⁺-selective ionophore. Note that our aim is to derive reliable trends of changes in the free energy of metal substitution (i.e., ranking the metal cations with respect to their relative affinity to salinomycin), rather than evaluating the absolute free energies of metal exchange. Such an approach has been successfully applied to various metal-containing systems such as monensin [35], macrocyclic cage molecules [37], model ion channels [38,39] and metalloproteins [40,41].

3.2. DFT/PCM Calculations

The Gaussian 09 suite of programs [42] was employed in performing the calculations. All the metal constructs were fully optimized in the gas phase at the B3LYP/6-31+G(d,p) level of theory producing the respective electronic energies, E_elect. The SDD basis set and the effective core potential were used for the heavier cations from the group IA (Rb⁺ and Cs⁺) and IB (Ag⁺ and Au⁺). This combination of method/basis set was selected for the current computations as it proved reliable and performed very well in reproducing the experimental structures of a series of Na⁺/M⁺ complexes similar to salinomycin ionophore—monensin, studied in our recent work [35]. The vibrational frequencies of the metal complexes were evaluated at the same level of theory. No imaginary frequency was detected for any of the optimized constructs signifying a local minimum of the potential energy surface. The vibrational frequencies were employed in evaluating the thermal energies, including the zero-point energy, E_T, and entropy, S. These values were used in estimating the gas phase Gibbs free energies for Equation (1), ΔG¹, at room temperature, T = 298.15 K, according to Equation (4):

ΔG¹ = ΔE_elect + ΔE_T − TΔS

(4)

where ΔE_elect, ΔE_T, and ΔS are the respective differences between the products and reactants. The solvation effects were assessed by employing polarizable continuum model calculations utilizing SMD scheme [43] as implemented in Gaussian 09.

Single point calculations of each gas-phase optimized structure were conducted in each solvent. The differences between the gas phase and solution energies yielded the free energies of solvation of the complex, ΔG^ε_solv (Table 5 and Table 6). Notably, the experimentally determined solvation free energies for the metal cations in aqueous solution (Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, and Ag⁺) were used in the computations. For other cations (Cu⁺ and Au⁺) and solvents, the estimated values were used, where the ratios in the theoretically evaluated quantities along with the experimental values were used to determine the respective ΔG^ε_solv of [M⁺/Na⁺-solution] [35]. Of note, such a calculation protocol utilizing a thermodynamic cycle that employs the experimental solvation free energies of some of participating entities (where available), has proven quite dependable in reproducing the experimental thermodynamic data [44]. The basis set superposition error (BSSE) for this type of exchange reaction (Equation (1)) was found to be negligible [38] and, therefore, was not considered in the current calculations.

4. Conclusions

The quantum chemical study performed disclosed the nature of the monovalent metal ions as the major factor in governing the affinity of salinomycinate anion towards group IA and IB cations. The smaller ions possessing higher positive charge density and higher ligand preference (Li⁺ and Cu⁺) predominantly coordinate in an environment with lower dielectric properties. In polar solvents, Li⁺ ions are superior to Na⁺, while the heavier alkali metal ions and their group IB counterparts are weaker competitors. The internal cavity of salinomycinate ion appears to be very flexible accommodating one or two water molecules that occupy, in addition to salinomycin binding groups, the inner coordination sphere of the metal center. The latter is crucial for the complex geometry where the metal coordination number varies from two in SalAu-W₁/W₂ to eight in SalCs-W₁.

Author Contributions

Conceptualization, I.P. and T.D.; methodology and analysis, T.D. and D.C.; writing—original draft preparation, review and editing, I.P., T.D., D.C. and P.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research and APC were funded by the Bulgarian National Science Fund (Grant number KP-06-H29/3).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available from the authors upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not available.

References

Miyazaki, Y.; Shibuya, M.; Sugawara, H.; Kawaguchi, O.; Hirsoe, C. Salinomycin, a new polyether antibiotic. J. Antibiot. 1974, 27, 814–821. [Google Scholar] [CrossRef]
Chappel, L.R. The site of action of the anticoccidial salinomycin (Coxistac). J. Parasitol. 1979, 65, 137–143. [Google Scholar] [CrossRef] [PubMed]
Long, P.L.; Jeffers, T.K. Studies on the stage of action of ionophorous antibiotics against Eimeria. J. Parasitol. 1982, 68, 363–371. [Google Scholar] [CrossRef] [PubMed]
Coombs, G.H.; Müller, S. Recent advances in the search for new anti-coccidial drugs. Int. J. Parasitol. 2002, 32, 497–508. [Google Scholar] [CrossRef]
Mahmoudi, N.; de Julian-Ortiz, J.-V.; Ciceron, L.; Galvez, J.; Mazier, D.; Danis, M.; Derouin, F.; Garcıa-Domenech, R. Identification of new antimalarial drugs by linear discriminant analysis and topological virtual screening. J. Antimicrob. Chemother. 2006, 57, 489–497. [Google Scholar] [CrossRef] [PubMed]
D’Alessandro, S.; Corbett, Y.; Ilboudo, D.P.; Misiano, P.; Dahiya, N.; Abay, S.M.; Habluetzel, A.; Grande, R.; Gismondo, M.R.; Dechering, K.J.; et al. Salinomycin and other ionophores as a new class of antimalarial drugs with transmission-blocking activity. Antimicrob. Agents Chemother. 2015, 59, 5135–5144. [Google Scholar] [CrossRef] [Green Version]
Yuan, C.; Huang, X.; Zhai, R.; Ma, Y.; Xu, A.; Zhang, P.; Yang, Q. In Vitro antiviral activities of salinomycin on porcine epidemic diarrhea virus. Viruses 2021, 13, 580. [Google Scholar] [CrossRef]
Gupta, P.B.; Onder, T.T.; Jiang, G.; Tao, K.; Kuperwasser, C.; Weinberg, R.A.; Lander, E.S. Identification of selective inhibitors of cancer stem cells by high-throughput screening. Cell 2009, 138, 645–659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Lu, D.; Choi, M.Y.; Yu, J.; Castro, J.E.; Kipps, T.J.; Carson, D.A. Salinomycin inhibits Wnt signaling and selectively induces apoptosis in chronic lymphocytic leukemia cells. Proc. Natl. Acad. Sci. USA 2011, 108, 13253–13257. [Google Scholar] [CrossRef] [Green Version]
Tang, Q.L.; Zhao, Z.Q.; Li, J.C.; Liang, Y.; Yin, J.Q.; Zou, C.Y.; Xie, X.B.; Zeng, Y.X.; Shen, J.N.; Kang, T.; et al. Salinomycin inhibits osteosarcoma by targeting its tumor stem cells. Cancer Lett. 2011, 311, 113–121. [Google Scholar] [CrossRef]
Naujokat, C.; Steinhart, R. Salinomycin as a drug for targeting human cancer stem cells. J. Biomed. Biotechnol. 2012, 2012, 950658. [Google Scholar] [CrossRef] [PubMed]
Zhou, S.; Wang, F.; Wong, E.T.; Fonkem, E.; Hsieh, T.C.; Wu, J.M.; Wu, E. Salinomycin: A novel anti-cancer agent with known anti-coccidial activities. Curr. Med. Chem. 2013, 20, 4095–4101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Zhang, C.; Tian, Y.; Song, F.; Fu, C.; Han, B.; Wang, Y. Salinomycin inhibits the growth of colorectal carcinoma by targeting tumor stem cells. Oncol. Rep. 2015, 34, 2469–2476. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Najumudeen, A.K.; Jaiswal, A.; Lectez, B.; Oetken-Lindholm, C.; Guzman, C.; Siljamaki, E.; Posada, I.M.D.; Lacey, E.; Aittkallio, T.; Abankwa, D. Cancer stem cell drugs target K-ras signaling in a stemness context. Oncogene 2016, 35, 5248–5262. [Google Scholar] [CrossRef] [Green Version]
Dewangan, J.; Srivastava, S.; Rath, S.K. Salinomycin: A new paradigm in cancer therapy. Tumour. Biol. 2017, 39, 1010428317695035. [Google Scholar] [CrossRef] [Green Version]
Mai, T.T.; Hamai, A.; Hienzsch, A.; Caneque, T.; Muller, S.; Wicinski, J.; Cabaud, O.; Leroy, C.; David, A.; Acevedo, V.; et al. Salinomycin kills cancer stem cells by sequestering iron in lysosomes. Nat. Chem. 2017, 9, 1025–1033. [Google Scholar] [CrossRef] [Green Version]
Lee, H.G.; Shin, S.J.; Chung, H.W.; Kwon, S.H.; Cha, S.D.; Lee, J.E.; Cho, C.H. Salinomycin reduces stemness and induces apoptosis on human ovarian cancer stem cell. J. Gynecol. Oncol. 2017, 28, e14. [Google Scholar] [CrossRef] [Green Version]
Liu, L.; Wang, Q.; Mao, J.; Qin, T.; Sun, Y.; Yang, J.; Han, Y.; Li, L.; Li, Q. Salinomycin suppresses cancer cell stemness and attenuates TGF-beta-induced epithelial-mesenchymal transition of renal cell carcinoma cells. Chem. Biol. Interact. 2018, 296, 145–153. [Google Scholar] [CrossRef]
Zhou, J.; Liu, S.; Wang, Y.; Dai, W.; Zou, H.; Wang, S.; Zhang, J.; Pan, J. Salinomycin effectively eliminates cancer stem-like cells and obviates hepatic metastasis in uveal melanoma. Mol. Cancer 2019, 18, 159. [Google Scholar] [CrossRef]
Liu, Y.; Hao, Y.; Li, Y.; Zheng, Y.; Dai, J.; Zhong, F.; Wei, W.; Fang, Z. Salinomycin induces autophagic cell death in salinomycin-sensitive melanoma cells through inhibition of autophagic flux. Sci. Rep. 2020, 10, 18515. [Google Scholar] [CrossRef]
Li, K.; Pang, L.; Pan, X.; Fan, S.; Wang, X.; Wang, Q.; Dai, P.; Gao, W.; Gao, J. GE11 modified PLGA/TPGS nanoparticles targeting delivery of salinomycin to breast cancer cells. Technol. Cancer Res. Treat. 2021, 20, 15330338211004954. [Google Scholar] [CrossRef] [PubMed]
Wang, H.; Zhang, H.; Zhu, Y.; Wu, Z.; Cui, C.; Cai, F. Anticancer mechanisms of salinomycin in breast cancer and its clinical applications. Front. Oncol. 2021, 11, 654428. [Google Scholar] [CrossRef] [PubMed]
Pinkerton, M.; Stenrauf, L.K. Molecular structure of monovalent metal cation complexes of monensin. J. Mol. Biol. 1970, 49, 533–546. [Google Scholar] [CrossRef]
Lutz, W.K.; Winkler, F.K.; Dunitz, J.D. Crystal structure of the antibiotic monensin. Similarities and differences between free acid and metal complex. Helv. Chim. Acta 1971, 54, 1103–1108. [Google Scholar] [CrossRef]
Duax, W.L.; Smith, G.D.; Strong, P.D. Complexation of metal ions by monensin. Crystal and molecular structure of hydrated and anhydrous crystal forms of sodium monensin. J. Am. Chem. Soc. 1980, 102, 6725–6729. [Google Scholar] [CrossRef]
Walba, D.M.; Hermsmeier, M.; Haltiwanger, R.C.; Noordik, J.H. Crystal structures of monensin B lithium and silver salts. J. Org. Chem. 1986, 51, 245–247. [Google Scholar] [CrossRef]
Pangborn, W.; Duax, W.; Langs, D. The hydrated potassium complex of the ionophore monensin A. J. Am. Chem. Soc. 1987, 109, 2163–2165. [Google Scholar] [CrossRef]
Huczynski, A.; Ratajczak-Sitarz, M.; Katrusiak, A.; Brzezinski, B. Molecular structure of the 1:1 inclusion complex of monensin A sodium salt with acetonitrile. J. Mol. Struct. 2007, 832, 84–89. [Google Scholar] [CrossRef]
Huczynski, A.; Ratajczak-Sitarz, M.; Katrusiak, A.; Brzezinski, B. Molecular structure of the 1:1 inclusion complex of monensin A lithium salt with acetonitrile. J. Mol. Struct. 2007, 871, 92–97. [Google Scholar] [CrossRef]
Huczynski, A.; Ratajczak-Sitarz, M.; Katrusiak, A.; Brzezinski, B. Molecular structure of rubidium six-coordinated dihydrate complex with monensin A. J. Mol. Struct. 2008, 883, 224–229. [Google Scholar] [CrossRef]
Paulus, E.F.; Kurz, M.; Matter, H.; Vertesy, L. Solid-state and solution structure of the salinomycin-sodium complex: Stabilization of different conformers for an ionophore in different environments. J. Am. Chem. Soc. 1998, 120, 8209–8221. [Google Scholar] [CrossRef]
Mitani, M.; Yamanishi, T.; Miyazaki, Y. Salinomycin: A new monovalent cation ionophore. Biochem. Biophys Res. Comm. 1975, 66, 12311236. [Google Scholar] [CrossRef]
Shannon, R.D. Revised effective ionic radii and systematic studies in interatomic distances in halides and chalcogenides. Acta Cryst. 1976, A32, 751–767. [Google Scholar] [CrossRef]
Gagne, O.C.; Hawthorne, F.C. Empirical Lewis acid strengths for 135 cations bonded to oxygen. Acta Cryst. B 2017, 73, 956–961. [Google Scholar] [CrossRef] [PubMed]
Dudev, T.; Cheshmedzhieva, D.; Dimitrova, R.; Dorkov, P.; Pantcheva, I. Factors governing the competition between group IA and IB cations for monensin A: A DFT/PCM study. RSC Adv. 2020, 10, 5734–5741. [Google Scholar] [CrossRef]
Laguna, A. Modern Supramolecular Gold Chemistry: Gold-Metal Interactions and Applications, 1st ed.; WILEY-VCH Verlag GmbH & Co. KgaA: Weinheim, Germany, 2008. [Google Scholar]
Angelova, S.E.; Nikolova, V.K.; Dudev, T.M. Determinants of the host-guest interactions between α-, β- and γ-cyclodextrins and group IA, IIA and IIIA metal cations: A DFT/PCM study. Phys. Chem. Chem. Phys. 2017, 19, 15129–15136. [Google Scholar] [CrossRef]
Dudev, T.; Lim, C. Determinants of K⁺ vs. Na⁺ selectivity in potassium channels. J. Am. Chem. Soc. 2009, 131, 8092–8101. [Google Scholar] [CrossRef] [PubMed]
Dudev, T.; Lim, C. Ion selectivity strategies of sodium channel selectivity filters. Acc. Chem. Res. 2014, 47, 3580–3587. [Google Scholar] [CrossRef] [PubMed]
Dudev, T.; Cheshmedzhieva, D.; Doudeva, L. Competition between abiogenic Al³⁺ and native Mg²⁺, Fe²⁺ and Zn²⁺ ions in protein binding sites: Implications for aluminum toxicity. J. Mol. Model. 2018, 24, 55. [Google Scholar] [CrossRef]
Dudev, T.; Grauffel, C.; Lim, C. How Pb²⁺ binds and modulates properties of Ca²⁺-signaling proteins. Inorg. Chem. 2018, 57, 14798–14809. [Google Scholar] [CrossRef]
Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Montgomery, J.A., Jr.; Vreven, T.; Kudin, K.N.; Burant, J.C.; et al. Gaussian 09; Gaussian Inc.: Pittsburgh, PA, USA, 2009. [Google Scholar]
Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 2009, 113, 63786396. [Google Scholar] [CrossRef] [PubMed]
Kircheva, N.; Dudev, T. Novel insights into gallium’s mechanism of therapeutic action: A DFT/PCM study of the interaction between Ga³⁺ and ribonucleotide reductase substrates. J. Phys. Chem. B 2019, 123, 5444–5451. [Google Scholar] [CrossRef] [PubMed]

Scheme 1. Chemical structure of salinomycinic acid with the numbering sequence of oxygen atoms.

Figure 1. Crystal structure of sodium salinomycinate (SalNa) containing (a) one and (b) two water ligands [31]. Color scheme: C—grey, O—red, Na—purple (hydrogens are omitted for clarity).

Figure 2. B3LYP/6-31+G(d,p) fully optimized structures of salinomycin complexes with Li⁺, Na⁺, K⁺, Rb⁺ and Cs⁺ cations, containing one (W₁) water molecule. Color scheme: C—green, O—red, H—light grey, Li—magenta, Na—purple, K—blue, Rb—yellow, Cs—deep olive (carbon hydrogens are omitted for clarity).

Figure 3. B3LYP/6-31+G(d,p) fully optimized structures of salinomycin complexes with Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺ cations, containing two (W₂) water molecules. Color scheme: C—green, O—red, H—light grey, Li—magenta, Na—purple, K—blue, Rb—yellow, Cs—deep olive (carbon hydrogens are omitted for clarity).

Figure 4. B3LYP/6-31+G(d,p) fully optimized structures of salinomycin complexes with Cu⁺, Ag⁺, and Au⁺, containing (a) one (W₁) and (b) two (W₂) water molecules. Color scheme: C—green, O—red, H—light grey, Cu—orange, Ag—light blue, Au—gray (carbon hydrogens are omitted for clarity).

Table 1. Structural and thermodynamic characteristics of alkali and coinage metal cations.

Cation	Ionic Radius (r, Å) [33] ¹	Charge Density ² (e/Å³)	Lewis Acidity (Valence Units) [34]	Charge Transfer to the Metal ³ (e)
Li⁺	0.59 (IV)	1.16	0.215	0.286/0.260
Na⁺	1.00/1.02 (V/VI)	0.24/0.22	0.159	0.215/0.198
K⁺	1.38 (VI)	0.09	0.108	0.166/0.152
Rb⁺	1.52/1.56 (VI/VII)	0.07/0.06	0.099	0.142/0.132
Cs⁺	-/1.74 (VII/VIII)	-/0.04	0.084	0.129/0.119
Cu⁺	0.60 (IV)	1.10	0.400	0.360/0.338
Ag⁺	1.09/1.15 (V/VI)	0.18/0.16	0.191	0.300/0.281
Au⁺	-/1.37 (II/VI)	-/0.09	-	0.412/0.381

¹ The respective coordination number(s) of the metal is/are given in parentheses; ² Charge density = Ion charge/Ion volume = 1/(4/3 π r³); ³ Charge transfer from CH₃OCH₃ (model for an ether group) and CH₃OH (model for a hydroxyl group) evaluated at B3LYP/6-31+G(d,p) level [35].

Table 2. Gibbs free energies of M⁺/Na⁺ metal substitution (ΔG^ε, kcal/mol) for Equation (1), evaluated in media with different dielectric constant ε.

Reaction	ΔG^ε	Reaction	ΔG^ε
[Li⁺-solution] + [SalNa-W₁] → [SalLi-W₁] + [Na⁺-solution]	ΔG¹ = −21.6 ΔG² = −7.8 ΔG⁴ = −3.2 ΔG³² = 1.5 ΔG⁷⁸ = 1.4	[Li⁺-solution] + [SalNa-W₂] → [SalLi-W₂] + [Na⁺-solution]	ΔG¹ = −24.0 ΔG² = −10.6 ΔG⁴ = −6.5 ΔG³² = −1.9 ΔG⁷⁸ = −1.7
[K⁺-solution] + [SalNa-W₁] → [SalK-W₁] + [Na⁺-solution]	ΔG¹ = 19.8 ΔG² = 9.9 ΔG⁴ = 6.8 ΔG³² = 2.8 ΔG⁷⁸ = 2.8	[K⁺-solution] + [SalNa-W₂] → [SalK-W₂] + [Na⁺-solution]	ΔG¹ = 20.6 ΔG² = 10.0 ΔG⁴ = 6.3 ΔG³² = 2.7 ΔG⁷⁸ = 2.7
[Rb⁺-solution] + [SalNa-W₁] → [SalRb-W₁] + [Na⁺-solution]	ΔG¹ = 58.2 ΔG² = 46.2 ΔG⁴ = 41.6 ΔG³² = 34.3 ΔG⁷⁸ = 33.2	[Rb⁺-solution] + [SalNa-W₂] → [SalRb-W₂] + [Na⁺-solution]	ΔG¹ = 56.1 ΔG² = 43.5 ΔG⁴ = 38.2 ΔG³² = 32.8 ΔG⁷⁸ = 31.2
[Cs⁺-solution] + [SalNa-W₁] → [SalCs-W₁] + [Na⁺-solution]	ΔG¹ = 71.2 ΔG² = 55.7 ΔG⁴ = 49.7 ΔG³² = 41.4 ΔG⁷⁸ = 40.4	[Cs⁺-solution] + [SalNa-W₂] → [SalCs-W₂] + [Na⁺-solution]	ΔG¹ = 67.3 ΔG² = 51.1 ΔG⁴ = 44.7 ΔG³² = 38.4 ΔG⁷⁸ = 36.7
[Cu⁺-solution] + [SalNa-W₁] → [SalCu-W₁] + [Na⁺-solution]	ΔG¹ = −34.2 ΔG² = −5.8 ΔG⁴ = 4.3 ΔG³² = 9.2 ΔG⁷⁸ = 18.0	[Cu⁺-solution] + [SalNa-W₂] → [SalCu-W₂] + [Na⁺-solution]	ΔG¹ = −31.6 ΔG² = −2.9 ΔG⁴ = 7.2 ΔG³² = 12.0 ΔG⁷⁸ = 11.4
[Ag⁺-solution] + [SalNa-W₁] → [SalAg-W₁] + [Na⁺-solution]	ΔG¹ = 15.8 ΔG² = 25.0 ΔG⁴ = 28.2 ΔG³² = 29.9 ΔG⁷⁸ = 29.6	[Ag⁺-solution] + [SalNa-W₂] → [SalAg-W₂] + [Na⁺-solution]	ΔG¹ = 16.0 ΔG² = 25.3 ΔG⁴ = 28.4 ΔG³² = 30.1 ΔG⁷⁸ = 29.5
[Au⁺-solution] + [SalNa-W₁] → [SalAu-W₁] + [Na⁺-solution]	ΔG¹ = 5.2 ΔG² = 20.0 ΔG⁴ = 24.8 ΔG³² = 27.5 ΔG⁷⁸ = 26.8	[Au⁺-solution] + [SalNa-W₂] → [SalAu-W₂] + [Na⁺-solution]	ΔG¹ = 0.7 ΔG² = 14.5 ΔG⁴ = 18.5 ΔG³² = 20.1 ΔG⁷⁸ = 18.0

Table 3. Atom charges after the charge transfer (CT) from ligands to the metal cation in salinomycin complexes with one water molecule from Hirshfeld, CM5 and NBO population analyses (in e).

Ion	q^hirsh	q^CM5	q^NBO
Li⁺	0.1746	0.5007	0.8892
Na⁺	0.2600	0.6675	0.8993
K⁺	0.2773	0.7816	0.9112
Rb⁺	0.3023	0.8319	0.9284
Cs⁺	0.3293	0.9555	0.9454
Cu⁺	0.2780	0.4605	0.7379
Ag⁺	0.3224	0.4651	0.6711
Au⁺	0.2805	0.4202	0.5418

Table 4. Atom charges after the charge transfer (CT) from ligands to the metal cation in salinomycin complexes with two water molecules from Hirshfeld, CM5 and NBO population analyses (in e).

Ion	q^hirsh	q^CM5	q^NBO
Li⁺	0.1725	0.5074	0.9005
Na⁺	0.2474	0.6128	0.9121
K⁺	0.2797	0.7392	0.9277
Rb⁺	0.3263	0.7472	0.9395
Cs⁺	0.3546	0.8293	0.9519
Cu⁺	0.2411	0.4046	0.6690
Ag⁺	0.3182	0.4450	0.6687
Au⁺	0.2562	0.3946	0.4588

Table 5. Gibbs free energies of solvation evaluated from PCM calculations for salinomycin complexes containing one water molecule in solvents with different dielectric constants (in kcal/mol).

SalM-W₁	ΔG^solv
SalM-W₁	Cyclohexane	Diethylether	Methanol	Water
Li	−19.97	−27.48	−39.73	−24.86
Na	−19.10	−26.06	−37.70	−22.72
K	−18.49	−25.17	−36.98	−22.24
Rb	−19.29	−26.42	−40.78	−26.00
Cs	−19.42	−26.75	−40.97	−26.20
Cu	−19.73	−27.12	−44.50	−20.88
Ag	−19.49	−26.74	−39.99	−25.21
Au	−19.82	−27.36	−42.13	−27.63

ΔG^solv = E_solvent − E_gas-phase.

Table 6. Gibbs free energies of solvation evaluated from PCM calculations for salinomycin complexes containing two water molecules in solvents with different dielectric constants (in kcal/mol).

SalM-W₂	ΔG^solv
SalM-W₂	Cyclohexane	Diethylether	Methanol	Water
Li	−20.50	−28.79	−45.65	−30.76
Na	−19.23	−26.52	−43.59	−27.90
K	−19.31	−26.91	−43.83	−28.26
Rb	−20.05	−28.07	−46.11	−31.25
Cs	−20.17	−28.27	−45.98	−31.21
Cu	−19.53	−27.17	−50.19	−35.30
Ag	−19.51	−27.09	−45.93	−30.65
Au	−20.50	−29.01	−48.56	−34.44

ΔG^solv = E_solvent − E_gas-phase.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Dudev, T.; Cheshmedzhieva, D.; Dorkov, P.; Pantcheva, I. A DFT/PCM Study on the Affinity of Salinomycin to Bind Monovalent Metal Cations. Molecules 2022, 27, 532. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27020532

AMA Style

Dudev T, Cheshmedzhieva D, Dorkov P, Pantcheva I. A DFT/PCM Study on the Affinity of Salinomycin to Bind Monovalent Metal Cations. Molecules. 2022; 27(2):532. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27020532

Chicago/Turabian Style

Dudev, Todor, Diana Cheshmedzhieva, Peter Dorkov, and Ivayla Pantcheva. 2022. "A DFT/PCM Study on the Affinity of Salinomycin to Bind Monovalent Metal Cations" Molecules 27, no. 2: 532. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27020532

Article Menu

A DFT/PCM Study on the Affinity of Salinomycin to Bind Monovalent Metal Cations

Abstract

1. Introduction

2. Results and Discussion

2.1. Salinomycin Complexes with Alkali Metal Ions

2.2. Salinomycin Complexes with Coinage Metal Ions

2.3. Salinomycin vs. Monensin

3. Methods

3.1. Gibbs Free Energy Change for the Reaction Modeled

3.2. DFT/PCM Calculations

4. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

Sample Availability

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI