Coordination Chemistry of a Bis(Tetrazine) Tweezer: A Case of Host-Guest Behavior with Silver Salts

The carbon-carbon cross-coupling of phenyl s-tetrazine (Tz) units at their ortho-phenyl positions allows the formation of constrained bis(tetrazines) with original tweezer structures. In these compounds, the face-to-face positioning of the central tetrazine cores is reinforced by π-stacking of the electron-poor nitrogen-containing heteroaromatic moieties. The resulting tetra-aromatic structure can be used as a weak coordinating ligand with cationic silver. This coordination generates a set of bis(tetrazine)-silver(I) coordination complexes tolerating a large variety of counter anions of various geometries, namely, PF6−, BF4−, SbF6−, ClO4−, NTf2−, and OTf−. These compounds were characterized in the solid state by single-crystal X-ray diffraction (XRD) and diffuse reflectance spectroscopy, and in solution by 1H-NMR, mass spectrometry, electroanalysis, and UV-visible absorption spectrophotometry. The X-ray crystal structure of complexes {[Ag(3)][PF6]}∞ (4) and {[Ag(3)][SbF6]}∞ (6), where 3 is 3,3′-[(1,1′-biphenyl)-2,2′-diyl]-6,6′-bis(phenyl)-1,2,4,5-tetrazine, revealed the formation of 1D polymeric chains, characterized by an evolution to a large opening of the original tweezer and a coordination of silver(I) via two chelating nitrogen atom and some C=C π-interactions. Electrochemical and UV spectroscopic properties of the original tweezer and of the corresponding silver complexes are reported and compared. 1H-NMR titrations with AgNTf2 allowed the determination of the stoichiometry and apparent stability of two solution species, namely [Ag(3)]+ and [Ag(3)2]2+, that formed in CDCl3/CD3OD 2:1 v/v mixtures.


Introduction
The s-tetrazine (Tz) unit (Figure 1, 1) is an object of high interest in photophysics and for biomedical applications, owing to its electron-poor poly(hetero)aromatic nature that is also amenable to click chemistry [1][2][3]. The diversification of synthetic routes to conceive new tetrazine structures has become very useful to expand current applications [3][4][5][6][7][8][9]. We recently reported a copper-catalyzed homocoupling protocol, which allowed the efficient synthesis of constrained bis(tetrazines) with a unique bridge clamp structure (Figure 1, 3) [3]. These compounds have been studied in the solid state by X-ray diffraction (XRD, Figure 1a) and analyzed by DFT. The results have shown that their bridge clamp structure (Figure 1b) is due to the unexpected existence of stacking interactions between the [C 2 N 4 ]·[C 2 N 4 ] electronpoor central tetrazine cores (Figure 1c). We have established that London dispersion forces between the Tz moieties clearly stabilize and favor the clamp geometry. In their absence, the tweezer-like interconversion barrier between the most stable rotamers corresponding s-Tetrazine derivatives may be successfully used as ligands for transition metals [10,11]. In comparison, the coordination properties of bis(tetrazines) have been rarely studied [12][13][14][15]. Herein, we investigated the binding properties of bis(tetrazine) 3 with various metals, focusing on the Ag + cation, thus assessing the global stability of the tweezer structure and its degree of flexibility around the ortho-C-C bond. We found that 3 is a weak ligand of several transition metals in non-coordinating low polarity solvents. Interactions with silver(I) salts generate bis(tetrazine)-silver(I) coordination polymers, which tolerate a large variety of counter anions including PF6 − , BF4 − , SbF6 − , ClO4 − , NTf2 -, and OTf − . Accordingly, the XRD (6) revealed the formation of 1D polymeric chains characterized by a large opening of the original tweezer structure (with θ = 144.0(2) ° for 4 and 145.20(3) ° for 6, vs. θ = 54° in 3). Electrochemical behavior and UV spectroscopic characterization of the tweezer and its silver complexes are reported, and further studies in solution focused on the 1 H-NMR of the bis(tetrazine) 3 in comparison with the related mono(tetrazine) 1 and the formation of the bis(tetrazine)-silver(I) coordination polymers 4 and 6. Our study showed a typical tweezer-Ag + host-guest interaction behavior in solution, which could be further exploited for easy anion exchange.

Results and Discussion
Bis(tetrazine) 3 and its functionalized derivatives display in the solid state a centroidto-centroid separation of the central rings ranging between 3.3 and 3.5 Å (Figure 1) [3]. Given this structure and the presence of eight nitrogen atoms as potential coordinating donors, we were intrigued about the potential coordination properties of this new bridge clamp structure with transition metals. Would it be possible that a sandwich-like association reminding classical metallocenes and hemi-metallocenes could form? Moreover, does the framework retain sufficient conformational freedom to accommodate a large range of coordination modes with various metals? With the view of addressing these issues, we conducted comparative coordination chemistry studies of the bis(tetrazine) 3 in solution and in the solid state. s-Tetrazine derivatives may be successfully used as ligands for transition metals [10,11]. In comparison, the coordination properties of bis(tetrazines) have been rarely studied [12][13][14][15]. Herein, we investigated the binding properties of bis(tetrazine) 3 with various metals, focusing on the Ag + cation, thus assessing the global stability of the tweezer structure and its degree of flexibility around the ortho-C-C bond. We found that 3 is a weak ligand of several transition metals in non-coordinating low polarity solvents. Interactions with silver(I) salts generate bis(tetrazine)-silver(I) coordination polymers, which tolerate a large variety of counter anions including PF 6  . Electrochemical behavior and UV spectroscopic characterization of the tweezer and its silver complexes are reported, and further studies in solution focused on the 1 H-NMR of the bis(tetrazine) 3 in comparison with the related mono(tetrazine) 1 and the formation of the bis(tetrazine)silver(I) coordination polymers 4 and 6. Our study showed a typical tweezer-Ag + host-guest interaction behavior in solution, which could be further exploited for easy anion exchange.

