Next Article in Journal
Binding Studies of Caffeic and p-Coumaric Acid with α-Amylase: Multispectroscopic and Computational Approaches Deciphering the Effect on Advanced Glycation End Products (AGEs)
Next Article in Special Issue
Novel Sterically Crowded and Conformationally Constrained α-Aminophosphonates with a Near-Neutral pKa as Highly Accurate 31P NMR pH Probes. Application to Subtle pH Gradients Determination in Dictyostelium discoideum Cells
Previous Article in Journal
Sustainable Protocol for the Synthesis of 2′,3′-Dideoxynucleoside and 2′,3′-Didehydro-2′,3′-dideoxynucleoside Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Acid–Base Equilibrium and Self-Association in Relation to High Antitumor Activity of Selected Unsymmetrical Bisacridines Established by Extensive Chemometric Analysis

by
Michał Kosno
1,
Tomasz Laskowski
1,*,
Joanna E. Frackowiak
1,
Agnieszka Potęga
1,
Agnieszka Kurdyn
1,
Witold Andrałojć
2,
Julia Borzyszkowska-Bukowska
1,
Katarzyna Szwarc-Karabyka
3 and
Zofia Mazerska
1,*
1
Department of Pharmaceutical Technology and Biochemistry and BioMedTech Centre, Faculty of Chemistry, Gdańsk University of Technology, Gabriela Narutowicza Str. 11/12, 80-233 Gdańsk, Poland
2
Institute of Bioorganic Chemistry, Polish Academy of Sciences, Zygmunta Noskowskiego Str. 12/14, 61-704 Poznań, Poland
3
Nuclear Magnetic Resonance Laboratory, Gdańsk University of Technology, Gabriela Narutowicza Str. 11/12, 80-233 Gdańsk, Poland
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(13), 3995; https://doi.org/10.3390/molecules27133995
Submission received: 3 June 2022 / Revised: 18 June 2022 / Accepted: 19 June 2022 / Published: 21 June 2022
(This article belongs to the Special Issue NMR Spectroscopy in Drug Discovery Research)

Abstract

:
Unsymmetrical bisacridines (UAs) represent a novel class of anticancer agents previously synthesized by our group. Our recent studies have demonstrated their high antitumor potential against multiple cancer cell lines and human tumor xenografts in nude mice. At the cellular level, these compounds affected 3D cancer spheroid growth and their cellular uptake was selectively modulated by quantum dots. UAs were shown to undergo metabolic transformations in vitro and in tumor cells. However, the physicochemical properties of UAs, which could possibly affect their interactions with molecular targets, remain unknown. Therefore, we selected four highly active UAs for the assessment of physicochemical parameters under various pH conditions. We determined the compounds’ pKa dissociation constants as well as their potential to self-associate. Both parameters were determined by detailed and complex chemometric analysis of UV-Vis spectra supported by nuclear magnetic resonance (NMR) spectroscopy. The obtained results indicate that general molecular properties of UAs in aqueous media, including their protonation state, self-association ratio, and solubility, are strongly pH-dependent, particularly in the physiological pH range of 6 to 8. In conclusion, we describe the detailed physicochemical characteristics of UAs, which might contribute to their selectivity towards tumour cells as opposed to their effect on normal cells.

1. Introduction

Unsymmetrical bisacridines (UAs) are a new class of compounds developed by our research group; they are characterized by a unique structure consisting of two heterocyclic ring systems linked by an aminoalkyl chain [1,2]. UAs exhibited high cytotoxic activity against all fourteen tested cancer cell lines as well as antitumor activity against Walker 256 rat adenocarcinoma and ten human tumor xenografts in nude mice. Notably, the compounds which displayed the highest activity strongly inhibited pancreatic cancer cell lines [3]. Studies on the biological effects of these compounds demonstrated their ability to suppress 3D cancer spheroid growth [4]. Moreover, their anticancer activity was enhanced when bound to quaternary quantum dots, resulting in selective upregulation of their cellular uptake [5,6,7]. Our previous studies demonstrated that UAs undergo metabolic transformations which retain the compounds’ dimeric structure. These include products of phase I and phase II metabolism such as N-oxide derivatives in the aminoalkyl side chain through flavin monooxygenase (FMO) transformations as well as uridine UDP-glucuronosyltransferase (UGT)-mediated glucuronide conjugates at the hydroxyl group in the imidazoacridinone ring, respectively [8]. We identified products resulting from glutathione S-transferase-mediated glutathione conjugation [9].
Our recent studies indicated that certain structural characteristics of UA dimers (Figure 1) which distinguish them from acridine monomers (Figure 2) might reflect their unique physicochemical properties under physiological conditions [10,11,12]. Although dimer compounds indeed consist of an imidazoacridinone and 1-nitroacridine monomer connected by a flexible linker, their properties are by no means a simple sum of characteristics displayed by monomeric compounds. For instance, acridine monomers have been shown to directly interact with DNA, which was confirmed by up to a 20-degree increase in the delta Tm of the double-stranded DNA helix [13,14]. However, our preliminary studies revealed that UAs increased this parameter by merely 5–7 degrees, which suggests non-specific interactions between UAs and the DNA helix [1]. Instead, it has been hypothesized that the unique dimer structure of UAs may promote stabilization of DNA quadruplexes [1]. However, due to the large structural variability of G-quadruplexes [15,16,17,18], it is difficult to determine the nature of these interactions.
In this context, it is essential to investigate how pH values relevant to normal and tumor tissue might affect the physicochemical properties of these anticancer compounds. Healthy tissue maintains a pH between 7.2 and 7.4 in the intracellular and extracellular space. However, while the intracellular pH of cancer cells remains approximately 7.4, the pH of their extracellular microenvironment is decidedly more acidic and can vary from 5.5 to 7.0 [19,20,21]. This suggests significance of a lower pH in relation to the physicochemical properties of UAs, and therefore their possible selectivity towards cancer cells.
In our preliminary studies, we identified at least two distinct pKa values for UAs and their three different pH-dependent different forms [8]. The compounds’ protonation state would therefore affect their solubility in water and ultimately the nature of their possible interactions with various cellular components [22,23,24]. In particular, DNA G-quadruplexes may be the main target of anticancer UAs owing to negatively charged groups present in nucleic acids. Considering the above, it is crucial to fully understand variations in charge between pH-dependent structures of UAs, especially considering their potential to interact with DNA [25,26,27].
In the present work, we aimed to establish structures of all possible pH-dependent ionic forms of UAs, particularly at a pH relevant for tumor cells as compared to normal cells. In order to identify possible structural elements of UAs which might be crucial for interactions with DNA, we performed extensive chemometric analysis of UV-Vis spectra of four UAs selected on the basis of their high antitumor activity. The obtained data were supported by nuclear magnetic resonance (NMR) experiments and by the interpretation of TOCSY and ROESY 2D NMR spectra. This allowed us to establish protonation states and assess the aggregation propensity of the studied compounds (C-2028, C-2041, C-2045, and C-2053), the structures of which are shown in Figure 1. The purpose of the study was to gain insight into the relevance of physicochemical properties displayed by UAs in the context of their previously observed significant antitumor effect.

