Ignez Caracelli1, Julio Zukerman-Schpector2, Henrique J Traesel3, Paulo R Olivato3, Mukesh M Jotani4, Edward R T Tiekink5. 1. Departamento de Física, Universidade Federal de São Carlos, 13565-905 São Carlos, SP, Brazil. 2. Departamento de Química, Universidade Federal de São Carlos, 13565-905 São Carlos, SP, Brazil. 3. Instituto de Química, Universidade de São Paulo, 05508-000 São Paulo, SP, Brazil. 4. Department of Physics, Bhavan's Sheth R. A. College of Science, Ahmedabad, Gujarat 380001, India. 5. Research Centre for Crystalline Materials, School of Science and Technology, Sunway University, 47500 Bandar Sunway, Selangor Darul Ehsan, Malaysia.
Abstract
The title compound, C15H13BrO2S, comprises three different substituents bound to a central (and chiral) methine-C atom, i.e. (4-bromo-phen-yl)sulfanyl, benzaldehyde and meth-oxy residues: crystal symmetry generates a racemic mixture. A twist in the mol-ecule is evident about the methine-C-C(carbon-yl) bond as evidenced by the O-C-C-O torsion angle of -20.8 (7)°. The dihedral angle between the bromo-benzene and phenyl rings is 43.2 (2)°, with the former disposed to lie over the oxygen atoms. The most prominent feature of the packing is the formation of helical supra-molecular chains as a result of methyl- and methine-C-H⋯O(carbon-yl) inter-actions. The chains assemble into a three-dimensional architecture without directional inter-actions between them. The nature of the weak points of contacts has been probed by a combination of Hirshfeld surface analysis, non-covalent inter-action plots and inter-action energy calculations. These point to the importance of weaker H⋯H and C-H⋯C inter-actions in the consolidation of the structure.
The title compound, C15H13BrO2S, comprises three different substituents bound to a central (and chiral) methine-C atom, i.e. (4-bromo-phen-yl)sulfanyl, benzaldehyde and meth-oxy residues: crystal symmetry generates a racemic mixture. A twist in the mol-ecule is evident about the methine-C-C(carbon-yl) bond as evidenced by the O-C-C-O torsion angle of -20.8 (7)°. The dihedral angle between the bromo-benzene and phenyl rings is 43.2 (2)°, with the former disposed to lie over the oxygen atoms. The most prominent feature of the packing is the formation of helical supra-molecular chains as a result of methyl- and methine-C-H⋯O(carbon-yl) inter-actions. The chains assemble into a three-dimensional architecture without directional inter-actions between them. The nature of the weak points of contacts has been probed by a combination of Hirshfeld surface analysis, non-covalent inter-action plots and inter-action energy calculations. These point to the importance of weaker H⋯H and C-H⋯C inter-actions in the consolidation of the structure.
Recently, the crystal structure determination of the chloro analogue of the title compound was described (Caracelli et al., 2018 ▸). This was evaluated as a part of on-going studies into the conformational and electronic characteristics of various β-thiocarbonyl, β-bis-thiocarbonyl and β-thio-β-oxacarbonyl compounds, and their selenium counterparts, employing infrared spectroscopy, computational chemistry and X-ray crystallographic methods (Vinhato et al., 2013 ▸; Zukerman-Schpector et al., 2015 ▸; Caracelli et al., 2015 ▸; Traesel et al., 2018 ▸). In particular, the evaluation of the anti-inflammatory activity of what could be selective COX-2 inhibitors (Cerqueira et al., 2017 ▸) motivates these investigations, which are supported by molecular docking studies designed to ascertain the mechanism(s) of inhibition (Baptistini, 2015 ▸). Subsequently, crystals of the title bromo analogue (I) were obtained: the crystal structure is reported herein along with an analysis of the calculated Hirshfeld surfaces, non-covalent interaction plots (for selected interactions) as well as a computational chemistry study.
Structural commentary
The molecular structure of (I), Fig. 1 ▸, is isostructural with the previously described chloro analogue, (II) (Caracelli et al., 2018 ▸). Here, the central chiral methine-C8 atom is connected to (4-bromophenyl)sulfanyl, phenylethanone and methoxy groups. There is a twist in the ethanone residue as seen in the value of the O1—C8—C9—O2 torsion angle of −20.8 (7)°, with the oxygen atoms being approximately syn. The dihedral angle between the bromobenzene and phenyl rings is 43.2 (2)°, indicative of an inclined relative disposition. Globally, the bromobenzene ring is orientated towards the ethanone residue.
Figure 1
The molecular structure of (I), showing the atom-labelling scheme and displacement ellipsoids at the 25% probability level.
The geometric parameters in (I) can be compared with those of (II): the twist about the central C8—C9 bond is approximately the same in (II), i.e. the the O1—C8—C9—O2 torsion angle is 19.3 (7)°, as is the dihedral angle of 42.9 (2)° between the aromatic rings. The overlay diagram in Fig. 2 ▸ highlights the close similarity between the molecular structures of (I) and (II).
Figure 2
Overlay diagram of (I) (red image) and (II) (blue image).
Supramolecular features
The main feature of the molecular packing of (I) is the presence of C—H⋯O interactions where the carbonyl-O2 atom accepts two contacts from methyl-C7-H and methine-C8-H atoms derived from the same molecule to generate six-membered {⋯O⋯HCOCH} synthons, Table 1 ▸. The result is a supramolecular chain propagating along [001] with an helical topology (21 symmetry), Fig. 3 ▸(a). The chains pack without directional interactions between them, Fig. 3 ▸(b).
Table 1
Hydrogen-bond geometry (Å, °)
D—H⋯A
D—H
H⋯A
D⋯A
D—H⋯A
C7—H7A⋯O2i
0.96
2.47
3.296 (9)
144
C8—H8⋯O2i
0.98
2.44
3.331 (6)
150
Symmetry code: (i) .
Figure 3
Molecular packing in (I): (a) view of the helical supramolecular chain parallel to the c axis sustained by C—H⋯O interactions shown as orange dashed lines and (b) view of the unit-cell contents shown in projection down the c axis; one chain is highlighted in space-filling mode.
Hirshfeld surface analysis
The Hirshfeld surface calculations for (I) were performed in accord with protocols described recently (Tan et al., 2019 ▸) employing Crystal Explorer (Turner et al., 2017 ▸). Over and above the analysis of the important surface contacts in the crystal of (I), the results are compared with those for the recently determined isostructural chloro analogue (II) (Caracelli et al., 2018 ▸). The crystal of (I) has similar intermolecular C—H⋯O interactions (Table 1 ▸) and short interatomic H⋯H, C⋯H and C⋯C contacts (Table 2 ▸) as in isostructural (II), as detailed below.
Table 2
Summary of short interatomic contacts (Å) in (I) and (II)
Contact
Distance
Symmetry operation
(I)
H7B⋯H14
2.15
1 − x, − y, − + z
H7C⋯C6
2.74
1 − x, 2 − y, + z
H12⋯Br1
3.02
− x, −1 + y, + z
C6⋯C9
3.355 (8)
1 − x, 1 − y, − + z
(II)
H7B⋯H14
2.10
1 − x, − y, + z
H7B⋯C14
2.76
1 − x, − y, + z
H7C⋯C6
2.73
1 − x, 1 − y, + z
C6⋯C9
3.334 (9)
1 − x, − y, + z
Notes: (a) The interatomic distances are calculated in Crystal Explorer (Turner et al., 2017 ▸) whereby the X—H bond lengths are adjusted to their neutron values.
The intermolecular contacts in (I), Tables 1 ▸ and 2 ▸, are characterized as the pair of bright-red spots near the carbonyl-O2 atom, and each of the methyl-H7A and methine-H8 atoms on the Hirshfeld surfaces mapped over d
norm in the images of Fig. 4 ▸. Further, interactions are indicated by the faint-red spots near the methyl-H7B and H7C, phenyl-H14, bromobenzene-C6 and carbonyl-C9 atoms in Fig. 4 ▸. On the Hirshfeld surfaces mapped over the calculated electrostatic potential in the images of Fig. 5 ▸, the donors and acceptors of intermolecular interactions are viewed as blue and red regions around the participating atoms corresponding to positive and negative potentials, respectively. The environment around a reference molecule within the d
norm-mapped Hirshfeld surface highlighting the intermolecular C—H⋯O interactions and short interatomic H⋯H, C⋯H/H⋯C and C⋯C contacts is illustrated in Fig. 6 ▸.
Figure 4
Two views of the Hirshfeld surface for (I) mapped over d
norm in the range −0.084 to +1.422 arbitrary units.
Figure 5
Two views of the Hirshfeld surface for (I) mapped over the electrostatic potential in the range −0.074 to + 0.053 atomic units. The red and blue regions represent negative and positive electrostatic potentials, respectively.
Figure 6
A view of the Hirshfeld surface for (I) mapped over d
norm in the range −0.084 to +1.422 arbitrary units highlighting intermolecular C—H⋯O, C⋯C, H⋯H and C⋯H/H⋯C contacts by black, red, yellow and sky-blue dashed lines, respectively.
From the overall two-dimensional fingerprint plot in Fig. 7 ▸(a), and also those delineated into H⋯H, O⋯H/H⋯O, C⋯H/H⋯C, C⋯C and Br⋯H/H⋯Br contacts in Fig. 7 ▸(b)–(f), respectively, it is evident that the plots are basically identical in shape to those calculated for the chloro analogue (II) with only slight differences in the distribution of points (Caracelli et al., 2018 ▸). The percentage contributions from the different interatomic contacts to the Hirshfeld surfaces of (I) and (II) are summarized in Table 3 ▸; these values again highlight the similarities between (I) and (II).
Figure 7
(a) The full two-dimensional fingerprint plot for (I) and (b)–(f) those delineated into H⋯H, O⋯H/H⋯O, C⋯H/H⋯C, C⋯C and Br⋯H/H⋯Br contacts.
Table 3
Percentage contributions of interatomic contacts to the Hirshfeld surface for (I) and (II)
Contact
Percentage contribution
(I), X = Br
(II), X = Cl
H⋯H
39.3
39.1
O⋯H/H⋯O
11.0
10.7
C⋯H/H⋯C
23.2
23.0
X⋯H/H⋯X
12.8
13.3
S⋯H/H⋯S
4.4
4.3
X⋯S/S⋯X
2.1
2.3
X⋯O/O⋯X
2.1
2.1
C⋯O/O⋯C
1.5
1.5
C⋯X/X⋯C
1.5
1.8
C⋯S/S⋯C
1.2
1.1
C⋯C
0.6
0.6
The C—H⋯O contacts significant in the crystal of (I), Table 1 ▸, are represented as the pair of spikes at d
e + d
i ∼2.3 Å in the fingerprint plot delineated into O⋯H/H⋯O contacts, Fig. 7 ▸(c). The short interatomic H⋯H, C⋯H/H⋯C and C⋯C contacts (Table 2 ▸) are characterized as pair of beak-shape tips at d
e + d
i ∼2.1 Å, Fig. 7 ▸(b), and forceps at d
e + d
i ∼2.8 Å, Fig. 7 ▸(d), and vase-shaped distribution of points at d
e + d
i ∼3.3 Å, Fig. 7 ▸(e), in the respective delineated fingerprint plots. In addition to these contacts, the crystal also features short interatomic Br⋯H/H⋯Br contacts appearing as the pair of forceps-like tips at d
e + d
i ∼3.0 Å in the delineated fingerprint plot of Fig. 7 ▸(f). The small contribution from other remaining interatomic contacts summarized in Table 3 ▸ have a negligible effect on the packing.
Interaction energies
The pairwise interaction energies between the molecules within the crystal are calculated by the summation of four energy components comprising electrostatic (E
ele), polarization (E
pol), dispersion (E
dis) and exchange-repulsion (E
rep) (Turner et al., 2017 ▸). These energies were obtained by using the wave function calculated at the HF/STO-3G level theory for each of (I) and (II). The individual energy components as well as total interaction energy relative to reference molecule within the molecular cluster were calculated. Table 4 ▸ summarizes quantitatively the strength and nature of intermolecular interactions in the crystals of (I) and (II).
Table 4
Summary of interaction energies (kJ mol−1) calculated for (I) and (II)
Contact
R (Å)
Eele
Epol
Edis
Erep
Etot
(I)
C7—H7A⋯O2i +
C8—H8⋯O2i +
H7B⋯H14i +
C6⋯C9i
6.40
−20.0
−12.1
−53.2
34.0
−48.0
H7C⋯C6ii
8.75
−7.0
−1.2
−16.7
9.3
−15.4
H12⋯Br1ii
10.83
−4.1
−0.9
−12.9
6.4
−11.2
(II)
C7—H7A⋯O2iii +
C8—H8⋯O2iii +
H7B⋯H14iii +
C6⋯C9iii +
H7B⋯C14iii
6.13
−19.5
−11.8
−52.7
35.1
−46.6
H7C⋯C6iv
9.06
−6.6
−1.4
−14.5
8.2
−14.0
Notes: Symmetry operations: (i) 1 − x, 1 − y, − + z; (ii) 1 − x, 2 − y, + z; (iii) 1 − x, − y, + z; (iv) 1 − x, 1 − y, + z.
It is observed from the interaction energies calculated between the reference molecule and the symmetry-related molecules at R = 6.40 and 6.13 Å (where R is the separation of the centres of gravity of the molecules), respectively (Table 4 ▸), that the almost identical values of the electrostatic energy component are due to intermolecular C—H⋯O interactions whereas the dispersive components are dominant owing to the short interatomic contacts between the same molecules. The other short interatomic C⋯H/H⋯C contact between the methyl-H7C and phenyl-C6 atoms in (I) and (II), and the H12⋯Br1 contact in (I) have a major contribution from dispersion components.The magnitudes of intermolecular energies are represented graphically in the energy frameworks for (I) and (II) viewed down the c axes are shown in Fig. 8 ▸. Here, the supramolecular architecture of the crystals is represented as cylinders joining centroids of molecular pairs. The red, green and blue coloration represent the energy components E
ele, E
disp and E
tot, respectively. The radius of the cylinder is proportional to the magnitude of interaction energy which are adjusted to the same scale factor (3 kJ mol−1) within 4 × 4 × 4 unit cells. From the energy frameworks for (I) and (II) illustrated in Fig. 8 ▸, it is clearly evident that the supramolecular associations viewed down the c axis are identical, reflecting the isostructural relationship between (I) and (II).
Figure 8
A comparison of the energy frameworks, plotted with the same scale, composed of (a) electrostatic potential force, (b) dispersion force and (c) total energy for the molecules of (I), upper images, and (II), lower images, all viewed down the c-axis direction. same scale factor of 50 with a cut-off value of 3 kJ mol−1 within 4 x 4 x 4 unit cells.
Non-covalent interaction plots
The non-covalent interaction plot (NCIplot) analysis was used in the present study in order to confirm the attractive nature of some of the specified intermolecular contacts (Contreras-García et al., 2011 ▸). This method is based on the electron density and its derivatives allowing the visualization of the gradient isosurfaces. The colour-based isosurfaces correspond to the values of sign(λ2)ρ(r), where ρ is the electron density and λ2 is the second eigenvalue of the Hessian matrix of ρ (Johnson et al., 2010 ▸). The isosurfaces for the interactions between the carbonyl-O2 and each of the methyl-H7B and phenyl-H14 atoms, the H7B and H14 atoms, and the chlorobenzene-C6 and methyl-H7C atoms are shown in the upper views of Fig. 9 ▸(a)–(c), respectively. The green isosurface observed in each of these indicates a weakly attractive interaction as opposed to attractive (blue isosurface) or repulsive (red). The lower views of Fig. 9 ▸, where the plots of the RDG versus sign(λ2)ρ(r) are depicted, the non-covalent interaction peaks appear at density values equal or lower than 0.01 a.u., consistent with weakly attractive interactions.
Figure 9
Non-covalent interaction plots for intermolecular interactions between (a) each of the methyl-C7- and methine-C—H atoms and the carbonyl-O2 atom, (b) the methyl-H7B and phenyl-H14 atoms and (c) bromobenzene-C6 and methyl-H7C atoms.
Database survey
There are three literature structures related to (I), namely the already mentioned (II) (NIBTAW; Caracelli et al., 2018 ▸), the S-bound 4-methoxybenzene derivative [(III); JUPLOZ; Caracelli et al., 2015 ▸] and the S-bound 4-tolyl species [NOVGIQ; (IV); Zukerman-Schpector et al., 2015 ▸] derivatives. All four compounds crystallize in the orthorhombic space group Pca21 and are isostructural. The differences between the molecules of (I)–(IV) relates to the relative orientations of the S-bound methoxybenzene ring in (III). This comes about owing to a twist about the C8—S1 bond as manifested in the C4—S1—C8—C9 torsion angles of 57.1 (4), 57.3 (5), 46.6 (3) and 57.9 (3)° for (I)–(IV), respectively. This difference notwithstanding, the angles between the S-bound benzene rings and the phenyl rings in (I)–(IV) span a relatively narrow range of values, i.e. 43.2 (2), 42.9 (2), 40.11 (16) and 44.03 (16)°, respectively.
Synthesis and crystallization
Firstly, 4′-bromothiophenol (10.0 g, 52.9 mmol) was reacted with bromine (3.1 ml, 56.0 mmol) in dichloromethane (400 ml) on a hydrated silica gel support (50 g of SiO2 and water (30 ml) to give 4′-bromophenyl disulfide (8.0 g, yield 80%). A brown solid was obtained after filtration and evaporation without further purification (Ali & McDermott, 2002 ▸). Then, a solution of 2-methoxy acetophenone (Sigma–Aldrich; 1.0 ml, 7.3 mmol) in THF (25 ml), was added dropwise to a cooled (195 K) solution of diisopropylamine (1.1 ml, 8.0 mmol) and n-butyllithium (5.4 ml, 7.3 mmol) in THF (30 ml). After 30 mins, a solution of 4′-bromophenyl disulfide (2.8 g, 7.3 mmol) with hexamethylphosphoramide (HMPA) (1.3 ml, ca 7.3 mmol) dissolved in THF (35 ml) was added dropwise to the enolate solution (Zoretic & Soja, 1976 ▸). After stirring for 3 h, water (70 ml) was added at room temperature and extraction with diethyl ether ensued. The organic layer was then treated with a saturated solution of ammonium chloride until neutral pH was reached and then dried over anhydrous magnesium sulfate. A brown oil was obtained after evaporation of the solvent. Purification through flash chromatography with n-hexane was used in order to remove the non-polar reactant (disulfide), then with dry acetone to give a mixture of both acetophenones (product and reactant). Crystallization was performed by vapour diffusion of n-hexane into a chloroform solution held at 283 K to give the pure product (0.6 g, yield = 70%). Irregular colourless crystals suitable for X-ray diffraction of (I) were obtained by the same pathway. M.p. 357.0–357.5 K. 1H NMR (CDCl3, 500 MHz, δ ppm): 3.67 (s, 3H), 5.87 (s, 1H), 7.20–7.23 (m, 2H), 7.39–7.41 (m, 2H), 7.44–7.47 (m, 2H), 7.57–7.62 (m, 1H), 7.92–7.94 (m, 2H). 13C NMR (CDCl3, 125 MHz, δ p.p.m.): 190.16, 135.73, 134.18, 133.53, 132.13, 129.92, 128.81, 128.57, 123.41, 89.28, 56.10. Microanalysis calculated for C15H13BrO2S (%): C 53.42, H 3.89. Found (%): C 53.19, H 3.85. High-Resolution MS [M
+, M
+] calculated: 335.9820, 337.9799; found: 335.9797, 337.9778.
Refinement details
Crystal data, data collection and structure refinement details are summarized in Table 5 ▸. The carbon-bound H atoms were placed in calculated positions (C—H = 0.93–0.98 Å) and were included in the refinement in the riding-model approximation, with U
iso(H) set to 1.2–1.5U
eq(C). The absolute structure was determined based on differences in Friedel pairs included in the data set (Parsons et al., 2013 ▸).
Table 5
Experimental details
Crystal data
Chemical formula
C15H13BrO2S
Mr
337.21
Crystal system, space group
Orthorhombic, Pca21
Temperature (K)
293
a, b, c (Å)
18.0683 (13), 8.0190 (6), 9.8513 (5)
V (Å3)
1427.35 (16)
Z
4
Radiation type
Mo Kα
μ (mm−1)
3.02
Crystal size (mm)
0.47 × 0.20 × 0.14
Data collection
Diffractometer
Bruker APEXII CCD
Absorption correction
Multi-scan (SADABS; Sheldrick, 1996 ▸)
Tmin, Tmax
0.545, 0.745
No. of measured, independent and observed [I > 2σ(I)] reflections
6329, 2820, 1903
Rint
0.037
(sin θ/λ)max (Å−1)
0.625
Refinement
R[F2 > 2σ(F2)], wR(F2), S
0.035, 0.086, 0.90
No. of reflections
2820
No. of parameters
173
No. of restraints
1
H-atom treatment
H-atom parameters constrained
Δρmax, Δρmin (e Å−3)
0.24, −0.36
Absolute structure
Flack x determined using 702 quotients [(I+)−(I−)]/[(I+)+(I−)] (Parsons et al., 2013 ▸)
Absolute structure parameter
0.013 (11)
Computer programs: APEX2 and SAINT (Bruker, 2009 ▸), SIR2014 (Burla et al., 2015 ▸), SHELXL2014 (Sheldrick, 2015 ▸), ORTEP-3 for Windows (Farrugia, 2012 ▸), DIAMOND (Brandenburg, 2006 ▸), MarvinSketch (ChemAxon, 2010 ▸) and publCIF (Westrip, 2010 ▸).
Crystal structure: contains datablock(s) I, global. DOI: 10.1107/S2056989019006765/hb7825sup1.cifStructure factors: contains datablock(s) I. DOI: 10.1107/S2056989019006765/hb7825Isup2.hklClick here for additional data file.Supporting information file. DOI: 10.1107/S2056989019006765/hb7825Isup3.cmlCCDC reference: 1915470Additional supporting information: crystallographic information; 3D view; checkCIF report
Secondary atom site location: difference Fourier map
Least-squares matrix: full
Hydrogen site location: inferred from neighbouring sites
R[F2 > 2σ(F2)] = 0.035
H-atom parameters constrained
wR(F2) = 0.086
w = 1/[σ2(Fo2)] where P = (Fo2 + 2Fc2)/3
S = 0.90
(Δ/σ)max < 0.001
2820 reflections
Δρmax = 0.24 e Å−3
173 parameters
Δρmin = −0.36 e Å−3
1 restraint
Absolute structure: Flack x determined using 702 quotients
[(I+)-(I-)]/[(I+)+(I-)]
(Parsons et al., 2013)
Primary atom site location: structure-invariant direct methods
Absolute structure parameter: 0.013 (11)
Geometry. All esds (except the esd in the dihedral angle between two l.s. planes)
are estimated using the full covariance matrix. The cell esds are taken
into account individually in the estimation of esds in distances, angles
and torsion angles; correlations between esds in cell parameters are only
used when they are defined by crystal symmetry. An approximate (isotropic)
treatment of cell esds is used for estimating esds involving l.s. planes.
Authors: Julia Contreras-García; Erin R Johnson; Shahar Keinan; Robin Chaudret; Jean-Philip Piquemal; David N Beratan; Weitao Yang Journal: J Chem Theory Comput Date: 2011-03-08 Impact factor: 6.006
Authors: Elisângela Vinhato; Paulo R Olivato; Julio Zukerman-Schpector; Maurizio Dal Colle Journal: Spectrochim Acta A Mol Biomol Spectrosc Date: 2013-07-08 Impact factor: 4.098
Authors: Erin R Johnson; Shahar Keinan; Paula Mori-Sánchez; Julia Contreras-García; Aron J Cohen; Weitao Yang Journal: J Am Chem Soc Date: 2010-05-12 Impact factor: 15.419