Abstract

Breast cancer is one of the major types of cancer around the world, and early diagnosis is essential for successful treatment. New contrast agents (CAs), with reduced toxicology, are needed to improve diagnosis. One of the most promising Magnetic Resonance Imaging (MRI) CA is based on rhenium conjugated with a benzothiazole derivate (ReABT). In this sense, DFT has been used to evaluate the best methodology for calculating the hyperfine coupling constant (Aiso) of ReABT. Then, a thermodynamic analysis was performed to confirm the stability of the complex. Furthermore, a docking study of ReABT at the enzyme PI3K active site and Aiso calculations of ReABT in the enzyme environment were carried out. The best methodology for the Aiso calculation of ReABT was using the M06L functional, SARC-ZORA-TZVP (for Re) and TZVP (for all other atoms) basis set, relativistic Hamiltonian, and the CPCM solvation model with water as the solvent which confirm that the relativistic effects are important for calculating the Aiso values. In addition, thermodynamic analysis indicates that ReABT presents a higher stability and a lower toxicity than Gd-based CAs. The docking studies point out that ReABT interacts with amino acids residues of alanine, aspartate, and lysine from the PI3K active site. Considering the enzyme environment, Aiso values decrease significantly. These findings indicate that the CA candidate ReABT could be a good candidate for a new contrast agent.

1. Introduction

Female breast cancer is expected to be the major type of cancer, besides non-melanoma skin tumors, in Brazil between 2020 and 2022, affecting over 66 thousand cases per year according to the National Institute of Cancer José Alencar (INCA) [1]. Across the world, data from the American Cancer Society estimate that in 2020 more than 275,000 cases should be diagnosed in the United States and, according to the European Cancer Information System, ECIS, in the European Union more than 355,000 diagnoses are expected by 2020 [2, 3].

Early diagnosis is fundamental to the successful treatment of breast cancer. In this context, numerous techniques can help with diagnoses, such as mammography, PET, CT, SPECT, and MRI [4]. The last one has aroused great interest because of its high sensitivity for the identification and characterization of breast cancer and its success in early diagnosis [57].

This technique is based on the principles of NMR spectroscopy, where the signals come from the relaxation time of the hydrogen and oxygen nuclei after the application of a radio frequency [8]. Using more efficient contrast agents (CAs) leads to a better resolution of the images and can cause a change in the longitudinal and transversal relaxation times, T1 and T2, respectively, [5, 8, 9]. Hyperfine interactions, which are related to interactions between nuclear and electron spins, of paramagnetic species change relaxation times and can be obtained by EPR spectroscopy [9, 10].

The hyperfine coupling (A) depends on the contributions of isotropic (Aiso) and dipolar (Adip) hyperfine coupling constants (equation (1)) [9, 10]. The first one comes from an unpaired spin in an s orbital, and the second one results from a nonlocal and local dipolar constant. Nevertheless, the dipolar hyperfine coupling constant is close to zero when fluid solutions are considered [9, 10]. These parameters, since they affect relaxation times and therefore present activity in the EPR, are crucial for studying and designing new contrast agents for MRI [9, 10]:

Traditionally, CAs are metal-based structures like Gd, Pt, and Tc, for example. The study of new CA is important due to the need for increasingly efficient agents, providing better images from the relaxation time of water molecules [11, 12]. Among these metals, we can highlight Tc, highly applied and studied as a contrast agent, and Re, a metal that has great potential as CA because of its low toxicity [5, 1317].

Re compounds, especially those that have the fac-[Re (CO)3]+ fragment, have been extensively studied as a diagnostic agent for breast cancer. Tzanopoulou et al. presented this fragment conjugated to the phenylbenzothiazolic derivative, 2-(4-Aminophenyl) benzothiazole (ABT) (Figure 1), as an agent with potential application for MRI [15]. The ABT compound shows great biological activity as antibacterial, antifungal, antiallergic, anti-inflammatory, and, in our interest, anticancer and inhibiting some enzymes, such as the enzyme kinase [14, 18], especially the phosphoinositide 3-kinase (PI3K) enzyme, which is related with the growth of tumors, like breast cancer [19].

