Journal of Chemistry
Volume 2017 (2017), Article ID 8124323, 9 pages
https://doi.org/10.1155/2017/8124323

Research Article

On the Interactions of Fused Pyrazole Derivative with Selected Amino Acids: DFT Calculations

Kornelia Czaja,¹ Jacek Kujawski,¹ Elżbieta Jodłowska-Siewert,¹ Paulina Szulc,¹ Tomasz Ratajczak,² Dominika Krygier,³ Marcin K. Chmielewski,³ and Marek K. Bernard¹

¹Department of Organic Chemistry, Faculty of Pharmacy, Poznan University of Medical Sciences, Grunwaldzka 6 Street, 60-780 Poznań, Poland
²Centre of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland
³Institute of Bioorganic Chemistry, Polish Academy of Sciences, Z. Noskowskiego Str. 12/14, 61-704 Poznań, Poland

Correspondence should be addressed to Jacek Kujawski; jacekkuj@ump.edu.pl

Received 21 March 2017; Accepted 18 April 2017; Published 14 May 2017

Academic Editor: James W. Gauld

Copyright © 2017 Kornelia Czaja et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Due to the increasing prevalence of neoplasms, there is a permanent need for new selective cytostatic compounds. Anticancer drugs can act in different ways, affecting protein expression and synthesis, including disruption of signaling pathways within cells. Continuing our previous research aiming at elucidating the mechanism of pyrazole’s anticancer activity, we carried out in silico studies on the interactions of fused pyrazole derivative with alanine, lysine, glutamic acid, and methionine. The objective of the study is to improve our understanding of the possible interactions of pyrazole derivatives with the above-mentioned amino acids. For this purpose, we apply the DFT formalism (optimization using the B3LYP, CAM-B3LYP, PBE0, and M06L functionals) and interaction energy calculations (counterpoise corrected method based on the basis set superposition error, BSSE) together with QTAIM approach and estimation of the ¹H NMR chemical shifts of analyzed pyrazole derivative using different basis sets and DFT functionals in CPCM solvation model (and water used as a solvent).

1. Introduction

Currently, there is a growing demand for new cytostatic compounds that act selectively. Drugs with anticancer properties exert them through different modes of action [1]. These include disruption of DNA replication and thus inhibition of cell proliferation and inhibition of translation of target proteins, such as kinases that regulate signaling pathways within cells. Eventually, these compounds induce apoptosis of neoplastic cells. So far, numerous oncogenes have been described and many of them are dysregulated protein kinases. Decreasing the activity of upregulated kinases results in a stunted tumor growth [2–5]. This is why kinase inhibitors have attracted significant attention for targeted therapy of many cancers.

Heterocycles containing pyrazole are a part of a group of anticancer drugs on which we focused our attention. This type of anticancer compounds includes those which interact with active sites of kinases, for example, through binding with amino acids. Alanine (Ala), lysine (Lys), glutamic acid (Glu), and methionine (Met) existing in these cavities are especially important as they can interact with pyrazole derivatives [2–5].

According to the results obtained by our group previously, certain condensed pyrazole derivatives exert a cytostatic activity and induce apoptosis in the HT29 human colorectal neoplastic cells [6, 7]. In our subsequent papers, we hypothesized that they could disrupt the intracellular Mg²⁺ ions balance through the action of pyridinic nitrogen within the pyrazole moiety [8, 9]. An NMR experiment elucidating the interaction of chosen amino acids and the condensed pyrazole derivative showed that this interaction primarily happens between the pyrrolic nitrogen atom (NH) and carbon fragment B (Scheme 1) [7–9].

Scheme 1

Continuing our studies that aimed at elucidating the potential mode of pyrazole cytostatic activity, we conducted further theoretical studies investigating simulated interactions of fused derivative 7 (synthesis is given in Scheme 1, 1: 4-nitro-1H-pyrazole; 2: (CH₃)₂SO₄ [8]; 3: Ts-CH₂-Cl; [H]: reduction of the nitro group; [HNO₂]: diazotization and cyclization; A-E: protons of the final product considered for NMR investigations) with certain amino acids (Ala, Lys, Glu, and Met). The objective of the study is to improve our understanding of the tendency regarding mechanism of interactions between pyrazole derivatives and selected amino acids. We focused our attention on the electrostatic contacts, especially hydrogen bond formation between the azole and amino acids. To achieve that, we decided to choose the DFT formalism and interaction energy calculations with QTAIM approach and calculation of the ¹H NMR chemical shifts of compound 7 with several basis sets and DFT functionals in conductor-like polarizable continuum model (CPCM) of solvation, so that the influence of the solvent environment in cell is also accounted for. What is more, we aimed to determine the methods and basis sets that gave the closest results to the experimental shifts in the ¹H NMR spectrum and the most stable and preferable conformation of compound 7.