2. Results

2.1. Determination of pKa Values for Studied Compounds

pKa values for the studied unsymmetrical bisacridines were determined by analysis of their UV-Vis spectra under several pH conditions. It was found that the differences between the spectra of individual forms are rather inconsiderable. Therefore, a complex chemometric analysis of the obtained spectra was performed following the procedure described in Materials and Methods. Figure 3 depicts molar fractions of individual spectral forms for the studied UAs. The points represent the coordinates of the observations, which were calculated from experimental data using PCA. Solid curves depict the best fit of the theoretical model to the experimental data according to Equations (1)–(3), provided in Section 4.5.1.
Figure 3 presents the determined molar fractions of five independent spectral forms in relation to pH solution for C-2028, C-2045, and C-2053, and four for C-2041. The initial pKa values were identified as points at the intersection of relevant molar fractions and were then optimized by simplex optimization. The obtained pKa values and their standard deviations (SDs) were evaluated by leave-one-out cross validation and are summarized in Table 1.
A possible model of successive molecule deprotonation events in the studied pH range from 1.5 to 11 was proposed as described below. First, it was crucial to determine the pH-dependent structures of UA monomer analogues under the studied conditions. In this respect, pH-dependent forms of monomers were established using combined UV-Vis and chemometric methodology as well as advanced NMR spectroscopy. The structures of three analysed compounds: 5. 1-alkylamino-8-hydroxyimidazoacridinone (Symadex®, C-1311; Xanthus Pharmaceuticals, Cambridge, MA, United States), 6. 1-alkylaminoimidazoacridinone, and 7. 1-nitro-9-alkylaminoacridine (Ledakrin, C-283), representing monomer units of UAs, are presented in Figure 2. Analysis of compounds (5) and (6) provided information regarding dissociation constants associated with the imidazoacridinone ring and the hydroxyl group substituted at position 8 in this ring, whereas compound (7), with the 1-nitroacridine ring. pKa values determined for the described monomers, is displayed in Table 2.
These results allowed us to propose a sequence of proton dissociation for UAs. Analogous to the pKa of phenol, the pKa3 (9.675) is most likely associated with the hydroxyl group in the imidazoacridinone ring; therefore, it is expected that the pKa5 value for the dimeric C-2045 represents the hydroxyl group. Moreover, the structure of C-2045 and C-2053 is identical apart from the hydroxyl group at position 8 of the imidazoacridinone ring in C-2045. The highest pKa value determined for C-2053 was 8.235, suggesting that pKa5 is related to the hydroxyl group present in C-2045.
In the present work, we aimed to establish the structures of all possible pH-dependent ionic forms of UAs, particularly at a pH relevant for tumor cells as compared to normal cells. In order to identify possible structural elements of UAs which might be crucial for interactions with DNA, we performed extensive chemometric analysis of UV-Vis spectra of four UAs selected on the basis of their high antitumor activity. The obtained data were supported by nuclear magnetic resonance (NMR) experiments and by the interpretation of TOCSY and ROESY 2D NMR spectra. This allowed us to establish protonation states and assess the aggregation propensity of the studied compounds (C-2028, C-2041, C-2045, and C-2053), the structures of which are shown in Figure 1. The purpose of the study was to gain insight into the relevance of physicochemical properties displayed by UAs in the context of their previously observed significant antitumor effect.

2.2. UA Structure Determination by NMR Spectroscopy

