Curtin–Hammett-Driven Intramolecular Cyclization of Heteroenyne

Article Cite This: J. Org. Chem. 2018, 83, 1790−1796

pubs.acs.org/joc

Curtin−Hammett-Driven Intramolecular Cyclization of Heteroenyne− Allenes to Phenanthridine-Fused Quinazoliniminiums Pamela S. Filby and Sundeep Rayat* Department of Chemistry, Ball State University, Cooper Physical Science Building, Muncie, Indiana 47304-0445, United States S Supporting Information *

ABSTRACT: Intramolecular cyclization of the heteroenyne−allene 2-((biphenyl-2-ylimino)methyleneamino)benzonitrile 1 to phenanthridine-fused quinazoliniminium salt PQ in the presence of a Lewis acid at room temperature involves formation of two new bonds: a C− C bond and a C−N bond. In this article, density functional theory (B3LYP and M06-2X) was employed in conjunction with 6-311G* basis set to gain insights into the mechanism of this cyclization reaction. The solvent effects were considered using Polarizable Continuum Model with nitromethane as the solvent. Our calculations show that C−C bond formation precedes the C−N bond formation. Precisely, the mechanism involves initial protonation of 1 at Nα and Nβ of the carbodiimide to form rapidly equilibrating conformers of the tautomers 2a,b and 3a,b. The Curtin−Hammett principle is invoked to determine the course of the reaction from these protonated species, which predicts that the intramolecular Friedel−Crafts type N-acylation (C−C bond formation) occurs between the protonated carbodiimide and biphenyl ring of the isomer 3b to form phenanthridinium cation 6b via transition state TS3b6b. Once 6b is formed, it converts to its most stable tautomers 8R and 9a. Once again, the Curtin−Hammett principle suggests that intramolecular nucleophilic attack is preferred from the tautomer 8R, where phenanthridine N-atom (Nβ) attacks the protonated nitrile group (C−N bond formation) and results in the formation of intermediate 11 via TS8R11. 11 then tautomerizes to the most stable cation 13. The coordination of the latter with the chloride anion yields the phenanthridine-fused heterocyclic salt PQ with overall release of energy. The pathways involving protonation at the nitrile (Nγ) of 1 were found to be energetically unfavorable and thus insignificant to the mechanism of cyclization.

■

INTRODUCTION Fused heterocycles are represented in many natural and synthetic compounds that exhibit a wide range of medicinal applications.1 As a result, there is a high interest in the pursuit of unique ring-fused structures because of their potential to provide new classes of pharmaceutical agents. For instance, the phenanthridine ring system is featured in compounds that exhibit anticancer,2−4 antituberculosis,5,6 antitrypanosomiasis,7,8 antihepatitis,9,10 and pesticidal activities.11,12 Similarly, the 4(3H)-quinazoliniminium motif is an important pharmacophore with cholinesterase inhibitory,13−16 antihypertensive,17 antimicrobial,18 and antiproliferative activities.19,20 This scaffold also represents a key precursor to 4(3H)-quinazolinones, which have emerged as important targets in medicinal chemistry with a myriad of interesting biological activities such as antimicrobial,21 antioxidant,22 antimalarial,23,24 anticonvulsive,25 analgesic and anti-inflammatory,26,27 anticancer,28,29 and antihypertensive.30,31 The fusion of phenanthridine and quinazolinimine cores offers the opportunity to expand the library of biologically active structures and to facilitate drug discovery. However, prior to our work,32 no synthetic routes existed in the literature that united the phenanthridine and quinazolinimine ring systems in one molecular skeleton. We have developed an efficient and versatile protocol to obtain phenanthridine-fused quinazolini© 2018 American Chemical Society