2. Materials and Methods

We carried out density functional calculations and optimized the compound geometries at the DFT [10] level of theory using the Gaussian 09 D.01 program [11], B3LYP/6-31G(d,p) [12], (b) CAM-B3LYP/6-31G(d,p) [13], (c) B3LYP/6-311+G(d,p) [14, 15], (d) PBE0/6-31G(d,p) [16, 17], and (e) M06L/6-31G(d,p) [18] approaches in the CPCM of solvation (water molecules as solvent). The ideal gas, rigid rotor, and harmonic oscillator approximations were applied to calculate the vibrational frequencies and thermodynamic properties. We obtained no imaginary frequencies within the vibrational spectra and therefore confirmed the energy minimum for all conformers. We created the conformers through the rotations of C1-S1-C2 bond of azole and bonds connecting N-H moiety of azole with carboxyl or amine group of the particular amino acid in dihedral angle increments of 20°.

AIMAll Professional package for topological analysis was used for QTAIM calculations [19]. For this purpose, the “Proaim” basin integration approach was used using “superfine” interatomic surface mesh and “very high” outer angular quadrature. Thanks to the atomic (integrated) Laplacian values, we could monitor the accuracy of basin integrations. What is more, the Poincaré-Hopf rule was satisfied for all carried out calculations. Theρ and parameters at the bond critical points (BCPs) concerning dipole polarizability tensors were calculated based on the optimized structures (taking into consideration the already calculated interaction energy) using the AIM approach and applying B3LYP/6-31G(d,p) (adducts 7AlaB, 7LysB, 7GluB, and 7MetB), B3LYP/6-311+G(d,p) (adducts 7GluB6 and 7MetB6), PBE0/6-31G(d,p) (adducts 7LysP), and M06L/6-31G(d,p) (adduct 7AlaM) levels of theory, respectively. The components of molecular dipole polarizability tensors were calculated as the set of first derivatives of the dipole moment components with respect to the components of an applied electric field. We calculated the value of NMR shielding for proton () in TMS at the B3LYP/6-31G(d,p) level of theory (CPCM of solvation and water as solvent). The above-mentioned method was also used for calculation of TMS (reference compound) and the investigated 7–amino acid adducts. Then the values obtained for TMS allowed calculating the chemical shifts of 7 and 7–amino acid adducts, in accordance with the following equation: , with being chemical shift of -nuclei of 7–amino acid adducts and and being the calculated isotropic magnetic shielding tensors of TMS and 7–amino acid adducts, respectively [8, 20–22]. An average of the estimated chemical shifts was calculated for the methyl groups protons A or E and also these of C or D (Scheme 1). We calculated interaction energy of all optimized adducts applying counterpoise corrected method based on the basis set superposition error (BSSE) at B3LYP/6-31G(d,p) and B3LYP/6-311++G(2d,3p) levels of theory, respectively [23, 24].

The ChemCraft 1.7 software was utilized for the visualization of all optimized conformers and adducts [25]. The calculations were carried out using resources provided by the Polish Grid Infrastructure (PL-Grid) and Wrocław Center for Networking and Supercomputing (Bem cluster, WCSS Grant number 327/2014).

3. Results and Discussion

The analysis of azole 7 coordination with magnesium ions proved that the pyrrolic nitrogen NH with a labile exchangeable proton is the principal position in this interaction [8]. Moreover, the pH alterations are unimportant for small molecules that, like compound 7, contain practically one protonation site. Additionally, we prove that such compounds do not tautomerize (tautomerism involving protonation of the nitrogen atoms can be neglected) [7–9]. Based on this information, we decided to examine the interaction of the pyrrolic nitrogen with the following amino acids: alanine, lysine, glutamic acid, and methionine. The geometry of the azole–amino acids adducts was optimized using four different functionals, including hybrid, GGA, and meta-GGA functionals, that is, B3LYP, CAM-B3LYP, PBE0, and M06L. The data concerning the lowest (in bold) and highest interaction energies (CPCM of solvation and water as solvent, counterpoise corrected method including the basis set superposition BSSE error, Gaussian G09 D.01 software [11]) for the azole–amino acids adducts are shown in Table 1.

Table 1: Interaction energies between compound 7 andalanine (7Ala), lysine (7Lys), glutamic acid (7Glu), and methionine (7Met). Suffixes added to the names of adducts mean that the structure was optimized at functional: B, B-C, B-N, B-αN: B3LYP/6-31G(d,p); B6: B3LYP/6-311+G(d,p); C: CAM-B3LYP/6-31G(d,p); P, P-N, P-αN: PBE0/6-31G(d,p); M: M06L/6-31G(d,p). Cartesian coordinates of all adducts are gathered in the Supplementary Material.

