Molecular Recognition, Conformational Behavior ... - ACS Publications

Dec 21, 2017 - Department of Chemistry, Savitribai Phule Pune University, Pune 411 007, India. •S Supporting Information. ABSTRACT: In the present w...
1 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Article

Molecular Recognition, Conformational Behavior and Spectral Characteristics of Oxatub[4]arene Macrocycle Dipali N. Lande, and Shridhar P. Gejji J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b12472 • Publication Date (Web): 21 Dec 2017 Downloaded from http://pubs.acs.org on December 22, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Molecular Recognition, Conformational Behavior and Spectral Characteristics of Oxatub[4]arene Macrocycle Dipali N. Lande and Shridhar P. Gejji Department of Chemistry, Savitribai Phule Pune University, 411 007, India

Corresponding author: [email protected] Fax No.: +91-20-225691728 Telephone No.: +91-20 2560122

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract In the present work we analyze molecular recognition behavior of synthetic hydroxylated oxatub[4]arene (TA4) receptor toward the methyl viologen in different redox states. The supramolecular binding of methyl viologen guest toward TA4 macrocyclic scaffold has been studied employing the dispersion corrected ωB97X-D based density functional theory (DFT). The methyl viologen in dicationic and neutral forms revealed distinct features in electronic, 1H NMR and infrared spectra. Quantum theory of atoms in molecules (QTAIM) in conjunction with the noncovalent interaction reduced density gradient (NCI-RDG) in real space, have been used as tools to characterize the underlying host-guest binding.

2 ACS Paragon Plus Environment

Page 2 of 33

Page 3 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Artificial molecular receptors endowed with the distinct structural characteristics and intriguing host-guest properties are the primary workhorses in the domain of supramolecular chemistry.1-4 The emergence of newer macrocyclic hosts have accelerated the expansion of this field with their potential applications in interdisciplinary areas including catalysis,5 molecular machines,6 supramolecular polymers,7 and drug delivery.8 During a past decade a variety of arene based molecular scaffolds for example, resorcinarenes,9 calixarenes,10 cyclotriveratrylenes,11 pillararenes12 and the modified cavitands incorporating pyrrole,13 pyridine14 or imidazolium15 heterocyclic groups instead of phenyl moieties in their architecture, have been obtained. A majority of these macrocyclic arenes are composed of mono-benzene or substituted benzene moieties. Extension from mono-arene to larger aromatic arenes including biphenyl, naphthols render resulting hosts with unique structural attributes such as deep and wider cavity, enhanced π-electron density etc., offer unquestionable advantages and make them fascinating.16,17 Interestingly the fluorescent characteristics of naphthalene endows the molecular container with natural affinity for sensing.18 The Li and coworkers19 synthesized methylene bridged biphenyl cyclo-oligomers (biphen[n]arene) using the reaction between 4,4’-biphenolether and paraformaldehyde in the presence of a Lewis acid catalyst. The hollow π-electron rich cavity of biphen[n]arene showed a large affinity toward multiple cationic as well as electron deficient neutral guests. As far as naphthalene monomer based macrocycles are concerned, calix[n]naphthalene scaffolds comprised of methylene bridge unit were derived from direct cyclocondensation reaction between 1-naphthol and formaldehyde under the basic conditions by Georghiou and coworkers.20 Further modifications at the bridging position of the macrocycle by amide or allyl functionalities was expected to

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

improve the guest binding ability and their potential applications as fluorescent sensors.21,22 However, majority of such naphthalene-based molecular cavitand suffered from difficult regioisomeric products with ill-defined and shallow cavities, and scarcely displayed satisfactory guest binding.23,24 To overcome these difficulties recently Jiang et.al.25-28 reported flexible naphthalene based macrocycles viz., oxatub[n]arene (n=4, 5 and 6) and synthesized by the one pot reaction under high dilution condition in which the repeating 2,6-dihydroxynaphthalene units are connected by ether (CH2-O-CH2) linkage at the 1,5-positions which minimize the potential isomerism and avoid the self-occupation of the cavity as well as helps in maintaining the deeper cavities. Additionally, the flexible naphthalene moieties lead to multiple conformers for effective guest binding. These oxatub[4]arene synthetic receptors can accommodate a variety of guests with varying size ranging from 1,4-diazabicyclo[2.2.2]octane, viologen derivatives to C60 or C70 fullerenes. In addition to this the guest binding and solid-state self-assemblies formation of zorb[4]arenes have been reported by the same group.29 It has been observed that the uncomplexed state, the such macrocycle assumes either collapsed or self-included conformations. In continuation to this Jiang et al. 30-34 reported a several molecular tubes based on the rigid bis-naphthalene cleft. Subsequently to mimic the biological functionalities these authors introduced converging functional groups such as urea and thiourea into the bis-napthalene cleft cavity which is referred as endo-functionalized bisurea/thiourea molecular tubes.35,36 On the theoretical front recently Gejji et al.37,38 studied structural features of inclusion complexes between the neutral guests and endofunctionalized urea/thiourea hosts. It has been shown that syn and anti configurational isomers of the macrocyle exhibit distinct selectivity toward guest binding and the complexation of centrosymmetric guest with the centrosymmetric anti isomer being energetically favored over the syn one.

4 ACS Paragon Plus Environment

Page 4 of 33

Page 5 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Viologens comprised of two pyridinium rings joined at para positions and have been known for electron relay in the photo- and electro- chemical reactions. Their redox and electrochromic behavior garnered considerable interest across the scientific community.39-41 As a result of strong electron accepting and the one/two- electron reversible chemical and electrochemical reduction ability the viologen are known to exist as stable radical cation or in neutral forms.42 Its radical-cation finds potential applications in printing and display technologies.43 Besides this, the host-guest chemistry of viologens has attracted growing attention since the bipyridinium ions can be used for design of redox-active ionic liquid crystals;44-47 as building blocks in molecular machines to trigger the motion of a subunit following proper electrochemical stimulus.48-50 Remarkably, these guests exhibit characteristic charge transfer transition bands which can be used for anion recognition.51-53 In particular, the highly noxious methyl viologen cation (trade name paraquate i.e. N,N-dialkyl,4,4-bipyridinium), (MV2+) has widely been act as a probe for analyses interactions of DNA or zeolites.54,55 Moreover, MV2+ conduces stable inclusion host-guest