Rhenium compounds, especially those with 186Re and 188Re, present properties, such as β emission and half-life time equal to 3.8 and 0.7 d, respectively, that give them a potential for medicinal applications [11]. The fac-[Re(CO)3] core presents a low spin. This characteristic, allied with the characteristics previously presented, provides high interest in this family of compounds for use as radiopharmaceuticals [14]. This leads to a broad study of the application of these compounds in EPR spectroscopy and MRI [14] and the need to evaluate the solvent and environment effects on the EPR parameters, which module the action of this compound as a potential MRI contrast agent. Therefore, the aim of this work was to study solvent and relativistic effects on spectroscopic properties (Aiso values for 1H and 17O) of the complex Re(CO)3(NNO) conjugated with 2-(4′-aminophenyl) benzothiazole.

2. Materials and Methods

2.1. Structural Evaluation

The first step was to perform a geometry optimization calculation for the Re(CO)3(NNO) conjugated to 2-(4′-aminophenyl) benzothiazole, ReABT. This calculation was performed using the DFT method, using the ωB97X-D3 hybrid-GGA functional [20], using TZVP basis set for all atoms except Re and SARC-ZORA-TZVP [21] for rhenium, and considering the zero order regular approximation (ZORA) [22] on software ORCA 4.2.1 [23].

2.2. Spectroscopic Calculations

For the hyperfine coupling constant calculations (Aiso), the density functional theory was also used, taking the geometry of minimum energy got and applying different approaches to define which methodology works better in our compound. To do that, three different functionals were used, the GGA functional of Perdew-Burke-Ernzerhof, PBE [24], the hybrid-GGA functional PBE0, which mixes the PBE exchange energy with the Hartree-Fock exchange energy [25], and the meta-hybrid GGA functional Minessota 06-L, M06-L [26].

The effect of the solvent was also analyzed by theoretical calculations carried out under vacuum and in the presence of water as a solvent [27]. The solvation method used was the Conductor-like Polarizable Continuum Model, CPCM, since it considers the solvent as a polarizable conductor, so it can describe molecules that interact with polar solvents, like water [28].

So, in the first step, the solvent effect was analyzed. In the second step, relativistic effects were evaluated; to this end, calculations considering these effects in the Hamiltonian with ZORA, using the SARC-ZORA-TZVP basis set for rhenium atom and TZVP for all other atoms, and nonrelativistic calculations were performed with def2-TZVP and TZVP basis sets for rhenium and other atoms, respectively. These three steps will provide us with a methodology that fits better with the compound under study.

2.3. Thermodynamics Investigation and Molecular Docking Calculation

To understand whether the complex is thermodynamically stable, a frequency calculation was performed, using the M06L functional, SARC-ZORA-TZVP basis set for rhenium and TZVP for all other atoms and relativistic effects in the Hamiltonian with ZORA, in the ReABT complex, in the Re(CO)3+ core (Figure 2(c)), and in the tridentate ligand, (NNO) ABT (Figure 2(b)), both in the gas phase, to obtain the Δ, and in solution, to obtain ΔGsolv (ReABT), ΔGsolv (Re(CO)3), and ΔGsolv ((NNO) ABT) values. This analysis will provide the free energy of complexation of the complex (Figure 2(a) and equation (2)), allowing us to understand if the rhenium core is strongly bonded in the complex [29]:

The molecular docking of the complex was performed at the active site of the PI3K enzyme to obtain Aiso values for the compound [7, 17]. For this, the complex was docked inside the PI3K enzyme (Protein Data Bank (PDB) code 3QJZ [30]) using the Molegro Virtual Docking (MVD) software [31]. The conformation that most closely matches the active ligand of the PI3K enzyme, with the amino acid residues that perform hydrogen bonds (H-Bonds) in the PI3K active site, was selected to calculate the hyperfine coupling constant values using the methodology defined previously.

3. Results and Discussion