Results and Discussion
Bis(tetrazine) 3 and its functionalized derivatives display in the solid state a centroidto-centroid separation of the central rings ranging between 3.3 and 3.5 Å (Figure 1) [3]. Given this structure and the presence of eight nitrogen atoms as potential coordinating donors, we were intrigued about the potential coordination properties of this new bridge clamp structure with transition metals. Would it be possible that a sandwich-like association reminding classical metallocenes and hemi-metallocenes could form? Moreover, does the framework retain sufficient conformational freedom to accommodate a large range of coordination modes with various metals? With the view of addressing these issues, we conducted comparative coordination chemistry studies of the bis(tetrazine) 3 in solution and in the solid state.

Bis(Tetrazine) 3 for Metal Complexation
The silver(I) cation is a suitable candidate for metal complexation because it can adopt several coordination numbers and geometries. Moreover, as a soft cation Ag + is known to interact with tetrazine ligands [16,17]. The reaction of the bis(tetrazine) 3 with one equivalent of AgPF 6 in a CH 2 Cl 2 /MeOH mixture (2:1 v/v) proceeds rapidly, affording in 88% yield an air-stable, soluble red-orange complex, denoted [Ag(3)]PF 6 , (4 in Scheme 1). When the silver salt was used in excess compared to 3, an orange solid was formed that turned out to be insoluble in CH 2 Cl 2, suggesting the formation of polymeric species.

Bis(Tetrazine) 3 for Metal Complexation
The silver(I) cation is a suitable candidate for metal complexation because it can adopt several coordination numbers and geometries. Moreover, as a soft cation Ag + is known to interact with tetrazine ligands [16,17]. The reaction of the bis(tetrazine) 3 with one equivalent of AgPF6 in a CH2Cl2/MeOH mixture (2:1 v/v) proceeds rapidly, affording in 88% yield an air-stable, soluble red-orange complex, denoted [Ag(3)]PF6, (4 in Scheme 1). When the silver salt was used in excess compared to 3, an orange solid was formed that turned out to be insoluble in CH2Cl2, suggesting the formation of polymeric species. Scheme 1. Synthesis of cationic silver(I) complexes 4-9 obtained from the stoichiometric reaction of AgX salts (X = PF6 − , BF4 − , SbF6 − , ClO4 − , NTf2 − , and OTf − ) with the bis(tetrazine) 3.
The complexes 5-9 (Scheme 1) were isolated quantitatively upon reacting 3 with a variety of AgX salts (where X = BF4 − , SbF6 − , ClO4 − , NTf2 − , and OTf − , respectively) under similar conditions. In all cases, the use of an excess of Ag(I) salt led to the formation of poorly soluble aggregates that hampered detailed characterization in solution. The binding of the bis(tetrazine) 3 to other transition metal cations, such as Cu + and Pd 2+ , was also investigated. [Cu(CH3CN)4]BF4, [Pd(CH3CN)2Cl2], and [Pd(CH3CN)4](BF4)2 were used as metal sources and clearly appeared to readily react with 3. The complex denoted [Cu(3)]BF4 was analyzed in solution by NMR spectroscopy (details in SI), but the compounds resulting from the reaction with palladium (tentatively denoted [Pd(3)] + and [Pd(3)Cl] + according to mass-spectrometric analyses) were poorly soluble and could not be easily characterized by NMR spectroscopy in solution. However, HR-mass spectrometry could be conducted, and ESI-MS spectra exhibited peaks at m/z = 529, 607, and 572 a.m.u. with an isotopic distribution pattern consistent with [Cu(3)] + , [Pd(3)] + , and [Pd(3)Cl] + , respectively (details in SI). Conversely, our attempts to isolate and fully characterize those complexes formed by reacting 3 with the zerovalent precursors [Pd2(dba)3] and Cr(CO)6 failed, obviating the eventuality of stabilizing metallocene-like complexes by the controlled opening of tweezer 3.

