Pedro de Sena Murteira
Pinheiro
ab,
Daniel Alencar
Rodrigues
ac,
Marina Amaral
Alves
ac,
Luzineide Wanderley
Tinoco
bd,
Glaucio Braga
Ferreira
e,
Carlos Mauricio Rabello
de Sant’Anna
af and
Carlos Alberto Manssour
Fraga
*ab
aLaboratório de Avaliação e Síntese de Substâncias Bioativas (LASSBio), Instituto de Ciências Biomédicas, Universidade Federal do Rio de Janeiro, PO Box 68023, 21941-902, Rio de Janeiro, RJ, Brazil. E-mail: cmfraga@ccsdecania.ufrj.br
bPrograma de Pós-Graduação em Farmacologia e Química Medicinal, Instituto de Ciências Biomédicas, Universidade Federal do Rio de Janeiro, 21941-902, Rio de Janeiro, RJ, Brazil
cPrograma de Pós-Graduação em Química, Instituto de Química, Universidade Federal do Rio de Janeiro, 21941-909, Rio de Janeiro, RJ, Brazil
dLaboratório de Análise e Desenvolvimento de Inibidores Enzimáticos (LADIE), Instituto de Pesquisa em Produtos Naturais (IPPN), Universidade Federal do Rio de Janeiro, 21941-902, Rio de Janeiro, RJ, Brazil
ePrograma de Pós-Graduação em Química, Instituto de Química, Universidade Federal Fluminense, Niterói, Rio de Janeiro, 24020-141, Brazil
fDepartamento de Química, Instituto de Ciências Exatas, Universidade Federal Rural do Rio de Janeiro, 23970-000, Seropédica, RJ, Brazil
First published on 28th November 2017
Sigma-hole (σ-hole) bonds are interactions that are gaining special attention in medicinal chemistry. This type of interaction, initially assigned to the halogens (group 17 of the periodic table), has been extended to atoms of groups 14, 15 and 16. Sulfur atoms have been outstanding for describing these interactions at the intramolecular level (to induce conformational stability) and the intermolecular level (participating in molecular recognition of bioactive compounds by their respective targets). Thus, this work describes the theoretical and experimental characterization of a 1,4-N⋯S σ-hole intramolecular interaction in the N-acylhydrazone cardioactive prototype LASSBio-294 (1), which leads to conformational stabilization and has a direct influence on the molecular properties of this inotropic prototype compared to a negative control for the interaction, LASSBio-897 (2), which is the regioisomer at the thiophene ring. Our theoretical results were reached using the B3LYP/6-311G(d) level of theory, including analysis of conformational, orbital and electrostatic properties. We performed experimental studies using IR, Raman, UV and NMR spectroscopies, which corroborated our theoretical data, showing significant differences between LASSBio-294 (1) and LASSBio-897 (2) in relation to the bond strength of the groups involved in the N⋯S interaction (S–C and NC bonds), the energies of the orbitals associated with the S lone pair (Lp(S)) and the antibonding NC π bond (π*(NC)), as well as the 15N chemical shifts in both systems. Together, our results show how this unusual interaction can influence the molecular properties of some organic compounds.
The importance of the σ-hole interactions involving elements of group 17, also called halogen bonds, has been recognized and exploited in the rational design of drug candidates.3,5,6 Some interesting examples of this interaction are involved in the molecular recognition of bioactive compounds by their respective biotargets.5,7–9
Other atoms are also attracting interest for their σ-hole interactions, including sulfur, where σ-hole regions are present on both sides of the atom (opposite the C–S bonds) allowing both intramolecular and intermolecular interactions with Lewis bases.10,11 As can be observed in the work published by Burling and Goldstein,12 who performed an investigation in the Cambridge Structural Database (CSD) and identified 78 structures containing a 1,4-N⋯S relationship, among which 69% showed the two atoms separated by a distance smaller than the sum of their van der Waals radii.
These interactions, mediated by sulfur atoms, were recently reviewed by Beno and co-workers (2015),3 highlighting that these σ-hole intramolecular interactions can induce conformational stabilization of compounds, which can consequently influence the interaction mode of bioactive compounds with their respective targets.13–15
Murray and colleagues (2008)16 studied this interaction, and these researchers analyzed the influence of electron-withdrawing and electron-donating groups in the sulfur σ-hole regions theoretically. In their report, electron withdrawing groups linked to a thiazole ring, such as –CN and –Cl, increased the electrostatic potential of the σ-hole regions; presence of electron-donating groups, such as –NH2, decreased the electrostatic potential of these regions, promoting a direct effect on the interaction energy between the sulfur atoms and Lewis bases.
Similar reports describing interactions with other groups that can act as Lewis bases were also published. Markham and Bock (1997) carried out theoretical research on 1,4-O⋯S intramolecular interactions and found that these interactions are very sensitive to the electronegative character of the group linked to the sulfur atom.17 Also related to the O⋯S intramolecular interactions, Rittner and colleagues (2011) performed an investigative study on the conformational preferences of 5-substituted 2-acetyl-thiophenes using theoretical calculations and infrared spectroscopy.18 They determined a population of approximately 80% for the O,S-cis conformation, highlighting the importance of this type of interaction in fixing specific conformations.18 Moreover, Zhang and colleagues (2015) reported the characterization of O⋯S intermolecular interactions through the use of theoretical calculations and statistical analysis in derivatives presenting thiophene, thiazole and thiadiazole rings.11 They showed increased stability of the O⋯S intermolecular interaction when the sulfur atom is present in thiazole and thiadiazole rings, which may be a consequence of increasing both the electrostatic potential and the size of the σ-hole regions.11
Thus, the importance of σ-hole interactions involving sulfur atoms is just now coming to light. Our group is interested in these interactions in bioactive compounds such as the cardioactive compound LASSBio-294 (1) (Fig. 1).19 This N-acylhydrazone derivative was developed by our research group and shows potent inotropic effects in isolated hearts of Wistar rats20 and moderate vasodilator activity in the intact aortic rings of Wistar rats.21 It is also able to prevent exercise intolerance and cardiac remodeling in rats submitted to myocardial infarction.22,23 These data make this compound an interesting drug candidate for the treatment of chronic heart failure. Many efforts were performed to optimize LASSBio-294 (1), but none were successful. One of the most interesting examples confirms that optimization of LASSBio-294 (1) is difficult, namely, the analysis of its regioisomer at the thiophene ring, LASSBio-897 (2) (Fig. 1). LASSBio-897 (2) showed an improved vasodilatory activity, but no inotropic effect.24 Based on recent publications3,14,25,26 and data about LASSBio-294 (1) and LASSBio-897 (2), we believe that the 1,4-N⋯S relationship present in the structure of LASSBio-294 (1) is important for a N⋯S σ-hole intramolecular interaction and possibly for the range of pharmacological activities evidenced for this cardioactive prototype. Beyond that, N-acylhydrazones are known privileged structures that can be used in the pharmacological modulation of several biological targets.27–29
Thus, in this work, we investigated, for the first time, the possibility of a N⋯S σ-hole intramolecular interaction in the N-acylhydrazone derivative LASSBio-294 (1), through the combined use of molecular modeling and spectroscopic techniques.
We first performed theoretical calculations to evaluate the intramolecular interaction and then performed experimental evaluations to validate the theoretical data. At all stages, we used LASSBio-897 (2) as a negative control for this intramolecular interaction, due to the great structural similarity between both compounds; the sulfur position in the thiophene ring of LASSBio-897 (2) makes the N⋯S interaction impossible.
To validate the conformational theoretical data, we conducted NOESY-2D analyses, where it was possible to observe spatial interactions between the hydrogens in LASSBio-294 (1) (A = δ 8.62, 7.43; B = δ 11.63, 8.62; C = δ 11.63, 7.47; D = δ 11.63, 7.43) and LASSBio-897 (2) (A = δ 8.50, 7.51; B = δ 11.61, 8.50; C = δ 11.61, 7.53; D = δ 11.63, 7.44), which corroborated our theoretical data (Fig. 2E and F). The interaction A in LASSBio-294 (1) is related to the conformation of the thiophene ring in relation to the imine moiety, which demonstrated the stability of the syn conformer, as predicted by the PES scan analysis. This NOESY-2D analysis showed the proximity of the nitrogen and sulfur atoms in DMSO for LASSBio-294 (1), that is another indication of the presence of a σ-hole 1,4-N⋯S intramolecular interaction. In addition, it was possible to identify two conformers in relation to the benzodioxole for LASSBio-294 (1) and LASSBio-897 (2) (Fig. 2E and F), which according to our theoretical data, it is probably due to the low energy rotational barrier (∼3 kcal mol−1) and to the small energy difference between these two conformers (∼0.3 kcal mol−1) (Fig. 2B).
σ-hole interactions are described in the literature both (a) from an electrostatic point of view, using Coulomb's Law,2,43 where the nitrogen atom carries a partial negative charge and the sulfur atom carries a partial positive charge; and (b) from the orbital point of view, as an interaction between an isolated electron pair, for example the electron lone pair of an imine nitrogen, and one of the σ* orbitals of the sulfur–carbon bonds.3,4,18
According to the Mulliken, NBO and CHELPG charges (Table 1), the sulfur atom can behave as an electrophile in both LASSBio-294 (1) and LASSBio-897 (2), with a more positive partial charge in LASSBio-294 (1) according to the Mulliken and NBO charges, which is a consequence of the greater delocalization of the perpendicular sulfur lone pair through the thiophene ring. The charge data also showed that the imine nitrogen can act as a nucleophile in both systems. However, the Mulliken charges showed that the imine nitrogen of LASSBio-294 (1) is more nucleophilic than LASSBio-897 (2), while the NBO charges indicated that the imine nitrogen has the same charge in both systems and the CHELPG charges indicated that the imine nitrogen of LASSBio-294 (1) is less nucleophilic than LASSBio-897 (2). According to the 15N chemical shift values detected by (1H–15N) HMBC experiments the imine nitrogen of LASSBio-294 (1) is more shielded than in LASSBio-897 (2) (Table 1), which is only corroborated by the Mulliken charges. In addition, chemical shift values indicated that the amide nitrogen atoms present the same charge in both systems, which was confirmed by the Mulliken's charge calculations, but it was not confirmed by NBO and CHELPG charges (Table 1). In general, regardless of the type of charge used, the theoretical data indicated that in LASSBio-294 (1) the imine nitrogen atom has a negative partial charge and the sulfur atom has a positive partial charge, which reinforces the electrostatic point of view of the interaction between these two atoms. This is also corroborated by the chemical shielding effect of the imine nitrogen of LASSBio-294 (1) in relation to LASSBio-897 (2), since this effect is similar to other types of interactions, such as intramolecular hydrogen bonds, which has a strong electrostatic character.44,45 These data further corroborate our theoretical and experimental observations.
Imine 15N chemical shift (ppm) | −67.4 | −63.0 |
Amide 15N chemical shift (ppm) | −212.2 | −212.0 |
Imine N Mulliken charges | −0.134 | −0.126 |
Imine N NBO charges | −0.227 | −0.227 |
Imine N CHELPG charges | −0.200 | −0.224 |
Amide N Mulliken charges | −0.492 | −0.492 |
Amide N NBO charges | −0.421 | −0.423 |
Amide N CHELPG charges | −0.289 | −0.237 |
S Mulliken charges | 0.326 | 0.300 |
S NBO charges | 0.444 | 0.405 |
S CHELPG charges | 0.023 | 0.037 |
From an orbital point of view, the Mulliken bond order analysis showed a direct interaction between the nitrogen and sulfur atoms (0.060) for LASSBio-294 (1), indicative of a weak interaction between these atoms. In addition, it was possible to identify an electronic transition area through the analysis of the electrostatic potential map of LASSBio-294 (1) (Fig. 3). We also performed NBO analysis to observe interactions between the Lp(N) orbital and the σ*(S–C) (Lp(N) → σ*(S–C)) in LASSBio-294 (1). However, an interaction between these orbitals was not observed, suggesting that in LASSBio-294 (1) this interaction energy is lower than 0.5 kcal mol−1. It is expected that in LASSBio-294 (1) the interaction has more electrostatic character, possibly because of the 1,4-N⋯S relationship, which prevents effective orbital overlap.3
Fig. 3 Molecular electrostatic potential map of LASSBio-294 (1) (B3LYP/6-311G(d)) at the 0.002 a.u. contour level. The map are in the range of +100 to −100 kJ mol−1. |
We also observed the orbital energy differences of the HOMO, classified as sulfur lone pair – Lp(S), and the LUMO, the π*(NC) orbital in both systems, through the UV spectra calculations and NBO analysis. In LASSBio-294 (1), an increase in the Lp(S) orbital energy was observed compared to the Lp(S) orbital of LASSBio-897 (2) (Fig. 4B and C), which was probably due to the field effect caused by the spatial approximation of the imine nitrogen atom.
In addition, due to the increased energy of the HOMO (Lp(S)) orbital in LASSBio-294 (1), there is greater delocalization of this electron pair toward the thiophene ring, as can be seen in Table 2, in which the Lp(S) orbital of LASSBio-294 (1) interacts with the π* adjacent orbitals of the thiophene ring with ∼1.8 kcal mol−1 higher in relation to LASSBio-897 (2). Consequently, there is also a stronger resonance between the imine moiety and the thiophene ring of LASSBio-294 (1) (∼2 kcal mol−1 higher) in comparison to LASSBio-897 (2) (Table 2). With this stronger resonance, it would be expected a lower energy of the LUMO orbital for LASSBio-294 (1) in relation to LASSBio-897 (2), due to the greater electron donation to the π*(NC) orbital, which was also corroborated by the theoretical data (Fig. 4B and C). Consequently, there is a lower energy gap between the HOMO and LUMO orbitals in LASSBio-294 (1) in comparison to LASSBio-897 (2), being expected a longer λMAX absorption in the UV spectra analysis for LASSBio-294 (1). These theoretical analysis were corroborated by our experimental data, where a bathochromic shift of the λMAX of LASSBio-294 (1) was observed in relation to the λMAX of LASSBio-897 (2), highlighting the ΔE(HOMO–LUMO) difference of both systems (Fig. 4A). Our theoretical UV spectra calculations confirmed that the λMAX absorption bands are related to HOMO–LUMO transitions in both LASSBio-294 (1) (oscillator strength = 0.8242) and LASSBio-897 (2) (oscillator strength = 0.8146) (see ESI,† Section 4 and Table S17).
Although our theoretical analysis indicated that the 1,4-N⋯S σ-hole intramolecular interaction in LASSBio-294 (1) has higher electrostatic character, it was possible to experimentally observe the presence of a relevant orbital character in this interaction using IR and Raman spectroscopic analysis. Using theoretical spectra calculations, it was possible to perform a total energy distribution analysis36 to determine the main bands related to the vibrations of the S–C and NC bonds (see ESI,† Section 4, Tables S15 and S16).46 Significant differences were observed in the main stretching wavenumbers of the S–C and NC bonds in LASSBio-294 (1) and LASSBio-897 (2) (Fig. 5). The S–C stretching wavenumbers in LASSBio-294 (1) appeared at 758 cm−1, whereas in LASSBio-897 (2) it appeared at 848 cm−1. In relation to the imine bonds (NC) it is possible to observe that in LASSBio-294 (1), the stretching wavenumber appeared at 1592 cm−1 while in LASSBio-897 (2) it appeared at 1606 cm−1, in agreement with the observations of NBO analysis for the resonance between the thiophene ring and the imine moiety of the N-acylhydrazone subunit (Table 2). The lower stretching wavenumber of the NC bond in LASSBio-294 (1) is experimentally indicative of the stronger resonance, as there is higher electronic donation to the π*(NC) orbital (∼3.3 kcal mol−1, Table 2), decreasing the bond strength of the imine bond in LASSBio-294 (1) in comparison to LASSBio-897 (2).
Fig. 5 Data from IR (red) and Raman (blue) spectroscopy in the region of 1700–1450/880–600 cm−1. (A) Data obtained for LASSBio-294 (1). (B) Data obtained for LASSBio-897 (2). |
Taking into account the theoretical and experimental data presented in this work, it is possible to confirm that there is an N⋯S intramolecular interaction in the N-acylhydrazone derivative LASSBio-294 (1), which causes a conformational restriction and has relevant consequences in its molecular properties when compared to its regioisomer, LASSBio-897 (2). This observation has implications in the medicinal chemistry community, since this type of interaction is only beginning to be understood and there are some other interesting examples in which N-acylhydrazone derivatives substituted with a 2-thiophene ring at the imine side were used and presented several pharmacological activities, such as antiviral,47 antibacterial,48,49 analgesic and anti-inflammatory.50 Beyond that, the study of the molecular properties of N-acylhydrazone derivatives associated with the N⋯S intramolecular interaction can help in the understanding of the photoswitch properties of this privileged structure.51 To our knowledge, this is the first report of a N⋯S interaction involving the privileged structure N-acylhydrazone and considering that this molecular fragment has been extensively explored in the pharmacological modulation of several biological targets,28,29,52 we hope that this work will open new doors for the rational use of this privileged structure, such as exploring conformationally restricted compounds through a N⋯S σ-hole interaction.3,13–15
To explore the analysis performed here, we applied the same studies in two similar systems with N-methylation of the amide nitrogen of LASSBio-294 (1) and LASSBio-897 (2), i.e., LASSBio-785 (4)40,53,54 and LASSBio-1289 (5),55–57 respectively (Fig. 6). The N-methylation of these compounds causes significant conformational changes,40,58 but did not affect the occurrence of the 1,4-N⋯S σ-hole intramolecular interaction in LASSBio-785 (4), as significant differences were also observed in the bond strength of the S–C and NC bonds, including a bathochromic shift of the λMAX in the UV spectra and the 15N chemical shifts in relation to LASSBio-1289 (5) (see ESI,† Section 3, Fig. S9–S13 and Tables S12–S14). These analyses were applied in two other systems, a positive control of these interactions already established in the literature 2-(thiophen-2-yl)pyridine (3)37,41 (used in our theoretical validation) and its respective negative control 2-(thiophen-3-yl)pyridine (6) (Fig. 6), which validated the results for the N-acylhydrazone derivatives (see ESI,† Section 2, Fig. S4–S8 and Tables S9–S11).
To our knowledge, this is the first time that σ-hole interactions have been investigated simultaneously using different spectroscopic techniques and theoretical evaluations and a negative control, which strengthens our results. This is also the first description of N⋯S σ-hole intramolecular interactions in N-acylhydrazone derivatives, which is very important because this versatile privileged structure can be used in the modulation of several biological targets. We hope that our results will assist with further characterizations of this type of interaction.
Footnote |
† Electronic supplementary information (ESI) available: Validation process of the theoretical methodology for the description of the 1,4-N⋯S σ-hole intramolecular interaction. Theoretical and experimental characterization of the 1,4-N⋯S σ-hole intramolecular interaction for the compounds 2-(thiophen-2-yl)pyridine (3), 2-(thiophen-3-yl)pyridine (6), LASSBio-785 (4) and LASSBio-1289 (5). See DOI: 10.1039/c7nj03543h |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2018 |