minium salts PQ at room temperature that involved one-pot Lewis acid-assisted cascade/tandem cyclization of heteroenyne−allenes, specifically the 2-((biphenyl-2-ylimino)methyleneamino)benzonitriles 1 (Scheme 1). The reaction is believed to involve in situ generation of a Bronsted acid (e.g., HX, where X = Cl− or BF4−) from Lewis acid (MX) and trace water. A series of PQ functionalized at various positions on the scaffold were prepared in this way from the corresponding heteroenyne−allenes 1 in moderate to good yields.32 The cyclization of 1 to PQ involves (a) C−N bond formation and (b) a C−C bond formation. The sequence in which C−C and C−N bond formations take place was not clearly evident. On the basis of a previous study from our research group,33 it was initially hypothesized that the reaction may involve the formation of intermediate I (Scheme 1) via the addition of the Bronsted acid HX on the carbodiimide moiety of 1 followed by nucleophilic addition of the Nβ on the protonated nitrile group (C−N bond formation). We envisaged that I would subsequently undergo intramolecular Friedel− Crafts-type N-acylation to form PQ (C−C bond formation). To explore this possibility, I was independently synthesized, Received: October 18, 2017 Published: January 19, 2018 1790

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796

Article

The Journal of Organic Chemistry

Scheme 1. Phenanthridine-Fused Quinazoliniminium Salt PQ from Heteroenyne−Allene 1 in the Presence of a Bronsted Acid (where X = Cl−, BF4−)

Figure 1. Protonation of heteroenyne−allene 1 in the presence of HCl. Energies are computed at B3LYP/6-311G* and include the free chloride ion. Optimized structures of 2a,b, 3a,b, and 4a are also shown.

and when treated with Lewis acid, it failed to cyclize to PQ.32 This indicated that the formation of PQ may not involve the intermediate I. In this article, we report a systematic computational investigation of several mechanistic avenues that could be invoked involving the cyclization of 1 to the fused heterocyclic salt PQ. Our primary goal was to identify critical reactive intermediates and the nature of transition states involved in the formation of the final product that may subsequently aid in substrate selection in our future work or enable us to modify or improve experimental conditions to increase the efficiency and scope of this reaction. The identification of key intermediates may also suggest opportunities to access other chemical frameworks of high importance using these reaction conditions.

■

potential energy surface were obtained using density functional theory (B3LYP35,36 and M06-2X37) with the employment of the 6-311G*38,39 basis set. Optimizations were followed by vibrational analyses to ensure an imaginary frequency of zero for the ground states and one for the transition states. Each transition state connecting the corresponding reactant and product was confirmed by following the Intrinsic Reaction Coordinate40−43 (IRC) path in both directions. Solvent effects on the reaction energies and activation barriers were considered by performing single-point Polarizable Continuum Model (PCM)44 calculations for nitromethane (ε = 38.2) on optimized gasphase geometries. Thermochemical data obtained in gas phase were used with the corresponding single-point solvation energies to calculate Gibbs free energies in the presence of solvent (ΔG(solv)) at both the levels of theory. The discussion below is based on ΔG(solv) for the level indicated. The conformational entropy effect was expected to be small45 and thus not considered in our computations.

■

RESULTS AND DISCUSSION Reaction Energetics at B3LYP Level of Theory. The various possible pathways for the cyclization of heteroenyne−

COMPUTATIONAL METHODS

All calculations were performed using the Gaussian 09 suite of programs.34 Optimized geometries of the stationary points on the 1791

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796

Article

The Journal of Organic Chemistry allene 1 to phenanthridine-fused heterocyclic salt PQ in the presence of HCl as Bronsted acid are discussed below. Protonation of Heteroenyne−Allene 1. In the presence of HCl, there are three possible sites for protonation of 1, at Nα, Nβ, and Nγ, to form tautomers 2, 3, and 4, respectively (Figure 1). Protonation at Nα and Nβ may produce different conformers of 2 and 3, respectively. We consider two synclinal conformers of 2 with respect to the protonated C−Nα bond, in which the NCN fragment displayed a dihedral angle of 23° (2a) and −44.5° (2b) with the plane of the benzonitrile ring, respectively. Similarly, two conformers are considered for 3. In 3a, ring A of the biphenyl is oriented nearly perpendicular with respect to the carbodiimide, while in 3b, the ring A is in the same plane as the carbodiimide. In the case of protonation at the nitrile (at Nγ), the optimized structure of the tautomer 4a is also shown in Figure 1. The highest occupied molecular orbital (HOMO) of heteroenyne−allene 1 was analyzed to determine the preferred sites of protonation (Figure 2). We note that the molecular