The results provided in this paper were divided into three major parts. The first one is the validation of the methodology, where the structural results were analyzed and the best theoretical methodology to calculate Aiso was defined. The second part will discuss the thermodynamics of the complex to evaluate the free energy of complexation. In the last part of the results, the methodology defined previously was applied to analyze how the Aiso values of the ReABT complex changed when the complex was docked at the active site of PI3K.

3.1. Validation of the Methodology

The geometry optimization of complex ReABT was performed using DFT calculation with the functional ωB97X-D3, considering relativistic effects by using ZORA, with TZVP basis sets for H, C, O, N, and S atoms and SARC-ZORA-TZVP basis sets for Re. The optimized geometry is presented in Figure 3, with the main bond lengths highlighted. These bond lengths are presented in Table 1. As, to our knowledge, there are no crystallographic data for ReABT, we use the experimental results obtained by Machura et al. for the complex Re(CO)3Cl conjugated to 2-(2-aminophenyl) benzothiazole (ReAPBT) [13]. This complex was chosen since both compounds present the fac-[Re(CO)3]+ core and coordinated in the rhenium, with one equatorial and one axial Re-N bond.

In order to evaluate if whether the theoretical values for the bonds of our complex (Table 1) are in agreement with experimental data and structural parameters from our theoretical findings were compared to Machura’s work [13]. It should be kept in mind, however, that the equatorial bond between Re-O in our complex was replaced by a Re-Cl bond in the compound studied by Machura and co-workers [13], so this change can affect the other bonds. The main difference is in the Re-Nax bond, with 0.114 Å longer in the theoretical result. The C-O bonds are shorter, ranging from 0.005 to 0.009 Å. Differences in bond lengths have a medium value of 0.026 Å.

To improve structural analysis, we compared the bond angles among theoretical and experimental values obtained from the work of Machura et al. (Table 2). The differences in bond angles, presented in Table 2, range from 0.25°, for the Ceq-Re-Ceq bond and 5.07° for the Ceq1-Re-Cax bond when compared with the experimental data of Machura and co-workers’ paper [13]. These results of bond lengths and bond angles show that our result of geometry optimization is in agreement with the experimental results.

Since the geometry of the complex was optimized, hyperfine coupling constant calculations, Aiso, were carried out. Equations (3) and (4) define the relaxation times (T1 and T2), i.e., the parameter that came from the interaction of the paramagnetic substances with water. This interaction causes a large drop in the relaxation times. This relationship was one of the most important parameters of the NMR relaxation for the study of MRI probes [5, 9, 32, 33]:

These two equations show the dependence of the relaxation time constant (T1 and T2) of the electron spin (S), the electronic and proton factors ( and ), the Bohr magneton (β), the distance between the nucleus and the ion (r), the Larmor frequencies for proton and electron spins (ω1 and ω2), the correlation times τC and τS, corresponding to the molecular rotational correlation time and the internal rotational correlation time, respectively, and the hyperfine coupling constant (A), which is described by equation (1) [5, 9, 32, 33].

First, the solvent effect has been evaluated on the Aiso values. For that, theoretical calculations considering ReABT in the gas phase and solution with relativistic effects were performed [34]. Our theoretical findings are presented in Table 3.

Table 3 shows that the differences among Aiso results in solution are small. The larger differences are between the gas phase and solution, for oxygen, about 1.084 MHz, using the M06L functional, 0.543 MHz, for the PBE functional, and 0.083 MHz, when the functional used is PBE0. For hydrogen atoms, these differences are 0.281, 0.106, and 0.009 MHz, for M06L, PBE, and PBE0 functionals, respectively. This result shows that the complex in water, as occurring in the biological environment and human body, can decrease the hyperfine coupling constant values when compared to the result in the gas phase.

A large difference between the Aiso results of Re can be observed, except for the PBE0 functional, where the difference is about 0.365 MHz, which is smaller in the solvent. For the PBE and M06L functionals, the decrease was approximately 26.948 and 60.347 MHz, respectively, when the solvation model was applied.

