Harshita Kumari*a,
Wei G. Wycoffb,
Collin M. Mayhanb,
Steven R. Klinec,
Joshua R. Sob,
Carol A. Deakyne*b,
John E. Adams*b and
Jerry L. Atwood*b
aJames L. Winkle College of Pharmacy, University of Cincinnati, OH 45267, USA. E-mail: kumariha@ucmail.uc.edu
bDepartment of Chemistry, University of Missouri-Columbia, Columbia, MO 65211, USA. E-mail: AtwoodJ@umsystem.edu; DeakyneC@umsystem.edu; AdamsJE@umsystem.edu
cNIST Center for Neutron Research, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, MD 20899, USA
First published on 15th January 2021
The structural stability and solution geometry of zinc-seamed-C-propylpyrogallol[4]arene dimers has been studied in solution using in situ neutron scattering and 2D-DOSY NMR methods. In comparison with the structures of the analogous copper-/nickel-seamed dimeric entities, the spherical geometry of the PgC3Zn species (R = 9.4 Å; diffusion coefficient = 1.05 × 10−10 m2 s−1) is larger due to the presence of ligands at the periphery in solution. This enhanced radius in solution due to ligation is also consistent with the findings of model molecular dynamics simulations of the zinc-seamed dimers.
One group of such assemblies comprises the metal-seamed pyrogallol[4]arene nanocapsules (Fig. 1).1,7 Although these nanocapsules have now been synthesized with a selection of alkyl groups, metals, guests and axial ligands, only a few have been examined in solution.8–15 Our motivation in studying solution-phase architectures is to determine whether the solution structures conform to the solid-state geometries. Many nanocapsules, specifically transition metal (Cu/Ni)-seamed hexamers and dimers and hydrogen-bonded dimers, retain their solid-state geometry in solution.10,13 In contrast, Ga- and Ga/Zn-seamed pyrogallol[4]arene hexamers rearrange from spheres in the solid state to toroids in solution.12 Similarly, the hydrogen-bonded PgC1⊂ferrocene/PgC6⊂pyrene nanotubes rearrange to dimers in solution.14,15 Solution-phase studies also have allowed us to discover novel structures of species that were either difficult to crystallize or expected to be unstable in solution.9,11 For example, the ellipsoidal and tubular architectures of the PgC17Cu and PgC1Fe nanoassemblies, respectively, have been studied exclusively in solution.9,11
The current study and analysis focus on the solution structure of the original dimeric pyrogallol[4]arene-based nanocapsules: [Zn8(C-propylpyrogallol[4]arene)2(pyridine)7(DMSO)1⊂pyridine], with propyl R-groups, a pyridine guest and pyridine/DMSO axial ligands (Fig. 1C and D, propyl R-groups and guest pyridine molecule are removed for clarity).16 The zinc dimer is composed of two pyrogallol[4]arene subunits that are coordinated by the ring of eight Zn+2 ions across the equatorial belt. The metal atoms displace 16 of the 24 hydroxyl protons of the pyrogallols to seam this spherical framework. Each of the zinc centers is ligated by four phenoxy groups, two from each pyrogallol, along with an axial DMSO or pyridine ligand, that gives a square pyramidal configuration to the metal center. The central phenoxy group from each pyrogallol is bridged between two adjacent zinc centers, whereas the two exterior phenoxy groups are each ligated to a single zinc center and hydrogen-bonded to a phenoxy group on the neighboring pyrogallol subunit, making the overall assembly neutral. In the transition from pentacoordinated Zn in the solid state to tetracoordinated Zn in the gas phase, as a result of the MALDI-TOF analysis, this local Zn–O binding geometry is essentially preserved.16–18 Of course, thermal fluctuations do occur, although, as we shall show via molecular dynamics simulations, the overall shape remains roughly spherical.
The PgC3Zn dimeric nanocapsule was synthesized by adding 1 equiv. of PgC3 (0.01 mol L−1), 4 equiv. of zinc(II) nitrate (0.04 mol L−1), and 14 equiv. of pyridine in DMSO solvent. The solution was heated, sonicated, and left overnight to obtain yellow crystals of the dimer. The unit cell was confirmed and the crystals were filtered and dissolved in d6-DMSO at mass fractions of 1% and 5% for SANS measurements. Preferential deuteration of solvent for SANS studies enhances the coherent scattering and improves the contrast between the solute and the solvent. The samples were left overnight to ensure saturation without precipitation, and the SANS measurement was conducted at 25 °C on the NG3 30m SANS instrument at the NIST Center for Neutron Research in Gaithersburg, MD.19 The collected scattering data was reduced and analyzed using IgorPro20 incorporating the instrumental q-resolution function.
The measurements on the dimeric sample at mass fractions of both 1% and 5% were conducted to ensure sufficient scattering and to detect a concentration dependence of the capsular structure. The scattering length densities (SLDs) of the solute and solvent were calculated and held fixed for data analyses. The SLDs of the solvent (d6-DMSO) and solute (PgC3Zn dimer) were fixed at 5.28 × 10−6 Å−2 and 1.37 × 10−6 Å−2, respectively.
The reduced SANS data was fitted globally, with the two data sets corresponding to different mass fractions being modelled simultaneously to obtain the best fit. The data was fitted to various cylindrical, ellipsoidal and core–shell models to detect the occurrence of structural transformation. It was quickly found, however, that the entity does not have a core–shell, ellipsoidal or cylindrical geometry (ESI†). Previous studies on copper and nickel nanocapsules have indicated the presence of both dimers and hexamers in solution;13 hence, it was important to investigate the possible presence of multiple species with similar or different architectures. For this purpose, the data was globally fitted to the bimodal Schulz sphere and Schulz sphere models (ESI†). The calculated SLDs for the hexamer (1.56 × 10−6 Å−2), dimer (1.37 × 10−6 Å−2) and solvent (5.28 × 10−6 Å−2) were held fixed for global fitting. An unphysical volume fraction was obtained for the second spherical species in the bimodal Schulz sphere global fit, clearly indicating that the PgC3Zn sample has a unimodal distribution. In fact, the data for the PgC3Zn was fit best to a polydisperse sphere (Schulz sphere) model21,22 with a chi-sqrt value of 1.156 and a radius of ≈9.4 Å, a result indicating the presence of only one spherical assembly in solution (ESI;† Fig. 1 and 2).
In comparison with the radii of the PgC3Cu and PgC3Ni dimers, the radius obtained for the PgC3Zn dimer is ≈2.4 Å larger, and the size/diameter of the PgC3Zn dimer (≈18.8 Å) is intermediate between that of a PgC3Cu/Ni dimer (≈14 Å) and that of a PgC3Cu/Ni hexamer (≈20 Å). This comparison indicates that PgC3Cu and PgC3Ni most likely lose their ligands in solution. As this was the case, we examined the dimensions of the solid-state PgC3Zn dimer (XRD). The solid-state radius of the PgC3Zn dimer is ≈9.4 Å and ≈7 Å with and without external pyridine/DMSO ligands, respectively.16 The geometric dimensions obtained for the PgC3Zn dimer in the solution phase are nearly identical to those obtained in the solid state. Thus, SANS data analysis for the solvent-solubilized PgC3Zn dimer not only confirms its spherical shape and structural stability in solution but also reveals the preference for a square pyramidal metal bonding geometry, which differs from the distorted square planar geometry found in other metal-seamed (PgC3Cu, PgC3Ni) dimers in solution.
These radii obtained for the PgC3Zn dimer are also commensurate with those calculated from molecular dynamics simulations of PgC0Zn dimers with and without 8 ligated molecules. These mixed quantum-classical calculations were performed using the AMBER program suite,23 with a tight-binding DFT scheme (DFTB3)24 being used to model the Zn, O, and bridging H atoms as well as the atoms of the ligand. All other atoms of the dimer were treated classically using the gaff2 force field.25 The simulation results reported here reflect the dynamics of room-temperature gas-phase species. In the absence of the ligands, the average dimer radius was found to be 4.9 ± 0.1 Å, whereas with 8 ligated ammonia molecules the average radius was 7.1 ± 0.1 Å, with 8 ligated pyridine molecules the average radius was found to be 9.8 ± 0.1 Å, and with 8 ligated DMSO molecules the average radius was found to be 8.4 ± 0.7 Å (the uncertainties here are standard deviations of the average radii as measured from the instantaneous center of mass of the 8 Zn atoms. That is, they reflect the time-dependent thermal fluctuations of the dimer capsule rather than uncertainties in the calculated forces and simulation methodology). The larger uncertainty associated with the calculated average radius for the DMSO-ligated species here derives from the nonlinear bonding at the oxygen and nonplanar bonding at the sulfur in DMSO and the rotation about the oxygen–sulfur bond. In the above ligated species, the calculated radius of the roughly spherical dimer cage itself also differs from that of the non-ligated dimer, being larger by about 0.1 Å in each case.
The solution-phase behavior of PgC4Zn was further investigated using diffusion-ordered NMR spectroscopy (DOSY). Specifically, we investigated the encapsulation of pyridine within zinc dimers.
Note that we are reporting here the first DOSY NMR study of the metal-seamed pyrogallol[4]arene dimeric framework. The diffusion NMR studies reported by Cohen and co-workers were on hydrogen-bonded pyrogallol[4]arene, octahydroxypyridine[4]arene26,27 and resorcin[4]arene nanocapsules in CDCl3.28 Atwood et al. have more recently reported DOSY NMR studies on PgC4Ga hexamers, PgC4GaZn toroids and hydrogen-bonded PgC6⊂pyrene nanoassemblies in d6-acetone.29
Capsules of PgC4Zn⊂pyridine for the NMR studies were synthesized by mixing methanolic solutions of C-butylpyrogallol[4]arene, zinc nitrate and pyridine in the ratio of 1:4:14. Sonication yielded a precipitate of zinc dimers. A 2% solution of PgC4Zn⊂pyridine was prepared in d6-DMSO for the NMR measurements. Notably, we used a different tail (butyl) here, versus the previously studied C-propyl tail, as the tail length has been shown to not impact the self-assembled metal-coordinated core architecture of the nanocapsule.
The diffusion coefficients of the phenoxy hydroxyls, alkyl moieties and bound pyridine are similar (Fig. 3), indicating the structural integrity of the PgC4Zn dimer in DMSO solvent. Overall, the macrocycle/capsule and bound guest pyridine have a diffusion coefficient of 1.04 × 10−10 m2 s−1. The integration of bound versus free pyridine is 1:8 suggesting that the external pyridine ligand has exchanged with DMSO. The diffusion coefficients of the “free” pyridine ligand, solvent DMSO, methanol and water are 6.79 × 10−10, 7.00 × 10−10, 8.71 × 10−10 and 9.08 × 10−10 m2 s−1, respectively. The diffusion coefficient of the bound pyridine is about 6.5 times lower than that of a free pyridine.
Fig. 3 The DOSY NMR spectra of the [PgC4Zn⊂(pyridine)] dimer. The Y-axis on the DOSY plot is displayed as logD. |
Overall, the combined neutron scattering, diffusion-ordered NMR spectroscopy and molecular dynamics results indicate (a) the structural stability of the dimer in solution; (b) guest encapsulation; (c) fast external ligand exchange with solvent, which can only be captured through SANS; and (d) an increase in the effective radius of the zinc dimer due to ligation at zinc centers showing a preference for square pyramidal bonding.
The robustness of dimers demonstrated in this study justifies the use of these frameworks for guest encapsulation. Furthermore, the presence of exterior ligands affords valuable insight into the nature of guest-ligand communication. Future studies thus will focus on studying these robust frameworks and on variations in their structural properties in response to varying the encapsulated guest.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ra10053f |
This journal is © The Royal Society of Chemistry 2021 |