The obtained pKa results, which are summarized in Table 1, revealed the total number of protonation forms possible for each studied UA compound. Although comparative UV-Vis spectra analyses of UAs and their monomeric analogues provided a general idea of pH-dependent proton dissociation sequences, in certain cases additional experiments were required in order to paint the full picture. The missing dots were connected upon NMR studies of UAs at four selected pH values. The lowest and highest pH conditions were chosen to study the fully protonated and fully deprotonated UA species, respectively.
In order to enable spectroscopic observations of labile protons of UA molecules, the relevant NMR experiments, i.e., one-dimensional and 2D TOCSY and ROESY spectra, were recorded at 5 °C. At these conditions, the exchange ratio of amino, imino, and hydroxyl protons was significantly lowered due to the increasing proton-dissociation energetic barrier. Unfortunately, we were unable to record standard NOESY experiments for UAs, as these compounds exhibit molecular masses just below 1 kDa, which corresponds to the phase-shift region of the nuclear Overhauser effect.
Our 1H NMR analysis of UAs (Figure 4) at pH~1 revealed signals (1) in the 11–12 ppm region and (2) in the 9.5–11 ppm region, both of which disappeared at pH = 7. These protons were identified upon unambiguous correlations, observed in ROESY and TOCSY spectra (data not shown), as: (1) an imino proton bound to the nitrogen atom of the 1-nitroacridine ring system and (2) an amino proton of the tertiary nitrogen in the linking chain, respectively. Hence, pKa2 values of 6.139, 5.843, 6.120, and 5.997 for C-2028, C-2041, C-2045, and C-2053, respectively, were assigned to the nitrogen atom of the 1-nitroacridine ring. Meanwhile, the amino group at position 9 of the 1-nitro-9-aminoacridine ring is not protonated and exists mainly as an imine at a pH lower than the pKa value of the aromatic acridine nitrogen. Dissociation of this proton results in a negative charge of the entire aromatic system, which in consequence causes the imine group at position 9 to be more susceptible to hydrolysis (data not shown).
Continuing the above, the next pKa3 values, ranging from 7.158 to 7.522, were related to the proton dissociating from the tertiary nitrogen atom in the linking chain (previously referred to as point 2), whereas pKa4 values equal to 8.235 and 8.546 were assigned to the nitrogen atom in the imidazole part of the imidazoacridinone ring. It should be noted that the proton attached to the imidazole ring was not visible in the NMR spectra. These results demonstrate that the higher pKa4 values for C-2028 and C-2053, rather than being associated with the tertiary amino group of the side chain, are actually related to the nitrogen atom in the imidazole part of the imidazoacridinone ring. Moreover, pKa5 equal to 9.137 should relate to dissociation of the proton of the 8-hydroxyl group of C-2045, which is in line with the standard pKa values of phenol-type molecules. A summary of this scheme is presented in Figure 5, illustrating the sequence of proton dissociation in UAs described above.

2.3. Self-Association of UAs

Several hundred UA spectra collected in selected pH and temperature conditions were used to study the self-association of UAs. Subsequently, chemometric analysis of the obtained spectra was performed according to the procedure described in Materials and Methods. Figure 6 depicts molar fractions of all studied compounds at pH 2.5. Points corresponding to molar fractions of independent spectral forms were determined on the basis of experimental data and were derived from coordinates of the observations obtained via PCA. Solid curves represent the best fit of the theoretical model to the experimental data.
The theoretical model of self-association was determined by numerical fitting. Three different aggregation models were considered: dimerization, unlimited aggregation, and micellar aggregation. The obtained self-association constants were evaluated via leave-one-out cross validation.
Figure 6 shows that at 35 °C and pH 2.5 aggregation proceeds similarly for C-2028, C-2045, and C-2053, while self-association of C-2041 is observed to a lesser extent. In this case, dimerization is the best model for all UAs. The results slightly differ for pH = 3.5 (Figure 7); while self-association remains a one-step process, the determined molar fractions no longer fit a simple aggregation model and unlimited self-association seems to be the most suitable model.
The cross-points of the individual molar fractions shifted towards lower concentration values, suggesting that the data fit a different self-association model, while the self-association constants clearly have lower values.
Moreover, in Figure 8 we demonstrate that pH equal to or greater than 6.5 increases the complexity of UA self-association, as it progresses in at least two stages. Aggregation of C-2041 at pH 6.5 is relatively less complicated than for the remaining compounds, as it occurs in a single step. While self-association analysis was additionally carried out at 4 °C, it failed to provide a consistent aggregation model due to the strong propensity for aggregation based on two forms. This was supported by elevated values of the self-association constant. All determined self-association constants including their standard deviations are listed in Table 3. These findings support the view that self-association of UAs is strongly pH-dependent, as it is directly correlated with increasing pH values.

3. Discussion

