Skip to main content


Front. Chem., 28 April 2021
Sec. Polymer Chemistry
Volume 9 - 2021 |

Iterative Exponential Growth of Oxygen-Linked Aromatic Polymers Driven by Nucleophilic Aromatic Substitution Reactions

www.frontiersin.orgTyler J. Jaynes1 www.frontiersin.orgMona Sharafi1 www.frontiersin.orgJoseph P. Campbell1 www.frontiersin.orgJessica Bocanegra1 www.frontiersin.orgKyle T. McKay1 www.frontiersin.orgKassondra Little1 www.frontiersin.orgReilly Osadchey Brown1 www.frontiersin.orgDanielle L. Gray2 www.frontiersin.orgToby J. Woods2 www.frontiersin.orgJianing Li1 www.frontiersin.orgSeverin T. Schneebeli1*
  • 1Department of Chemistry, University of Vermont, Burlington, VT, United States
  • 2George L. Clark X-Ray Facility and 3M Materials Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL, United States

This work presents the first transition metal-free synthesis of oxygen-linked aromatic polymers by integrating iterative exponential polymer growth (IEG) with nucleophilic aromatic substitution (SNAr) reactions. Our approach applies methyl sulfones as the leaving groups, which eliminate the need for a transition metal catalyst, while also providing flexibility in functionality and configuration of the building blocks used. As indicated by 1) 1H-1H NOESY NMR spectroscopy, 2) single-crystal X-ray crystallography, and 3) density functional theory (DFT) calculations, the unimolecular polymers obtained are folded by nonclassical hydrogen bonds formed between the oxygens of the electron-rich aromatic rings and the positively polarized C–H bonds of the electron-poor pyrimidine functions. Our results not only introduce a transition metal-free synthetic methodology to access precision polymers but also demonstrate how interactions between relatively small, neutral aromatic units in the polymers can be utilized as new supramolecular interaction pairs to control the folding of precision macromolecules.


The backbones of conjugated and heteroatom-linked aromatic polymers tend to possess fewer conformational degrees of freedom than polymers with more flexible aliphatic or partially aliphatic backbones. This reduced amount of conformational freedom can help enhance the folding of aromatic polymers, to advance a variety of useful properties such as selective supramolecular recognition (Goodman et al., 2007; Liu et al., 2015; Otte, 2016; Schneebeli et al., 2016; Adhikari et al., 2017; Xie et al., 2020), selective catalysis (Rajappan et al., 2020; Sharafi et al., 2020), and self-assembly (Cole et al., 2017; Greene et al., 2017; Bonduelle, 2018; Ong and Swager, 2018; Delawder et al., 2019; Zhao et al., 2019). However, while the precise chemical structures, lengths, and sequences of such macromolecules (Dobscha et al., 2019; Gerthoffer et al., 2020; Zhao et al., 2020) dictate their folding, and with it their functionalities and physical properties (Chen et al., 2015; Hanlon et al., 2017), it remains challenging to synthesize polyaromatic structures with precise lengths and/or sequences as unimolecular entities. One of the most efficient ways to precisely control the length and sequence of synthetic polymers is by iteratively coupling (Jones et al., 1997) polymer strands together in a convergent/divergent fashion (Hawker et al., 1997; Read et al., 2000; Grayson and Frechet, 2001; Li et al., 2005; Liess et al., 2006; Binauld et al., 2011). This methodology (Figure 1) is generally referred to as iterative exponential growth (IEG) (Barnes et al., 2015; Leibfarth et al., 2015).


FIGURE 1. Comparison of previously reported transition metal (TM)-catalyzed IEG methods to our TM-free strategy for the construction of well-defined polyaromatic precision macromolecules. Our approach eliminates the need for transition metal catalysts by introducing methyl sulfones as efficient leaving groups for SNAr-based IEG processes.

IEG requires a dormant monomer that includes two distinct functionalities, which are orthogonally deprotected/activated for coupling, to afford two reactive species that can be coupled together selectively. These two reactive species are then combined in a chemo-selective manner, resulting in a dormant dimer, which contains, again, protected/masked functionalities identical to the original dormant monomer (Figure 1). Repeating this simple protocol results in the creation of architecturally defined structures with an exponential gain in polymer length. The IEG method has been applied extensively to synthesize flexible linear polymers (Barnes et al., 2015; Huang et al., 2017)—for example, via copper-catalyzed azide–alkyne click (CuAAC) reactions. However, it still remains a challenge to efficiently utilize arylation chemistry with an IEG-reaction framework, as is required in order to access conjugated or heteroatom-linked polyaromatic polymers with IEG methodology. In particular, while transition metal–catalyzed arylation reactions have been utilized previously to synthesize conjugated and/or heteroatom-linked polyaromatic macromolecules (Figure 1), (Zhang et al., 1992; Louie and Hartwig, 1998; Sadighi et al., 1998; Liess et al., 2006), there are (to the best of our knowledge) currently no transition metal-free reactions available to achieve exponential polymer growth for conjugated or heteroatom-linked polyaromatic structures. We are now able to meet this synthetic challenge by marrying iterative exponential polymer growth with nucleophilic aromatic substitution (SNAr) reactions (Beugelmans et al., 1994; Blaziak et al., 2016; Landovsky et al., 2019) with a unique masking/unmasking technique. Unmasking was achieved by simply oxidizing aromatic sulfides to sulfones—a mild, simple, and scalable reaction well suited for IEG applications. At the same time, the IEG coupling reactions also proceeded under mild conditions (mild heating to ∼50°C), which represents an advantage over transition metal–catalyzed alternatives [which tend to require stronger heating (Sadighi et al., 1998)].