Figure 3. Energy diagram computed at B3LYP/6-311G* for C−C bond formation from protonated carbodiimide species 2a,b and 3a,b. Energies are relative to the reactants (1 + HCl) and thus include the free chloride ion.

Figure 2. HOMO of heteroenyne−allene 1 at B3LYP/6-311G*.

orbital component is mostly localized on the carbodiimide moiety (NCN) and on the biphenyl ring attached to Nβ. There is a little contribution of the molecular orbital from the CN group. This implies that Nα and Nβ of 1 are more basic than Nγ, and therefore protonation at the former is expected to be thermochemically more favorable than the latter. This hypothesis is also supported by literature precedence that reports carbodiimides to exhibit proton affinities higher than that of nitriles.46 A preference for protonation at carbodiimide was also noted in our previous mechanistic study.47 Analogously, our current data indicate that protonation of 1 at Nα to form conformers 2a and 2b (+22.3 and +21.7 kcal/ mol) and at Nβ to form the conformers 3a and 3b (+23.2 and +24.8 kcal/mol) is more favorable than protonation at Nγ to give 4a (+30.5 kcal/mol) as shown in Figure 1. There is also not a significant difference in the ΔG(solv) values for the protonation of Nα and Nβ, even though the frontier molecular orbital suggests that Nβ is slightly more basic than Nα and hence more likely to get protonated (Figure 2). Chemistry from Protonated Carbodiimide. C−C Bond Formation. On the basis of the ΔG(solv) values, we predict that protonation of 1 occurs at both Nα and Nβ, leading to the formation of more stable tautomers 2b and 3a, respectively, which are expected to rapidly interconvert (Figure 3). 2b also undergoes fast conversion to its conformer 2a which is only 0.6 kcal/mol higher in energy. The intramolecular Friedel−Craftstype N-acylation involving the protonated carbodiimide and the ring B of 2a and 3a leads to the formation of syn arenium ions 5a and 6a via transition states TS2a5a and TS3a6a, respectively (Figure 4). This process is slightly endergonic for 2a → 5a

Figure 4. Transition states TS2a5a, TS2b5b, TS3a6a, and TS3b6b for the C−C bond formation (B3LYP/6-311G*).

(+0.7 kcal/mol), while exergonic for 3a → 6a (−3.0 kcal/mol), and occurs with an activation barrier of 15.8 and 17.8 kcal/mol, respectively (Figure 3). Note that in 5a and 6a, the benzonitrile group is syn to the CN and C−N bond of the phenanthridine ring system. Similarly, the conversion of 3a to its slightly higher energy conformer 3b (+1.6 kcal/mol) is expected to be fast, and intramolecular Friedel−Crafts-type N-acylation from 2b and 3b leads to the formation the anti arenium ions 5b and 6b via transition states TS2b5b and TS3b6b (Figures 3 and 4). The 1792

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796

Article

The Journal of Organic Chemistry

Figure 5. Energies of the various tautomers computed at B3LYP/6-311G*. Structures in blue could not be obtained as stationary points. All energies are in kcal/mol and are relative to the reactants (1 + HCl) and thus include the free chloride ion.

Figure 6. Energy diagram for C−N bond formation computed at B3LYP/6-311G*. Reaction energies are relative to the reactants (1 + HCl) and thus include the free chloride ion.

bond formation from 2b → 5b and 3b → 6b is endergonic by +5.4 and +0.8 kcal/mol, respectively. The transitions states