With respect to the aim of this study, our findings outlined in Table 1 indicate that the four UA antitumor derivatives (Figure 1) selected for this study exist as five (for C-2028, C-2041, and C-2053) or six (for C-2045) independent forms in the pH range from 1.5 to 11. Notably, at physiological pH (between 6 and 8) there are three to four individual protonation forms, each most likely capable of distinct interactions with DNA and therefore of diversely affecting cellular processes at the molecular level. Moreover, our studies were prompted by the current understanding that the pH value maintained in the extracellular environment of cancer cells is lower (5.5 to 7.0) than that in healthy tissue (7.2 to 7.4). This suggests that certain protonation states of the studied compounds might be selective towards cancer cells.
It was challenging to assign a particular pKa value acquired by UV-Vis analysis, to the respective (de)protonating atom using of NMR. We were able to obtain preliminary assignments through an approach based on UA structural analogues. NMR analysis of the monomers Symadex® (C-1311), Ledakrin (C-283), and 1-alkylaminoimidazoacridinone (Figure 3) provided evidence that pKa values of approximately 2 and 9 were associated with the imidazoacridione ring, whereas 6 was related to the 1-nitroacridine ring.
NMR analysis confirmed the validity of the preliminary pKa assignments obtained through analysis of UA structural analogues. Although detection of all exchangeable protons by NMR was not possible due to different exchange ratios with solvent molecules, this technique allowed us to recognize exchangeable protons related to the alkyl part of UAs. Moreover, NMR enabled us to determine the sequence of proton detachment. Our approach, incorporating the several strategies described above, allowed us to clearly define the exact structures of the individual UA protonation states, which are presented in Figure 5.
Furthermore, we were able to demonstrate that the investigated unsymmetrical bisacridines are prone to self-association. UV-Vis spectra extensively analysed by chemometric techniques revealed that all studied UAs underwent aggregation, which was strongly pH-dependent. Principal component analysis demonstrated the presence of at least two sources of variability in the spectra of various UA concentrations collected at several pH values. Based on the obtained results, we propose dimerization as the best fitted model of self-association in acidic media. Determining the nature of the aggregation processes at physiological pH proved to be more challenging. Presumably, it involves a dimer as the first order aggregate followed by a second order aggregate of greater complexity. The complicated nature of the aggregation process may result from distinct structural forms related to the described protonation states.
In conclusion, the results presented here can provide an insightful background for further studies including metabolic transformations of UAs and their possible interactions with nucleic acids and/or proteins considering pH differences between normal and tumor tissues. Proton exchange significantly rearranges the distribution of electron density in a compound, which in turn strongly influences its ability to self-associate, potentially impacting the anticancer activity of these compounds.

4. Materials and Methods

4.1. Chemicals

Unsymmetrical bisacridines: 9′-{N-[(imidazo [4,5,1-de]-acridin-6-on-5-yl)aminopropyl]-N-methylaminopropylamino}-10-nitroacridine (C-2028, 1); 1-[3-(imidazo [4,5,1-de]-acridin-6-on-5-yl)aminopropyl]-4-[3′-(10-nitroacridin-1-yl)-aminopropyl]piperazine (C-2041, 2); 9′-{N-[(8-hydroxyimidazo [4,5,1-de]-acridin-6-on-5-yl)aminopropyl]-N-methylaminopropylamino}-40-methyl-10-nitroacridine (C-2045, 3); and 9′-{N-[(imidazo [4,5,1-de]-acridin-6-on-5-yl)aminopropyl]-N-methylaminopropylamino}-40-methyl-10-nitroacridine (C-2053, 4) were synthesized and purified in our laboratory according to the method described in the relevant patents. HEPES and PIPES buffers, potassium cacodylate, sodium acetate, sodium hydroxide, and sodium tetraborate were purchased from Merck (Merck Millipore, Darmstadt, Germany). Deionized water (conductivity 0.056 μS·cm−1, Milli-Q water purification system; Merck Millipore) was used in all experiments.

4.2. General Procedures

All solvents and chemicals used required no further purification. Stock UA solutions (10−2 M) were prepared in deionized water and then diluted to a 10−3 M or 10−4 M concentration if necessary.

4.3. pKa Determination

Hydrochloric acid or buffer (2 × 10−1 M) solutions of UAs (2 × 10−5 M) were prepared with the following pH values: 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6.1, 6.3, 6.5, 6.7, 6.9, 7.1, 7.3, 7.5, 7.7, 7.9, 8.1, 9, 10 and 11. Based on their buffering capacity, we used acetate, PIPES, HEPES, and borate buffer in the pH range from 4 to 11 and hydrochloric acid in the pH range from 1.5 to 3.5. UV-Vis spectra of the samples were recorded in a quartz cuvette with an optical path length of 1 cm.

4.4. Aggregation Study

UAs in the following concentrations: 10−6, 2 × 10−6, 5 × 10−6, 10−5, 2 × 10−5, 5 × 10−5, 10−4, 2 × 10−4, 5 × 10−4, and 10−3 M were prepared as HEPES buffer solutions (10−1 M) with pH values of 2.5, 6.5, 7.4, and 8.5. We investigated the propensity of UAs to undergo aggregation under the conditions described above at a temperature of 37 °C, as required for enzymatic transformations of UAs. UV-Vis spectra were collected at various temperatures (4 °C, 22 °C, and 37 °C) with the optical path length of a quartz cuvette varying from 0.1 cm to 1 cm.

4.5. Chemometric Analysis: General Procedures

All absorbance values were expressed as molar extinction coefficients. Spectral data were subjected to appropriate pre-processing involving correction of signal distortion and elimination of outliers. Subsequently, data were organized into a matrix and subjected to numerical decomposition into eigenvectors (Principal Component Analysis, PCA). Significant principal components were determined through a linear map and residual spectral analysis. Based on the retained number of significant principal components (n), PCA was repeated to extract the (n − 1) dimensional hyperplane from the n-dimensional space. For hyperplanes with dimensions greater than or equal to three, the coordinates of the observations (V) were set in an n-dimensional simplex. The corners of these entities were manually determined and then optimised using the Nelder–Mead method. The contributions of the principal components were assessed using the V and then optimized using the penalty function. Molar fractions were obtained based on the contributions of individual principal components and calculated using the coefficients of the plane equations. A relevant plane was identified using the linear regression method. For two-dimensional space, the initial and final points along the principal component were identified using the Nelder–Mead method. All steps of the performed chemometric analysis are described in detail in the Supplementary Materials.

4.5.1. Chemometric Analysis for pKa Determination