Solid-State Characterization of Silver(I) Complexes
Single crystals of complexes 4 or 6 were obtained by slow evaporation of their solution in dichloromethane and in a chloroform/pentane mixture, respectively. The crystallographic structure of those complexes (Figures 2 and 3) was solved in the orthorhombic Pnna and Pbcm space groups, for 4 and 6, respectively. In these crystal structures, the bis(tetrazine) 3 adopts an open shape where the two tetrazine rings are interconnected through the coordination of a silver(I) cation by the doublets of N1 nitrogen atoms. This arrangement formed a rare nine-membered metallacycle. In the solid state, the bis(tetrazine) 3 displayed a torsion angle between the two Tz moieties equal to 51.1(3) °. Upon silver binding, the corresponding C2-C7-C7 i -C2 i torsion angle opens up to 144.0(2) ° for 4, and the analogous C2-C3-C3 i -C2 i torsion angle is consistently found equal to 145.2(3) ° for 6. Accordingly, the distance between the two C2N4 centroids (Ct) in the Tz changed from 3.320(1) Å in 3 to 6.6846(9) Å in 4 and 6.6895(16) Å in 6 after silver coordination. In the complexes, the Tz cores in the monomer unit are roughly perpendicular to each other, with angles equal to 76.73 ° for 4 and 78.05 ° for 6, and each of them is parallel to the tetrazine ring of a neighboring monomer (see Figures 2 and 3). Complexes 4 and 6 exhibit . Conversely, our attempts to isolate and fully characterize those complexes formed by reacting 3 with the zerovalent precursors [Pd 2 (dba) 3 ] and Cr(CO) 6 failed, obviating the eventuality of stabilizing metallocene-like complexes by the controlled opening of tweezer 3.
1 H-NMR spectra of 3 were also recorded in a series of solvents of various polarities, including toluene-d 8 , CDCl 3 , THF-d 8 , (CD 3 ) 2 CO, CD 3 CN, and (CD 3 ) 2 SO. The same spectral profile was observed in CDCl 3 and CD 2 Cl 2 , but noticeable splitting and shifts were observed with the other solvents ( Figure S1). More polar solvents such as CD 3 CN or DMSO-d 6 -but also toluene-showed prominent shifts for most of the resonances, suggesting solvent-induced structural changes in 3. These could be related either to the change of the equilibrium value of the twist-angle θ, since intramolecular π-π interactions are less favored in more polar media, or some solvent-π intercalation in toluene. Accordingly, the molecular tweezer was predicted to possibly open up, i.e., widening of the θ angle, in such dissociating solvents.  1 H-NMR spectra of 3 were also recorded in a series of solvents of various polarities, including toluene-d8, CDCl3, THF-d8, (CD3)2CO, CD3CN, and (CD3)2SO. The same spectral profile was observed in CDCl3 and CD2Cl2, but noticeable splitting and shifts were observed with the other solvents ( Figure S1). More polar solvents such as CD3CN or DMSO-d6 -but also toluene-showed prominent shifts for most of the resonances, suggesting solvent-induced structural changes in 3. These could be related either to the change of the equilibrium value of the twist-angle θ, since intramolecular π-π interactions are less favored in more polar media, or some solvent-π intercalation in toluene. Accordingly, the molecular tweezer was predicted to possibly open up, i.e., widening of the θ angle, in such dissociating solvents.
While the 1 H-NMR spectrum of the uncoordinated compound 3 in CDCl3 exhibits seven overlapping signals (H1,5, H2,4, H3, H6, H7, H8, and H9) at room temperature in agreement with its symmetry, this set of resonances is split in CDCl3 for the corresponding bis(tetrazine) silver(I) coordination complexes 4-9 ( Figure 6). This symmetry disruption is a good indication for the Ag + uptake by the bis(tetrazine) ligand 3 in solution. The coordination mode in solution could be derived from the single-crystal structure of 4 obtained by XRD (see above section). The complete signal attribution is reported in Table 1, including the J spin-spin coupling constants (SSCCs). While the 1 H-NMR spectrum of the uncoordinated compound 3 in CDCl 3 exhibits seven overlapping signals (H 1,5 , H 2,4 , H 3 , H 6 , H 7 , H 8 , and H 9 ) at room temperature in agreement with its symmetry, this set of resonances is split in CDCl 3 for the corresponding bis(tetrazine) silver(I) coordination complexes 4-9 ( Figure 6). This symmetry disruption is a good indication for the Ag + uptake by the bis(tetrazine) ligand 3 in solution. The coordination mode in solution could be derived from the single-crystal structure of 4 obtained by XRD (see above section). The complete signal attribution is reported in Table 1, including the J spin-spin coupling constants (SSCCs).