complexes

with

cyclodextrin,56

cucurbit[n]uril,57

calix[n]arene,58

pillar[n]arene59 and other macrocycles. Thus, it is of interest to investigate the complexation of viologen with recently reported naphthalene based supramolecular scaffold by mean of density functional theory protocol with the purpose of designing redox-switchable molecule. In this spirit, the present endeavor emphasis on the molecular insights for the inclusion of methyl viologen in its different redox state to the smart artificial hydroxylated oxatub[4]arene (TA4) receptor. We explore the distinct behavior of cationic and neutral guest upon encapsulation within such host-cavity. Theoretical investigations also enlightening on how cationic and neutral binding reflect in their spectral characteristics. A synopsis on the interplay of diverse non-covalent

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 33

interactions underlying such host-guest binding has been presented using the NCI-RDG and QTAIM methods. The methodology adopted in this work is outlined below.

Computational Details Optimization of oxatub[4]arene (hosts) conformer, methyl viologen in its neutral and dicationic forms (guests) and their inclusion complexes were performed using the Gaussian 09 suit of program60 and Gauss View 5.0 utilized as graphic interface.61 The ωB97X-D hybrid exchange correlation functional proposed by Head-Gordon and coworkers62 was employed in conjunction with the internally stored 6-311+G(d,p) basis set. This level of theory is adequate for a faithful representation of underlying supramolecular interactions and describe well the electron distribution in host-guest chemistry.63-67 The tationary point structures of the host, guest and their complexes were confirmed to be the local minima on the multivariate potential energy surfaces from vibrational frequency calculations. The complexation of TA4 cavitands can be understood from the charge distribution within its cavity which can be characterized in terms of the molecular electrostatic potential (MESP) topography.68-70 To elucidate the nature and strength of underlying molecular interactions accompanying the guest encapsulation QTAIM analyses were carried out with the help of AIMAll-2000 software.71,72 Furthermore the non-localized dispersion and repulsive nonbonded interactions can be assured employing the NCI-RDG method as a tool.73 The RDG defined through s =



|∇ |

( )/ /

within the

QTAIM approach was calculated using the Multiwfn program74 and envisaged through the use of visual molecular dynamics software.75 Further the chemical shifts (δH) in 1H NMR spectra were derived by subtracting the nuclear magnetic shielding tensors of protons in the individual host, guest and their complexes from those in the tetramethylsilane used as the reference by the gauge-independent atomic orbital (GIAO) method.76 Effect of solvation (with the dichloromethane as solvent) on the structure and 1H NMR chemical 6 ACS Paragon Plus Environment

Page 7 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

shifts was simulated employing the self-consistent reaction field (SCRF) theory incorporating the polarization continuum model (PCM).77

Results and Discussion The naphthalene based TA4 supramolecular scaffold is made up of 2,6dihydroxynaphthalene units linked by ether bridges (CH2-O-CH2) at 1,5- position. The flipping of naphthalene monomer relative to linkages engenders different conformers with all the naphthalene monomers pointing along the same direction denoted as zigzag (I); one naphthalene ring flipped in a different direction referred as the disrupted zigzag (II); two neighboring naphthalene moieties flipped in different directions called as 1,2alternate (III) conformer and lastly two opposing naphthalene monomers flipped in different directions called as 1,3-alternate (IV) conformer. Stationary point structures thus derived were optimized in the framework of ωB97x-D/6-311+G(d,p) theory, which are displayed in Figure 1. All these conformers possess well defined cavities of tubular shape and varying dimensions. The cavity sizes are estimated by knowing the separation(s) of the radially opposite oxygens atoms in the naphthalene monomers. In an effort to understand the binding modes of hitherto conformers of TA4 the MESP topography was used. It has earlier been recognized that three dimensional (3D) MESP maps give a direct perspective of the segments in the molecular system along with its reactivity and affinity toward different guests. The MESP, V(r) generated by a molecule at a point r, is given by 

V(r) = 



ZA ρ(r )   −  d r |r − R  | |r − r  |

The positive potential produced by N nuclei with nuclear charges (ZA), situated at (RA) has been signified by first term, whilst the second term denotes the negative potential engendered by continuous electron density, (r). MESP attains positive as well as 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