Dissociation constants (pKa) were determined based on molar fractions using an extended version of the Henderson–Hasselbalch equation (Equations (1)–(3)):
γ = i = 1 4 K i i = 1 4 H i
u 1 = 1 1 + γ  
u n = u 1 · γ · n = { 1 , 4 }  
K—dissociation constant,
un—molar fractions.
The initial pKa values were identified as intersection points of the relevant molar fractions and then optimised using the Nelder–Mead method [28] (consult Supplementary Materials for more details). The accuracy of the calculated molar fractions was estimated using leave-one-out cross-validation (LOOCV). This sophisticated methodology was employed to determine the pKa values of multiprotic compounds in the 2–11 pH range.

4.5.2. Chemometric Analysis for Self-Association Study

Self-association constants (K), for two-dimensional space, were identified by fitting the determined molar fractions into an appropriate theoretical model (dimerization), then adjusted using simplex optimisation. For three-dimensional space, K values were determined as reciprocals of the logarithm of concentration corresponding to the intersection of the molar fractions. All calculations were performed in the R environment (R Foundation for Statistical Computing, version 4.0.2, Vienna, Austria).

4.6. NMR Studies of UA Protonation Forms

In order to track the resonances of the labile protons of the studied UAs, relevant solutions (10−3 M) were prepared in the following buffers or acids (10−2 M): hydrochloric or sulphuric acid (pH = 1), cacodylate buffer (pH values of 5, 7 or 8), and borate buffer (pH = 11). Solutions contained 10% D2O, while those analysed by TOCSY and ROESY were prepared in pure D2O. Subsequently, sets of 1D and 2D NMR (TOCSY and ROESY) were recorded by a 700 MHz Avance III spectrometer (Bruker) equipped with a QCl CryoProbe (in collaboration with the Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznań, Poland). All spectra were processed using TopSpin software (Bruker GmbH, version 4.1.3, Karlsruhe, Germany). The TOCSY mix time was set to 80 ms, whereas the ROESY mix time was set to 300 ms.
Furthermore, a complete assignment of the non-exchangeable protons of the UAs was performed on the basis of an additional set of DQF-COSY, TOCSY, and ROESY spectra. These experiments were recorded at pH = 4.5 (acetate buffer) and a temperature of 25 °C using samples of compounds 14 at concentrations of 10−3 M dissolved in pure D2O. The TOCSY mix time was set to 80 ms, whereas the ROESY mix time was set to 300 ms. The resulting assignments as well as fragments of selected 2D NMR spectra are included in the Supplementary Materials (Table S1, Figures S6–S10).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27133995/s1, Figure S1: Spectra for analysis of the C-2045 protonation forms; Figure S2: Initial spectra for aggregation study of C-2045 (left), centered spectra (right); Figure S3: Spectra of the C-2045 components (pKa determination); Figure S4: Spectra of the C-2045 components (aggregation study); Figure S5: Quaternary plot for C-2041 in orthogonal projection (left) and perspective projection (right); Figure S6: The general structure of UA’s with proton codenames, used in NMR assignments; Figure S7: Aromatic region of the DQF-COSY spectrum of C-2041; Figure S8: Fragment of the ROESY spectrum of C-2041, depicting dipolar couplings between the aromatic protons and the linker sidechain; Figure S9: Aromatic region of the TOCSY spectrum of C-2045; Figure S10: Fragment of the ROESY spectrum of C-2045, depicting dipolar couplings between the aromatic protons and the aminoalkyl linker. Also, diagnostic Me2′-H3′ ROE could be observed; Table S1: Chemical shifts (δ) of resonance signals of non-exchangeable protons of UA’s, recorded at pH = 4.5 in D2O. Multiplet codenames: s—singlet, d—doublet, t—triplet, m—unresolved multiplet. All chemical shifts were referenced to residual water signal at δ = 4.780 ppm. Protons were referenced as shown at Figure S6.

Author Contributions

Conceptualization, Z.M.; Data curation, M.K., T.L. and J.E.F.; Formal analysis, M.K. and J.B.-B.; Funding acquisition, Z.M.; Investigation, M.K., T.L., J.E.F., A.K., W.A., J.B.-B. and K.S.-K.; Project administration, Z.M.; Software, M.K. and J.B.-B.; Supervision, Z.M.; Writing—original draft, M.K., T.L., J.E.F. and A.P.; Writing—review and editing, Z.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by a grant from the Polish National Science Centre, no. 2019/33/B/NZ7/02534.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analysed during this study are included in this published article and its supplementary material files.

Acknowledgments