Attribution
Chemical Shift (ppm in CD 2 Cl 2 ) SSCCs (Hz) Figure 6. COSY 1 H-NMR (400 MHz, 298 K) spectrum of complex 4 in CD2Cl2. As shown in Figure 7 for the complexes 4-8, the H9 proton resonance is clearly the most affected one by the coordination, as it moves to stronger fields, with chemical shift changes of Δδ = 0.41 (4, PF6 -), 0.48 (5, BF4 -), 0.55 (6, SbF6 -), 0.65 (7, ClO4 -), and 0.70 (8, ppm. In addition to the strong effect of silver coordination to the o-C-H protons, the nature of the counter anion also has a significant impact on the shift in a given solvent that is consistent with a strong ion-pairing in solution.  Figure S2, AgPF6 case). However, Ag + being a very labile cation, it is difficult to conclude whether or not these ML2 and M2L2 species exist in solution, or if they might be formed in the MS gas phase by a rapid displacement of the equilibria during the drying of the spray. No evidence for any higher nuclearity species (oligomers) could be gained by ESI-MS.
The UV-Vis diffuse reflectance spectra of microcrystalline compounds diluted in dry BaSO 4 were acquired with a "praying mantis" (Figure S4). The spectral shape of a solid sample of 3 ( Figure S4) resembles that found in CH 2 Cl 2 solution, albeit both band maxima λ max = 301 and 568 nm displayed a bathochromic shift with respect to solution data, while the UV π-π* band is flanked at low energy by a broad shoulder (λ sh ca. 390 nm). The latter is assigned to an n-π* transition, with a HOMO → LUMO+2 character. The amplitude of the bathochromic shift of both UV and vis bands is the largest in methanol, with that is defined as ∆λ = λ solidλ solution = 14 and 26 nm, respectively (Table S1). bis(tetrazine) 3. A weak bathochromic shift is observed for the band around 292 nm, and a hypsochromic shift for the band around 552 nm, relative to 3, although the displacement does not exceed 4 nm ( Figure S5, Table S2). In the solid state, both electronic transitions occur at lower energy with respect to the parent compound 3, which translates into 10 and 6 nm shifts of the UV and vis bands, respectively. Moreover, the 390 nm shoulder observed in the diffuse reflectance spectrum of 3 undergoes likewise a red shift upon silver binding ( Figure S6, Table S3).

Electrochemical Properties
Cyclic voltammetry (CV) was employed to study the electrochemical behavior of 3,6diphenyl-1,2,4,5-tetrazine 1, the bis(tetrazine) 3 and the complex {[Ag(3)][BF4]}∞ (5). The CV experiments were carried out in CH2Cl2 containing 0.1 M tetra(n-ethyl)ammonium tetrafluoroborate (TEABF4), see Figure 9, Figure S7, and Table S4, as supporting electrolyte. The voltammogram of 1 displays a single reversible reduction wave (E1/2 = −0.85 V/SCE) corresponding to the formation of the tetrazine radical anion [22,23], whereas bis(tetrazine) 3 undergoes two successive one-electron reductions at E1/2 = −0.83 and −0.93 V/SCE. According to theoretical approaches [24], molecules undergoing two successive reversible one-electron transfers centered on chemically equivalent, and fully independent redox sites (such as, for example, bis(ferrocenyl)porphyrins) [25] give rise to a 35.6 mV shift between each successive formal half-wave potentials. The resulting CV curve for each one-electron transfer should display the same 60 mV peak-to-peak separation [24,26]. The UV-Vis spectra of complexes 4-8 were collected likewise in solution and in the solid state. In solution, the complexes 4-8 absorb in the same wavelength region as the bis(tetrazine) 3. A weak bathochromic shift is observed for the band around 292 nm, and a hypsochromic shift for the band around 552 nm, relative to 3, although the displacement does not exceed 4 nm ( Figure S5, Table S2). In the solid state, both electronic transitions occur at lower energy with respect to the parent compound 3, which translates into 10 and 6 nm shifts of the UV and vis bands, respectively. Moreover, the 390 nm shoulder observed in the diffuse reflectance spectrum of 3 undergoes likewise a red shift upon silver binding ( Figure S6, Table S3).

Electrochemical Properties
Cyclic voltammetry (CV) was employed to study the electrochemical behavior of 3,6-diphenyl-1,2,4,5-tetrazine 1, the bis(tetrazine) 3 and the complex {[Ag(3)][BF 4 ]} ∞ (5). The CV experiments were carried out in CH 2 Cl 2 containing 0.1 M tetra(n-ethyl)ammonium tetrafluoroborate (TEABF 4 ), see Figure 9, Figure S7, and Table S4, as supporting electrolyte. The voltammogram of 1 displays a single reversible reduction wave (E 1/2 = −0.85 V/SCE) corresponding to the formation of the tetrazine radical anion [22,23], whereas bis(tetrazine) 3 undergoes two successive one-electron reductions at E 1/2 = −0.83 and −0.93 V/SCE. According to theoretical approaches [24], molecules undergoing two successive reversible one-electron transfers centered on chemically equivalent, and fully independent redox sites (such as, for example, bis(ferrocenyl)porphyrins) [25] give rise to a 35.6 mV shift between each successive formal half-wave potentials. The resulting CV curve for each one-electron transfer should display the same 60 mV peak-to-peak separation [24,26]. Herein, for 3, the potential gap between the successive one-electron reduction reactions is equal to 100 mV. This feature indicates that the redox sites, albeit chemically equivalent, are clearly not electronically independent from each other. This is attributed to a significant electronic communication that occurs between the two monomeric tetrazine moieties. The exact pathway for this electronic communication still remains unclear because disruption of aromaticity is expected to be rather significant when both heteroaromatic rings depart from a stacked, parallel arrangement, as observed in the solid state upon complex formation. The Ag(I) complex 5 is irreversibly reduced at E pc = 0.70 V/SCE. At this potential, silver(I) is reduced into Ag(0), which plates out at the electrode surface. The second reduction peak at E pc ca. −0.84 V/SCE is close to the value found for 3, which confirms the dissociation of the Ag(I) complex. On the backward scan, a clear anodic stripping peak is observed at E pa = 1.08 V/SCE for the dissolution of silver metal that was deposited on the electrode surface during the forward scan (Figure 7).
Herein, for 3, the potential gap between the successive one-electron reduction reactions is equal to 100 mV. This feature indicates that the redox sites, albeit chemically equivalent, are clearly not electronically independent from each other. This is attributed to a significant electronic communication that occurs between the two monomeric tetrazine moieties. The exact pathway for this electronic communication still remains unclear because disruption of aromaticity is expected to be rather significant when both heteroaromatic rings depart from a stacked, parallel arrangement, as observed in the solid state upon complex formation. The Ag(I) complex 5 is irreversibly reduced at Epc = 0.70 V/SCE. At this potential, silver(I) is reduced into Ag(0), which plates out at the electrode surface. The second reduction peak at Epc ca. −0.84 V/SCE is close to the value found for 3, which confirms the dissociation of the Ag(I) complex. On the backward scan, a clear anodic stripping peak is observed at Epa = 1.08 V/SCE for the dissolution of silver metal that was deposited on the electrode surface during the forward scan (Figure 7).