Adduct
Interaction energy (kcal mol⁻¹)

B3LYP/6-31G(d,p)
B3LYP/6-311++G(2d,3p)

7AlaB
−9.36
−8.77

7AlaB6
−9.24
−8.67

7AlaC
−9.28
−8.68

7AlaP
−9.33
−8.74

7AlaM
−8.56
−8.07

7AlaB-N
−10.85
−8.73

7LysB
−5.11
−4.73

7LysB6
−14.88
−14.73

7LysC
−15.70
−14.81

7LysP
−15.96
−15.06

7LysM
−15.46
−14.57

7LysP-N
1.05
1.41

7LysP-αN
−8.69
−7.20

7GluB
−15.88
−14.85

7GluB6
−10.90
−9.97

7GluC
−15.53
−14.46

7GluP
−15.88
−14.84

7GluM
−14.30
−13.34

7GluB-αN
−9.70
−7.71

7GluB-C
−9.94
−9.14

7MetB
−13.69
−11.67

7MetB6
−8.63
−8.39

7MetC
−13.19
−11.11

7MetP
−13.19
−11.11

7MetM
−3.08
−2.37

7MetB-αN
−8.42
−7.79

Interaction energy between two atoms or molecules A and B can be generally calculated to be a result of subtraction of the energies of components A and B from the energy of the fused AB structure:

The index reflects the geometry of the fused structure, whereas the index reflects that of separate and . However, this equation often results in incorrect estimates, particularly when hydrogen bonding or dispersion interactions are present within the fused structures. The error associated with dispersion interactions can be lessened by using smaller basis sets, which make the structure more stable than and separately using the basis set superposition error (BSSE). This result is achieved thanks to the expansion of the wavefunction over fewer basis functions for and molecules compared to the structure [23]. Using very big basis sets with the BSSE is impossible for the majority of investigated compounds. Therefore, counterpoise method (CP) can be applied, which can estimate the BSSE size. Then the structure is described in the same way, but and have basis sets which are of equal size, the same as that of the basis set for the structure. Thus, interaction energy is calculated using CP in the following equation:

The indices in (2) show that absolute basis sets for the calculations concerning all the components are the same and are equal to those of the structure [23]. Hence, it can be concluded that such methodology significantly better describes the interaction phenomenon within the analyzed herein adducts.

The energy of interactions within the azole 7–alanine adduct 7AlaB, which was optimized at the B3LYP/6-31G(d,p) level of theory, was –9.36 when calculated at the B3LYP/6-31G(d,p) or –8.77 kcal mol⁻¹ when calculated at the B3LYP/6-311++G(2d,3p) level of theory. It resulted mainly from one strong hydrogen contact O-H⋯N ( Å,  Å, and θ = 160.1°; Figure 1), although the distance between pyrrolic proton (N2-H3 in Figure 1) of azole and hydroxyl group of alanine was 3.433 Å. Moreover, a potential weak hydrogen bond binding the carbonyl oxygen and phenylic proton (d = 2.236 Å, Figure 1) could contribute to some extent to the overall interaction energy.

Figure 1: The optimized 7–alanine adduct 7AlaB (B3LYP/6-31G(d,p)).