The authors wish to thank Jan Mazerski (Department of Pharmaceutical Technology and Biochemistry, Faculty of Chemistry, Gdańsk University of Technology, Gdańsk, Poland) for his support and substantial help in chemometric analyses. We thank Ewa Paluszkiewicz (Department of Pharmaceutical Technology and Biochemistry, Faculty of Chemistry, Gdańsk University of Technology, Gdańsk, Poland) for synthesis of UAs as described in relevant patents.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Konopa, J.; Horowska, B.; Paluszkiewicz, E.; Borowa-Mazgaj, B.; Augustin, E.; Skwarska, A.; Mazerska, Z. Unsymmetrical Bisacridines with Antitumor Activity and Use Thereof. European Patent EP 3070078 A1, 10 October 2017. [Google Scholar]
  2. Konopa, J.; Horowska, B.; Paluszkiewicz, E.; Borowa-Mazgaj, B.; Augustin, E.; Skwarska, A.; Mazerska, Z. Asymmetric Bis-Acridines with Antitumour Activity and Their Uses. United States Patent US10202349B2, 12 February 2019. [Google Scholar]
  3. Paluszkiewicz, E.; Horowska, B.; Borowa-Mazgaj, B.; Peszyńska-Sularz, G.; Paradziej-Łukowicz, J.; Augustin, E.; Konopa, J.; Mazerska, Z. Design, synthesis and high antitumor potential of new unsymmetrical bisacridine derivatives towards human solid tumors, specifically pancreatic cancers and their unique ability to stabilize DNA G-quadruplexes. Eur. J. Med. Chem. 2020, 204, 112599. [Google Scholar] [CrossRef] [PubMed]
  4. Kulesza, J.; Pawłowska, M.; Augustin, E. The Influence of Antitumor Unsymmetrical Bisacridines on 3D Cancer Spheroids Growth and Viability. Molecules 2021, 26, 6262. [Google Scholar] [CrossRef] [PubMed]
  5. Pilch, J.; Matysiak-Brynda, E.; Kowalczyk, A.; Bujak, P.; Mazerska, Z.; Nowicka, A.M.; Augustin, E. New Unsymmetrical Bisacridine Derivatives Noncovalently Attached to Quaternary Quantum Dots Improve Cancer Therapy by Enhancing Cytotoxicity toward Cancer Cells and Protecting Normal Cells. ACS Appl. Mater. Interfaces 2020, 12, 17276–17289. [Google Scholar] [CrossRef] [PubMed]
  6. Pilch, J.; Kowalik, P.; Bujak, P.; Nowicka, A.M.; Augustin, E. Quantum Dots as a Good Carriers of Unsymmetrical Bisacridines for Modulating Cellular Uptake and the Biological Response in Lung and Colon Cancer Cells. Nanomaterials 2021, 11, 462. [Google Scholar] [CrossRef] [PubMed]
  7. Pilch, J.; Kowalik, P.; Kowalczyk, A.; Bujak, P.; Kasprzak, A.; Paluszkiewicz, E.; Augustin, E.; Nowicka, A.M. Foliate-Targeting Quantum Dots-β-Cyclodextrin Nanocarrier for Efficient Delivery of Unsymmetrical Bisacridines to Lung and Prostate Cancer Cells. Int. J. Mol. Sci. 2021, 23, 1261. [Google Scholar] [CrossRef] [PubMed]
  8. Mieszkowska, A.; Nowicka, A.; Kowalczyk, A.; Potęga, A.; Pawłowska, M.; Kosno, M.; Augustin, E.; Mazerska, Z. Metabolic Profiles of New Unsymmetrical Bisacridine Antitumor Agents in Electrochemical and Enzymatic Noncellular Systems and in Tumor Cells. Pharmaceuticals 2021, 14, 317. [Google Scholar] [CrossRef]
  9. Potęga, A.; Kosno, M.; Mazerska, Z. Novel insights into conjugation of antitumor-active unsymmetrical bisacridine C-2028 with glutathione: Characteristics of non-enzymatic and glutathione S-transferase-mediated reactions. J. Pharm. Anal. 2021, 11, 791–798. [Google Scholar] [CrossRef]
  10. WHO Chronicle. International Non-proprietary Names for Pharmaceutical Substances. Suppl. WHO Chron. 1976, 30, 1–18. [Google Scholar]
  11. Wiśniewska, A.; Chrapkowska, A.; Kot-Wasik, A.; Konopa, J.; Mazerska, Z. Metabolic transformations of antitumor imidazoacridinone, C-1311, with microsomal fractions of rat and human liver. Acta Biochim. Pol. 2007, 54, 831–838. [Google Scholar] [CrossRef] [Green Version]
  12. Gniazdowski, M.; Szmigiero, L. Nitracrine and its congeners-An overview. Gen. Pharmacol. 1995, 26, 473–481. [Google Scholar] [CrossRef]
  13. Laskowski, T.; Borzyszkowska, J.; Grynda, J.; Mazerski, J. C-1311 (Symadex), a potential anti-cancer drug, intercalates into DNA between A and G moieties. NMR-derived and MD-refined stereostructure of the d(GAGGCCTC) 2:C-1311 complex. J. Mol. Struct. 2017, 1141, 357–367. [Google Scholar] [CrossRef]
  14. Laskowski, T.; Czub, J.; Sowiński, P.; Mazerski, J. Intercalation complex of imidazoacridinone C-1311, a potential anticancer drug, with DNA helix d(CGATCG)2: Stereostructural studies by 2D NMR spectroscopy. J. Biomol. Struct. Dyn. 2016, 34, 653–663. [Google Scholar] [CrossRef] [PubMed]
  15. Bochman, M.L.; Paeschke, K.; Zakian, V.A. DNA secondary structures: Stability and function of G-quadruplex structures. Nat. Rev. Genet. 2012, 13, 770–780. [Google Scholar] [CrossRef] [Green Version]
  16. Kolesnikova, S.; Curtis, E.A. Structure and Function of Multimeric G-Quadruplexes. Molecules 2019, 24, 3074. [Google Scholar] [CrossRef] [Green Version]
  17. Chen, Y.; Yang, D. Sequence, Stability, Structure of G-Quadruplexes and Their Drug Interactions. Curr. Protoc. Nucleic Acid Chem. 2012, 50, 17.5.1–17.5.17. [Google Scholar] [CrossRef]
  18. Phan, A.T.; Luu, K.N.; Patel, D.J. Different loop arrangements of intramolecular human telomeric (3+1) G-quadruplexes in K+ solution. Nucleic Acids Res. 2006, 34, 5715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. White, K.A.; Grillo-Hill, B.K.; Barber, D.L. Cancer cell behaviors mediated by dysregulated pH dynamics at a glance. J. Cell Sci. 2017, 130, 663–669. [Google Scholar] [CrossRef] [Green Version]
  20. Galenkamp, K.M.O.; Sosicka, P.; Jung, M.; Recouvreux, M.V.; Zhang, Y.; Moldenhauer, M.R.; Brandi, G.; Freeze, H.H.; Commisso, C. Golgi acidification by NHE7 regulates cytosolic pH homeostasis in pancreatic cancer cells. Cancer Discov. 2020, 10, 823–835. [Google Scholar] [CrossRef] [Green Version]
  21. Xu, J.; Li, Y.; Li, Z.; Shao, W.; Song, J.; Wei, J. Acidic tumor microenvironment promotes pancreatic cancer through miR-451a/MEF2D axis. J. Oncol. 2022, 12, ID3966386. [Google Scholar] [CrossRef]
  22. Petukh, M.; Stefl, S.; Alexov, E. The Role of Protonation States in Ligand-Receptor Recognition and Binding. Curr. Pharm. 2013, 19, 4182. [Google Scholar] [CrossRef] [Green Version]
  23. Onufriev, A.V.; Alexov, E. Protonation and pK changes in protein-ligand binding. Q. Rev. Biophys. 2013, 46, 181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Longo, G.S.; Pérez-Chávez, N.A.; Szleifer, I. How protonation modulates the interaction between proteins and pH-responsive hydrogel films. Curr. Opin. Colloid Interface Sci. 2019, 41, 27–39. [Google Scholar] [CrossRef]
  25. Nagesh, N.; Krishnaiah, A.; Dhople, V.M.; Sundaram, C.S.; Jagannadham, M.V. Noncovalent interaction of G-quadruplex DNA with acridine at low concentration monitored by MALDI-TOF mass spectrometry. Nucleosides Nucleotides Nucleic Acids 2007, 26, 303–315. [Google Scholar] [CrossRef] [PubMed]
  26. Ferreira, R.; Artali, R.; Benoit, A.; Gargallo, R.; Eritja, R.; Ferguson, D.M.; Sham, Y.Y.; Mazzini, S. Structure and Stability of Human Telomeric G-Quadruplex with Preclinical 9-Amino Acridines. PLoS ONE 2013, 8, e57701. [Google Scholar] [CrossRef]
  27. Gowan, S.M.; Heald, R.; Stevens, M.F.; Kelland, L.R. Potent inhibition of telomerase by small-molecule pentacyclic acridines capable of interacting with G-quadruplexes. Mol. Pharmacol. 2001, 60, 981–988. [Google Scholar] [CrossRef]
  28. Nelder, J.A.; Mead, R. A Simplex Method for Function Minimization. Comput. J. 1965, 7, 308–313. [Google Scholar] [CrossRef]