Host-Guest Behavior of Tweezer 3 with Silver Salts
The coordination of the cation Ag + to tweezer 3 was confirmed both in the solid state and in solution. Considering the strong ion-pairing effect revealed by NMR spectroscopy for complexes 4-8 and the detection by mass spectrometry of dinuclear silver species, we investigated by 1 H-NMR the speciation of Ag + in the presence of tweezer 3; this in order to establish the stoichiometry of the formed complexes and their stability. We selected a binary solvent mixture (CDCl3/CD3OD 2:1 v/v) for solubility reasons, as it enabled us to dissolve the silver salt precursor in a fairly large concentration range (2.5-34.9 mM), but also the tweezer 3 (in 25.6 mM), and the resulting coordination complexes.
Upon successive additions of an AgX salt solution to a solution of 3, all 1 H-NMR resonances undergo progressive shifts, suggesting a fast exchange between the free and

Host-Guest Behavior of Tweezer 3 with Silver Salts
The coordination of the cation Ag + to tweezer 3 was confirmed both in the solid state and in solution. Considering the strong ion-pairing effect revealed by NMR spectroscopy for complexes 4-8 and the detection by mass spectrometry of dinuclear silver species, we investigated by 1 H-NMR the speciation of Ag + in the presence of tweezer 3; this in order to establish the stoichiometry of the formed complexes and their stability. We selected a binary solvent mixture (CDCl 3 /CD 3 OD 2:1 v/v) for solubility reasons, as it enabled us to dissolve the silver salt precursor in a fairly large concentration range (2.5-34.9 mM), but also the tweezer 3 (in 25.6 mM), and the resulting coordination complexes.
Upon successive additions of an AgX salt solution to a solution of 3, all 1 H-NMR resonances undergo progressive shifts, suggesting a fast exchange between the free and complexed Tz at the typical NMR time-scale ( Figure 10, Table S5, AgNTf 2 case). As previously mentioned, H 9 protons, which are located in α-position with respect to the central C-C bond connecting the two phenyl-tetrazine moieties, experienced the largest upfield shift during the titration until reaching a limiting chemical shift value after addition of 3.0 equivalents of Ag + . Accordingly, we selected that particular signal to monitor the metal uptake process. Notably, solution equilibrium studies were restricted to the sole AgNTf 2 /3 system for which no precipitation was observed up to the addition of 3.0 equivalents of Ag + . Indeed, the occurrence of precipitates upon addition of 0.6 to 1.0 equivalent of the various silver salts AgPF 6 , AgBF 4 , AgSbF 6 , AgClO 4 , and Ag(OTf) hampered a quantitative investigation. The formation of precipitates that turned out to be insoluble in all common solvents suggests that the soluble molecular complexes formed at the early stage of the titration progressively evolved into oligomeric or polymeric species. This is a common feature in supramolecular chemistry for labile host-guest assemblies of weak stability [27][28][29][30][31][32].
complexed Tz at the typical NMR time-scale ( Figure 10, Table S5, AgNTf2 case). As previously mentioned, H9 protons, which are located in α-position with respect to the central C-C bond connecting the two phenyl-tetrazine moieties, experienced the largest up-field shift during the titration until reaching a limiting chemical shift value after addition of 3.0 equivalents of Ag + . Accordingly, we selected that particular signal to monitor the metal uptake process. Notably, solution equilibrium studies were restricted to the sole AgNTf2/3 system for which no precipitation was observed up to the addition of 3.0 equivalents of Ag + . Indeed, the occurrence of precipitates upon addition of 0.6 to 1.0 equivalent of the various silver salts AgPF6, AgBF4, AgSbF6, AgClO4, and Ag(OTf) hampered a quantitative investigation. The formation of precipitates that turned out to be insoluble in all common solvents suggests that the soluble molecular complexes formed at the early stage of the titration progressively evolved into oligomeric or polymeric species. This is a common feature in supramolecular chemistry for labile host-guest assemblies of weak stability [27][28][29][30][31][32]. To obtain a first insight into the stoichiometry of the formed species, an additional titration following the continuous variation method was performed [33][34][35][36][37][38]. According to the Job plot shown in Figure 11 [27][28][29][30][31][32], the maximum of the variation of the chemical shift of H9 protons is reached for x = [3]tot/([Ag]tot + [3]tot) = 0.5. This suggests a 1:1 Ag/3 stoichiometry for the prevailing complex in solution. However, the formation of additional species (for example, the 1:2 or 2:1 complexes) possibly present in lower amounts is often difficult to detect by this method, mainly because the position of the maximum and the overall shape of the graph are only slightly affected in those cases [39-41]. To obtain a first insight into the stoichiometry of the formed species, an additional titration following the continuous variation method was performed [33][34][35][36][37][38]. According to the Job plot shown in Figure 11 [27][28][29][30][31][32], the maximum of the variation of the chemical shift of H 9 protons is reached for x = [3] tot /([Ag] tot + [3] tot ) = 0.5. This suggests a 1:1 Ag/3 stoichiometry for the prevailing complex in solution. However, the formation of additional species (for example, the 1:2 or 2:1 complexes) possibly present in lower amounts is often difficult to detect by this method, mainly because the position of the maximum and the overall shape of the graph are only slightly affected in those cases [39][40][41].
Therefore, the chemical shifts of protons H 9 (δ obs ), which were determined upon incremental addition of 20 aliquots of an AgNTf 2 solution (up to 3 equivalents) to a 25.6 mM solution of 3 in CDCl 3 /CD 3 OD 2:1 v/v (Figure 10), were adjusted by nonlinear least-squares with the HypNMR 2008 program (Figure 12a) [42] for the set of complexation equilibria we propose below (Equation (1)).
Since the binding/dissociation kinetics is much faster than the NMR observation time, the experimental chemical shift can be expressed according to Equation (2), where αml, δml, and βml stand for the molar fraction, the intrinsic chemical shift, and the global stability constant of the MmLl species, respectively, while m and l correspond to the stoichiometric coefficients (by convention, β10 = β01 = 1). The HypNMR software iteratively refines the equilibrium constants βml for a given chemical model and calculates the intrinsic chemical shift δml for each species incorporating one or more molecules of ligand by solving at each iteration and for each analyzed sample the mass-balance equations.  Therefore, the chemical shifts of protons H9 (δobs), which were determined upon incremental addition of 20 aliquots of an AgNTf2 solution (up to 3 equivalents) to a 25.6 mM solution of 3 in CDCl3/CD3OD 2:1 v/v (Figure 10), were adjusted by nonlinear least-squares with the HypNMR 2008 program (Figure 12a) [42] for the set of complexation equilibria we propose below (Equation (1)).
Since the binding/dissociation kinetics is much faster than the NMR observation time, the experimental chemical shift can be expressed according to Equation (2), where αml, δml, and βml stand for the molar fraction, the intrinsic chemical shift, and the global stability constant of the MmLl species, respectively, while m and l correspond to the stoichiometric coefficients (by convention, β10 = β01 = 1). The HypNMR software iteratively refines the equilibrium constants βml for a given chemical model and calculates the intrinsic chemical shift δml for each species incorporating one or more molecules of ligand by solving at each iteration and for each analyzed sample the mass-balance equations.  Since the binding/dissociation kinetics is much faster than the NMR observation time, the experimental chemical shift can be expressed according to Equation (2), where α ml , δ ml , and β ml stand for the molar fraction, the intrinsic chemical shift, and the global stability constant of the M m L l species, respectively, while m and l correspond to the stoichiometric coefficients (by convention, β 10 = β 01 = 1). The HypNMR software iteratively refines the equilibrium constants β ml for a given chemical model and calculates the intrinsic chemical shift δ ml for each species incorporating one or more molecules of ligand by solving at each iteration and for each analyzed sample the mass-balance equations.
A fairly satisfactory fit (dashed line in Figure 12a) with a root mean square deviation (RMSD) of 0.012 ppm was obtained by adjusting the data with a single equilibrium involving the 1:1 complex (m = l = 1) as suggested by the Job plot. The value of the refined stability constant is log β 11 = 2.4 ± 0.1. The chemical shift calculated for 3 (δ 01 = 7.530 ± 0.006 ppm) deviates by 0.03 ppm from the experimental value measured for a pure sample (δ 01 = 7.560 ppm), while that refined for [Ag(3)] + is δ 11 = 7.064 ± 0.006 ppm. Accordingly, the 1:1 complex was formed only at 70% at the end of the titration, which corresponded to the addition of 3 equivalents of silver salt.
A second equilibrium was introduced in the chemical model that improve the goodness of fit. The refinement diverged when the binuclear [Ag 2 (3)] 2+ complex was considered. However, a significant improvement (RMSD = 0.003 ppm) could be achieved by including the [Ag(3) 2 ] + complex in the model (solid line in Figure 12a). The latter species fully compensates for the lack of fit (up to 0.03 ppm) seen for all data points corresponding to metalover-ligand molar concentration ratios below 1 equivalent with the one-equilibrium model. The optimization procedure converged for the following parameters: log β 11 = 1.72 ± 0.01, log β 12 = 3.30 ± 0.01, δ 01 = 7.561 ± 0.003 ppm, δ 11 = 6.875 ± 0.003 ppm, and δ 12 = 6.981 ± 0.003 ppm. A Job plot was simulated for this chemical model and the same concentrations as those used experimentally. The maximum of the concave curve was found at x = 0.55 instead of the theoretical value of 0.5 expected for the sole formation of a 1:1 complex. As already mentioned, the slight shape dissymmetry introduced by the [Ag(3) 2 ] + species makes its detection by the continuous variation method difficult [40,41]. This two-equilibria model is also supported by the intrinsic chemical shift refined for bis(tetrazine) 3, which perfectly matches with the experimental value. Moreover, the fact that the optimized δ 12 value (6.98 ppm) for [Ag(3) 2 ] + is ranging in-between those obtained for the free tweezer 3 (7.56 ppm) and the corresponding 1:1 silver complex (6.88 ppm) provides further support. Accordingly, an intermediate conformation of both tweezer molecules in [Ag(3) 2 ] + , located somewhere between the closed form of 3 stabilized by [C 2 N 4 ]···[C 2 N 4 ] intramolecular π-π interactions [3], and the more open structure found in [Ag(3)] + , might be postulated.
The low value computed for the first (log β 11 = 1.72 ± 0.01) and second (log β 12 = 3.30 ± 0.01) global equilibrium constants clearly indicates that the bis(tetrazine) 3 forms extremely weak complexes with the soft Ag + cation in the selected binary solvent mixture (CDCl 3 /CD 3 OD 2:1 v/v). The low affinity of 3 is best illustrated by the distribution diagram corresponding to the titration conditions (Figure 12b).
Although such low values of equilibrium constants have to be interpreted with some caution, the difference ∆ 12 = log K 11 − log K 12 = 0.14 ± 0.04 is lower than the statistical factor (log 4 ≈ 0.60) expected for a stepwise ligand uptake by a bis-coordinated silver center [43], as it was found in the crystal structure of {[Ag(3)][SbF 6 ]·CHCl 3 } ∞ . Therefore, the cooperative formation of [Ag(3) 2 ] + in the early stages of the titration, when the ligand is still in excess over Ag + , rules out strong steric interactions between the molecules of bis(tetrazine), as their occurrence would contribute to an increase in ∆ 12 . The analysis of non-covalent interactions in [Ag(3) 2 ] + also indicates that the dispersion interaction between the two ligands may compensate for the steric repulsion ( Figure S8). It is, however, likely that the counter anion contributes to stabilizing the species in solution by an ion-pairing mechanism (electrostatic outer-sphere interactions) in weakly polar media, as suggested by the different titration profiles (chemical shift variations) for the H 9,9' proton pair that could be observed for the other considered silver salts.