benzonitrile group in 5b and 6b is anti to the CN bond and C−N bond of the phenanthridine ring, respectively. The C−C 1793

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796

Article

The Journal of Organic Chemistry

Figure 7. Tautomerization of 11 and 12 to the final cation 13 and its coordination with chloride to form the product PQ (B3LYP/6-311G*). All energies are relative to the reactants (1 + HCl) and include the free chloride ion (except the salt PQ).

to form the C−N bond. Similarly, the C−N bond formation from the high energy tautomer 10R to 13 was also investigated, which occurred via rotation to 10a followed by cyclization. All of these pathways involve transition states that are considerably higher in energy on the potential energy surface compared to the pathway from tautomer 8R to 11. For instance, transition states TS5a11 and TS5b11 are +24.6 kcal/mol and +26.2 kcal/ mol, and TS10R13 is +62.2 kcal/mol higher than TS8R11 on the potential energy surface and are therefore not expected to play any significant role (see Supporting Information). Chemistry from Protonated Nitrile Group. Because the protonation of heteroenyne−allene 1 at Nγ to form 4a is less favored (Figure 1) compared to protonation at Nα, and Nβ, any further chemistry emanating from this species is expected to have a minor contribution. Our exploration of the potential energy surface reveals that the formation of the arenium ion 15 from 4a occurs in a stepwise fashion (Figure 8). If the C−C

TS2b5b and TS3b6b are found to be +20.6 and +11.5 kcal/mol higher in energy with respect to the corresponding 2b and 3b. The Curtin−Hammett principle is employed to predict the reaction pathway. This principle states that the product ratio (e.g., 5a,b and 6a,b) from rapidly interconverting isomers (e.g., 2a,b and 3a,b) is derived by the relative energies of the transition states, (e.g., TS2a5a, TS2b5b, TS3a6a and TS3b6b) and is not significantly affected by the relative energies of the interconverting isomers.48,49 Therefore, we propose that tautomer 6b is predominantly formed, as the energy of TS3b6b is the lowest on the potential energy surface (Figure 3). Tautomerizations on the Potential Energy Surface and C−N Bond Formation. Once 6b is formed, it has the option of undergoing rotational change to form 6a or tautomerization/ conformational change to 5a,b, 7a,b, and 9a,b (Figure 5). The structures 5a,b, and 7a,b can further tautomerize to 8a,b and 10a,b, respectively as shown. We were able to obtain minimum structures for 9a and 9b. However, we were unable to obtain 7a,b, 8a,b, and 10a,b. Instead, we obtained the rotamers 7R, 8R, and 10R, featuring the benzonitrile group syn to the CN of the phenanthridine ring analogous to 5a and 6a; however, its CN group is oriented opposite to the CN of the phenanthridine ring (Figure 5). The tautomer 9a is predicted to be most stable. The structures 7R and 10R are considerably higher in energy (+46.8 and +37.9 kcal/mol, respectively) from 9a, owing to their highly unstable quinoid-type structure. The tautomers 8R and 9b are only 0.7 and 1.8 kcal/mol less stable than 9a. Thus, we predict that 8R and 9a,b predominate on the potential energy surface. Our search for transition states for C− N bond formation from these tautomers revealed that 8R undergoes rotational change to 8a which then undergoes intramolecular nucleophilic attack from Nβ to its protonated CN to yield 11 via TS8R11 with an activation barrier of 5.5 kcal/mol (Figures 6). The cyclization of 8R → 11 is highly exergonic by −24.8 kcal/mol. Similarly, 9b undergoes E to Z isomerism with respect to its CNα bond to form slightly more stable 9a via TS9b9a (Ea = +2.4 kcal/mol) (see Supporting Information) which then ring closes to yield 12 via transition state TS9a12 (Ea = +15.8 kcal/mol) (Figure 6). This process is endergonic by +14.2 kcal/mol. Because 8R and 9a are rapidly interconverting isomers, once again the Curtin− Hammett principle48,49 suggests that the pathway 8R → 11 will be favored because transition state TS8R11 is −9.6 kcal/mol lower in energy compared to TS9a12 (Figure 6). Once 11 is formed, it tautomerizes to its most stable tautomer 13 (−6.4 kcal/mol more stable). The coordination of a chloride anion to 13 gives the phenanthridine-fused heterocyclic salt PQ with −8.5 kcal/mol free energy decrease (Figure 7). The transition states for the C−N bond formation from tautomers 5a,b to 11 were also located. These pathways occur via conversion to tautomers 8a,b which is spontaneously followed by nucleophilic attack from Nβ to the protonated CN