A similar geometry was found for the same optimized adduct involving alanine and compound 7 but at the M06L/6-31G(d,p) level of theory (adduct 7AlaM, Figure  S2, in Supplementary Material available online at https://doi.org/10.1155/2017/8124323). The calculated energy equaled −8.56 or −8.75 kcal·mol⁻¹ depending on the B3LYP/6-31G(d,p) or B3LYP/6-311++G(2d,3p) approach, respectively. The distance d between the carboxylic OH group and pyridinic nitrogen was 1.906 Å ( Å and θ = 173.0°). As above, a weak contact involving carbonyl oxygen and phenylic ortho hydrogen was also observed (d = 2.146 Å). In this case, the distance between pyrrolic proton of azole and hydroxyl group of alanine was 3.796 Å.

Although the interactions between amino acid carboxylic group and heterocyclic nitrogen are significantly stronger than the interactions involving amino group of amino acid and heterocyclic nitrogen [20], the adduct 7–Ala was additionally optimized regarding the interaction of the alanine amino group and pyrazole nitrogen atoms (adduct 7AlaB-N, B3LYP/6-31G(d,p) level of theory; Figure  S3, Supplementary Material).

The estimated interaction energy [B3LYP/6-31G(d,p), counterpoise corrected method] was –10.85 kcal/mol. However, this apparently lower energy increased to −8.73 kcal·mol⁻¹ after calculation at the B3LYP/6-311++G(2d,3p) level of theory (counterpoise corrected method). The weaker interaction between the alanine amino group and pyrazole 7 is due to one detected hydrogen bond N-H⋯NH₂ ( Å,  Å, and θ = 161.2°).

The more negative energy (strongest interaction), that is, –15.96 kcal·mol⁻¹ [interaction energy calculated for adduct optimized at the B3LYP/6-31G(d,p) level of theory] or −15.06 kcal·mol⁻¹ [interaction energy calculated for adduct optimized at the B3LYP/6-311++G(2d,3p) level of theory], respectively, concerning the interaction between azole 7 and lysine (7–Lys, Table 1), was observed for adduct 7LysP (Figure 2, geometry optimized at the PBE0/6-31G(d,p) level of theory). This energy originated from the formation of two strong hydrogen bonds: N-H⋯O=C ( Å,  Å, and θ = 147.6°) and O-H⋯N ( Å,  Å, and θ = 169.4°).

Figure 2: The optimized 7–lysine adduct 7LysP (PBE0/6-31G(d,p)).

Instead, the interaction energy for adduct 7LysB (Figure  S4, Supplementary Material; optimized at the B3LYP/6-31G(d,p) level of theory) was merely −5.11 kcal·mol⁻¹ and resulted from the only one strong hydrogen bond being formed: N-H⋯O=C ( Å,  Å, and θ = 147.7°).

The interaction energies between the lysine amine group and azole 7 for optimized adduct 7LysP-N at the PBE0/6-31G(d,p) level of theory (Figure  S5, Supplementary Material) were similar and amounted to a low value of ca 1 kcal·mol⁻¹. This effect was connected with the presence of one relatively weak hydrogen bond: N-H⋯NH₂ ( Å,  Å, and θ = 163.6°).

The interaction between the amine group of lysine and azole 7 pyrrolic nitrogen (adduct 7LysP- N, Figure  S6, Supplementary Material) optimized at the PBE0/6-31G(d,p) level of theory resulted in a strong hydrogen contact formation, N-H⋯NH₂ ( Å,  Å, and θ = 167.7°), with interaction energy of −8.69 [obtained at the B3LYP/6-31G(d,p) level of theory] or −7.20 kcal·mol⁻¹ [obtained at the B3LYP/6-311++G(2d,3p) level of theory].

The above discussion led to a conclusion that adduct 7LysP (optimized at the PBE0/6-31G(d,p) level of theory) involving azole 7 and lysine has the strongest linking due to the interaction between the amino acid carboxylic group and pyrrolic nitrogen of compound 7.

The application of B3LYP/6-31G(d,p) (adduct 7GluB, Figure 3) or PBE0/6-31G(d,p) (adduct 7GluP) approach for optimization of adduct 7–Glu led to rotamers of almost identical energies of −15.88 or −14.85 kcal·mol⁻¹ using B3LYP/6-31G(d,p) or B3LYP/6-311++G(2d,3p) approaches (Table 1). The most significant contribution to the overall interaction energy of adducts 7GluB and 7GluP came from two strong hydrogen bonds, namely, N-H⋯O=C ( Å,  Å, and θ = 144.5°) and O-H⋯N ( Å,  Å, and θ = 169.7°). Conversely, the interaction energy for adduct 7GluB6 (Figure  S7, Supplementary Material; geometry optimized at the B3LYP/6-311+G(d,p) level of theory) was notably lower, that is, −10.9 kcal·mol⁻¹ (energy determined at the B3LYP/6-31G(d,p) level of theory, Table 1). Here, the most important contribution originated from one strong hydrogen bond: N-H⋯O=C ( Å,  Å, and θ = 148.5°). Adduct 7GluB-C (optimized at the B3LYP/6-31G(d,p) level of theory, Figure  S9, Supplementary Material) regarding interaction of azole with peripheral carboxyl group of Glu had similar interaction energy values between −9 and −10 kcal·mol⁻¹ (Table 1) which derived from one strong hydrogen bond: N-H⋯O=C ( Å,  Å, and θ = 175.3°).

Figure 3: The optimized 7–glutamic acid adduct 7GluB (B3LYP/6-31G(d,p)).

Analogously to the above-described complex of compound 7 with lysine, the interaction energy involving the glutamic acid α amino group and pyrazole nitrogen [adduct 7GluB- N, B3LYP/6-31G(d,p) level of theory, Figure  S8, Supplementary Material] was significantly lower, namely, −9.7 [B3LYP/6-31G(d,p) level of theory] or −7.71 [B3LYP/6-311++G(2d,3p) approach] kcal mol⁻¹. Its relatively low value came from the presence of only one hydrogen contact: N-H⋯NH₂ ( Å,  Å, and θ = 163.8°).

The highest (considering absolute value) interaction energy (Table 1) for the complex of azole 7 with methionine (7–Met) was obtained for adduct 7MetB optimized at the B3LYP/6-31G(d,p) level of theory and amounted to −13.69 [B3LYP/6-31G(d,p) level of theory] or −11.67 [B3LYP/6-311++G(2d,3p) level of theory] kcal mol⁻¹. This structure was stabilized by two strong hydrogen bonds involving the pyrrolic nitrogen atom (N2-H3 in Figure 4) of azole 7 and amine group (N5 in Figure 4) of amino acid or the hydroxyl function of methionine and sulfonyl group of condensed azole 7, namely, N-H⋯NH₂ ( Å,  Å, and θ = 148.6°) or O-H⋯O₂S ( Å,  Å, and θ = 144.6°).

Figure 4: The optimized 7–methionine adduct 7MetB (B3LYP/6-31G(d,p)).

Another structure, 7MetB6, of the adduct 7Met was obtained using B3LYP/6-311+G(d,p) level of theory (Figure  S10, Supplementary Material). The interaction energy calculated at B3LYP/6-31G(d,p) or B3LYP/6-311++G(2d,3p) level of theory was lower than that described in the preceding paragraph, that is, −8.63 or −8.39 kcal mol⁻¹, respectively, and resulted mainly from one hydrogen contact between the methionine carboxylic group and sulfone functionality: O-H⋯O₂S ( Å,  Å, and θ = 172.3°).

The optimization of the adduct 7MetB- N at the B3LYP/6-31G(d,p) level of theory (Figure  S11, Supplementary Material) resulted in a relatively low interaction energy, that is, −8.42 [interaction energy calculated at the B3LYP/6-31G(d,p) level of theory] or −7.79 kcal mol⁻¹ [interaction energy calculated at the B3LYP/6-311++G(2d,3p) level of theory]. The interaction involved the methionine α amine group and the condensed pyrazole system led to two relatively weak hydrogen bonds: N-H⋯NH₂ ( Å,  Å, and θ = 106.9°) or NH₂⋯O₂S ( Å,  Å, and θ = 142.2°).

3.1. AIM Analysis

The application of AIM method includes determination of type of bond (ionic or covalent) and its properties [26] also when it comes to hydrogen bonding. In our AIM calculations, we focused on the interactions of the amino acids carboxylic group with condensed azole 7. Particular attention was devoted to adducts with the highest and lowest interaction energies.

As a principle, there is a relation between the parameter at the bond critical point (BCP) and the strength of the hydrogen bond. The fact whether the Laplacian at the BCP is positive or negative determines whether the hydrogen bond dominates the covalent () or electrostatic () interactions. The values for the electron density topological parameters and (related with a physical basis via local statement of the virial theorem for bonding interactions [26]) for the analysed adducts are given in Tables 2 and  S1–S3 (Supplementary Material). It is noteworthy that the AIM parameters in the literature were correlated with the strength of analyzed bond [27]; however, in cited paper only one bond within the structure of the analyte was considered, whereas in our models we have detected even two hydrogen contacts within the structure of the particular azole–amino acid rotamer being analyzed.

Table 2: Electron densities and Laplacian fields [eÅ] determined for the 7Lys adducts with a hydrogen bond; nd: no data.

Parameter
Adduct
Hydrogen bond

N-H⋯O=C
O-H⋯N

7LysB
0.032433
0.048718

7LysB
0.096056
0.108884

7LysP
0.035598
nd

7LysP
0.103572
nd

In case of all analyzed adducts, the absolute value of bond ellipticity () was ca 0.01–0.04, which indicates cylindrical, directional character of investigated contacts.

The values were positive but relatively small; thus they indicate that the adducts were connected through strong hydrogen bonds of electrostatic character. Moreover, the hydrogen bonds of adducts azole 7–amino acid with lower (more negative) interaction energies have higher values of electron densities (Table 1). This demonstrates that adducts 7AlaB, 7LysB6, LysM, and 7MetB are connected by stronger hydrogen bonds in comparison with adducts 7AlaM, 7LysP, 7GluB6, and 7MetB6, respectively.

3.2. NMR Analysis

The mechanism of interactions of metal ions and organic compounds, such as synthesized drugs and proteins, is of utmost significance from the biological point of view. Various analytical methods can be applied to investigate them, including ¹H NMR spectroscopy. Noteworthy is the fact that NMR spectroscopy has found its use for the determination of interactions of some ligands with proteins and allows assessing quantitatively these interactions if the equilibrium constants and changes in chemical shifts are taken into account [8, 20, 21, 26, 28–31]. As a result, we decided to focus only on ¹H NMR spectroscopy in our study. Albeit compound 7 has four nitrogen atoms within its structure, the low abundance of the ¹⁵N isotope in comparison with the ¹H isotope as well as the significant signals broadening due to the large quadrupole moment of ¹⁴N renders the nitrogen NMR spectroscopy impractical [8, 30, 31]. Considering the low utility of ¹³C NMR spectroscopy for the investigations of interactions between metal ions and small molecule organic compounds [8, 9, 29] as well as the low natural abundance of ¹³C, we decided to use ¹H NMR for the studies depicted in our paper. The theoretical ¹H NMR spectrum was generated for all azole–amino acid adducts using the gauge-including atomic orbital (GIAO) method [22] as implemented in Gaussian G09 D.01. In order to maintain clarity, we considered four adducts with the lowest (more negative) energy, previously optimized in CPCM of solvation (water as solvent) at the B3LYP/6-31G(d,p) (Tables  S4, S9, S13, and S18; Supplementary Material), B3LYP/6-311+G(d,p) (Tables  S5, S10, S14, and S19; Supplementary Material), CAM-B3LYP/6-31G(d,p) (Tables  S6, S11, S15, and S20; Supplementary Material), PBE0/6-31G(d,p) (Table 1 and Tables  S7, S16, and S21; Supplementary Material), and M06L/6-31G(d,p) (Tables  S8, S12, S17, and S22; Supplementary Material) level of theory.

The theoretical values of chemical shifts correlate well with the experimental ones for pyrazole derivative 7, which is shown by low mean absolute deviation (MAD) [8]. These low MAD values for the theoretical ¹H NMR spectra of the adducts azole 7–amino acids (Tables 3 as well as S4, S16, and S18; Supplementary Material) correspond to the method and basis used for optimization of the adducts 7AlaB (B3LYP/6-31G(d,p) approach), 7LysP (PBE0/6-31G(d,p) approach), 7GluB (B3LYP/6-31G(d,p) approach), and 7MetB (B3LYP/6-31G(d,p) approach) with the highest (more negative) interaction energies calculated at the B3LYP/6-311++G(2d,3p) level of theory (Table 1) and emphasize the importance of hydrogen bonds within these adducts.

Table 3: Experimental () and calculated chemical shifts (rotamers ) for compound 7 within the 7Lys adduct; absolute errors (), average absolute error (), and relative percentage errors (); calculated NMR shielding for proton  ppm for TMS (PBE0/6-31G(d,p)/GIAO); MAD = 0.27.

Locant

7LysP

A
2.37
2.30
2.29
2.32
2.32
0.08
0.09
0.06
0.06
0.07
3

B
7.65
7.31
7.32
7.30
7.31
0.34
0.33
0.35
0.34
0.34
4

NH
13.67
13.16
13.16
13.17
13.11
0.51
0.51
0.50
0.56
0.52
4

C
7.85
8.26
8.24
8.26
8.28
0.41
0.39
0.41
0.42
0.41
5

D
7.43
7.49
7.48
7.51
7.51
0.05
0.05
0.08
0.08
0.07
1

E
4.08
4.32
4.33
4.31
4.31
0.24
0.26
0.23
0.24
0.24
6

From theoretical point of view, fascinating conclusions can be drawn from the calculations involving spin-spin coupling constants between the nuclei of hydrogen bond donor and acceptor as calculated at the B3LYP/6-311++G(2d,3p) level of theory. The conclusions are related with the adducts with the highest (more negative) and lowest (more positive) interaction energies (Figures  S12–S19 and Table  S23, Supplementary Material). For the interactions of azole 7 and alanine, the coupling constants involving atoms connected by hydrogen bonds O-H⋯N have values of ca −2 Hz. The increase in length of this bond in adducts 7AlaB and 7AlaM leads to an increase in both the torsion angle of the hydrogen bond donor and negative value of the dihedral angle . Conversely, the increase in length of C-H⋯O=C bond (adducts 7AlaM and 7AlaB) resulted in an insignificant decrease of the angle and more negative value of the dihedral angle as well as a decrease in coupling constant from 2 to 1 Hz. It should be added that the hydrogen contacts in adduct 7AlaB are characterized by greater values of electron density calculated by the QTAIM method (Table  S1, Supplementary Material).

The N-H⋯O=C hydrogen contact in azole 7–lysine adduct (adducts 7LysB and 7LysP) can be characterized by comparable values of the bond length, torsion angle, and coupling constant for hydrogen bond acceptor–donor (ca 5 Hz), as well as electron density (Table 2). However, the dihedral angle τ has a negative value for adduct 7LysB.

The higher interaction energy of adduct 7LysP in comparison with complex 7LysB is probably due to the already noted smaller distance of the hydroxyl group (carboxylic function) of lysine from the pyrazole N2 nitrogen atom, as well as the second hydrogen contact O-H⋯N.

Comparison of the interaction energy of 7Met adducts 7MetB and 7MetB6 reveals that the higher energy of adduct 7MetB is connected with the presence of a second hydrogen bond O-H⋯O₂S, despite the fact that this contact has lower electron densityρ and is longer compared to adduct 7MetB6 (Tables  S3 and S23, Supplementary Material). However, this is associated with a smaller value of torsion angle and dihedral angle as well as coupling constant that decreases to ca 4 Hz (adduct 7MetB) in comparison with 7 Hz (adduct 7MetB6).

The interaction of azole 7 with glutamic acid leads to adducts 7GluB and 7GluB6 that are linked through a N-H⋯O=C hydrogen bond. This hydrogen contact has similar values of the bond length, torsion angle , and coupling constant (ca 4 Hz, Table  S23, Supplementary Material), although the dihedral angle (O-H⋯N) in adduct 7GluB is shorter and has somewhat higher value of torsion angle and lower coupling constant, that is, −1.18 Hz versus −9.38 Hz, in comparison with the same bond parameters in adduct 7GluB6. It should be added that the electron densities in the paths of both hydrogen bonds between the nuclear attractors of adduct 7GluB attain higher values.

4. Conclusions

The paper has focused attention on the interactions of the cytotoxic condensed derivative of pyrazole 7 with selected amino acids present in the receptor cavities of human protein kinases. The most significant interactions involve the amino acids carboxylic groups with azole 7 pyrrolic nitrogen. Results of interaction energy calculations among the 7–amino acid adducts regarding method and basis set used for initial optimization protocol for each of the analyzed adducts are gathered in Table 1. From this point of view, we observed that for 7Ala, 7Glu, and 7Met adducts the optimization at the B3LYP/6-31G(d,p) level of theory was undoubtedly a reasonable choice. However, for the 7Lys adduct, the PB0 functional resulted with optimized structure, within which more significant interactions were detected.

We have shown that the low MAD values for the theoretical ¹H NMR spectra of the adducts azole 7–amino acids correspond to the method and basis used for optimization of the adducts 7AlaB, 7LysP, 7GluB, and 7MetB6 with the highest interaction energies.

The results presented in the previous sections show that, for adducts with similar geometric parameters and coupling constants, the electron densitiesρ of the paths of bonds between nuclear attractors (adducts 7Lys) have a decisive influence on the interaction strength. For the remaining adducts, that is, 7Ala, 7Glu, and 7Met, a second hydrogen bond significantly affects the interaction. Furthermore, an increase in the hydrogen bond length in these adducts is accompanied by a decrease in the value of donor torsion angle and coupling constant.

Data Access

Cartesian coordinates and visualization of all formed adducts and results of NMR and QTAIM calculations are available via the Internet or upon request from the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The calculations were carried out using resources provided by the Polish Grid Infrastructure (PL-Grid) and Wrocław Center for Networking and Supercomputing (Bem cluster, WCSS Grant no. 327/2014). The authors are grateful to Dr. Jarosław J. Panek (Faculty of Chemistry, University of Wrocław, Poland) and Dr. Dariusz Kędziera (Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Poland) for help in the NMR and DFT calculations and for manuscript correction.

References

C. Avendaño and J. C. Menendéz, Medicinal Chemistry of Anticancer Drugs, Elsevier, Amsterdam, The Netherlands, 2nd edition, 2015.

S. S. Abd El-Karim, M. M. Anwar, N. A. Mohamed, T. Nasr, and S. A. Elseginy, “Design, synthesis, biological evaluation and molecular docking studies of novel benzofuran-pyrazole derivatives as anticancer agents,” Bioorganic Chemistry, vol. 63, pp. 1–12, 2015.

W. Aman, J. Lee, M. Kim, S. Yang, H. Jung, and J.-M. Hah, “Discovery of highly selective CRAF inhibitors, 3-carboxamido-2H-indazole-6-arylamide: in silico FBLD design, synthesis and evaluation,” Bioorganic & Medicinal Chemistry Letters, vol. 26, pp. 1188–1192, 2016.

L. Ye, X. Ou, Y. Tian et al., “Indazoles as potential c-met inhibitors: design, synthesis and molecular docking studies,” European Journal of Medicinal Chemistry, vol. 65, pp. 112–118, 2013.

H.-P. Buchstaller, H.-M. Eggenweiler, C. Sirrenberg et al., “Fragment-based discovery of hydroxy-indazole-carboxamides as novel small molecule inhibitors of Hsp90,” Bioorganic and Medicinal Chemistry Letters, vol. 22, no. 13, pp. 4396–4403, 2012.

E. Totoń, E. Ignatowicz, M. K. Bernard, J. Kujawski, and M. Rybczyńska, “Evaluation of apoptotic activity of new condensed pyrazole derivatives,” Journal of Physiology and Pharmacology, vol. 64, no. 1, pp. 115–123, 2013.

J. Kujawski, K. Czaja, M. Doskocz et al., “On the interactions of indazole derivative with nucleosides - toward modeling the cytotoxic activity mechanism,” Computational and Theoretical Chemistry, vol. 1059, pp. 45–50, 2015.

K. Czaja, J. Kujawski, U. Girreser, J. J. Panek, M. Doskocz, and M. K. Bernarda, “Possible interactions between fused pyrazole derivative and magnesium ions - NMR experiments and theoretical calculations,” Arkivoc, vol. 4, pp. 22–43, 2016.

J. Kujawski, M. Doskocz, H. Popielarska et al., “Interactions between indazole derivative and magnesium cations-NMR investigations and theoretical calculations,” Journal of Molecular Structure, vol. 1047, pp. 292–301, 2013.

P. Hohenberg and W. Kohn, “Inhomogeneous electron gas,” Physical Review, vol. 136, pp. B864–B871, 1964.

M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09 D.01. Revision A.1, Gaussian, Inc, Wallingford, Conn, USA, 2009.

A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.

T. Yanai, D. Tew, and N. A. Handy, “A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP),” Chemical Physics Letters, vol. 393, pp. 51–57, 2004.