Despite its great importance, to our knowledge, this is the first investigation of the hyperfine coupling constant values 185|187Re for the ReABT complex. To compare our result, we use the result obtained by Lunsford et al. which points to the Aiso value for 185|187Re in a Re(CO)3 core as 38 MHz [35]. Thus, our findings point out that the calculation using the Minessota 06-L functional in the solution agrees with the result obtained in the literature.

Keeping in mind that although applying the relativistic effects in the Hamiltonian and the basis set leads to an improvement in the accuracy of the calculation, it also leads to an increase in computational demand. With this, we perform a nonrelativistic calculation to investigate the impact of the relativistic effects on the Aiso values of the spectroscopic probe ReABT (Figure 1). Since the previous results show that the solvent effect was important, all calculations were performed, including the solvent effect. These results are presented in Table 4 and, as expected, show poor results for nonrelativistic calculations of 185|187Re hyperfine coupling constant values. So, our findings put in evidence that the relativistic effects are needed for calculating more accurate Aiso values of ReABT.

The results presented in Table 4 point out that, as expected, the Aiso value decreases drastically for rhenium, being 9.179, 24.886, and 36.305 MHz smaller for PBE0, PBE, and M06L functionals, respectively, when the relativistic effect was not included in the calculations. For the oxygen and hydrogen atoms, the Aiso value increased when the relativistic effect was not present. For hydrogen, the increase was 1.180, 0.480, and 0.777 MHz for PBE0, PBE, and M06L functionals, respectively, and for oxygen, the increase was 5.465, 1.929, and 14.344 MHz for PBE0, PBE, and M06L functionals, respectively.

From these Aiso theoretical findings, it is possible to define the best theoretical strategy, combination of functional, solvent model, and relativistic effect. For the spectroscopic probe ReABT, the best combination is using the M06L functional considering relativistic effects on zeroth order regular approximation, with the SARC-ZORA-TZVP basis set for Re and solvent effects evaluated at the CPCM level.

3.2. Thermodynamics Evaluation for the Complex Stability: Toxicologic Properties

Currently, most contrast agents (CAs) used for MRI are Gd-based compounds. They perform well for the improved images provided by MR. However, in the early 2000s, a life-threatening disease, Nephrogenic Systemic Fibrosis, was associated with these CAs, since they release free Gd ions. Since then, the search for new CAs with low toxicity and high performance has intensified [36, 37].

Rhenium-based CAs, such as ReABT, were studied due to their potential as spectroscopy probes and low toxicity [15]. It is important to understand the behavior of the complex in solution, i.e., if the ReABT complex will release free Re (CO)3+ ions in solution. For that, the thermodynamic cycle (Figure 2) has been used to obtain the Gibbs free energy of complexation. Although a study by Breitz et al. shows that compounds with rhenium atoms, despite their accumulation in the kidneys, do not influence renal activity [38], it is important to understand the stability of ReABT complex in solution.

The result obtained for the Gibbs free energy of complexation, −744.37 kJ/mol, at the M06-L level, with the ZORA and SARC-ZORA-TZVP basis sets for rhenium and TZVP for all other atoms, is strongly negative and indicates that the complex is stable in solution. Experimental results in the literature indicate that Gd-DOTA, one of the most common commercial MRI probes, has a Gibbs free energy of complexation around −60 kJ/mol [39]. We are quite aware of the limitations of direct comparison between theoretical and experimental result for the Gibbs free energy of complexation, but as shown in other publications, at least for small sized complexes, the theoretical strategy does not introduce significant changes in the mechanism and energetics of complexation reactions [40, 41]. Therefore, the mechanism of the reaction, as well as the relative energies, can be well established with this model.

Therefore, it can be understood that the ReABT complex is significantly more stable when compared to Gd-DOTA. In this way, the complex ReABT shows low toxicity due to its high stability, which will theoretically not show the dissociation of Re ions. Now, our findings have pointed out that the complex shows high stability in solution; it is possible to move on to the next step of this paper, where the hyperfine coupling constant values for ReABT when docked at the active site of the enzyme will be evaluated.