Figure 8. High energy pathways from protonated nitrile 4a computed at B3LYP/6-311G*. Energies are relative to the reactants (1 + HCl) and thus include free chloride ion (except for the salt PQ).

bond formation (Friedel−Crafts-type N-acylation) occurs prior to C−N bond formation (electrocyclization) via transition state TS4a15 (Ea = +30.5 kcal/mol), intermediate 7a is involved. However, if C−N bond formation (intramolecular nucleophilic attack) occurs before C−C bond formation (Friedel−Crafts type N-acylation) via transition state TS4a15′ (Ea = +32.1 kcal/ mol), then intermediate 14a is involved. The process 4 → 15 is exergonic by −5.9 kcal/mol. We note that 15 is a tautomer of 11, 12, and 13, albeit very unstable. Thus, once 15 forms, it would undergo proton migration to form the most stable tautomer 13 with a free energy decrease of −50.9 kcal/mol. The latter upon coordination to chloride anion would yield the final product PQ, as shown in Figure 7. Both 7a and 14a are not viable intermediates on the potential energy surface. These high energy structures form during the cyclization of 4a → 15, 1794

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796

Article

The Journal of Organic Chemistry and spontaneously undergo ring closure. The pathway 4a → 15 is energetically demanding as the high energy points involving transition states TS4a15 or TS4a15′ are significantly higher on the potential energy surface (at +61.0 kcal/mol and +62.6 kcal/ mol relative to the reactants, see Supporting Information) compared to high energy points on the pathway from protonated carbodiimide involving transition states TS3b6b and TS8R11 (at +36.3 kcal/mol and +10.3 kcal/mol, respectively relative to the reactants). Therefore, we deduce that cyclization 1 → PQ does not occur from the protonated nitrile intermediate 4a. Reaction Energetics at M06-2X Level of Theory. The major pathways involved in the cyclization of 1 to PQ were also investigated using dispersion corrected functional M06-2X with the employment of 6-311G* basis set. Although M06-2X energies were lower than B3LYP (see Supporting Information), the overall conclusions remain the same at both the levels of theory. Briefly, the protonation at nitrile group of 1 is found to be less favorable than protonation at the carbodiimide, as tautomer 4a (+21.1 kcal/mol) is significantly higher in energy than 2a (+12.3 kcal/mol) and 2b (+13.4 kcal/mol) as well as 3a (+12.1 kcal/mol) and 3b (+11.5 kcal/mol) at M06-2X/6311G* (Figure S21). The C−C bond formation is predicted to take place from 3b via transition state TS3b6b (+26.0 kcal/mol) to form 6b in accordance with Curtin−Hammett principle, as the other transition states TS2a5a (+27.8 kcal/mol), TS2b5b (+30.0 kcal/mol), and TS3a6a (+32.8 kcal/mol) are higher in energy (Figure S22). Once 6b is formed, it tautomerizes to the more stable tautomers 8R (−6.5 kcal/mol), 9a (−9.2 kcal/ mol), and 9b (−5.7 kcal/mol), and 8R (−6.5 kcal/mol). Analogous to B3LYP, the tautomers 7R (+41.9 kcal/mol) and 10R (+29.1 kcal/mol) are found to be significantly higher in energy at M06-2X/6-311G* (Figure S23). 9b undergoes E−Z isomerism to more stable 9a via TS9b9a (−0.7 kcal/mol). Again, the Curtin−Hammett principle suggests that C−N bond formation takes place from 8R via transition state TS8R11 (−3.4 kcal/mol) because it is lower in energy compared to TS9a12 (+6.5 kcal/mol) (Figure S24). Once 11 (−35.8 kcal/ mol) is formed, it tautomerizes to more stable 13 (−41.7 kcal/ mol) which upon coordination with chloride anion forms the salt PQ (−41.8 kcal/mol) (Figure S25).