Our new synthetic technology for precision polymer growth helps resolve a number of concerns, which exist with traditional, transition metal–catalyzed cross-coupling methodology (Sun and Shi, 2014; Li et al., 2020). These issues include, but are not limited to 1) the inherent toxicity of many transition metal catalysts, which goes hand in hand with the need to accomplish/demonstrate complete removal of any residual transition metal, especially for healthcare-related applications (which can be difficult to achieve with polymers containing many Lewis basic groups like the ones presented in this article), 2) the relatively high cost of the required transition metal catalysts and complex ligands, which are often needed for high-yielding cross-coupling reactions suitable for polymer growth. Transition metal-free cross-coupling protocols like the one presented in this article have the potential to overcome these fundamental challenges associated with transition metal–catalyzed polymer growth.

Results and Discussion

To implement our sulfide oxidation–driven IEG coupling strategy, we first tested a variety of sulfone leaving groups for their ability to couple with phenolate anions in our IEG scheme (for related work, see: Guan et al., 2015). We discovered that simple methyl sulfones (Figure 1) are best suited for this purpose; as with sulfone substituents larger than methyl (e.g., phenyl), we tended to get lower IEG coupling yields. For our synthetic scheme, we decided to protect the phenols with methyl groups (which can readily be removed in high yield with BBr3). However, it is very likely that instead of just methyl groups, alternative phenol protection groups (e.g., benzyl) can also be also be utilized in the future to further enhance the functional group tolerance of our IEG coupling scheme. To start the IEG couplings, we first condensed 3-methoxyphenol (1) with 4-chloro-6-(methylthio)pyrimidine (2) in quantitative yield under standard SNAr conditions with K2CO3 as the base (Van Rossom et al., 2012). The resulting dimeric product 4-(3-methoxyphenoxy)-6-(methylthio)pyrimidine (3) was then split in half for the next IEG deprotection/activation/coupling cycle. The first aliquot was subjected to BBr3 for methoxyl deprotection, while the second aliquot was activated via oxidization to the methyl sulfone with mCPBA. Following these deprotection/activation steps, we were able to perform an SNAr-based IEG coupling, again in the presence of K2CO3 as the base. Ultimately, by iteratively implementing (Scheme 1) our new, transition metal-free IEG methodology, we were able to access macromolecules with up to 16 perfectly alternating AB units—where A and B represent resorcinol and pyrimidine rings—for the first time in a fully unimolecular fashion.


SCHEME 1. Synthesis of unimolecular polymers with alternating resorcinol and pyrimidine units with SNAr-based IEG couplings, based on masking/unmasking via methyl sulfide to methyl sulfone oxidation. The general reaction conditions for deprotection/activation and IEG coupling are the following: (i) Coupling: K2CO3, DMF, 60°C, 12 h, room temperature. (ii) Deprotection of terminal methoxyl groups to generate nucleophilic phenols: BBr3, CH2Cl2, –78°C to room temperature, 12 h. (iii) Activation of methyl sulfides to generate electrophilic methyl sulfones: mCPBA, EtOAc, room temperature, 24 h.

We hypothesized that the electron-deficient pyrimidine units in our unimolecular oligomers might be able to promote folding (Scheme 1, inset) via hydrogen bonding of and/or π-stacking with the electron-rich resorcinol units. To investigate this hypothesis, we grew single crystals suitable for single-crystal X-ray diffraction analysis of compound 3 by slow vapor diffusion of hexanes into ethyl acetate solutions of 3. The packing observed (Figures 2A,B) upon analysis of the diffraction pattern of a single crystal of compound 3 clearly demonstrates the ability of the positively polarized -N=CH-N= hydrogens in the pyrimidine rings to form intermolecular hydrogen bonds with the methoxyl group of another molecule of 3 in the solid state. There are other hydrogen bonding interactions that are weaker, but still significant, that also contribute (Supplementary Table S1) to the supramolecular packing of this structure. To estimate the strength of the [pyrimidine–CHOCH3] hydrogen bonds observed in the crystal structure of 3, we optimized a model dimer—1,3-dimethoxyphenol in complex with 4-methoxy-6-(methylthio)pyrimidine—with density functional theory at the B3LYP (Lee et al., 1988; Becke, 1993)/6-31G** level of theory. Next, we used the noncovalent interactions (NCI) code (Johnson et al., 2010) implemented in the Jaguar (Bochevarov et al., 2013) software package to find all major attractive supramolecular interactions present in our model dimer, characterized as critical points of the electron density. The NCI critical points of the DFT-optimized structure of the model dimer (Figure 2C) clearly show the presence of a [pyrimidine–CHOCH3] hydrogen bond, with a strength of ∼5.8 kcal/mol.