negative values with the complementary of the above terms and emerge with the effective electron rich regions in the molecular system. To examine the chemical reactivity V(r) was computed on the molecular surface. We consider 0.001 au ( ≡ !"# ) contour of the isoelectronic density that encloses more than 95% of the electronic charge and generally located 0.2–0.5 Å beyond the van der Waals radii of the atoms.78,79 Projection of these isosurfaces along the molecular plane for the TA4 conformers are displayed in Figure 2. The variation in MESP mapped are depicted using the coding scheme showing: red the most negative potential regions appropriate for the electrophilic attack whereas the blue signifying most positive electrostatic potential imply spatial regions which are more susceptible for nucleophilic attack. The prevalence of green further corresponds to potential halfway between the two extremes red and blue. As shown, the electrostatic potential intensifications follows the order red < orange < yellow < green < blue. As can be seen from Figure 2 MESP maps in zigzag, disrupted zigzag, 1,2 alternate and 1,3 alternate TA[4] conformers exhibit a trend which is qualitatively similar. The oxygen atoms from hydroxyl group and those from ether linkages are rendered with the relatively large negative electrostatic potential (dark red), accordingly the oxygen lone pairs are available for facilitate intermolecular interactions with the guest. In addition to this, the delocalized π-electron clouds of naphthalene rings renders large negative valued MESP (red) continuum within the inner surface of the conformer than the outer surface (reddish yellow). It may conjecture that the endo-binding of guest is more favored over the exo-binding. The most reactive protons from the hydroxyl functionalities endowed with the maximum brunt of positive charge and are largely electron deficient. Further the positive MESP was observed near the aromatic protons, the ortho position protons showing large electron deficient character than those at the meta position while ether linkage protons (-CH2) reveal a zero valued electrostatic potential. It can be speculated

8 ACS Paragon Plus Environment

Page 8 of 33

Page 9 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that all TA4 conformers are electron-rich and capable of accommodating cationic or electron deficient neutral guests in their cavities. In this digest, we choose methyl viologen as a proto type guest in its reduced (MV0) and oxidized (MV2+) which are portrayed in Figure 3. MESP mapped onto density isosurfaces are shown along with. The MV0 displays the negative potential near viologen core while the MV2+ reveals large positive potential (dark blue). In order to achieve allosteric binding of such guests with TA4 conformers, the closer contact of the complementary (highest positive and negative) regions in the MESP has been crucial to promotes strong electrostatic attractions. Taking cue from these considerations different host-guest complexes was generated by electrostatic docking of the electron deficient sites around the MESP critical points in the macrocycle. In simple words the lock and key mechanism based on the complimentary characteristics of electron distribution governs the complexation. Optimized structures of complexes between MV0 as well as MV2+ guest and I to IV conformers of TA4 host are illustrated in Figure S1 to S2 along with the relative stabilization energies given in parentheses. The interaction energies are calculated from the supramolecular approach, by subtracting the sum of energies of the isolated TA4 conformer and individual guest from that of the complex. The inclusion complexes of neutral and dicationic form of methyl viologen with 1,2-alternate (III) conformer of the TA4 reveals stronger binding than other conformers owing to matching of shape and size. These inferences are in consonance with the single crystal X-ray structure experiments on the MV2+ complexes of octabutyloxatub[4]arene reported by Jiang and coworkers.25 Hence it may be surmise that the substitution of the alkyl chain at -OH functionality does not affect the host-guest binding patterns. To scrutinize the thermodynamic parameters of the inclusion phenomenon Gibbs free energy (∆G) parameters were computed which turned out to be < 0, indicating formation of the inclusion complex is spontaneous and

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

exergonic. Noteworthy enough the complexation of 1,2-alternate conformer binds to the charged species more favorably (~2.5 times) than the reduced viologen which can be attributed to attractive cation···π interactions those stem from the polarization of the electronic cloud in TA4 macrocycle. Structural parameters of the host, guests and their energy minimal complexes from the ωB97x-D/6-311+G(d,p) theory are summarized in Table 1. The isolated cationic and neutral methyl viologen guests distinguished based on the number of electrons; the MV2+ has positive charged on two quaternary nitrogen atoms. Optimization of free MV0 reveal classical quinoid type structure where both intra-ring C1-C2 and inter-ring C3-C3’ bonds have double bond character with the pyridine rings being coplanar (dihedral ∠C2C3-C3’-C2’=180" ) to each other. The methyl (CH3) functionality shows deviation of 15+ from the molecular skeleton which suggests the nitrogen lone pair is not fully involved in π-conjugation. Since the hybridization of nitrogen centers are modified from sp2 to sp3 during reduction process (accommodation of two electron) which subsequently increase the pyramidal feature of nitrogen atom. On the other hand the central C3-C3’ bond joining bipyridinum rings turns to be 1.481 Å which is consistent with the single bond character between sp2 carbon atoms. Here two non-planar pyridine moieties deviate by 42" from the planarity as also observed in the earlier experiments. Encapsulation of viologen within the tubular TA4 cavity brings about structural perturbations and reflect in the structural parameters of individual host and guest. The double bonds in conjugated framework of encapsulated MV0 are elongated up to 0.003 Å contrary to the alternate C1N1 and C2-C3 single bonds which are shortened; the concomitant loss of p character relative to its free analogue was thus inferred. As may be noticed the two pyridine rings in the MV2+@TA4 complex tend to attain planarity (deviation up to 22" ). The penetration of

10 ACS Paragon Plus Environment

Page 10 of 33

Page 11 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

MV2+ within the TA4 receptor brings about contraction of its cavity which explains its stronger binding than MV0 guest. In an attempt to shed some light on the supramolecular interactions between hostguest complexes Bader69 QTAIM methodology was used. Numerical data of the topological parameters viz., electron density (ρ), Laplacian of electron density (∇ ), kinetic and potential electron energy densities (G and V), and ellipticity parameters (.) for supramolecular interactions are collected in Table 3. The multiple C-H···π interactions are evident from the presence of the (3,-1) bond critical point (BCP) in between the bond path linking the methyl viologen hydrogen atoms with one or more carbon atoms of naphthalene moieties of the host, with the ρ being in the range 0.0049-0.0092 au for both complexes. The resultant ∇2ρ parameters for the same are positive evincing closed-shell type interactions. An appearance of BCP (ρ being in 0.0058 to 0.0074 au) for intermolecular C···C bond path between the viologen and naphthalene rings indicating attractive through-space interactions between the stacked π−π* orbitals in the complexes. Additionally the C-H···O, O-H···O contact are ascertained from the ρ values reported in the table.

Besides, the presence of intramolecular hydrogen-hydrogen

bonding interactions between the ortho-hydrogen atoms from adjacent pyridines in the MV0@TA4 complex explain its planar structure. The additional cation···/ interactions render large stability to the MV2+@TA4 complex. It has been observed here the total electron energy density (H) composed of counterbalancing kinetic (G > 0) and potential (V < 0) electron energy densities contributions. Furthermore, the ellipticity parameters (Ɛ=

0

0

− 1) at BCP were calculated which turn out be larger for C-H···π interactions owing

to the π-electrons serve as proton acceptors. NCI-RDG methods are utilized to capture details of the weak noncovalent interactions. Figure 3 displays graphical representation of the two dimensional (2D) 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

stiletto-heel shape scatter plots (s verses sign (1 )) and 3D spatial visualization of NCI isosurfaces and for the MV0@TA4 and MV2+@TA4 complexes. The sign (1 ) parameters are color-mapped onto the s-isosurfaces with two density-cutoffs being # and 2 and utilized to set the color scale.69 A red-green-blue scale color coding scheme has widely been preferred with the red representing  2 (repulsive) and blue for # (attractive). A direct visualization of 3D plot of MV0@TA4 complex impart the hydrogen bonding interactions which appear as round, localized, greenish blue NCI domain, in the region where hydrogen and oxygen atoms interact with each other and in the 2D plot shows a peak at sign (1 ) ≈ −0.03 au. The rest of weak intermolecular interactions for instance, C-H···π ones are represented by wider band around zero density (due to annihilation of the density gradient) region in 2D plot while a brownish green extended sheet like isosurfaces was observed in the 3D plots. The high ρ critical value (≈ 0.022 au) domain corresponds to destabilized steric clashes in closed shell ring interactions which emerge with the bright red cigar-shaped surfaces elongated along the direction of decreasing density up in the central region of aromatic moiety. Moreover the almond-shaped bicolored isosurfaces manifests between the hydrogen atoms of bipyridinium rings indicating dihydrogen interactions. The attractive component observed for these interactions are counterbalanced by destabilization arising from the induced ring closure. We find that the 2D scatter plot of cationic complex is not qualitatively different from the neutral complex. As demonstrated here the weak interactions between the host and guest shed light on unmistakable compatibility between QTAIM and NCI-RDG. To gain in-depth insights into the host-guest interactions with an emphasis on electronic spectra at molecular level, the time-dependent DFT (TDDFT) was adopted. Most representative molecular frontier orbitals are illustrated in Figure 6. The isolated MV0 and MV2+ guests show strong absorption bands near 345 nm (f= 1.00) and 279 nm

12 ACS Paragon Plus Environment

Page 12 of 33

Page 13 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(f=0.72) respectively, assigned to HOMO → LUMO transition whereas the TA4 scaffold reveal the corresponding band near ~312 nm. A new intermolecular charge transfer band was noticed in the visible region (625 nm) that arises from π (HOMO) → π∗ (LUMO) transition was observed for the MV0@TA4 complex. Noteworthy enough, such weak absorption band serve as an indicator of the host-guest binding. From the visual inspection of the frontier orbitals of the complex it become discernible that the interactions between MV0 and TA4 are of typical donor (guest) –acceptor (host) interactions. These inferences are significant and suggest the formation of complex is driven by intermolecular charge transfer interactions. A qualitatively different picture emerges for the charged (MV2+@TA4) complex in which the HOMO resides over the πorbitals of the host with polarization being induced from the guest. On the other hand the LUMO of the complex is rendered with large antibonding π-character and resides exclusively over the guest. Distinct binding behavior of cationic and neutral species toward molecular receptor can be monitored by the 1H NMR experiments. The GIAO derived δH values of isolated guest, molecular cavitand and their complexes were computed with the use of the dichloromethane as solvent which are compared in Table 3. The methyl viologen protons are broadly classified as: pyridine protons (H1 and H2) and methylene protons (H3) as displayed in Figure 3. The hierarchy of the calculated δH values turns out to be: H1> H2 > H3 for the isolated oxidized and reduced forms of methyl viologen. In comparison of MV2+ protons viz., H1 (9.94ppm), H2 (9.35 ppm) and H3 (5.49 ppm) are more de-shielded than those of MV0 guest. Upon complexation the protons exhibit upfield signals implying complete embedment of the guest within the host cavity. The host protons are nearly insensitive to complexation (cf. Table 3).

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Encapsulation of the guest within the host cavity bring about the structural changes which reflect in the frequency shifts of their characteristic normal vibrations. The vibrational spectra of the free host, guest and their complex from theory in the form of the molar absorption coefficient (or molar absorptivity in units of 0.1 m2 mol−1) plotted versus the frequency (in cm−1) of neutral and charged complexes are illustrated in Figure S1 of the Supporting Information. The isolated MV0 shows a band near 3240 cm-1 arises from the –CH stretches of pyridine rings participating in C-H···π interactions. Upon complexation this vibration engenders the frequency up-shift to ~3250 cm-1 (blue shifted hydrogen bond) with a diminutive intensity. The C···C contacts between viologen and naphthalene moieties show signature in the C=C stretching of the MV0@TA4 complex and accordingly the wavenumber lowering of ~10 cm-1 is predicted. These arguments further can be extended to MV2+@TA4 complex as well. Moreover the direction of frequency shifts and intensity patterns of normal vibrations in their infrared spectra serve as a “fingerprint” for characterizing molecular interactions underlying the host-guest binding.

Conclusions In the present venture we report systematic analysis of hydroxylated oxatub[4]arene receptor employing the ωB97x-D based density functional theory. The by flipping of naphthalene monomers in the macrocyle yields zigzag, disrupted zigzag, 1,2alternate and 1,3-alternate conformers, which possess well-defined tubular shape cavities with varying dimensions. The oxidized and reduced forms of methyl viologen favors encapsulation within the 1,2-alternate conformer which is attributed to its perfect complementary size and shape, and stabilized via multiple non-covalent interactions of viologen protons with naphthalene rings and ether linked oxygens as well. It has been shown that MV2+ binds strongly to TA4 receptor which has been explained from the additional cation···π interactions. The MV0@TA[4] represents the donor-acceptor

14 ACS Paragon Plus Environment

Page 14 of 33

Page 15 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

complex and accompanied by the intermolecular charge transfer interactions which manifests as the ~625 nm band in its electronic spectra. An encapsulation of the guest within host cavity further shows up-field perturbations of viologen (guest) protons in the 1H

NMR spectra. Besides the direction of shifts of -CH and -C=C characteristic stretchings

in the calculated vibrational spectra signify the formation of inclusion complex.

In conspectus, the molecular recognition of the methyl viologen guest by the oxatub[4]arene synthetic receptor discussed here should give impetus for design and modeling of supramolecular assemblies based on the naphthalene based receptors. Acknowledgements SPG acknowledges support from the Research Project (37(2)/14/11/2015-BRNS) from the Board of Research in Nuclear Sciences (BRNS), India. DNL is thankful to Savitribai Phule Pune University for the award of research fellowship through the University of Potential excellence scheme from the University Grants Commission, New Delhi, India. Authors thank the Center for Development of Advanced Computing (CDAC), Pune for providing National Param Supercomputing Facility.

References 1. Gloe, K. Macrocyclic Chemistry: Current Trends and Future Perspectives, ed. Springer, Dordrecht, 2005. 2. Davis, F.; Higson, S. Macrocycles: Construction, Chemistry and Nanotechnology Applications, Wiley, Chichester, 2011. 3. Lehn, J.-M. Supramolecular Chemistry: Concepts and Perspectives, Wiley-VCH, Weinheim, 1995. 4. Steed, J. W.; Atwood, J. L. Supramolecular Chemistry, Wiley, Chichester, 2nd edn, 2009.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

5. Jiang, J.; Meng, Y.; Zhang, L.; Liu, M. H. Self-Assembled Single-Walled Metal-Helical Nanotube (M-HN): Creation of Efficient Supramolecular Catalysts for Asymmetric Reaction. J. Am. Chem. Soc., 2016, 138, 15629–15635. 6. Mandl, C.; König, B. Chemistry in Motion—Unidirectional Rotating Molecular Motors. Angew. Chem. Int. Ed., 2004, 43, 1622–1624. 7. Li, C. Pillararene-Based Supramolecular Polymers: from Molecular Recognition to Polymeric Aggregates. Chem. Commun., 2014, 50, 12420–12433. 8. Sun, T.; Zhang, Y. S.; Pang, B.; Hyun, D. C.; Yang, M.; Xia, Y. Engineered nanoparticles for drug delivery in cancer therapy. Angew. Chem. Int. Ed., 2014, 53, 12320-12364. 9. Oshovsky, G. V.; Reinhoudt D. N.; Verboom, W. Triple-ion interactions for the construction of supramolecular capsules. J. Am. Chem. Soc., 2006, 128, 5270-5278. 10. Gutsche, C. D. Calixarenes, Monographs. In Supramolecular Chemistry; Stoddart, J. F., Ed.; The Royal Society of Chemistry: Cambridge, 1989. 11. Hardie, M. Recent advances in the Chemistry of Cyclotriveratrylene. J. Chem. Soc. Rev., 2010, 39, 516-527. 12. Ogoshi, T.; Kanai, S.; Fujinami, S.; Yamagishi T.; Nakamoto, Y. Para-Bridged Symmetrical Pillar[5]arenes: their Lewis Acid Catalyzed Synthesis and Host–guest Property. J. Am. Chem. Soc., 2008, 130, 5022-5023. 13. Schmidtchen, F. P. Surprises in the Energetics of Host− guest Anion Binding to Calix [4]pyrrole. Org. Lett., 2002, 4, 431-434. 14. Wang, M.-X.; Zhang X.-H.; Zheng, Q.-Y. Synthesis, Structure, and [60] Fullerene Complexation Properties of Azacalix[m]arene[n]pyridines. Angew. Chem. Int. Ed., 2004, 43, 838-842.

16 ACS Paragon Plus Environment

Page 16 of 33

Page 17 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

15. Chun, Y.; Singh, N. J.; Hwang, I.-C.; Lee, J. W.; Yu S. U.; Kim, K. S. Calix[n]imidazolium as a New Class of Positively Charged Homo-calix Compounds. Nat. Commun. , 2013, 4, 1797. 16. Ma, J.; Meng, Q.; Hu, X.; Li, B.; Ma, S.; Hu, B.; Li, J.; Jia, X.; Li, C. Synthesis of a Watersoluble

Carboxylatobiphen[4]arene

and

its

Selective

Complexation

toward

Acetylcholine. Org. lett., 2016, 18, 5740-5743. 17. Poh, B.-L.; Lim C. S.; Khoo, K. S. A Water-soluble Cyclic Tetramer from Reacting Chromotropic acid with Formaldehyde. Tetrahedron Lett., 1989, 30, 1005-1008. 18. Poh,

B.-L.;

Lim,

C.

S.

Complexations

of

Amines

with

Water-soluble

Cyclotetrachromotropylene. Tetrahedron, 1990, 46, 3651-3658. 19. Chen, H.; Fan, J.; Hu, X.; Ma, J.; Wang, S.; Li, J.; Yu, Y.; Jia, X.; Li, C. Biphen[n]arenes. Chem. Sci., 2015, 6, 197–202. 20. Georghiou, P. E.; Valluru,

G.; Schneider, C.; Liang,

S.; Woolridge, K.; Mulla, K.;

Adronovc, A.; Zhao, Y. Dispersion of Single-Walled Carbon Nanotubes into Aqueous Solutions using Poh's Cyclotetrachromo-Tropylene (CTCT). RSC Adv., 2014, 4, 3161431617. 21. Shorthill, B. J.; Avetta, C. T.; Glass, T. E. Shape-Selective Sensing of Lipids in Aqueous Solution by a Designed Fluorescent Molecular Tube. J. Am. Chem. Soc., 2004, 126, 12732–12733. 22. Avetta, C. T.; Shorthill, B. J.; Ren, C.; Glass, T. E. Molecular Tubes for Lipid Sensing: Tube Conformations Control Analyte Selectivity and Fluorescent Response. J. Org. Chem., 2012, 77, 851–857. 23. AlHujran, T. A.; Dawe L. N.; Georghiou, P. E. Synthesis of Functionalized Acenaphthenes and a New Class of Homooxacalixarenes. Org. Lett., 2012, 14, 3530– 3533. 17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

24. Andreetti, G. D.; Boehmer, V.; Jordon, J. G.; Tabatabai, M.; Ugozzoli, F.; Vogt W.; Wolff, A. Dissymmetric calix [4] arenes with C4-and C2-symmetry. Synthesis, X-ray Structures, Conformational Fixation, and Proton NMR Spectroscopic Studies. J. Org. Chem., 1993, 58, 4023–4032. 25. Jia, F.; He, Z.; Yang, L.-P.; Pan, Z.-S.; Yi, M.; Jiang, R.-W.; Jiang, W. Oxatub[4]arene: A Smart Macrocyclic Receptor with Multiple Interconvertible Cavities. Chem. Sci. 2015, 6, 6731–6738. 26. Jia, F.; Wang, H.-Y.; Li, D.-H.; Yang, L.-P.; Jiang, W. Oxatub[4]arene: A Molecular “transformer” Capable of Hosting a Wide Range of Organic Cations. Chem. Commun. 2016, 52, 5666-5669. 27. Jia, F.; Li, D.-H.; Yang, T.-L.; Yang, L.-P.; Dang, L.; Jiang, W. Oxatub[5,6]arene: Synthesis, Conformational Analysis, and the Recognition of C60 and C70. Chem. Commun., 2017, 53, 336–339. 28. Jia, F.; Yang, L. P.; Li, D. -H.; Jiang, W., Electronic Substituent Effects of Guests on the Conformational Network and Binding Behavior of Oxatub [4] arene. J. Org. Chem., 2017, 82, 10444-10449. 29. Yang, L.-P.; Jia, F.; Zhou, Q.-H.; Pan, F.; Sun, J.-N.; Rissanen, K.; Chung, L. W.; Jiang, W. Guest-Induced Folding and Self-Assembly of Conformationally Adaptive Macrocycles into Nanosheets and Nanotubes. Chem. Eur. J., 2017, 23, 1516–1520. 30. Huang, G. B.; Wang, S. H.; Ke, H.; Yang, L.P.; Jiang, W. Selective Recognition of Highly Hydrophilic Molecules in Water by Endo-Functionalized Molecular Tubes. J. Am. Chem. Soc., 2016, 138, 14550-14553. 31. Wang, L. L.; Chen, Z.; Liu, W.E.; Ke, H.; Wang, S.H.; Jiang, W. Molecular Recognition and Chirality Sensing of Epoxides in Water Using Endo-Functionalized Molecular Tubes. J. Am. Chem. Soc., 2017, 139, 8436–8439.

18 ACS Paragon Plus Environment

Page 18 of 33

Page 19 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

32. Ma, Y. L.; Ke, H.; Valkonen, A.; Rissanen, K. ; Jiang, W. Achieving Strong Positive Cooperativity through Activating Weak Non-Covalent Interactions. Angew. Chem. Int. Ed. 2017, 56, 1 – 6. 33. He, Z.; Ye, G.; Jiang, W. Imine Macrocycle with a Deep Cavity: Guest-Selected Formation of Syn/anti Configuration and Guest-Controlled Reconfiguration. Chem. Eur. J., 2015, 21, 3005–3012. 34. He, Z.; Yang, X.; Jiang, W. Synthesis, Solid-State Structures, and Molecular Recognition of Chiral Molecular Tweezer and Related Structures Based on a Rigid Bis-Naphthalene Cleft. Org. Lett., 2015, 17, 3880–3883. 35. Huang, G.; He, Z.; Cai, C.-X.; Pan, F.; Yang, D.; Rissanen, K.; Jiang, W. Bis-Urea Macrocycles with a Deep Cavity. Chem. Commun. 2015, 51, 15490–15493. 36. Huang, G.; Valkonen, A.; Rissanen, K.; Jiang, W. Endo-Functionalized Molecular Tubes: Selective Encapsulation of Neutral Molecules in Non-Polar Media. Chem. Commun., 2016, 52, 9078-9081. 37. Shewale, M. N.; Lande, D. N.; Gejji, S. P. Density Functional Investigations on the Selective Binding of an Endo-Functionalized Bis-Urea Macrocycle. J. Phys. Chem. A, 2017, 121, 288–297. 38. Lande, D. N.; Shewale, M. N.; Gejji, S. P. Host-guest Interactions Accompanying the Encapsulation

of

1,4-diazabicyclo[2.2.2]octane

within

endo-Functionalized

Macrocycles. J. Phys. Chem. A, 2017, 121, 288-297. 39. Monk, P. M. The Viologens, John Wiley & Sons, Chichester (UK), 1998. 40. Bird, C. L. Electrochemistry of the Viologens. Chem. Soc. Rev., 1981, 10, 49-82. 41. Sliwa, W.; Bachowska, B.; Zelichowicz, N. Chemistry of Viologens. Heterocycles, 1991, 32, 2241-2273.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

42. Kim, S. M.; Jang, J. H.; Kim, K. K.; Park, H. K.; Bae, J. J.; Yu, W. J.; Lee, I. H.; Kim, G.; Loc, D. D.; Kim, U. J.; Lee, E. H.; Shin, H. J.; Choi J. Y.; Lee, Y. H. Reduction-controlled Viologen in Bisolvent as an Environmentally stable n-type Dopant for Carbon Nanotubes. J. Am. Chem. Soc., 2009, 131, 327-331. 43. Moller, M.; Asaei, S.; Corr, D.; Ryan, M.; Walder, L., Switchable Electrochromic images based on a Combined top–down bottom–up Approach. Adv. Mater., 2004, 16, 1558– 1562. 44. Tanabe, K.; Yasuda, T.; Yoshio, M.; Kato, T. Viologen-based Redox-active Ionic liquid Crystals forming Columnar Phases. Org. Lett., 2007, 9, 4271-4274. 45. Causin V.; Saielli, G. Effect of Asymmetric Substitution on the Mesomorphic Behaviour of low-melting Viologen Salts of Bis (trifluoromethanesulfonyl) amide. J. Mater. Chem., 2009, 19, 9153-9162. 46. Causin, V.; Saielli, G. Effect of a Structural Modification of the Bipyridinium Core on the Phase Behaviour Of Viologen-based Bistriflimide Salts. J. Mol. Liq., 2009, 145, 4147. 47. Bhowmik, P. K.; Han, H. S.; Cebe, J. J.; Burchett, R. A.; Acharya, B.; Kumar, S. Ambient Temperature Thermotropic Liquid Crystalline Viologen Bis(triflimide) Salts. Liq. Cryst., 2003, 30, 1433-1440. 48. Kay, E. R.; Leigh D. A.; Zerbetto, F. Synthetic Molecular Motors and Mechanical Machines. Angew. Chem., Int. Ed., 2007, 46, 72-191. 49. Arduini, A.; Bussolati, R.; Credi, A.; Pochini, A.; Secchi, A.; Silvi, S.; Venturi, M. Rotaxanes with a Calix[6]Arene Wheel and Axles of different length. Synthesis, Characterization, and Photophysical and Electrochemical Properties. Tetrahedron, 2008, 64, 82798286.

20 ACS Paragon Plus Environment

Page 20 of 33

Page 21 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

50. Arduini, A.; Ciesa, F.; Fragassi, M.; Pochini, A.; Secchi, A. Selective Synthesis of Two Constitutionally Isomeric Oriented Calix [6] arene-Based Rotaxanes. Angew. Chem., Int. Ed., 2005, 44, 278-281. 51. Credi, A.; Dumas, S.; Silvi, S.; Venturi, M.; Arduini, A.; Pochini, A.; Secchi, A. Viologencalix[6]arene Pseudorotaxanes. Ion-pair Recognition and threading/dethreading molecular motions. J. Org. Chem., 2004, 69, 5881-5887. 52. Kamata, K.; Kawai, T.; Iyoda, Anion-controlled Redox Process in a Cross-linked Polyviologen film toward Electrochemical Anion Recognition. T. Langmuir, 2001, 17, 155-163. 53. Saielli, G. Ion-pairing of Octyl Viologen Diiodide in low-polar Solvents: An Experimental and Computational Study. J. Phys. Chem. A, 2008, 112, 7987-7995. 54. Brun, A. M.; Harriman, A. Photochemistry of Intercalated Quaternary Diazaaromatic Salts. J. Am. Chem. Soc., 1991, 113, 8153–8159. 55. Clennan, E. L.; Viologen Embedded Zeolites. Coord. Chem. Rev., 2004, 248, 477–492. 56. Mirzoian, A.; Kaifer, A. E. Reactive Pseudorotaxanes: Inclusion Complexation of Reduced

Viologens

by

the

Hosts

β-Cyclodextrin

and

Heptakis

(2,

6-di-o-Methyl)-β-Cyclodextrin. Chem. – A Eur. J., 1997, 3, 1052–1058. 57. Xiao, X.; Tao, Z.; Xue, S.-F.; Zhu, Q.-J.; Zhang, J.-X.; Lawrance, G. A.; Raguse B.; Wei, G. Interaction between Cucurbit[8]uril and Viologen Derivatives. J. Incl. Phenom. Macrocycl. Chem., 2008, 61, 131–138. 58. Guo, D.-S.; Wang, L.-H.; Liu, Y. Highly Effective Binding of Methyl Viologen Dication and its Radical Cation by p-sulfonatocalix[4,5]arenes. J. Org. Chem., 2007, 72, 7775–7778. 59. Li, C.; Xu, Q.; Li, J.; Yao, F.; Jia, X. Complex Interactions of Pillar[5]arene with Paraquats and bis (pyridinium) Derivatives. Org. Biomol. Chem., 2010, 8, 1568–1576.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

60. 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. et al. Gaussian, Inc., Wallingford CT,2009. 61. Dennington, R.; Keith, T.; Milliam, J. Semichem Inc., Shawnee Mission, KS, 2009. 62. Chai, J.-D.; Head-Gordon; M. Long-range Corrected Hybrid Density Functionals with damped atom–atom Dispersion Corrections, Phys. Chem. Chem. Phys., 2008, 10, 66156620. 63. Remya, K.; Suresh, C. H. Which Density Functional is Close to CCSD Accuracy to Describe Geometry and Interaction Energy of Small Noncovalent Dimers? A Benchmark Study Using Gaussian 09. J. Comput. Chem., 2013, 34, 1341–1353. 64. Liu, Y.; Zhao, J.; Li, F.; Chen, Z. Appropriate Description of Intermolecular Interactions in the Methane Hydrates: An Assessment of DFT Methods. J. Comput. Chem., 2013, 2, 121–131. 65. Rao, S. S.; Lande, D. N.; Gejji, S. P. Density Functional Theory Investigations on Binding and Spectral Features of Complexes of Ferrocenyl Derivatives with Cucurbit [7]Uril. J. Mol. Liq., 2016, 216, 298–308. 66. Das, R.; Chattaraj. P. K. Host–Guest Interactions in ExBox4+. ChemPhysChem., 2014, 15, 4108-4116. 67. McLean, A. D.; Chandler, G. S. Contracted Gaussian-basis sets for molecular calculations. 1. 2nd row atoms, Z=11-18, J. Chem. Phys., 1980, 72, 5639-48. 68. Murray, J. S.; Shields, Z. P.; Seybold, P. G.; Politzer, P. Intuitive and Counterintuitive Noncovalent Interactions of Aromatic π Regions with the Hydrogen and the Nitrogen of HCN. J. Comp. Sci., 2015, 10, 209–216 69. Politzer, P.; Murray, J. S.; Clark, T. Mathematical Modeling and Physical Reality in Noncovalent Interactions. J. Mol. Model., 2015, 21, 52.

22 ACS Paragon Plus Environment

Page 22 of 33

Page 23 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

70. Murry, J. S.; Seminario, J. M.; Politzer. P. A Computational Study of the Structures and Electrostatic Potentials of Some Azines and Nitroazines. J. Mol. Struct., (TheoChem) 1989, 87, 95-108. 71. Bader, R. F. W., Atoms in Molecule Oxford Science Publication: Oxford U. K. 1990. 72. AIMAll (Version 14.11.23), Todd A. Keith, TK Gristmill Software, Overland Park KS, USA, 2014 (aim.tkgristmill.com). 73. Johnson, E. R.; Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A. J.; Yang, W. Revealing Noncovalent Interactions. J. Am. Chem. Soc., 2010, 132, 6498–6506. 74. Lu, T.; Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem., 2012, 33 (5), 580–592. 75. Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics. J. Mol. Graphics, 1996, 14, 33−38. 76. Wolinski, K.; Hinton, J. F.; Pulay, P. Efficient Implementation of the Gauge-Independent Atomic Orbital Method for NMR Chemical Shift Calculations. J. Am. Chem. Soc., 1990, 112, 8251–8260. 77. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999−3093. 78. Bader, R. F. W.; Carroll, M.T.; Cheeseman, J. R.; Chang, C. Properties of Atoms in Molecules: Atomic Volumes. J. Am. Chem. Soc., 1987, 109, 7968–7979. 79. Murray, J. S.; Politzer, P. "Molecular surfaces, van der Waals radii and Electrostatic Potentials in Relation to Noncovalent Interactions. Croat. Chem. Acta., 2009, 82, 267– 275.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 33

Figures:

Top view

Side view

Zigzag (I)

Disrupted zigzag (II)

1,2-alternate (III)

1,3-alternate (IV)

Figure 1. Different conformers of parent hydroxylated oxatub[4]arene (TA4). 24 ACS Paragon Plus Environment

Page 25 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Zigzag (I)

Disrupted zigzag (II)

1,2-alternate (III)

1,3-alternate (IV)

Figure 2. Electron density isosurfaces (0.001 au) overlaid with MESP (from +0.06 au to -0.06 au) in TA4 conformers.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Reduced methylviologen (MV0)

Oxidized methylviologen (MV2+) Figure 3. Optimized geometries and electron density mapped isosurfaces of methylviologen guests.

26 ACS Paragon Plus Environment

Page 26 of 33

Page 27 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Top view

Side View

MV0@TA4

MV2+@TA4 Figure 4. Energy minimal optimized geometries of TA4 complexes.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

MV0@TA4

MV2+@TA4

Figure 5. Two dimensional (2D) stiletto-heel shape scatter plots (s verses sign (1 )) and three dimensional (3D) spatial visualization of NCI isosurfaces and for MV0@TA4 and MV2+@TA4 complexes.

28 ACS Paragon Plus Environment

Page 28 of 33

Page 29 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

HOMO

LUMO

TA4

MV0

MV0@TA4

MV2+

MV2+@TA4

Figure 6. Most representative frontier molecular orbitals.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 33

Table 1. Selected structural parameters (bond-distances (R) in Å and dihedral angles (θ) in °) of MV0, MV2+ and their inclusion complexes. MV0

MV2+

MV0@TA4

MV2+@TA4

R(N-C1)

1.386

1.346

1.379

1.344

R(C1-C2)

1.342

1.375

1.343

1.375

R(C2-C3)

1.462

1.395

1.458

1.396

R(C3-C3’)

1.367

1.481

1.367

1.485

R(C4-N)

1.445

1.476

1.448

1.475

R(N-N)

7.124

7.011

7.081

6.973

0

42

0

22

165

179

176

179

θ(C2-C3-C3’-C2) θ (C1-C2-N-C4)

30 ACS Paragon Plus Environment

Page 31 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 2. QTAIM parameters corresponding to noncovalent interactions present in TA4 inclusion complexes.

Complex MV0@TA4

MV2+@TA4

Non-covalent Interactions C2-H2···O C2’-H2’···O C1-H1···π C2-H2···π C1’-H1’···π C2’-H2’···π H2-H2’ π ···π C2-H2···O C2’-H2’···O C1-H1···π C2-H2···π C1’-H1’···π C2’-H2’···π π ···π cation···π

ρ

∇2ρ

G

V

Ɛ

0.0084 0.0091 0.0075 0.0092 0.0079 0.0070 0.0121 0.0074 0.0144 0.0073 0.0050 0.0077 0.0065 0.0088 0.0062 0.0061

0.0281 0.0280 0.0215 0.0282 0.0230 0.0222 0.0461 0.0213 0.0478 0.0233 0.0145 0.0214 0.0184 0.0270 0.0180 0.0235

0.0061 0.0062 0.0046 0.0060 0.0048 0.0047 0.0093 0.0044 0.0103 0.0050 0.0030 0.0044 0.0038 0.0056 0.0036 0.0049

-0.0053 -0.0055 -0.0037 -0.0049 -0.0039 -0.0037 -0.0071 -0.0032 -0.0087 -0.0042 -0.0025 -0.0036 -0.0030 -0.0046 -0.0027 -0.0041

1.1105 0.0681 0.2803 0.5106 2.2631 6.0888 0.8561 0.9668 0.1046 0.0239 0.9793 1.2646 1.4249 2.3528 0.0824 0.0049

31 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 33

Table 3. 1H NMR chemical shifts (ppm) of host, guests and their complexes. δH

TA4

Ha Hb Hc Hd H1 H2 H3

8.08 9.32 5.19 6.24

MV0

6.85 6.45 3.69

MV2+

MV0@TA4

MV2+@TA4

9.94 9.35 5.49

8.12 9.35 5.38 6.27 4.50 3.63 2.79

8.09 9.24 5.53 6.33 6.53 7.71 4.64

32 ACS Paragon Plus Environment

Page 33 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

325x117mm (96 x 96 DPI)

ACS Paragon Plus Environment