intermediate 4a is unlikely at room temperature due to the high energetics computed for this pathway. In summary, the mechanistic insights gained herein suggest that this chemistry could be exploited to construct unique fused heterocycles by careful choice of the aromatic substituent on Nβ of 1 or the electrophilic units on the scaffold. To our knowledge, this is the first report of Friedel−Crafts-type reaction employing protonated carbodiimides as N-acylating agents and thus may provide a viable route to synthesize N,Ndisubstituted amidines, which are a prominent class of nitrogencontaining functional groups exhibiting a wide range of biological and pharmaceutical profiles50−54 as well as material applications.55,56 This work is under investigation and will be reported in the future.

■

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.joc.7b02643. Proposed mechanistic pathway; optimized structures of the reactants, intermediates, transition states, and products; total energies and vibrational data of all optimized structures, and single-point PCM energies; high energy pathways computed at B3LYP/6-311G*; all energy diagrams computed at M06-2X/6-311G*; optimized Cartesian coordinates of relevant structures (PDF)

■

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Sundeep Rayat: 0000-0002-0071-3110 Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS The authors acknowledge the College of Sciences and Humanities Beowulf Cluster at Ball State University funded by the Office of Sponsored Program for computational time and support. S. Rayat is also grateful for the Senior Research Grant from the Indiana Academy of Science and the ASPiRE Junior Faculty Award from Ball State University for the partial support of this work.

■

CONCLUSION In this article, density functional theory has been employed to investigate the mechanism of tandem/cascade cyclization of heteroenyne allene 1 to the phenanthridine-fused quinazoliniminium salt PQ through the formation of C−C and C−N bonds. Our work reveals that C−C bond formation (Friedel− Crafts type N-acylation) precedes the C−N bond formation (nucleophilic attack), and the reaction involves initial protonation of 1 at the carbodiimide unit to form 2a,b and 3a,b. The mechanism is driven by the Curtin−Hammett principle which suggests that C−C bond formation involves the pathway 1 → 2b → 3a → 3b → TS3b6b → 6b. Once 6b is formed, it tautomerizes to its most stable tautomers 8R and 9a,b. Again, the Curtin−Hammett principle guides the C−N bond formation from 8R → 11 via transition state TS8R11. Once 11 is formed, it tautomerizes to the more stable structure 13. The capture of Cl anion by the latter gives the salt PQ. Overall, the formation of PQ is predicted to occur with a −34.9 and −41.8 kcal/mol free energy decrease from the reactants 1 and HCl at B3LYP and M06-2X levels, respectively. Our investigation reveals that cyclization from the protonated nitrile

■

REFERENCES

(1) Bucher, G.; Mahajan, A. A.; Schmittel, M. J. Org. Chem. 2008, 73, 8815. (2) Elsayed, M. S. A.; Zeller, M.; Cushman, M. Synth. Commun. 2016, 46, 1902. (3) Zhan, G.; Qu, X.; Liu, J.; Tong, Q.; Zhou, J.; Sun, B.; Yao, G. Sci. Rep. 2016, 6, 33990. (4) Tang, Z.-H.; Cao, W.-X.; Wang, Z.-Y.; Lu, J.-H.; Liu, B.; Chen, X.; Lu, J.-J. Redox Biol. 2017, 12, 367. (5) Cinar, M. E.; Vavilala, C.; Fan, J.; Schmittel, M. Org. Biomol. Chem. 2011, 9, 3776. (6) Cohen, A. J.; Mori-Sanchez, P.; Yang, W. Chem. Rev. 2012, 112, 289. (7) Newton, B. A.; Gilbert, R. J.. Agrochemicals: Fate and Food in the Environment; Proceedings of a Symposium, Rome, June 7−11, 1982; International Atomic Energy Agency: Vienna, 1982; p 255. (8) Boibessot, I.; Turner, C. M. R.; Watson, D. G.; Goldie, E.; Connel, G.; McIntosh, A.; Grant, M. H.; Skellern, G. G. Acta Trop. 2002, 84, 219. 1795

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796