General Synthetic Conditions
The bis(tetrazine) 3 was prepared following our previously reported procedure [3]. All reagents and precursors (AgPF 6 , AgBF 4 , AgSbF 6 , AgClO 4 , AgNTf 2 , and Ag(OTf)) were purchased from commercial suppliers and used without purifications. All reactions were performed under an atmosphere of dry argon in Schlenk tubes, glass vials, or in a microwave reaction vessel. Microwave heating was carried out using a CEM Discover microwave reactor equipped with an infrared temperature sensor. The microwave-assisted reactions were run in closed vessels under magnetic stirring. All solvents were dried and deoxygenated prior to use.

Physical Methods
All NMR spectra ( 1 H, 13 C, 19 F, and 31 P) were recorded on Bruker Avance III spectrometers operating at 400, 500, or 600 MHz in solution with the deuterated solvent as the lock and the residual solvent or TMS as the internal reference. Chemical shifts are reported in ppm relative to deuterated solvent, and coupling constants J are given in Hz. High-resolution Mass Spectra (HRMS) were obtained on a Thermo LTQ-Orbitrap XL with an ESI source. Flash chromatography was performed on silica gel (230-400 mesh).

Titration Procedures
1 H-NMR titrations were performed using deuterated solvents. All ligand and metalsalt solutions used for titrations were prepared in volumetric flasks by dissolving weighted amounts (±0.1 mg) of each reagent in the mixed solvent CDCl 3 /CD 3 OH 2:1 v/v. Additions were made with 100 µL Hamilton micro-syringes.

