- 1College of Chemistry and Molecular Engineering, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing, China
- 2ReadCrystal Technology Co., Jiangsu, China
The crystalline sponge method has shown to be a novel strategy for the structure determination of noncrystalline, oily, or trace amount of a compound. A target compound was absorbed and oriented orderly in the pregrown porous crystal for x-ray diffraction analysis. However, the diffusion in the micron-sized crystals is rather difficult. Lots of trial-and-error experiments are needed to optimize the guest-soaking process and to improve data quality. Nanocrystals are better in diffusion, yet it could not conduct a single crystal x-ray diffraction (SCXRD) analysis. Three-dimensional electron diffraction (3D-ED) is a powerful diffraction tool for the structure determination of small crystals. In this work, we successfully carried out the crystalline sponge method by 3D-ED technique using {(ZnI2)3-[2,4,6-tris(4-pyridyl)-1,3,5-triazine]2·x(guest)}n (1-Guest) porous complex nanocrystals. On account of the better diffuse ability of nanocrystals, the time needed for solvent exchange and guest soaking protocols are shortened 50-fold faster versus the original protocol. The crystal structure of the crystalline sponge incorporated with three different guests was fully resolved using a merged dataset. The structure model was identical to previously reported ones using x-ray, showing that the accuracy of the 3D-ED was comparable with SCXRD. The refinement results can also give the precise occupancy of the guest molecule soaked in the porous crystal. This work not only provides a new data collection strategy for crystalline sponge method but also demonstrates the potential of 3D-ED techniques to study host-guest interaction by solving the fine structure of porous material.
Introduction
A molecular structure of a substance determines its intrinsic properties, such as reactivity, stereochemistry selectivity, etc. Chemists have developed nuclear magnetic resonance spectroscopy (NMR) and mass spectroscopy (MS) to provide essential structural information, but there are still chances to assign the structure incorrectly. Single-crystal x-ray diffraction (SCXRD) is the standard method that provides a direct three-dimensional structure at the atomic level. However, the major limitation of SCXRD is that it requires the target compound to crystallize into a single crystal with sufficient quality and size (>5 μm × 5 μm × 5 μm).
To address the challenge of structure determination on noncrystalline samples, Fujita and coworkers reported the “crystalline sponge” (CS) method (Inokuma et al., 2013; Hoshino et al., 2016). The CSs are porous crystalline materials that can absorb molecules in the solution and enable them to orient orderly within the pore. By soaking the molecule of interest into the pregrown CS single crystal, SCXRD can be performed on the host–guest complex, and therefore, the structure of the molecule of interest and the host can be elucidate simultaneously. The most commonly used CS is [(ZnI2)3-(tpt)2·x(Guest)]n porous complex [noted as 1-Guest, = 2,4,6-tris(4-pyridyl)-1,3,5-triazine]. The CS method has successfully determined the structure of natural products (Urban et al., 2016), metabolites (Zigon et al., 2015), and pharmaceutical compounds (Taniguchi et al., 2020). This method can also determine the absolute configuration of the organic compound of interest (Sairenji et al., 2017; Yan et al., 2017). However, the CS method is facing two major drawbacks. The yield of the porous host single crystal suitable for the CS method is quite low. More importantly, the critical guest-soaking step still needs trial-and-error experiments to optimize (Hoshino et al., 2016). Unoptimized soaking conditions, such as soaking time, temperature, and evaporation rate, may lead to crystal cracking and lowering guest diffusion, thus decreasing the chance of success structure determination on the molecule of interest.
Three-dimensional electron diffraction (3D-ED) is an emerging technique that collects single-crystal diffraction data from a nanocrystal (Gemmi et al., 2019; Gemmi and Lanza, 2019). It has demonstrated its powerful ability on determining complex structures including zeolites (Su et al., 2014), metal–organic frameworks (MOFs) (Huang et al., 2021), covalent organic frameworks (COFs) (Huang et al., 2021), and proteins (Nannenga and Gonen, 2019; Nguyen and Gonen, 2020). Combining a continuous rotation collection protocol, a noise-free direct electron detector, and a cryogenic sample holder, a single 3D-ED dataset could be recorded in 3 min with minimized beam damage (Nannenga et al., 2014; Nannenga and Gonen, 2019). Although many unknown structures were studied by 3D-ED, the host–guest interaction studies by 3D-ED were less reported. Owing to the vacuum environment and beam damage, the guest molecules more easily escaped from the host crystal. Until now, only few studies reported the identification of guests in the porous materials (Wang et al., 2018), or ligand binding in proteins (Clark et al., 2021).
In this work, we demonstrate the crystalline sponge method using the 3D-ED technique with nanosized CS crystal. Using a merged dataset with high completeness and resolution, the position of the framework and the guest molecule can be resolved ab initio and have comparable accuracy with a traditional x-ray diffraction method.
Materials and Methods
Materials
2,4,6-Tris(4-pyridyl)-1,3,5-triazine was purchased from the Energy Chemical Co. ZnI2 and guaiazulene were purchased from the Shanghai Aladdin Biochemical Technology Co., Ltd. All solvents used were purchased from Energy Chemical Co.
Host Synthesis and Guest Inclusion Procedures
The [(ZnI2)3-(tpt)2·x(Guest)]n porous complex [tpt = 2,4,6-tris(4-pyridyl)-1,3,5-triazine] (donated as 1-Guest) host was synthesized following the protocol of Fujita (Hoshino et al., 2016). The crystalline sponge was prepared by layering a methanol solution of ZnI2 (14.4 mg in 1.5 ml) on a nitrobenzene solution of (9.5 mg in 6 ml). The solution was left for 7 days for crystallization. The powder-like 1-Nitrobezene at the bottom of the bottle was transferred to a new glass vial. It was then washed 2–3 times cyclohexane using the same volume of the synthesis liquid. The powder-like crystals were heated in cyclohexane at 50°C for 2 h and stored in the same liquid after it cooled to room temperature. The solvent exchange process was monitored by IR spectroscopy (Supplementary Figure S1). To include guest, about 1 mg of 1-Cyclohexane crystal was added into a capped glass vial. The excess solvent was removed carefully. Then 1 ml of guaiazulene cyclohexane solution (concentration ∼1 mg/ml) was added in the vial. The pierced capped vial was heated at 50°C for 12 h and transferred to 4°C for storage. This sample was noted as 1-Guaiazulene.
3D-ED Data Collection and Processing
The crystal suspension was drop-casted onto a copper grid (R1.2/1.3, QUANTIFOIL). After the solvent was almost evaporated, the grid was plunged into liquid nitrogen rapidly. The grid was then transferred to the Fischone 2550 cryo holder and TEM at liquid nitrogen temperature. The cRED data were collected on a JEOL 2100-plus TEM equipped with MerlinEM direct electron detector under 200 kV acceleration voltage and installed with Heimdall data collection software (software developed by the ReadCrystal Tech Co.). The tilting range depends on the location of the crystals on the grid. Each frame was collected with exposure time of 1 s, resulting in a 1° wedge per frame. The data were visualized with program REDp (Wan et al., 2013) and processed using XDS (Kabsch, 2010) with the aid of with the aid of Coeus (software developed by ReadCrystal Tech Co.) for batch processing and merging.
Results and Discussion
The sample for the CS method was prepared by the following steps: 1) synthesizing host, 2) exchanging solvent, and 3) soaking the molecule of interest. The pores of the as-made CS crystal were filled the nitrobenzene. Due to the high affinity between nitrobenzene and the host, the target molecule cannot be absorbed efficiently into the as-made sample. Thus, it is necessary to exchange the nitrobenzene to a more inert solvent, typically cyclohexane, to facilitate the inclusion of the target molecule. After solvent exchange, the target compound can be absorbed more efficiently into the host.
As a proof-of-concept experiment, we chose to elucidate the structure change in the three phases of the CS method: 1) as synthesized, 2) solvent exchanged, and 3) guest included. The crystalline sponge [(ZnI2)3-(tpt)2·x(Guest)]n (noted as 1-Guest) was synthesized following the protocol of Fujita (Hoshino et al., 2016), using nitrobenzene as the solvent. The as-made sample was donated as 1-Nitrobenzene. Instead of the large crystals used in the traditional CS method, we are more interested about the powder-like crystals at the bottom of the synthesis liquid. Crystals in a size of few hundred nanometers are better for a guest molecule to defuse and also suitable for 3D-ED data collection. To exchange the solvent in the pore from nitrobenzene (high affinity to host) to cyclohexane (less affinity) for better guest-soaking result, the powder-like 1-Nitrobenzene was washed with cyclohexane several times and soaked in it at 50°C for 2 h (donated as 1-Cyclohexane). The IR spectroscopy signal at 1,346 cm−1, assigned to nitrobenzene, completely disappeared (Supplementary Figure S1), indicating that the nitrobenzene within the pore was exchanged by cyclohexane completely. The exchange process made a great improvement in time versus a 7-day-exchange time in the original method using micron-sized crystal. We chose guaiazulene as our target molecule because it was widely tested from different groups and different light sources [Mo (Inokuma et al., 2013), Cu (Hoshino et al., 2016), and synchrotron (Ramadhar et al., 2015)]. The guest was included by slow-evaporation method. 1-Cyclohexane was soaked in the cyclohexane solution of guaiazulene (1 mg/ml) in a pierced capped vial at 50°C for 12 h and then transferred to 4°C for storage. This guest-soaked sample was noted as 1-Guaiazulene. The guest inclusion time (12 h) needed in this work was about fourfold faster than that in the original protocol (2 days) (Hoshino et al., 2016). In addition, these parameters could be further optimized to achieve faster guest exchange.
The as-made (1-Nitrobenzene), solvent-exchanged (1-Cyclohexane), and guest-included (1-Guaiazulene) samples were transferred to TEM, and 3D-ED data under cryogenic temperature were collected to prevent beam damage and the removal/shifting of guest molecules. Based on our experience, the size of the crystal ranging from 500 nm to 1 μm is optimum for 3D-ED data collection. Crystals in this size are large enough to have strong diffractions, yet not too thick for electron to penetrate. The tilting range depends on the crystals’ location on the grid. Each frame was collected with an exposure time of 1 s, resulting in a 1° wedge per frame. To keep diffraction points well-separated, 40 cm of camera length was selected. Thus, the resolution was limited by the detector edge to 1 Å. However, some high-resolution points were still observable at the corner of the detector. The 3D-ED data were processed and merged by XDS (Kabsch, 2010). As shown in Figure 1, the reconstructed 3D reciprocal lattices suggested that 1-Nitrobenzene crystallized in
FIGURE 1. TEM images of (A) 1-Nitrobenzene, (C) 1-Cyclohexane, and (E) 1-Guaiazulene and the 3D reciprocal lattice of (B) 1-Nitrobenzene, (D) 1-Cyclohexane, and (F) 1-Guaiazulene.
TABLE 1. The crystallographic statistics and refinement detail of 1-Nitrobenzene, 1-Cyclohexane, and 1-Guaiazulene.
FIGURE 2. The initial structure model of (A) 1-Nitrobenzene, (B) 1-Cyclohexane, and (C) 1-Guaiazulene obtained from SHELXT.
The final R1 value converged to 0.2265, 0.2192, and 0.2276 for 1-Nitrobenzene, 1-Cyclohexane, and 1-Guaiazulene, respectively. The refined model is shown in Figure 3. The pore size and cell volume of 1-Nitrobenzene are much smaller than the other samples due to the π–π interaction between two complex chains. As shown in Figure 4, the two ligands are offset stacked with their centroid distance of 3.663(7) Å. One nitrobenzene molecule is found in 1-Nitrobenzene, with its occupancy of 0.81(3). After solvent exchange, the framework is considerably expanded. No obvious residue electron density peaks in 1-Cyclohexane could assign to cyclohexane, and the R1 value is comparable with the others. Thus, further SQUEEZE (Spek, 2015) procedure is not needed. For 1-Guaiazulene, the missing part of the target molecule could be found on the difference Fourier map after few rounds of refinement. The low occupancy of guaiazulene, 0.465(15), leads to a slightly inaccurate bond length/angle and slightly distorted conformation. It is worth noting that these three samples have a large fraction of empty void where the residue peaks could not assign to any solvent molecules. The SCXRD CS method often suffers from disordered solvent molecules. It blurs the electron density distribution and complicates the structure refinement. However, in this case, the high-vacuum environment in the TEM chamber may extract the solvent molecules that have a weaker interaction with the host. Without the interference of a disordered solvent, target molecules with low occupancy can still be recognized.
FIGURE 3. The structure model, guest molecule conformation, and Fo density map (drawn in 0.5e-/Å3 level) of (A) 1-Nitrobenzene, (B) 1-Cyclohexane, and (C) 1-Guaiazulene.
If we compare the structure solved by 3D-ED with respect to the reported ones, as shown in Figure 5, the model for 1-Nitrobenzene solved by 3D-ED and that by x-ray (CCDC 187830) (Biradha and Fujita, 2002) are almost identical. For 1-Guaiazulene, the framework position is slightly distorted and the guest molecule location and number were totally different to the synchrotron-resolved one (Ramadhar et al., 2015). We believe this difference originated from the guest-soaking process. The [(ZnI2)3-(tpt)2]n framework is flexible to the response of the guest. Therefore, the guaiazulene-soaking condition will influence the fine structure, i.e., host conformation, number, and the position of the guests, dramatically. The guest location and occupancy in the literature-reported ones are also different (Inokuma et al., 2013; Ramadhar et al., 2015; Hoshino et al., 2016). This result also demonstrates the accuracy of the 3D-ED in the structure determination routine.
FIGURE 5. The comparison of the structure model in this work and the reported one for (A) [(ZnI2)3-(tpt)2·2(C6H5NO2)]n (1-Nitrobenzene) and CCDC 187830 (Biradha and Fujita, 2002) and (B) [(ZnI2)3-(tpt)2·x(C15H18)]n (1-Guaiazulene) and CCDC 1007929 (Ramadhar et al., 2015).
Conclusion
In conclusion, the crystalline sponge method was successfully carried out by the 3D-ED technique using [(ZnI2)3-(tpt)2·x(Guest)]n porous complex nanocrystals. Utilizing the better diffuse ability of small crystals, the time needed for solvent-exchange and guest-soaking protocols were shortened to 2 and 12 h respectively, about 50-fold faster versus the original protocol. It should be noticed that these conditions can be further optimized. The structures of 1-Nitrobenzene, 1-Cyclohexane, and 1-Guaiazulene were fully resolved using a merged dataset. The structure model was identical to previously reported ones using x-ray, showing that the accuracy of the 3D-ED was comparable with SCXRD. The occupancy of the guest molecule, nitrobenzene and guaiazulene, were 0.81(3) and 0.465(15), respectively. This work not only provides a new data collection strategy for crystalline sponge method but also demonstrates the potential of the 3D-ED techniques to study host–guest interaction by solving fine structures of porous materials.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material. The final structures (coordinates, reflections, and structure factors) of 1-Nitrobenzene, 1-Cyclohexane, and 1-Guaiazulene have been deposited in the Cambridge Crystallographic Data Center (CCDC). The CCDC numbers of the samples are 2124118, 2124119, 2124120.
Author Contributions
PC conceived the idea and designed and analyzed the experiments. YL, CZ, and FH collected the 3D-ED data and helped in analyzing the data. PC, LL, and JS cowrote the manuscript. All authors contributed to the discussion of the data and writing of the paper.
Funding
This study was supported by the Ministry of Science and Technology of China (grant number 2020YFA0210700) and National Natural Science Foundation of China (grant number 22125102, 21871009, 21527803).
Conflict of Interest
YL, CZ, FH, and LL are employed by the company ReadCrystal Technology Co., a contract research organization (CRO) that provides MicroED and electron microscopy services to the pharmaceutical/biotechnology industries and academic institutions. JL is the scientific consultant of ReadCrystal Technology Co.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank ReadCrystal Technology Co., for providing the TEM and 3D-ED (MicroED) platform.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.821927/full#supplementary-material
References
Biradha, K., and Fujita, M. (2002). A Springlike 3D-Coordination Network that Shrinks or Swells in a Crystal-to-Crystal Manner upon Guest Removal or Readsorption. Angew. Chem. Int. Ed. 41 (18), 3392–3395. doi:10.1002/1521-3773(20020916)41:18<3392:aid-anie3392>3.0.co;2-v
Clark, L. J., Bu, G., Nannenga, B. L., and Gonen, T. (2021). MicroED for the Study of Protein-Ligand Interactions and the Potential for Drug Discovery. Nat. Rev. Chem. 5, 853–858. doi:10.1038/s41570-021-00332-y
Gemmi, M., and Lanza, A. E. (2019). 3D Electron Diffraction Techniques. Acta Crystallogr. Sect B 75 (Pt 4), 495–504. doi:10.1107/S2052520619007510
Gemmi, M., Mugnaioli, E., Gorelik, T. E., Kolb, U., Palatinus, L., Boullay, P., et al. (2019). 3D Electron Diffraction: The Nanocrystallography Revolution. ACS Cent. Sci. 5 (8), 1315–1329. doi:10.1021/acscentsci.9b00394
Hoshino, M., Khutia, A., Xing, H., Inokuma, Y., and Fujita, M. (2016). The Crystalline Sponge Method Updated. Int. Union Crystallogr. J. 3 (Pt 2), 139–151. doi:10.1107/S2052252515024379
Huang, Z., Grape, E. S., Li, J., Inge, A. K., and Zou, X. (2021). 3D Electron Diffraction as an Important Technique for Structure Elucidation of Metal-Organic Frameworks and Covalent Organic Frameworks. Coord. Chem. Rev. 427, 213583. doi:10.1016/j.ccr.2020.213583
Inokuma, Y., Yoshioka, S., Ariyoshi, J., Arai, T., Hitora, Y., Takada, K., et al. (2013). X-Ray Analysis on the Nanogram to Microgram Scale Using Porous Complexes. Nature 495 (7442), 461–466. doi:10.1038/nature11990
Kabsch, W. (2010). XDS. Acta Crystallogr. D Biol. Crystallogr. 66 (Pt 2), 125–132. doi:10.1107/S0907444909047337
Nannenga, B. L., and Gonen, T. (2019). The Cryo-EM Method Microcrystal Electron Diffraction (MicroED). Nat. Methods 16 (5), 369–379. doi:10.1038/s41592-019-0395-x
Nannenga, B. L., Shi, D., Leslie, A. G. W., and Gonen, T. (2014). High-Resolution Structure Determination by Continuous-Rotation Data Collection in MicroED. Nat. Methods 11 (9), 927–930. doi:10.1038/nmeth.3043
Nguyen, C., and Gonen, T. (2020). Beyond Protein Structure Determination with MicroED. Curr. Opin. Struct. Biol. 64, 51–58. doi:10.1016/j.sbi.2020.05.018
Ramadhar, T. R., Zheng, S.-L., Chen, Y.-S., and Clardy, J. (2015). Analysis of Rapidly Synthesized Guest-Filled Porous Complexes with Synchrotron Radiation: Practical Guidelines for the Crystalline Sponge Method. Acta Cryst. A Found Adv. 71 (1), 46–58. doi:10.1107/S2053273314019573
Sairenji, S., Kikuchi, T., Abozeid, M. A., Takizawa, S., Sasai, H., Ando, Y., et al. (2017). Determination of the Absolute Configuration of Compounds Bearing Chiral Quaternary Carbon Centers Using the Crystalline Sponge Method. Chem. Sci. 8 (7), 5132–5136. doi:10.1039/C7SC01524K
Sheldrick, G. M. (2015a). Crystal Structure Refinement withSHELXL. Acta Crystallogr. C Struct. Chem 71 (Pt 1), 3–8. doi:10.1107/S2053229614024218
Sheldrick, G. M. (2015b). SHELXT- Integrated Space-Group and Crystal-Structure Determination. Acta Crystallogr. A Found Adv. 71 (Pt 1), 3–8. doi:10.1107/S2053273314026370
Spek, A. L. (2015). PLATONSQUEEZE: A Tool for the Calculation of the Disordered Solvent Contribution to the Calculated Structure Factors. Acta Crystallogr. C Struct. Chem. 71 (1), 9–18. doi:10.1107/S2053229614024929
Su, J., Kapaca, E., Liu, L., Georgieva, V., Wan, W., Sun, J., et al. (2014). Structure Analysis of Zeolites by Rotation Electron Diffraction (RED). Microporous Mesoporous Mater. 189, 115–125. doi:10.1016/j.micromeso.2013.10.014
Taniguchi, Y., Matsumoto, R., and Kadota, T. (2020). An Expansion of Crystalline Sponge X‐Ray Analysis to Elucidate the Molecular Structure of Reactive Compounds via Ion Pair Formation. Chem. Eur. J. 26 (68), 15799–15803. doi:10.1002/chem.202002504
Urban, S., Brkljača, R., Hoshino, M., Lee, S., and Fujita, M. (2016). Determination of the Absolute Configuration of the Pseudo‐Symmetric Natural Product Elatenyne by the Crystalline Sponge Method. Angew. Chem. Int. Ed. 55 (8), 2678–2682. doi:10.1002/anie.201509761
Wan, W., Sun, J., Su, J., Hovmöller, S., and Zou, X. (2013). Three-Dimensional Rotation Electron Diffraction: softwareREDfor Automated Data Collection and Data Processing. J. Appl. Crystallogr. 46 (Pt 6), 1863–1873. doi:10.1107/S0021889813027714
Wang, B., Rhauderwiek, T., Inge, A. K., Xu, H., Yang, T., Huang, Z., et al. (2018). A Porous Cobalt Tetraphosphonate Metal-Organic Framework: Accurate Structure and Guest Molecule Location Determined by Continuous-Rotation Electron Diffraction. Chem. Eur. J. 24 (66), 17429–17433. doi:10.1002/chem.201804133
Yan, K., Dubey, R., Arai, T., Inokuma, Y., and Fujita, M. (2017). Chiral Crystalline Sponges for the Absolute Structure Determination of Chiral Guests. J. Am. Chem. Soc. 139 (33), 11341–11344. doi:10.1021/jacs.7b06607
Keywords: crystalline sponge, electron diffraction, host-guest interaction, porous materials, structure solution
Citation: Chen P, Liu Y, Zhang C, Huang F, Liu L and Sun J (2022) Crystalline Sponge Method by Three-Dimensional Electron Diffraction. Front. Mol. Biosci. 8:821927. doi: 10.3389/fmolb.2021.821927
Received: 25 November 2021; Accepted: 22 December 2021;
Published: 07 February 2022.
Edited by:
Hongyi Xu, Stockholm University, SwedenReviewed by:
Brent L. Nannenga, Arizona State University, United StatesAndrew Kentaro Inge, Stockholm University, Sweden
Copyright © 2022 Chen, Liu, Zhang, Huang, Liu and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Junliang Sun, anVubGlhbmcuc3VuQHBrdS5lZHUuY24=