Article

The Journal of Organic Chemistry

(37) Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215. (38) Frisch, M. J.; Pople, J. A.; Binkley, J. S. J. Chem. Phys. 1984, 80, 3265. (39) Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A. J. Chem. Phys. 1980, 72, 650. (40) Fukui, K. Acc. Chem. Res. 1981, 14, 363. (41) Hratchian, H. P.; Schlegel, H. B. Theory and Applications of Computational Chemistry: The First Forty Years; Elsevier: New York, 2005; p 195 10.1016/B978-044451719-7/50053-6. (42) Hratchian, H. P.; Schlegel, H. B. J. Chem. Theory Comput. 2005, 1, 61. (43) Gonzalez, C.; Schlegel, H. B. J. Chem. Phys. 1989, 90, 2154. (44) Miertus, S.; Scrocco, E.; Tomasi, J. Chem. Phys. 1981, 55, 117. (45) Viskolcz, B.; Fejer, S. N.; Csizmadia, I. G. J. Phys. Chem. A 2006, 110, 3808. (46) Elrod, M. J. Int. J. Mass Spectrom. 2003, 228, 91. (47) Filby, P. S.; Rayat, S. J. Org. Chem. 2015, 80, 4299. (48) Cramer, C. J. Essentials of Computational Chemistry Theories and Models; John Wiley & Sons Ltd.: West Sussex, England, 2004. (49) Anslyn, E. V.; Dougherty, D. A. Modern Physical Organic Chemistry; University Science Books: Sausalito, CA, 2005. (50) Greenhill, J. V.; Lue, P. Prog. Med. Chem. 1993, 30, 203. (51) Kalisiak, J.; Ralph, E. C.; Zhang, J.; Cashman, J. R. J. Med. Chem. 2011, 54, 3319. (52) Ma, Z.; Rao, L.; Bierbach, U. J. Med. Chem. 2009, 52, 3424. (53) Mathews, T. P.; Kennedy, A. J.; Kharel, Y.; Kennedy, P. C.; Nicoara, O.; Sunkara, M.; Morris, A. J.; Wamhoff, B. R.; Lynch, K. R.; MacDonald, T. L. J. Med. Chem. 2010, 53, 2766. (54) Calas, M.; Ouattara, M.; Piquet, G.; Ziora, Z.; Bordat, Y.; Ancelin, M. L.; Escale, R.; Vial, H. J. Med. Chem. 2007, 50, 6307. (55) Beh, E. S.; Tong, L.; Gordon, R. G. Organometallics 2017, 36, 1453. (56) Cho, W.; Cho, H.; Lee, C. S.; Lee, B. Y.; Moon, B.; Kang, J. Organometallics 2014, 33, 1617.