FIGURE 2. (A, B) Single-crystal X-ray structure of 4-(3-methoxyphenoxy)-6-(methylthio)pyrimidine (3). (C) DFT-optimized structure (B3LYP/6-31G** level) of a model hydrogen-bonded dimer, 1,3-dimethoxyphenol in complex with 4-methoxy-6-(methylthio)pyrimidine. NCI critical points of the electron density (calculated at the B3LYP/6-31G** level of theory) are illustrated with blue spheres. Color code: C: gray; H: white; N: blue; O: red; S: yellow.

Having established that the pyrimidine rings are attracted to the resorcinol units in the solid state via [CHO] hydrogen bonding, we next embarked on investigating the ability of the polymer 6—the longest unimolecular, alternating pyrimidine/resorcinol polymer synthesized to date—to fold, driven by [CHO] hydrogen bonding interactions. To demonstrate folding, we recorded a 1H-1H NOESY NMR spectrum of the hexadecamer 6 in CDCl3 at 298 K. The 1H-1H NOESY NMR spectrum displays (Figure 3A) a NOE cross peak between the proton resonance x of the terminal methoxyl group and the -NCHN- proton resonances (b1b7, which all overlap at the same resonance frequency at 8.48 ppm). Notably, this NOE cross peak is absent in the 1H-1H NOESY NMR spectrum of the dimer 3 (Supplementary Figure S1), which demonstrates that folding and/or supramolecular aggregation (and not just simply the fact that the proton corresponding to resonance b7 is present in a pyrimidine ring next to the terminal methoxyl group) is responsible for the observed NOE cross peak. Furthermore, a 1H DOSY NMR spectrum of 6 (Supplementary Figure S10) confirmed that the compound does not aggregate significantly in solution, and thus the NOE cross peak observed must be stemming from intramolecular folding of the polymer chains. A folded model of the hexadecamer 6 (optimized with DFT at the B3LYP-D3/6-31G** level) is shown in Figure 3B. The structure shown represents the lowest energy conformation obtained from a MacroModel conformational search (carried out without constraints and the OPLS3e force field (Harder et al., 2016) of the hexadecamer 6. While other potential, folded structures exist, the observed NOE cross peak is consistent with the DFT-optimized structure shown in Figure 2B, which shows a [CHO] hydrogen bonding contact between a pyrimidine CH functionality and the terminal methoxyl group in 6. For property characterization, we also performed differential scanning calorimetry (DSC) of 6 (see Supplementary Figure S11 for the DSC thermogram). As 6 is still a relatively low molecular weight compound and fully unimolecular, we did not observe any phase transitions in the temperature range investigated (30–335°C).


FIGURE 3. (A) Partial 1H-1H NOESY NMR (500 MHz, CDCl3, 298 K) of the hexadecamer 6. The key NOE cross peak, which is consistent with polymer folding, is highlighted with an arrow. (B) The corresponding model of a potential folded structure of compound 6. The structure shown represents the lowest energy conformation found with a MacroModel conformational search (OPLS3e force field), which was then refined with a DFT optimization at the B3LYP-D3/6-31G** level of theory.

Overall, our results demonstrate that methyl sulfones can be effective leaving groups for transition metal-free, iterative exponential growth processes. Furthermore, the electrophilic pyrimidine cores present in our polymers not only help further enhance the yields for the SNAr-based IEG coupling reactions but also provide polarized C–H bonds, capable of directing the folding of the unimolecular alternating pyrimidine–resorcinol polymers presented in this work.

Materials and Methods

General Methods and Materials of Synthesis

All commercially available starting materials were purchased from Sigma-Aldrich, Fisher Scientific, or Oakwood Chemical. Unless noted otherwise, all reagents were used as received without further purification. When needed, dichloromethane (CH2Cl2) and dimethylformamide (DMF) were dried using a Glass Contour solvent purification system by SG Water United States, LLC. Removal of solvents was accomplished on a Büchi R-210 rotary evaporator and further concentration was attained under a Fisher Scientific Maxima C-Plus vacuum line. Column chromatography was performed manually with Sorbent grade 60 silica with a mesh size between 230 and 400 using a forced flow of indicated solvents or automatically with a Teledyne CombiFlash® Rf+ chromatography system.

All 1H NMR spectra (see: Supplementary Figures S2, S4, S6, S8) and 13C (1H) NMR spectra (see: Supplementary Figures S3, S5, S7, S9) were recorded on a Bruker ARX 500 (125 MHz) spectrometer, and all 1H-1H NOESY NMR were recorded at 298 K on a Varian Unity Inova 500 (500 MHz) spectrometer. Samples for NMR spectroscopy were dissolved in CDCl3, and the spectra were referenced to the residual solvent peak (CDCl3: 7.26 ppm for 1H and 77.16 ppm for 13C (1H) NMR; or to tetramethylsilane [TMS, 0.00 ppm for 1H and 13C (1H) NMR]) as the internal standard. Chemical shift values were recorded in parts per million (ppm). Data are reported as follows: chemical shift, multiplicity (s = singlet, d = doublet, t = triplet, q = quartet, m = multiplet, and dd = doublet of doublets), coupling constants (Hz), and number of protons. 1H-1H NOESY NMR spectra were acquired with a NOE mixing time of 600 ms. The datasets were processed with MestReNova v10.0.2-15,465 using Bernstein polynomial fits (with a polynomial order of 3) for baseline corrections. The 1H DOSY NMR spectra were acquired on a Varian Unity Inova 500 (500 MHz) spectrometer, equipped with a HCN probe with Z-axis gradients, and a Highland Technologies L700 gradient amplifier. The standard DOSY Varian pulse program “Dbppste” was used, with a stimulated echo sequence and bipolar gradient pulse pairs. All experiments were acquired at 25°C, and DOSY spectra were processed/analyzed using Agilent’s VnmrJ (version 4.2) software, employing the discrete approach for the inverse Laplace transform in the diffusion dimension. High resolution mass spectra were obtained on a Waters XEVO G2-XS QTof in positive ESI mode. Differential scanning calorimetry (DSC) was performed on a PerkinElmer Pyris 1 differential scanning calorimeter on aluminum plates.

General Synthetic Procedure for the IEG Deprotection/Activation and Coupling Steps

Methoxyl Deprotection: For deprotection of the terminal methoxyl groups, dormant oligomers 3, 4, and 5 (0.224 mmol) were dissolved in anhydrous CH2Cl2 (4 ml) under an atmosphere of dry nitrogen. The reaction mixtures were then cooled to –78°C, and a 1.0 M solution of BBr3 in CH2Cl2 (0.896 mmol of BBr3) was added dropwise via syringe into the stirring solution. Next, the reaction mixtures were allowed to warm to room temperature over the course of ∼1 h, stirred at room temperature for an additional 12 h, cooled to 0°C, and finally quenched by adding MeOH (1 ml) and ice chips until fuming ceased. The mixtures were then added to a separatory funnel, and the organic layers were removed. The aqueous layers were extracted with CH2Cl2 (3 × 5 ml), and the combined organic layers were dried over anhydrous Na2SO4, filtered, and concentrated under reduced pressure to afford the deprotected phenols in 88–98% yield (see also Scheme 1). The free phenols were carried forward for the IEG coupling steps without further purification.

Activation of the Electrophilic Coupling Sites via Sulfide Oxidation: Dormant oligomers 3, 4, and 5 (0.23 mmol) were dissolved in anhydrous EtOAc (4 ml), and a solution of 70% meta-chloroperoxybenzoic acid in water (0.92 mmol of mCPBA) was added. After stirring for 16 h at room temperature, the reaction mixtures were added to a separatory funnel, the organic layers were washed with a 0.5 N aqueous NaOH solution (5 × 4 ml) and brine (2 × 3 ml). Finally, the organic layers were dried over anhydrous Na2SO4, filtered, and concentrated under reduced pressure to afford the electrophilic methyl sulfones in 93–95% yield (see also Scheme 1). The activated sulfones were carried forward for the IEG coupling steps without further purification.

Coupling of Active IEG Units: The methyl sulfone oligomers (0.26 mmol) and K2CO3 (0.55 mmol) were added to a flame-dried, 5 ml round-bottomed flask under a nitrogen atmosphere, and anhydrous DMF (4 ml) was added. Next, the oligomers with the free phenolic ends (0.18 mmol) were added to the reaction mixtures at room temperature. Afterward, the reaction mixtures were warmed to 40°C and stirred at 40°C for 12 h. Next, the reaction mixtures were diluted with a 0.5 N aqueous HCl solution (5 ml), transferred to a separatory funnel, and extracted with EtOAc (3 × 5 ml). The combined organic phases were washed with brine (3 × 3 ml) and with a 0.5 N aqueous NaOH solution (2 × 3 ml) to ensure that any potentially remaining phenol was removed. Finally, the combined organic layers were dried over anhydrous Na2SO4, filtered, and concentrated under reduced pressure. The crude products were purified via flash column chromatography over silica gel (eluent: 0–15% ethyl acetate in hexanes) to afford the coupled IEG products 4, 5, and 6 in a 48–98% yield (see also Scheme 1). For the hexadecameric product 6, the yield was determined from a quantitative 1H NMR spectrum obtained from the crude reaction mixture (in the presence of dimethyl sulfone as the internal standard for 1H NMR integration) since a significant amount of material was lost during silica column chromatography purification as a result of the reduced solubility of this relatively large oligomer in the purification solvent. In general, we found that reduced solubility in hexanes of the IEG polymers with increasing lengths presented a challenge for purification of the compounds via chromatographic methods, which resulted in variable yields for the purified products.

Single-Crystal X-Ray Crystallography

Single crystals of 3 were grown by slow vapor diffusion of hexanes into ethyl acetate solutions of 3. Intensity data were collected on a Bruker D8 Venture kappa diffractometer equipped with a Photon II detector. An Iμs microfocus source provided the Mo Kα radiation (λ = 0.71073 A) that was monochromated with multilayer mirrors. The collection, cell refinement, and integration of intensity data were carried out with the APEX3 software (Bruker, 2018). A multi-scan absorption correction was performed with SADABS (Krause et al., 2015). The initial structure solution was solved with the intrinsic phasing methods SHELXT (Sheldrick, 2015b) and refined with the full-matrix least-squares SHELXL (Sheldrick, 2015a) program. The 1 3 5 and -2 5 4 reflections were omitted from the final refinement due to being partially obscured by the beam stop in some orientations.

Crystal data for 3: C12H12N2O2S, Mr = 248.30, crystal size 0.374 × 0.363 × 0.254 mm3, monoclinic, space group P21/n, a = 7.9829(5), b = 14.2529(8), c = 20.6830(11) Å, α = 90°, β = 91.442(2)°, γ = 90°, V = 2352.6(2) Å3, Z = 8, ρcalcd = 1.402 mg m−3, T = 100 (2) K, R1[F2 > 2σ(F2)] = 0.0288, and wR2 = 0.0790.

Density Functional Theory

All structures were optimized with the Jaguar software package at the B3LYP-D3/6-31G** and B3LYP-D3/6-31G** levels (for the NCI critical point calculations) of theory. NCI critical points of the electron density were calculated with the Jaguar software package at the B3LYP-D3/6-31G** level.


This work demonstrates a unique IEG technique, which employs nucleophilic aromatic substitution (SNAr) reactions as a tool to fabricate architecturally defined oxygen-linked ABAB polymers with a metal-free coupling strategy. The monomers are designed to exhibit both nucleophilic and electrophilic character, either of which can be accessed selectively via deprotection of methoxyl groups and via oxidation of aromatic sulfides. Integrating SNAr reactions with iterative exponential growth eliminates the need for expensive transition metal catalysts, which have been required so far in order to carry out iterative exponential growth of conjugated and/or unimolecular, heteroatom-linked aromatic polymers. Our transition metal-free IEG methodology allowed us to access heteroatom-linked, unimolecular aromatic polymers with up to 16 alternating resorcinol and pyrimidine units, enabling us to investigate, for the first time, the [CHO] hydrogen bond–driven folding of such polymers with single-crystal X-ray crystallography, 1H-1H NOESY NMR spectroscopy, and density functional theory (DFT). Overall, our new IEG methodology advances fundamental capacity to access new unimolecular polymers needed to investigate the diverse physical and electronic properties of heteroatom-linked polyaromatic systems. We are currently applying our new transition metal-free IEG methodology to longer precision polymers with extended solubilizing chains.

Data Availability Statement

The crystal data for the single crystal X-ray structure reported in this paper have been deposited in the Cambridge Crystallographic Data Centre (CCDC: 2036161).

Author Contributions

TJ, MS, JC, and SS wrote the article. SS obtained funding for the project. TJ, MS, JC, JB, KM, KL, and RB conducted the experiments. DG and TW solved the crystal structure reported in the article. All authors helped analyze the data and revise the article.


The IEG synthesis of the polymers was supported by the National Science Foundation (Grants CHE-1609137 and CHE-1848444 awarded to SS) as well as by UVM Chemistry, SURF, and Biochemistry undergraduate research fellowship grants awarded to TJ, KL, and ROB. JL was supported by an NSF grant (CHE-1945394) and KM was partially supported by an NIH R01 award (R01GM129431 to JL). The folding studies of the sequence-defined IEG polymers were supported by the Army Research Office (Grant 71015-CH-YIP awarded to SS). The UVM Mass Spectrometry facilities were supported by the National Institutes of Health (Grants S10-OD018126 and P30-GM118228).

Conflict of Interest

The 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.


We thank B. O’Rourke for high-resolution mass spectrometry and M. DiPinto for help with DSC.

Supplementary Material

The Supplementary Material for this article can be found online at:


Adhikari, B., Yamada, Y., Yamauchi, M., Wakita, K., Lin, X., Aratsu, K., et al. (2017). Light-induced unfolding and refolding of supramolecular polymer nanofibres. Nat. Commun. 8, 15254. doi:10.1038/ncomms15254

PubMed Abstract | CrossRef Full Text | Google Scholar

Barnes, J. C., Ehrlich, D. J., Gao, A. X., Leibfarth, F. A., Jiang, Y., Zhou, E., et al. (2015). Iterative exponential growth of stereo- and sequence-controlled polymers. Nat. Chem. 7, 810–815. doi:10.1038/nchem.2346

PubMed Abstract | CrossRef Full Text | Google Scholar

Becke, A. D. (1993). Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. doi:10.1063/1.464913

CrossRef Full Text | Google Scholar

Beugelmans, R., Singh, G. P., Bois-Choussy, M., Chastanet, J., and Zhu, J. (1994). SNAr-based macrocyclization: an application to the synthesis of vancomycin family models. J. Org. Chem. 59, 5535–5542. doi:10.1021/jo00098a010

CrossRef Full Text | Google Scholar

Binauld, S., Damiron, D., Connal, L. A., Hawker, C. J., and Drockenmuller, E. (2011). Precise synthesis of molecularly defined oligomers and polymers by orthogonal iterative divergent/convergent approaches. Macromol. Rapid Commun. 32, 147–168. doi:10.1002/marc.201000548

PubMed Abstract | CrossRef Full Text | Google Scholar

Blaziak, K., Danikiewicz, W., and Makosza, M. (2016). How does nucleophilic aromatic substitution really proceed in nitroarenes? computational prediction and experimental verification. J. Am. Chem. Soc. 138, 7276–7281. doi:10.1021/jacs.5b13365

PubMed Abstract | CrossRef Full Text | Google Scholar

Bochevarov, A. D., Harder, E., Hughes, T. F., Greenwood, J. R., Braden, D. A., Philipp, D. M., et al. (2013). Jaguar: a high-performance quantum chemistry software program with strengths in life and materials sciences. Int. J. Quantum Chem. 113, 2110–2142. doi:10.1002/qua.24481

CrossRef Full Text | Google Scholar

Bonduelle, C. (2018). Secondary structures of synthetic polypeptide polymers. Polym. Chem. 9, 1517–1529. doi:10.1039/c7py01725a

CrossRef Full Text | Google Scholar

Bruker, (2018). APEX3. Madison, WI: Bruker AXS, Inc.

Chen, L., Wang, H., Zhang, D. W., Zhou, Y., and Li, Z. T. (2015). Quadruple switching of pleated foldamers of tetrathiafulvalene-bipyridinium alternating dynamic covalent polymers. Angew. Chem. Int. Ed. Engl. 54, 4028–4031. doi:10.1002/anie.201410757

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, J. P., Hanlon, A. M., Rodriguez, K. J., and Berda, E. B. (2017). Protein-like structure and activity in synthetic polymers. J. Polym. Sci. Part A: Polym. Chem. 55, 191–206. doi:10.1002/pola.28378

CrossRef Full Text | Google Scholar

Delawder, A. O., Natraj, A., Colley, N. D., Saak, T., Greene, A. F., and Barnes, J. C. (2019). Synthesis, self-assembly, and photomechanical actuator performance of a sequence-defined polyviologen crosslinker. Supramol. Chem. 31, 523–531. doi:10.1080/10610278.2019.1632453

CrossRef Full Text | Google Scholar

Dobscha, J. R., Castillo, H. D., Li, Y., Fadler, R. E., Taylor, R. D., Brown, A. A., et al. (2019). Sequence-defined macrocycles for understanding and controlling the build-up of hierarchical order in self-assembled 2D arrays. J. Am. Chem. Soc. 141, 17588–17600. doi:10.1021/jacs.9b06410

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerthoffer, M. C., Wu, S., Chen, B., Wang, T., Huss, S., Oburn, S. M., et al. (2020). ‘Sacrificial’ supramolecular assembly and pressure-induced polymerization: toward sequence-defined functionalized nanothreads. Chem. Sci. 11, 11419–11424. doi:10.1039/D0SC03904G

CrossRef Full Text | Google Scholar

Goodman, C. M., Choi, S., Shandler, S., and Degrado, W. F. (2007). Foldamers as versatile frameworks for the design and evolution of function. Nat. Chem. Biol. 3, 252–262. doi:10.1038/nchembio876

PubMed Abstract | CrossRef Full Text | Google Scholar

Grayson, S. M., and Fréchet, J. M. (2001). Convergent dendrons and dendrimers: from synthesis to applications. Chem. Rev. 101, 3819–3868. doi:10.1021/cr990116h

PubMed Abstract | CrossRef Full Text | Google Scholar

Greene, A. F., Danielson, M. K., Delawder, A. O., Liles, K. P., Li, X., Natraj, A., et al. (2017). Redox-responsive artificial molecular muscles: reversible radical-based self-assembly for actuating hydrogels. Chem. Mater. 29, 9498–9508. doi:10.1021/acs.chemmater.7b03635

CrossRef Full Text | Google Scholar

Guan, Y., Wang, C., Wang, D., Dang, G., Chen, C., Zhou, H., et al. (2015). Methylsulfone as a leaving group for synthesis of hyperbranched poly(arylene pyrimidine ether)s by nucleophilic aromatic substitution. RSC Adv. 5, 12821–12823. doi:10.1039/c5ra00634a

CrossRef Full Text | Google Scholar

Hanlon, A. M., Martin, I., Bright, E. R., Chouinard, J., Rodriguez, K. J., Patenotte, G. E., et al. (2017). Exploring structural effects in single-chain “folding” mediated by intramolecular thermal Diels-Alder chemistry. Polym. Chem. 8, 5120–5128. doi:10.1039/c7py00320j

CrossRef Full Text | Google Scholar

Harder, E., Damm, W., Maple, J., Wu, C., Reboul, M., Xiang, J. Y., et al. (2016). OPLS3: a force field providing broad coverage of drug-like small molecules and proteins. J. Chem. Theory Comput. 12, 281–296. doi:10.1021/acs.jctc.5b00864

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawker, C. J., Malmström, E. E., Frank, C. W., and Kampf, J. P. (1997). Exact linear analogs of dendritic polyether macromolecules: design, synthesis, and unique properties. J. Am. Chem. Soc. 119, 9903–9904. doi:10.1021/ja972027x

CrossRef Full Text | Google Scholar

Huang, Z., Zhao, J., Wang, Z., Meng, F., Ding, K., Pan, X., et al. (2017). Combining orthogonal chain-end deprotections and thiol-maleimide Michael coupling: engineering discrete oligomers by an iterative growth strategy. Angew. Chem. Int. Ed. Engl. 56, 13612–13617. doi:10.1002/anie.201706522

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, E. R., Keinan, S., Mori-Sánchez, P., Contreras-García, J., Cohen, A. J., and Yang, W. (2010). Revealing noncovalent interactions. J. Am. Chem. Soc. 132, 6498–6506. doi:10.1021/ja100936w

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, L., Schumm, J. S., and Tour, J. M. (1997). Rapid solution and solid phase syntheses of oligo(1,4-phenylene ethynylene)s with thioester termini: molecular scale wires with alligator clips. derivation of iterative reaction efficiencies on a polymer support. J. Org. Chem. 62, 1388–1410. doi:10.1021/jo962336q

CrossRef Full Text | Google Scholar

Krause, L., Herbst-Irmer, R., Sheldrick, G. M., and Stalke, D. (2015). Comparison of silver and molybdenum microfocus X-ray sources for single-crystal structure determination. J. Appl. Crystallogr. 48, 3–10. doi:10.1107/S1600576714022985

PubMed Abstract | CrossRef Full Text | Google Scholar

Landovský, T., Eigner, V., Babor, M., Tichotová, M., Dvořáková, H., and Lhoták, P. (2019). Regioselective SNAr reaction of the phenoxathiin-based thiacalixarene as a route to a novel macrocyclic skeleton. Chem. Commun. (Camb) 56, 78–81. doi:10.1039/c9cc08335a

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C., Yang, W., and Parr, R. G. (1988). Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B Condens. Matter 37, 785–789. doi:10.1103/physrevb.37.785

PubMed Abstract | CrossRef Full Text | Google Scholar

Leibfarth, F. A., Johnson, J. A., and Jamison, T. F. (2015). Scalable synthesis of sequence-defined, unimolecular macromolecules by Flow-IEG. Proc. Natl. Acad. Sci. U.S.A. 112, 10617–10622. doi:10.1073/pnas.1508599112

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Nieves-Quinones, Y., Zhang, H., Liang, Q., Su, S., Liu, Q., et al. (2020). Transition-metal-free formal cross-coupling of aryl methyl sulfoxides and alcohols via nucleophilic activation of C-S bond. Nat. Commun. 11, 2890. doi:10.1038/s41467-020-16713-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Wang, X., and Wang, F. (2005). A novel in situ deprotection/coupling and iterative divergent/convergent strategy for the synthesis of oligo(1,4-phenyleneethynylene)s. Tetrahedron Lett. 46, 8971–8973. doi:10.1016/j.tetlet.2005.10.113

CrossRef Full Text | Google Scholar

Liess, P., Hensel, V., and Schlüter, A.-D. (2006). Oligophenylene rods: a repetitive approach. Eur. J. Org. Chem. 1996, 1037–1040. doi:10.1002/jlac.199619960703

CrossRef Full Text | Google Scholar

Liu, X., Weinert, Z. J., Sharafi, M., Liao, C., Li, J., and Schneebeli, S. T. (2015). Regulating molecular recognition with C-shaped strips attained by chirality-assisted synthesis. Angew. Chem. Int. Ed. Engl. 54, 12772–12776. doi:10.1002/anie.201506793

PubMed Abstract | CrossRef Full Text | Google Scholar

Louie, J., and Hartwig, J. F. (1998). The largest discrete oligo(m-aniline). An Exponential growth strategy using palladium-catalyzed amination of aryl sulfonates. Macromolecules 31, 6737–6739. doi:10.1021/ma981073u

CrossRef Full Text | Google Scholar

Ong, W. J., and Swager, T. M. (2018). Dynamic self-correcting nucleophilic aromatic substitution. Nat. Chem. 10, 1023–1030. doi:10.1038/s41557-018-0122-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Otte, M. (2016). Size-selective molecular flasks. ACS Catal. 6, 6491–6510. doi:10.1021/acscatal.6b01776

CrossRef Full Text | Google Scholar

Rajappan, S. C., Mccarthy, D. R., Campbell, J. P., Ferrell, J. B., Sharafi, M., Ambrozaite, O., et al. (2020). Selective monofunctionalization enabled by reaction-history-dependent communication in catalytic rotaxanes. Angew. Chem. Int. Ed. Engl. 59, 16668. doi:10.1002/anie.202006305

PubMed Abstract | CrossRef Full Text | Google Scholar

Read, M. W., Escobedo, J. O., Willis, D. M., Beck, P. A., and Strongin, R. M. (2000). Convenient iterative synthesis of an octameric tetracarboxylate-functionalized oligophenylene rod with divergent end groups. Org. Lett. 2, 3201–3204. doi:10.1021/ol000204k

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadighi, J. P., Singer, R. A., and Buchwald, S. L. (1998). Palladium-catalyzed synthesis of monodisperse, controlled-length, and functionalized oligoanilines. J. Am. Chem. Soc. 120, 4960–4976. doi:10.1021/ja980052c

CrossRef Full Text | Google Scholar

Schneebeli, S., Li, J., Sharafi, M., Weinert, Z., Cohen, I., Liao, C., et al. (2016). Controlled self-assembly inside C-shaped polyaromatic strips. Synlett 27, 2145–2149. doi:10.1055/s-0035-1561479

CrossRef Full Text | Google Scholar

Sharafi, M., Mckay, K. T., Ivancic, M., Mccarthy, D. R., Dudkina, N., Murphy, K. E., et al. (2020). Size-selective catalytic polymer acylation with a molecular tetrahedron. Chem 6, 1469–1494. doi:10.1016/j.chempr.2020.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheldrick, G. M. (2015a). Crystal structure refinement with SHELXL. Acta Crystallogr. C Struct. Chem. 71, 3–8. doi:10.1107/S2053229614024218

CrossRef Full Text | Google Scholar

Sheldrick, G. M. (2015b). SHELXT-integrated space-group and crystal-structure determination. Acta Crystallogr. A. Found Adv. 71, 3–8. doi:10.1107/S2053273314026370

CrossRef Full Text | Google Scholar

Sun, C. L., and Shi, Z. J. (2014). Transition-metal-free coupling reactions. Chem. Rev. 114, 9219–9280. doi:10.1021/cr400274j

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Rossom, W., Caers, J., Robeyns, K., Van Meervelt, L., Maes, W., and Dehaen, W. (2012). (Thio)ureido anion receptors based on a 1,3-alternate oxacalix[2]arene[2]pyrimidine scaffold. J. Org. Chem. 77, 2791–2797. doi:10.1021/jo300004p

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, J., Li, X., Wang, S., Li, A., Jiang, L., and Zhu, K. (2020). Heteroatom-bridged molecular belts as containers. Nat. Commun. 11, 3348. doi:10.1038/s41467-020-17134-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Moore, J. S., Xu, Z., and Aguirre, R. A. (1992). Nanoarchitectures. 1. Controlled synthesis of phenylacetylene sequences. J. Am. Chem. Soc. 114, 2273–2274. doi:10.1021/ja00032a060

CrossRef Full Text | Google Scholar

Zhao, W., Qiao, B., Tropp, J., Pink, M., Azoulay, J. D., and Flood, A. H. (2019). Linear supramolecular polymers driven by anion-anion dimerization of difunctional phosphonate monomers inside cyanostar macrocycles. J. Am. Chem. Soc. 141, 4980–4989. doi:10.1021/jacs.9b00248

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, W., Tropp, J., Qiao, B., Pink, M., Azoulay, J. D., and Flood, A. H. (2020). Tunable adhesion from stoichiometry-controlled and sequence-defined supramolecular polymers emerges hierarchically from cyanostar-stabilized anion-anion linkages. J. Am. Chem. Soc. 142, 2579–2591. doi:10.1021/jacs.9b12645

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: iterative convergent/divergent polymer synthesis, SNAr reactions, nonclassical hydrogen bonding, nuclear magnetic resonance spectroscopy, polymer folding, iterative exponential polymer growth, transition metal-free coupling

Citation: Jaynes TJ, Sharafi M, Campbell JP, Bocanegra J, McKay KT, Little K, Osadchey Brown R, Gray DL, Woods TJ, Li J and Schneebeli ST (2021) Iterative Exponential Growth of Oxygen-Linked Aromatic Polymers Driven by Nucleophilic Aromatic Substitution Reactions. Front. Chem. 9:620017. doi: 10.3389/fchem.2021.620017

Received: 21 October 2020; Accepted: 11 February 2021;
Published: 28 April 2021.

Edited by:

Jie He, University of Connecticut, United States

Reviewed by:

Mina Mazzeo, University of Salerno, Italy
Lingyao Meng, University of New Mexico, United States

Copyright © 2021 Jaynes, Sharafi, Campbell, Bocanegra, McKay, Little, Osadchey Brown, Gray, Woods, Li and Schneebeli. 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: Severin T. Schneebeli,

These authors have contributed equally to this work