This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.
Article http://pubs.acs.org/journal/acsodf
Inductive or Field Substituent Effect? Quantum Chemical Modeling of Interactions in 1‑Monosubstituted Bicyclooctane Derivatives Halina Szatylowicz,*,† Tomasz Siodła,‡ and Tadeusz M. Krygowski*,§ †
Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89b, 61-614 Poznań, Poland § Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland ‡
S Supporting Information *
ABSTRACT: Inductive substituent constants were obtained for systems lacking the resonance effect. The application of the charge of the substituent active region concept to study the substituent effect in 1-X-substituted bicyclooctane derivatives (B3LYP/6-311++G** calculations, X = NMe2, NH2, OH, OMe, CH3, H, F, Cl, CF3, CN, CHO, COMe, CONH2, COOH, NO2, NO) has revealed inductive interactions, which are through bonds.
■
space or through bonds? The first would be a field substituent effect (SE), whereas the other, an inductive SE. In the first case, the effect should be described by electrostatic rules.7,8 Some attempts were carried out by the use of the isolated molecule approach.9,10 Exner and Bohm1,11 performed a very detailed study using both these strategies and arrived at the conclusion that transmission of the SE through space, at least of noncharged functional groups, is of “little validity”. The alternative mechanism of interaction, that is, the transmission of the effect through bonds, “cannot estimate the SE explicitly; it describes only attenuation within the molecule”.1 This means that each bond weakens the effect by a constant ratio. In summary, according to the results of even the most detailed analysis of various theories of SEs, it is impossible to validate these approaches on the basis of experimental facts.1 Recently, several quantum chemical models have been developed that are of substantial importance for the discussion of SEs in terms of physically defined concepts. It should also be noted that particular physical characteristics are compared or correlated with substituent constants. One of the frequently applied methodologies made use of various defined electrostatic potentials, as follows: on particular atoms, at atoms of reaction sites, or in other defined sites of molecules.12−16 Additionally and importantly, application of the molecular electrostatic potential topography allowed us to appraise the through-bond
INTRODUCTION 1,4-Disubstituted bicyclo[2.2.2]octane (BCO) derivatives are key molecular systems for the estimation of so-called inductive/ field substituent constants.1 The concept of inductive substituent constants (σind) was introduced by Roberts and Moreland.2 For this purpose, ionization constants measured in a 50% (by volume) ethanol−water solution for 4-substituted BCO carboxylic acids (Scheme 1) were applied. Inductive Scheme 1. 4-X-Bicyclo[2.2.2]octane-1-carboxylic Acids (4-XBCO-COOH); X = H, OH, CO2C2H5, Br, CN
substituent constants obtained for these systems characterize mutual interactions between two relatively distant functional groups acting without the possibility of a resonance effect. Later, for this purpose, 4-substituted quinuclidine3,4 and gasphase acidities of 4-substituted BCO carboxylic acids5 were used. Many other attempts for the estimation of substituent constants (denoted as: σind, σI, σF, or F) have appeared, and in most cases, they have led to similar numerical results. This is well illustrated in Table 2 in the review by Hansch et al.6 Along with the development of research methods and accumulation of new results, a question concerning the mechanisms of interactions between the functional groups in 1,4-disubstituted BCO derivatives has arisen. Is the effect transmitted through © 2017 American Chemical Society
Received: April 1, 2017 Accepted: April 21, 2017 Published: May 1, 2017 1746
DOI: 10.1021/acsomega.7b00390 ACS Omega 2017, 2, 1746−1749
ACS Omega
Article
and through-space SEs.17 Implementation of the energy decomposition analysis18 allowed us to document that the pielectron energy of SE can be correlated with the substituent constants.19 Another very important issue is the energetic characteristic termed SE stabilization energy (SESE), obtained by the use of the isodesmic or homodesmotic reaction approach.20−22 In many cases, it was demonstrated that the SESE values correlated well with the substituent constants.23−25 Quite recently, it has been shown that the sum of charges at the substituent and ipso carbon atom, termed cSAR (acronym derived from the charge of the substituent active region),23−27 is also well correlated with the Hammett substituent constants. Moreover, it should be noted that usually the values of cSAR(X), calculated by means of various atomic charge assessments, are mutually well correlated.28 The motivation for our report is to provide a new perspective on the inductive/field effect in monosubstituted derivatives of bicycloctane (1-X-BCO) derivatives. For this purpose, the cSAR approach will be applied and the cSAR(X) values will be confronted with the substituent constants.
expected from the theory of the inductive effect.1 It appears that this result may be accepted as a strong argument for the inductive mechanism of the SE in BCO systems. Additionally, it should be mentioned that the through-space interactions seem to be negligible in our case because the CH2 and CH groups have a very small dipole moment and hence very low electrostatic energy. Actually, both kinds of interactions, via bonds and through space, are described by the substituent constants, F.6 Figure 2
■
RESULTS AND DISCUSSION To study the nature of the SE in 1-X-BCO derivatives, the cSAR parameter has been used to characterize all fragments of the studied systems (Scheme 2). Therefore, cSAR(X),
Figure 2. Dependence of cSAR(X) on F constants, taken from ref 6, for 1-X-BCO derivatives (red, σp < 0; blue, σp > 0); for all data, y = −0.218x + 0.024 and R2 = 0.825.
Scheme 2. 1-X-Bicyclo[2.2.2]octane Derivatives (1-X-BCO), X = NMe2, NH2, OH, OMe, CH3, H, F, Cl, CF3, CN, CHO, COMe, CONH2, COOH, NO2, NO
shows the relationship between cSAR(X) and F values for two groups of substituents: (i) with electron-donating properties (σp < 0, X = Me, OMe, OH, NH2, NMe2) and (ii) with electron-accepting properties (positive σp values, X = NO, NO2, CN, CF3, COMe, COOH, CHO, CONH2, Cl, F). In both cases, cSAR(X) values decrease with an increase in F constants. However, the electron-attracting ability of the substituents of the first group is significantly weaker than that of the substituents of the second group. The next, new result is associated with the application of cSAR(X) as a measure of the electron-donating/attracting properties of X. It has already been shown that the values of cSAR(X) depend on the kind of system to which X is attached,23,24 and this kind of SE is termed reverse SE.28 The obtained cSAR(X) values for 1-X-BCO and monosubstituted benzene (X-Ph) derivatives are listed in Table 1, together with the differences (Δ) between the values estimated for the X-Ph and 1-X-BCO systems, whereas Figure 3 illustrates these differences by means of linear regression. Its slope indicates that the SE in phenyl derivatives is ∼1.7 times stronger than that in the BCO derivatives. This is due to the difference in the nature of interactions. In the BCO series, only the inductive effect works, whereas in the phenyl series, both inductive and resonance effects are in use. Additionally, according to the Taft approach,29 the separation of the resonance and inductive effects of a substituent, Δ, may be considered as a measure of the resonance effect. Indeed, when Δ is plotted against the resonance substituent constants, R, the image is as shown in Figure 4. The high determination coefficient supports this assumption.
cSAR(CH2) at the 2 and 3 positions, and cSAR(CH) at the 4 position have been obtained; in the case of cSAR(CH2), the mean values of all three CH2 groups at the C2 and C3 positions have also been taken into account. Figure 1 presents the results
Figure 1. Dependence of cSAR(CH2) in positions 2−4 on cSAR(X) for 1-X-BCO derivatives.
■
and regression lines, wherein cSAR(CH2) in positions 2 and 3 and cSAR(CH) in position 4 of 1-X-BCO are plotted against cSAR(X). In all three cases, the regression lines have high determination coefficients, (R2 > 0.906); hence, the slopes of these lines are reliable data. As observed, their ratio was 0.19:0.12:0.06, which is near the ideal 3:2:1 ratio, as could be
CONCLUSIONS The estimated values of cSAR(CH2) and cSAR(CH) for groups in the 2, 3, and 4 positions of 1-X-BCO, influenced by substituent X, are attenuated in a regular ratio, which is 3:2:1. 1747
DOI: 10.1021/acsomega.7b00390 ACS Omega 2017, 2, 1746−1749
ACS Omega
Article
The SE in monosubstituted benzene derivatives is ∼1.7 times stronger than that in the 1-X-substituted BCO systems, that is, the aliphatic ones. Furthermore, the differences in cSAR(X) values for substituents in phenyl and BCO derivatives may be used as measures of the resonance effect of substituent X.
Table 1. cSAR(X) Values for Monosubstituted Benzene and BCO Derivatives cSAR(X) NO NO2 CN CF3 COMe COOH CHO CONH2 Cl F H Me OMe OH NH2 NMe2 range a
σpa
Fa
Ra
1-X-BCO
X-Ph
Δ
0.91 0.78 0.66 0.54 0.50 0.45 0.42 0.36 0.23 0.06 0.00 −0.17 −0.27 −0.37 −0.66 −0.83
0.49 0.65 0.51 0.38 0.33 0.34 0.33 0.26 0.42 0.45 0.03 0.01 0.29 0.33 0.08 0.15
0.42 0.13 0.15 0.16 0.17 0.11 0.09 0.10 −0.19 −0.39 0.00 −0.18 −0.56 −0.70 −0.74 −0.98
−0.074 −0.109 −0.131 −0.065 −0.041 −0.054 −0.066 −0.035 −0.078 −0.066 0.007 0.007 −0.015 −0.015 0.010 0.004 0.141
−0.132 −0.140 −0.139 −0.091 −0.069 −0.089 −0.101 −0.055 −0.050 −0.028 0.000 0.030 0.057 0.044 0.091 0.119 0.259
−0.058 −0.031 −0.008 −0.026 −0.028 −0.035 −0.035 −0.019 0.028 0.038 −0.007 0.024 0.072 0.059 0.081 0.115
■
COMPUTATIONAL METHODS An optimization, without any symmetry constraints, of all studied systems was performed using the B3LYP hybrid functional30,31 with the 6-311++G(d,p) basis set.32 The vibrational frequencies were calculated at the same level of theory to confirm that all calculated structures correspond to the minima on the potential energy surface. The cSAR parameter26,33 was calculated by summing the atomic charges of all atoms of group X and the charge at the ipso carbon atom to which X is attached cSAR(X) = q(X) + q(Cipso)
For assessments of the atomic charges, the Hirshfeld34 method was applied. All calculations were performed using the Gaussian09 program.35
■
Taken from ref 6.
ASSOCIATED CONTENT
S Supporting Information *
The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsomega.7b00390. Cartesian coordinates of equilibrium geometries of 1-XBCO and monosubstituted benzene derivatives (PDF)
■
AUTHOR INFORMATION
Corresponding Authors
*E-mail:
[email protected] (H.S.). *E-mail:
[email protected] (T.M.K.). ORCID
Halina Szatylowicz: 0000-0002-7034-6985
Figure 3. Dependence of cSAR(X) in X-Ph on cSAR(X) in 1-X-BCO derivatives.
Notes
The authors declare no competing financial interest.
■
ACKNOWLEDGMENTS We gratefully acknowledge the Interdisciplinary Center for Mathematical and Computational Modeling (Warsaw, Poland) for providing computer time and facilities. H.S. and T.M.K. thank the National Science Centre and Ministry of Science and Higher Education of Poland for supporting this work under grant no. UMO-2013/11/B/ST4/00531.
■ ■
DEDICATION Dedicated to our friend, professor Guenter Haefelinger, on the occasion of his 80th birthday. Figure 4. Dependence of Δ on resonance substituent constant R.
REFERENCES
(1) Exner, O.; Bohm, S. Theory of substituent effects: Recent advances. Curr. Org. Chem. 2006, 10, 763−778. (2) Roberts, J. D.; Moreland, W. T. Electrical effects of substituent groups in saturated systems. Reactivities of 4-substituted bicyclo[2.2.2]octane-1-carboxylic acids. J. Am. Chem. Soc. 1953, 75, 2167−2173. (3) Grob, C. A.; Schlageter, M. G. Polar Substituent Effects on Rates and Equilibria of 4-Substituted Quinuclidines. Preliminary communication. Helv. Chim. Acta 1974, 57, 509−511.
Hence, it seems justified to accept that F constants mostly represent the inductive effect. Electron-donating substituents, by means of Hammett σp < 0, are sometimes electron attracting, when determined in monosubstituted BCO derivatives. 1748
DOI: 10.1021/acsomega.7b00390 ACS Omega 2017, 2, 1746−1749
ACS Omega
Article
(4) Paleček, J.; Hlavatý, J. Azabicyclo compounds. XII. 4-Substituted quinuclidines as model compounds for the study of polar effects transfer. Collect. Czech. Chem. Commun. 1973, 38, 1985−2002. (5) Koppel, I. A.; Mishima, M.; Stock, L. M.; Taft, R. W.; Topsom, R. D. Acidities of 4-substituted benzoic, bicyclo[2.2.2]oct-1-yl and bicyclo[2.2.2]oct-2-enyl carboxylic acids. J. Phys. Org. Chem. 1993, 6, 685−689. (6) Hansch, C.; Leo, A.; Taft, R. W. A survey of Hammett substituent constants and resonance and field parameters. Chem. Rev. 1991, 91, 165−195. (7) Kirkwood, J. G.; Westheimer, F. H. The electrostatic influence of substituents on the dissociation constants of organic acids. I. J. Chem. Phys. 1938, 6, 506−512. (8) Westheimer, F. H.; Kirkwood, J. G. The electrostatic influence of substituents on the dissociation constants of organic acids. II. J. Chem. Phys. 1938, 6, 513−517. (9) Topsom, R. D. The isolated molecule approach. Theoretical studies of the inductive effect. J. Am. Chem. Soc. 1981, 103, 39−44. (10) Topsom, R. D. Contribution of theoretical chemistry to an understanding of electronic substituent effects. Acc. Chem. Res. 1983, 16, 292−298. (11) Exner, O.; Bohm, S. Electrostatic calculation of the substituent effect: An efficient test on isolated molecules. Chem. − Eur. J. 2003, 9, 4718−4723. (12) Gadre, S. R.; Suresh, C. H. Electronic perturbations of the aromatic nucleus: Hammett constants and electrostatic potential topography. J. Org. Chem. 1997, 62, 2625−2627. (13) Galabov, B.; Ilieva, S.; Schaefer, H. F., III An efficient computational approach for the evaluation of substituent constants. J. Org. Chem. 2006, 71, 6382−6387. (14) Galabov, B.; Ilieva, S.; Hadijeva, B.; Atanasov, Y.; Schaefer, H. F., III Predicting reactivities of organic molecules. Theoretical and experimental studies on the aminolysis of phenyl acetates. J. Phys. Chem. A 2008, 112, 6700−6707. (15) Sadlej-Sosnowska, N. Molecular similarity based on atomic electrostatic potential. J. Phys. Chem. A 2007, 111, 11134−11140. (16) Remya, G. S.; Suresh, C. H. Quantification and classification of substituent effects in organic chemistry: A theoretical molecular electrostatic potential study. Phys. Chem. Chem. Phys. 2016, 18, 20615−20626. (17) Sayyed, F. B.; Suresh, C. H.; Gadre, S. R. Appraisal of throughbond and through-space substituent effects via molecular electrostatic potential topography. J. Phys. Chem. A 2010, 114, 12330−12333. (18) te Velde, G.; Bickelhaupt, F. M.; Baerends, E. J.; van Gisbergen, S. J. A.; Fonseca Guerra, C.; Snijders, J. G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931−967. (19) Fernández, I.; Frenking, G. Correlation between Hammett substituent constants and directly calculated π-conjugation strength. J. Org. Chem. 2006, 71, 2251−2256. (20) Pross, A.; Radom, L.; Taft, R. W. Theoretical approach to substituent effects. Phenols and phenoxide ions. J. Org. Chem. 1980, 45, 818−826. (21) Hehre, W. J.; Radom, L.; Schleyer, P. vR.; Pople, J. A. Ab Initio Molecular Orbital Theory; John Wiley & Sons: New York, 1986. (22) Cyrański, M. K. Energetic aspects of cyclic pi-electron delocalization: Evaluation of the methods of estimating aromatic stabilization energies. Chem. Rev. 2005, 105, 3773−3811. (23) Siodła, T.; Oziminski, W. P.; Hoffmann, M.; Koroniak, H.; Krygowski, T. M. Toward a physical interpretation of substituent effects: The case of fluorine and trifluoromethyl groups. J. Org. Chem. 2014, 79, 7321−7331. (24) Krygowski, T. M.; Oziminski, W. P. Substituent effects in 1nitro-4-substituted bicyclo[2.2.2]octane derivatives: Inductive or field effects? J. Mol. Model. 2014, 20, 2352−2359. (25) Szatylowicz, H.; Siodla, T.; Stasyuk, O. A.; Krygowski, T. M. Towards physical interpretation of substituent effects: The case of meta- and para- substituted anilines. Phys. Chem. Chem. Phys. 2016, 18, 11711−11721.
(26) Sadlej-Sosnowska, N. Substituent active region - a gate for communication of substituent charge with the rest of a molecule: Monosubstituted benzenes. Chem. Phys. Lett. 2007, 447, 192−196. (27) Krygowski, T. M.; Sadlej-Sosnowska, N. Towards physical interpretation of Hammett constants: Charge transferred between active regions of substituents and a functional group. Struct. Chem. 2011, 22, 17−22. (28) Stasyuk, O. A.; Szatylowicz, H.; Fonseca Guerra, C.; Krygowski, T. M. Theoretical study of electron-attracting ability of the nitro group: classical and reverse substituent effects. Struct. Chem. 2015, 26, 905−913. (29) Taft, R. W. Sigma values from reactivities. J. Phys. Chem. 1960, 64, 1805−1815. (30) Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648−5652. (31) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785−789. (32) Ditchfield, R.; Hehre, W. J.; Pople, J. A. Self-consistent molecular-orbital methods. IX. An extended gaussian-type basis for molecular-orbital studies of organic molecules. J. Chem. Phys. 1971, 54, 724−728. (33) Sadlej-Sosnowska, N. On the way to physical interpretation of hammett constants: How substituent active space impacts on acidity and electron distribution in p-substituted benzoic acid molecules. Polish J. Chem. 2007, 81, 1123−1134. (34) Hirshfeld, F. L. Bonded-atom fragments for describing molecular charge densities. Theor. Chim. Acta 1977, 44, 129−138. (35) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, M. J.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09; Gaussian, Inc.: Wallingford, CT, 2009.
1749
DOI: 10.1021/acsomega.7b00390 ACS Omega 2017, 2, 1746−1749