Determination of the Stoichiometry of the Complexes
The stoichiometry of the complexes formed between the bis(tetrazine) and the metal salt was determined by applying the continuous variation method, also known as the Job method [27], to a 1 H-NMR titration. Two equimolar stock solutions of tetrazine ([3] tot = 8.6 mM) and of the metal salt ([Ag] tot = 8.6 mM) in the binary solvent mixture were mixed in different proportions varying from 0 to 100% in 10% steps, in order to keep the total concentration corresponding to the sum of the two-species (tetrazine + metal salt) constant at 8.6 mM. The Job plot is a graphical representation of the x∆δ product as a function of x = [3] tot /([Ag] tot + [3] tot ), where ∆δ = δ free ligand − δ obs .

Determination of the Association Constants
To determine the association constants between the tetrazine 3 and Ag + , a second 1 H-NMR titration was carried out by incremental addition of aliquots of a silver salt stock solution to a solution of 3 contained in a 5-mm NMR tube. The initial volume was 0.5 mL.
After each metal addition, the sample was mixed by vigorous shaking, and the 1 H-NMR spectrum was recorded on a 400 MHz spectrometer at a fixed temperature of 298 K. Each spectrum was carefully calibrated using the CD 3 OD signal as reference (δ = 3.31 ppm). The variations of the chemical shift of protons H 9,9' as a function of the total ligand and metal concentrations were processed with the nonlinear least-squares program HypNMR 2008 (Protonic Software) [42]. Since a unitary weighing scheme (ω i = 1) was selected, the returned RMSD of the fit can be directly compared to the measurement precision of chemical shifts (typically, σ ≈ 0.002 ppm). Hence, in the absence of systematic errors and for a correct chemical model (i.e., randomly distributed residuals), the expectation value of the RMSD should equal σ.