C. Lee, W. Yang, and R. G. Parr, “Development of the colle-salvetti correlation-energy formula into a functional of the electron density,” Physical Review, vol. 37, no. 2, pp. 785–789, 1988.

F. Eckert and A. Klamt, “Fast Solvent Screening via Quantum Chemistry: COSMO-RS Approach,” AIChE Journal, vol. 48, no. 2, pp. 369–385, 2002.

J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996.

J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple [Phys. Rev. Lett. 77, 3865 (1996)],” Physical Review Letters, vol. 78, no. 7, p. 1396, 1997.

Y. Zhao and D. G. Truhlar, The Journal of Chemical Physics, vol. 125, no. 194101, pp. 1–18, 2006.

T. A. Keith, AIMAll, Version 14.11.23, TK Gristmill Software, Overland Park, Kan, USA, 2014.

J. Kujawski, M. K. Bernard, E. Jodłowska, K. Czaja, and B. Drabińska, “On the interactions of leflunomide and teriflunomide within receptor cavity — NMR studies and energy calculations,” Journal of Molecular Modeling, vol. 21, p. 105, 2015.

J. Kujawski, K. Czaja, E. Jodłowska et al., “Structural and spectroscopic properties of econazole and sulconazole - Experimental and theoretical studies,” Journal of Molecular Structure, vol. 1119, pp. 250–258, 2016.