Figure 1. Chemical structures of the studied unsymmetrical bisacridines (UAs): C-2028 (1), C-2041 (2), C-2045 (3), and C-2053 (4).
Figure 1. Chemical structures of the studied unsymmetrical bisacridines (UAs): C-2028 (1), C-2041 (2), C-2045 (3), and C-2053 (4).
Molecules 27 03995 g001
Figure 2. Monomers of the studied unsymmetrical bisacridine structures: 1-alkylamino-8-hydroxyimidazoacridinone (Symadex®, C-1311, (5)), 1-alkylaminoimidazoacridinone (6), and 1-nitro-9-alkylaminoacridine (Ledakrin, C-283, (7)).
Figure 2. Monomers of the studied unsymmetrical bisacridine structures: 1-alkylamino-8-hydroxyimidazoacridinone (Symadex®, C-1311, (5)), 1-alkylaminoimidazoacridinone (6), and 1-nitro-9-alkylaminoacridine (Ledakrin, C-283, (7)).
Molecules 27 03995 g002
Figure 3. Dependence of UA molar fractions on pH values for each protonated form (points—optimal molar fractions; curves—the best fit of the theoretical model to the data). Each color corresponds to a unique spectral (i.e., protonation) form of a given molecule. The spectral forms possible for studied compounds are presented in Figure 5.
Figure 3. Dependence of UA molar fractions on pH values for each protonated form (points—optimal molar fractions; curves—the best fit of the theoretical model to the data). Each color corresponds to a unique spectral (i.e., protonation) form of a given molecule. The spectral forms possible for studied compounds are presented in Figure 5.
Molecules 27 03995 g003
Figure 4. 1H NMR spectra of UAs: C-2028 at pH = 1 (A) and pH = 7 (B), C-2041 at pH = 1 (C) and pH = 7 (D), C-2045 at pH = 1 (E) and pH = 7 (F), C-2053 at pH = 1 (G) and pH = 7 (H). Temperature = 5 °C; H2O/D2O collected in 10−2 M cacodylate buffer.
Figure 4. 1H NMR spectra of UAs: C-2028 at pH = 1 (A) and pH = 7 (B), C-2041 at pH = 1 (C) and pH = 7 (D), C-2045 at pH = 1 (E) and pH = 7 (F), C-2053 at pH = 1 (G) and pH = 7 (H). Temperature = 5 °C; H2O/D2O collected in 10−2 M cacodylate buffer.
Molecules 27 03995 g004
Figure 5. Proposed structures of individual pH-dependent forms of studied UAs (14) and the corresponding approximated pKa values. C–2045 (3, black structures) was used as an example for the last pKa = 9.137.
Figure 5. Proposed structures of individual pH-dependent forms of studied UAs (14) and the corresponding approximated pKa values. C–2045 (3, black structures) was used as an example for the last pKa = 9.137.
Molecules 27 03995 g005
Figure 6. UA molar fractions for UAs studied at pH = 2.5 and 35 °C. Points—optimal molar fractions; curves—the best fit of the theoretical model of dimerization to the data. Red color—monomer, blue color—dimer.
Figure 6. UA molar fractions for UAs studied at pH = 2.5 and 35 °C. Points—optimal molar fractions; curves—the best fit of the theoretical model of dimerization to the data. Red color—monomer, blue color—dimer.
Molecules 27 03995 g006
Figure 7. Molar fractions for UAs studied at pH = 3.5 and 35 °C. Points—optimal molar fractions; curves—the best fit of the theoretical model of dimerization to the data. Red color—monomer, blue color—dimer.
Figure 7. Molar fractions for UAs studied at pH = 3.5 and 35 °C. Points—optimal molar fractions; curves—the best fit of the theoretical model of dimerization to the data. Red color—monomer, blue color—dimer.
Molecules 27 03995 g007
Figure 8. Molar fractions of UAs at pH = 6.5, 7.4, 8.5 at 35 °C. Dotted curves—molar fractions of individual self-associates, solid curves—the best fit of the theoretical model. Red color—monomer, green color—first order aggregate (presumably a dimer), blue color—higher order aggregate.
Figure 8. Molar fractions of UAs at pH = 6.5, 7.4, 8.5 at 35 °C. Dotted curves—molar fractions of individual self-associates, solid curves—the best fit of the theoretical model. Red color—monomer, green color—first order aggregate (presumably a dimer), blue color—higher order aggregate.
Molecules 27 03995 g008
Table 1. pKa and standard deviation values, calculated using leave-one-out cross validation.
Table 1. pKa and standard deviation values, calculated using leave-one-out cross validation.
No.UAspKa1SDpKa2SDpKa3SDpKa4SDpKa5SD
1C-20282.1940.0056.1390.0017.3970.0058.5460.02--
2C-20412.4350.0045.8430.0057.1580.004----
3C-20452.3010.0036.1200.0027.3270.006--9.1370.04
4C-20532.2850.0055.9970.0027.5220.0038.2350.004--
Table 2. pKa and standard deviation values determined for the selected monomer structures presented in Figure 2.
Table 2. pKa and standard deviation values determined for the selected monomer structures presented in Figure 2.
No.pKa1SDpKa2SDpKa3SD
52.7680.0047.5970.0079.6750.02
62.6620.0077.6860.009--
7--6.2270.0007--
Table 3. Self-association constants (K) with SD at various pH values at 35 °C. Constants marked by ‘*’ were determined via leave-one-out cross validation.
Table 3. Self-association constants (K) with SD at various pH values at 35 °C. Constants marked by ‘*’ were determined via leave-one-out cross validation.
No.pH2.53.56.57.48.5
UAsKSDKSDKSDKSDKSD
1C-202810.89 *0.0226.94 *0.033~20 ~40, ~10 ~220,
~50
2C-20411.330.0450.77 *0.0233.25 *0.037~200, ~13 ~250,
~143
3C-204515.72 *0.026~63 ~125, ~40 ~200, ~40 ~125,
~10
4C-20539.14 *0.026~50 ~60, ~5 ~250, ~20 ~400,
~60
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kosno, M.; Laskowski, T.; Frackowiak, J.E.; Potęga, A.; Kurdyn, A.; Andrałojć, W.; Borzyszkowska-Bukowska, J.; Szwarc-Karabyka, K.; Mazerska, Z. Acid–Base Equilibrium and Self-Association in Relation to High Antitumor Activity of Selected Unsymmetrical Bisacridines Established by Extensive Chemometric Analysis. Molecules 2022, 27, 3995. https://doi.org/10.3390/molecules27133995