(9) Aoyama, H.; Sugita, K.; Nakamura, M.; Aoyama, A.; Salim, M. T. A.; Okamoto, M.; Baba, M.; Hashimoto, Y. Bioorg. Med. Chem. 2011, 19, 2675. (10) Chen, D.; Cai, J.; Yin, J.; Jiang, J.; Jing, C.; Zhu, Y.; Cheng, J.; Di, Y.; Zhang, Y.; Cao, M.; Li, S.; Peng, Z.; Hao, X. Future Med. Chem. 2015, 7, 561. (11) Miao, F.; Yang, X.-J.; Ma, Y.-N.; Zheng, F.; Song, X.-P.; Zhou, L. Chem. Pharm. Bull. 2012, 60, 1508. (12) Lv, P.; Chen, Y.; Shi, T.; Wu, X.; Li, Q. X.; Hua, R. Pestic. Biochem. Physiol. 2017, DOI: 10.1016/j.pestbp.2017.06.009. (13) Decker, M. Eur. J. Med. Chem. 2005, 40, 305. (14) Decker, M. J. Med. Chem. 2006, 49, 5411. (15) Decker, M.; Kraus, B.; Heilmann, J. Bioorg. Med. Chem. 2008, 16, 4252. (16) Decker, M.; Krauth, F.; Lehmann, J. Bioorg. Med. Chem. 2006, 14, 1966. (17) Patel, H. U.; Patel, R. S.; Patel, C. N. J. Appl. Pharm. Sci. 2013, 3, 171. (18) Prabhakar, V.; Jagadeesh, S.; Bhargavi, P.; Vani, C. D. Heterocycl. Lett. 2016, 6, 725. (19) Becerra, A.; Quintero, C.; Morales, V.; Valderrama, M.; Aguirre, A.; Faundez, M. A.; Rojas, R. S. Bioorg. Med. Chem. 2017, 25, 2681. (20) Perchellet, J.-P. H.; Waters, A. M.; Perchellet, E. M.; Naganaboina, V. K.; Chandra, K. L.; Desper, J.; Rayat, S. Anticancer Res. 2011, 31, 2083. (21) Chand, M.; Gupta, A.; Jain, S. C. Heterocycl. Lett. 2017, 7, 201. (22) Prashanth, M. K.; Revanasiddappa, H. D. Lett. Drug Des. Discovery 2014, 11, 712. (23) Birhan, Y. S.; Bekhit, A. A.; Hymete, A. BMC Res. Notes 2015, 8, 589. (24) Bule, M. H.; Haymete, A.; kefale, B. Indo Am. J. Pharm. Res. 2015, 5, 963. (25) Noureldin, N. A.; Kothayer, H.; Lashine, E.-S. M.; Baraka, M. M.; El-Eraky, W.; El Awdan, S. A. Arch. Pharm. 2017, 350, 1600332. (26) Abdel-Aziz, A. A. M.; Abou-Zeid, L. A.; ElTahir, K. E. H.; Mohamed, M. A.; Abu El-Enin, M. A.; El-Azab, A. S. Bioorg. Med. Chem. 2016, 24, 3818. (27) Dewangan, D.; Verma, V. S.; Nakhate, K. T.; Tripathi, D. K.; Kashyap, P.; Dhongade, H. Med. Chem. Res. 2016, 25, 2143. (28) Patel, H. M.; Pawara, R.; Ansari, A.; Noolvi, M.; Surana, S. Bioorg. Med. Chem. 2017, 25, 2713. (29) Sun, J.; Zhou, L.; Tan, G.; Li, S.; Wang, S.; Chen, H.; Li, X. Youji Huaxue 2017, 37, 455. (30) Levin, J. I.; Venkatesan, A. M.; Chan, P. S.; Bailey, T. K.; Vice, G.; Coupet, J. Bioorg. Med. Chem. Lett. 1994, 4, 1819. (31) Alagarsamy, V.; Pathak, U. S. Bioorg. Med. Chem. 2007, 15, 3457. (32) Alawode, O. E.; Naganaboina, V. K.; Liyanage, T.; Desper, J.; Rayat, S. Org. Lett. 2014, 16, 1494. (33) Naganaboina, V. K.; Chandra, K. L.; Desper, J.; Rayat, S. Org. Lett. 2011, 13, 3718. (34) Frisch, M. J.; 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, Ö .; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09; Wallingford, CT, 2009. (35) Becke, A. D. J. Chem. Phys. 1993, 98, 5648. (36) Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785. 1796

DOI: 10.1021/acs.joc.7b02643 J. Org. Chem. 2018, 83, 1790−1796