K. Wolinski, J. F. Hinton, and P. Pulay, “Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations,” Journal of the American Chemical Society, vol. 112, no. 23, pp. 8251–8260, 1990.

S. F. Boys and F. Bernardi, “The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors,” Molecular Physics, vol. 19, article 553, 1970.

A. Famulari, R. Specchio, M. Sironi, and M. Raimondi, “New basis set superposition error free ab initio MO-VB interaction potential: molecular-dynamics simulation of water at critical and supercritical conditions,” The Journal of Chemical Physics, vol. 108, article 3296, 3303 pages, 1998.

http://www.chemcraftprog.com.

R. F. W. Bader, Atoms in Molecules: A Quantum Theory, Clarendon, Oxford, UK, 1994.

S. J. Grabowski, “An estimation of strength of intramolecular hydrogen bonds — ab initio and AIM studies,” Journal of Molecular Structure, vol. 562, pp. 137–143, 2001.

J. Kujawski, K. Czaja, T. Ratajczak, E. Jodłowska, and M. K. Chmielewski, “Investigations on synperiplanar and antiperiplanar isomers of losartan: Theoretical and experimental NMR studies,” Molecules, vol. 20, no. 7, pp. 11875–11890, 2015.

L. Fielding, “NMR methods for the determination of protein–ligand dissociation constants,” Progress in Nuclear Magnetic Resonance Spectroscopy, vol. 51, article 219, 2007.

F. Wada, T. Kadonaga, and S. Matsuo, “15N abundance in nitrogen of naturally occuring substances and global assessment of denitrification from isotropic viewpoint,” Geochemical Journal, vol. 9, article 139, 1975.

E. Anders and N. Grevesse, “Abundances of the elements: meteoritic and solar,” Geochimica et Cosmochimica Acta, vol. 53, article 197, 1998.