AMA Style

Kosno M, Laskowski T, Frackowiak JE, Potęga A, Kurdyn A, Andrałojć W, Borzyszkowska-Bukowska J, Szwarc-Karabyka K, Mazerska Z. Acid–Base Equilibrium and Self-Association in Relation to High Antitumor Activity of Selected Unsymmetrical Bisacridines Established by Extensive Chemometric Analysis. Molecules. 2022; 27(13):3995. https://doi.org/10.3390/molecules27133995

Chicago/Turabian Style

Kosno, Michał, Tomasz Laskowski, Joanna E. Frackowiak, Agnieszka Potęga, Agnieszka Kurdyn, Witold Andrałojć, Julia Borzyszkowska-Bukowska, Katarzyna Szwarc-Karabyka, and Zofia Mazerska. 2022. "Acid–Base Equilibrium and Self-Association in Relation to High Antitumor Activity of Selected Unsymmetrical Bisacridines Established by Extensive Chemometric Analysis" Molecules 27, no. 13: 3995. https://doi.org/10.3390/molecules27133995

APA Style

Kosno, M., Laskowski, T., Frackowiak, J. E., Potęga, A., Kurdyn, A., Andrałojć, W., Borzyszkowska-Bukowska, J., Szwarc-Karabyka, K., & Mazerska, Z. (2022). Acid–Base Equilibrium and Self-Association in Relation to High Antitumor Activity of Selected Unsymmetrical Bisacridines Established by Extensive Chemometric Analysis. Molecules, 27(13), 3995. https://doi.org/10.3390/molecules27133995

Article Metrics

Back to TopTop