3.3. Molecular Interaction between ReABT and the Target Protein: Docking Studies

Early diagnosis is crucial for the successful treatment of any disease, especially breast cancer [4]. For that, the use of spectroscopic probes that can interact with specific molecular targets is important. In breast cancer, this molecular target is the enzyme phosphoinositide 3-kinase (PI3K), which is closely related to the growth of cancer cells in this tissue [19]. The compounds of the benzothiazole family are known inhibitors of the enzyme kinase family, especially the PI3K enzyme [7, 17, 42, 43]. In this context, the study of the interaction between the complex ReABT and the enzyme PI3K is important and was investigated in our work using docking calculations.

First, it is possible to check the position of the compound in the active site by taking the overlap and orientation according to the active ligand, N-{6-[(2-(methylsulfanyl) pyrimidin-4yl]- 1, 3-benzothiazol-2-yl} acetamide, a benzothiazole derived (Figure 4).

To continue the analysis of the docking calculation, we can investigate the hydrogen bonds that the ReABT performs inside the active site. In Figure 5, it is possible to observe that ReABT interacts by means of H-bonds, with amino acids residues Val882, Lys883, and Asp884. It is important to note that the active ligand of the PI3K enzyme interacts with the amino acid residue Val882. The intermolecular interaction energy of the active ligand was equal to −121.17 kcal·mol−1 and −163.62 kcal·mol−1 for the ReABT. These results pointed out that the complex ReABT could have better stability in the active site of the PI3K enzyme. Meanwhile, it should be noted that these results do not consider pharmacokinetic properties.

Now, the theoretical methodology defined in the hyperfine coupling constant study has been applied to analyze the change in the Aiso value when the ReABT complex is in the active site of the PI3K enzyme. The calculation was performed using the DTF, with M06L functional, ZORA in the Hamiltonian, SARC-ZORA-TZVP basis sets for Re, and TZVP for all other atoms. This calculation shows that the Aiso value for 187Re decreased from 36.311 MHz to 0.0571 MHZ, while for 1H the decrease was from 2.656 MHz to 0.0236 MHz and from 1.445 MHz to 0.0018 MHz for 17O. This dramatic change indicates that ReABT shows activity on EPR at the active site of PI3K demonstrating its potential for the development of a new spectroscopic probe.

4. Conclusions

In this work, we study a potential spectroscopic probe, the complex (CORe)3(NNO) conjugated to 2-(4′-aminophenyl) benzothiazole, to obtain the EPR activity of this complex in an enzymatic environment. From our data, the best calculation level for obtaining the hyperfine coupling constant (Aiso) for the complex ReABT is with the M06L functional, SARC-ZORA-TZVP basis set for Re and TZVP for all other atoms, considering the relativistic effects in the Hamiltonian with ZORA and using implicit solvation method CPCM. Thermodynamic analysis shows that the ReABT complex is stable and should not release free rhenium into biological systems. In this sense, the toxicological risks associated with the complex are considerably reduced.

The best theoretical strategy has been selected and applied to ReABT. This complex was docked in the PI3K active site, presenting electrostatic interactions with the amino acids residues of alanine, aspartate, and lysine. Our findings put in evidence, in this environment, a huge decrease in the Aiso values for 1H, 17O, and 187Re, showing that the ReABT complex has activity on EPR spectroscopy. This finding reveals that the complex Re (CO)3(NNO) conjugates to 2-(4′-aminophenyl) benzothiazole has the potential to lead to the development of a new spectroscopic probe for diagnosis cases of breast cancer.

Data Availability

Data are available on request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Brazilian financial agencies Conselho Nacional de Desenvolvimento Científico e Tecnologico (CNPq), Fundação de Amparo ao Ensino e Pesquisa de Minas Gerais (FAPEMIG), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), and Federal University of Lavras (UFLA) and LCC-UFLA (physical infrastructure and working space). This work was also supported by Long-Term Development Plan UHK. The authors also would like to thank the preprint publication of this research at ChemRxvi (10.26434/-2022-kkgk0).