UV-Visible Absorption and Diffuse Reflectance Spectrophotometry
UV-Vis absorption spectra in solution were recorded at room temperature with a Cary 50 (Varian) spectrophotometer using a quartz cell (Hellma) of 1 cm optical path length. Concentrations were about 10 -5 M (between 230 and 430 nm) or 10 -3 M (between 430 and 630 nm). Diffuse reflectance spectra of pure microcrystalline or amorphous compounds diluted in dry barium sulfate (avocado, precipitated, >99%) were acquired between 200 and 800 nm on a CARY 5000 (Agilent) UV-Vis-NIR spectrophotometer fitted with a Praying Mantis TM accessory (Harrick). The baseline was recorded with BaSO 4 . Corrected reflectance data were converted to f (R) values using the Kubelka-Munk function expressed as f (R) = (1 -R) 2 /2R, where R stands for the reflectance.

Cyclic Voltammetry Measurements
All manipulations were performed using Schlenk techniques in an oxygen-free atmosphere of dry argon at room temperature (θ = 20 ± 3 • C). The supporting electrolyte (TEABF 4 ) was degassed under vacuum before use and then dissolved in CH 2 Cl 2 to a concentration of 0.1 M. Voltammetric analyses were carried out in a standard three-electrode cell with an Autolab PGSTAT 302 N potentiostat, connected to an interfaced computer and controlled by the Electrochemistry Nova software (v. 1.11). The reference electrode was a KCl saturated calomel electrode (SCE) separated from the analyzed solution by a sintered-glass salt bridge filled with the background electrolyte solution. The auxiliary electrode was a platinum foil separated from the analyzed solution by a sintered-glass salt bridge filled with the background electrolyte solution. For all voltammetric measurements, the working electrode was a platinum (GC) electrode (Ø = 2 mm) that was polished with a diamond suspension. In these conditions (0.1 M TEABF 4 in CH 2 Cl 2 ), the formal half-wave potential for the ferrocenium/ferrocene couple was +0.46 V vs. SCE.

Computational Details
Quantum mechanics calculations were performed with the Gaussian16 software package [44]. Energy and forces were computed by density functional theory with the hybrid B3PW91 exchange-correlation functional. In addition, dispersion effects have been considered using the D3 dispersion correction suggested by Grimme et al. with Becke-Johnson damping [45]. This computational level will be denoted by B3PW91-D3. The solvent effects were modeled by means of the continuum model as implemented in Gaussian. Geometries were optimized and characterized with the aug-cc-pVTZ basis set [46]. The non-covalent interactions were visualized with the Independent Gradient Model [47,48] that can be seen as an extension of the NCI method of Yang et al. [49]. The UVvisible transitions were computed using the B3PW91-D3 functional with the 6-31+G(d,p) basis sets. The orbital transitions of selected excited states were characterized by using the natural transition orbital (NTO) method with an isovalue of 0.04 au [50].

Conclusions
New supramolecular polymeric architectures were achieved by self-assembly of silver(I) cations with the bis(tetrazine) 3. The intrinsic tweezer-like flexibility of 3 has been confirmed in self-assembly processes with transition metal ions: Ag I , Cu I , and Pd II . The resulting discrete coordination polymers were identified by the combination of physicochemical analyses in solution and XRD identification in the solid state. The supramolecular polymeric structure of the crystalline complexes {[Ag (3) 6 ]} ∞ were fully described, and although we could not further confirm by X-ray diffraction the structures of the other metal complexes, the combination of analytical characterization methods clearly indicated similar structures in solution and in the solid state for all the obtained host-guest silver compounds. The physical measurements (NMR, XRD, UV, and electrochemistry) evidenced weak interactions of the ligand with the metal cation, as well as within the polymeric backbone. These interactions were confirmed by the variations in chemical shift values, absorption maxima, electrochemical reduction potentials, as well as low values of association constants obtained by host-guest titrations. These trends were confirmed by theoretical calculations, which revealed weak non-covalent interactions in the new coordination polymers crystallized from the bis(tetrazine) 3.