Next Article in Journal
Interaction of Vinyl-Type Carbocations, C3H5+ and C4H7+ with Molecules of Water, Alcohols, and Acetone
Next Article in Special Issue
Molecular and Electronic Structures of Macrocyclic Compounds Formed at Template Synthesis in the M(II)—Thiocarbohydrazide—Diacetyl Triple Systems: A Quantum-Chemical Analysis by DFT Methods
Previous Article in Journal
Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways
Previous Article in Special Issue
Calix[6]arene-Based [3]Rotaxanes as Prototypes for the Template Synthesis of Molecular Capsules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles

by
Alexander J. Menke
,
Nicholas C. Henderson
,
Lola C. Kouretas
,
Anne N. Estenson
,
Benjamin G. Janesko
* and
Eric E. Simanek
*
Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, TX 76109, USA
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(3), 1144; https://doi.org/10.3390/molecules28031144
Submission received: 15 December 2022 / Revised: 11 January 2023 / Accepted: 15 January 2023 / Published: 23 January 2023

Abstract

:
In the absence of preorganization, macrocyclization reactions are often plagued by oligomeric and polymeric side products. Here, a network of hydrogen bonds was identified as the basis for quantitative yields of macrocycles derived from the dimerization of monomers. Oligomers and polymers were not observed. Macrocyclization, the result of the formation of two hydrazones, was hypothesized to proceed in two steps. After condensation to yield the monohydrazone, a network of hydrogen bonds formed to preorganize the terminal acetal and hydrazine groups for cyclization. Experimental evidence for preorganization derived from macrocycles and acyclic models. Solution NMR spectroscopy and single-crystal X-ray diffraction revealed that the macrocycles isolated from the cyclization reaction were protonated twice. These protons contributed to an intramolecular network of hydrogen bonds that engaged distant carbonyl groups to realize a long-range order. DFT calculations showed that this network of hydrogen bonds contributed 8.7 kcal/mol to stability. Acyclic models recapitulated this network in solution. Condensation of an acetal and a triazinyl hydrazine, which adopted a number of conformational isomers, yielded a hydrazone that adopted a favored rotamer conformation in solution. The critical hydrogen-bonded proton was also evident. DFT calculations of acyclic models showed that the rotamers were isoenergetic when deprotonated. Upon protonation, however, energies diverged with one low-energy rotamer adopting the conformation observed in the macrocycle. This conformation anchored the network of hydrogen bonds of the intermediate. Computation revealed that the hydrogen-bonded network in the acyclic intermediate contributed up to 14 kcal/mol of stability and preorganized the acetal and hydrazine for cyclization.

1. Introduction

The cyclic, undecapeptide cyclosporin, is an archetype for the design of bioactive macrocyclic molecules. Comprising a ring of 33 atoms, cyclosporin serves as a foundation for building out the so-called BRo5 (beyond Lipinski’s Rule of 5) principles [1,2,3] that define drug-like characteristics for large molecules [4]. These design principles address the roles of conformational flexibility (so-called chameleonicity) [5], N-alkylation [6], and solubility [7,8]. Cyclosporin also serves to anchor new computational modeling strategies [9]. Given the current interests in these topics, cyclosporin’s historic role in shaping synthetic strategy receives less attention. It served as a model for evaluating methodology as well as identifying sites for cyclization [10]. In addition, and more germane to this work, the synthesis of cyclosporin points to the role that hydrogen bonding networks can play in retrosynthetic analysis. Although executed at high dilution (10−3–10−4 M), the high yield (62%) obtained for cyclization was attributed to preorganization of the linear peptide into a product-like conformation.
In facilitating macrocycle synthesis, preorganization serves as a powerful strategy for influencing product distributions [11,12]. Oftentimes, this distribution can be defined as discrete versus polymeric products. An early example of the use of hydrogen bonding to promote macrocyclization over polymerization involved the condensation of a dialdehyde with a diamine. Intramolecular hydrogen bonds preorganized the aldehyde groups in a cyclic arrangement. In their absence, only a polymeric material was observed [13]. A similar observation was made with the reaction of a diamine with an activated diacid. Hydrogen bonding has also directed dimerization by olefin metathesis [14].
Preorganization based on a network of hydrogen bonds appears to be relevant in our 24-atom macrocycles (Scheme 1). Treatment of the bifunctional monomers with trifluoroacetic acid leads to quantitative dimerization. Upon the acid-catalyzed deprotection, the hydrazine reacts with the acetal to yield the first hydrazone, our so-called “hypothetical preorganized intermediate”. We posited that protonation of the triazines contributes to a network of hydrogen bonds that templates the cyclization reaction by localizing the acetal and hydrazine groups to promote dimer formation. Experimental and computational evidence both supported this thesis. Indeed, the lack of evidence for oligomeric or polymeric species was surprising because such products are commonly observed in dynamic covalent chemistry [15,16,17,18,19].
Experimental evidence for templating derived from (1) the quantitative nature of the reaction across a variety of substrates, (2) the observation of the proton template and (3) network of hydrogen bonds in the products, and (4) the reproduction of the templating effect within acyclic models. Computation provided (1) an evaluation of ground-state energies of isomeric models, both neutral and protonated, that explained the site of protonation and shift in isomer distributions and (2) the energetic contribution of hydrogen bonding to stability. Each topic is elaborated upon in the following paragraphs. The molecules probed in this study appear in Chart 1.

2. Results and Discussion

Macrocyclization is quantitative. The current method employed for macrocyclization proceeded in quantitative yields to the limits of detection by NMR spectroscopy (>95%). That is, upon evaporation of the solvent, no workup was necessary. However, even before the adoption of the current method, the yields for cyclization were high. In the initial report, macrocycles were subjected to purification by centrifugation and chromatography [20]. Yields averaged 90%. In a second study of the ring size, centrifugation yielded a pellet that was washed with dichloromethane [21]. Yields ranged from 83%–100%. The observation of a preponderance of dimer over oligomeric or polymeric products was consistent with templating.
Protonation of the triazine ring was consistently observed in the macrocyclic products. The reaction conditions employed for macrocyclization, 50% trifluoroacetic acid in dichloromethane, provided the source of protons as well as the trifluoroacetate counterion. Protonation of the triazine ring was consistently observed in the macrocycles prepared to date in both solution and the solid state. Protonation occurred at the site opposite the auxiliary group, dimethylamine or morpholine, as illustrated in Scheme 1. Figure 1 shows the crystal structure of a valine macrocycle. The proton was not assigned. It was identified based on the observed electron density.
The hydrogen bonding network was present in the macrocycle. Figure 1 also reveals the network of hydrogen bonds in the products matched that proposed in the preorganized intermediate. The amide carbonyl of one subunit engaged in hydrogen bonding with the other subunit. The carbonyl formed two hydrogen bonds, one with the protonated triazine and the other with α-NH of the amino acid.
Hydrogen bonding was also apparent in the 1H NMR spectra of these macrocycles. Figure 2 shows the appearance of well-resolved resonances in the downfield “fingerprint” region of the NMR spectrum. Both the chemical shift and shape were indicative of hydrogen bonding interactions. While the intermediates showed numerous broad resonances between 6 and 9 ppm, these shifted to sharp signals appearing between 7 and 13 ppm for the macrocycles [22]. Analysis of 2-D ROESY spectra established that the solution structure matched that observed in the solid state. The threonine and isoleucine macrocycles also adopted the same folded conformation in solution with the same intramolecular network of hydrogen bonds.
Rotamers were observed in the NMR spectra of intermediates but not in the macrocycle. Substituted triazines existed in an equilibrium of rotational isomers, the so-called rotamers, due to a hindered rotation about the triazine-N bond [23,24,25]. Experimental and computational models placed the barrier of rotation at 15–16 kcal/mol [23]. These rotamers are illustrated in Chart 2 using an ethylamide intermediate. The geometry is indicated by the naming convention (I–IV), which is conserved throughout this report. Evidence for rotamers appeared in the 1H NMR spectrum of the carboxylic acid intermediates and the monomer, manifesting as multiple resonances in the NMR spectrum (see the Supplementary Materials for examples). For example, the α-proton of the amino acid often produced three of the four independent resonances at ~4.5 ppm, each corresponding to a different rotamer. The exchangeable NH protons region between 6 ppm and 9 ppm showed a large number of resonances as well.
In stark contrast, however, only a single rotamer, I, was observed in the macrocycles. This situation could not be dismissed readily as the result of constraints imposed by cyclization. Computational models using well-tempered metadynamics suggested that these macrocycles can adopt multiple conformations [26]. The physical manipulation of Corey–Pauling–Koltun models corroborated these conclusions. Reports of the use of substituted triazines, melamines, in molecular recognition led to a second explanation: the rotamer equilibrium can be influenced by noncovalent interactions. From Figure 2, it is evident that protonation was sufficient to bias the rotamer equilibrium. Computation offered insight into these energetic costs.
Computational models of ground-state energies showed that the rotamers were isoenergetic when neutral. DFT calculations were used to predict relative ground-state energies of models, shown in Figure 3. These models were intended to reproduce the dimethylamine auxiliary group (with fidelity), the amino acid (with methylamine that apes the α-NH and α-carbon), and the hydrazone (limited to the hydrazone of acetaldehyde). The negligible differences in the ground-state energies of the rotamers were consistent with the expectation that bias to one species could be realized with weak, noncovalent interactions—including one or more hydrogen bonds.
Computational models revealed energetic preferences for specific rotamers upon protonation. To elaborate the models to account for protonation, isodesmic reactions were used, as illustrated in Scheme 2.
Figure 4 shows the relative energy values (in kcal/mol) for protonation at different sites for each of the four rotamers of the model. Only the trans-hydrazone was modeled. The cis-hydrazone was not observed experimentally. To validate the computational model, we analyzed the protonation events at the possible sites of triaminotriazines (melamines) bearing from 1 to 6 N-methyl groups and compared the predicted energies to measured pKa values [27,28].
The data revealed a few trends. First, protonation of the hydrazone was significantly disfavored (red) over other sites, with an average energy of 9.7 kcal/mol above the minimum. Second, the most favorable site of protonation (blue) was on a triazine nitrogen adjacent to the hydrazone and oriented such that the lone pair of the hydrazone was available to contribute electron density. The calculated energies for these sites across the rotamer population ranged from the 0.0 to 0.9 kcal/mol above the minimum. We considered these rotamers to be isoenergetic at this level of theory. This site agreed with chemical intuition and with the solution-state and solid-state structures as well. The site opposite the hydrazone (green) also represented a favorable site for protonation, with an average energy of 1.6 kcal/mol above the minimum. The final site for protonation (orange) put the proton in proximity of the hydrazone but without the advantage of the hydrazone’s lone pair of electrons. These energies ranged from 2.7–3.2 kcal/mol above the minimum, a difference that was significant compared to the lowest-energy sites (blue).
Experimental data derived from solution and solid-phase structure analyses of macrocycles using NMR and X-ray crystallography, respectively, showed that only rotamer I was observed. This protonation event occurred at 0.9 kcal/mol higher energy than the minimum. Given the negligible energetic penalty, we turned our attention to whether the hydrogen bonding network provided a greater amount of stabilizing energy.
Increasing the number of hydrogen bonds increasingly stabilized the preorganized intermediate. Using M06-2X/6-31+G(d) DFT calculations and a continuum model for solvent, several optimizations of the extended monohydrazone converged to partially folded states. Figure 5 shows representative structures. The low-energy structure of the neutral intermediate was stabilized by a single hydrogen bond, conveying a 6.4 kcal/mol benefit over the extended state. Enthalpic penalties appeared to preclude the formation of two hydrogen bonds in the lowest-energy structure.
Two different monoprotonated intermediates were modeled. Both contributed more stabilizing energy than the neutral molecule. Protonating the distal triazine afforded two hydrogen bonds and 10.2 kcal/mol of stabilization. Protonating the triazine closer to the hydrazone afforded two hydrogen bonds and 12.4 kcal/mol of stabilization. The doubly protonated intermediate yielded four hydrogen bonds and the greatest stabilization energy, 14.2 kcal/mol. A plot of folding energy vs. number of protons was reasonably linear (R2 = 0.90) with a slope of 3.9 kcal/mol of folding energy for each proton added. The data supported the belief that templating and dimer selectivity should be favored at low pH.
Computational models of the products revealed significant stabilization energy from hydrogen bonding. Solution-phase studies of these macrocycles corroborated the hydrogen bonding network illustrated in Scheme 1. To determine the energetic contributions of these hydrogen bonds in the product macrocycles, we modeled neutral, monoprotonated, and doubly protonated macrocycles using a solvent continuum. To simplify calculations, the sidechain was truncated to a single methyl group (an alanine equivalent). The calculated structures are shown in Figure 6.
In contrast to the glycine macrocycle, which oriented the carbonyl groups outward (out) to form bridging contacts between unit cells in the crystal structure, the valine and isoleucine macrocycles oriented the carbonyl groups inward (in) to form four productive hydrogen bonding interactions with the protonated triazines. In solution, all macrocycles favored the inward conformation. In the neutral molecule, orienting the carbonyls inward created unfavorable interactions, leading to a penalty of 4.7 kcal/mol when compared with the outward orientation. When doubly protonated, however, an inward orientation of the carbonyls was favored over the outward orientation by 8.7 kcal/mol. The energy difference between the carbonyl orientations for a monoprotonated macrocycle was within 1 kcal/mol at this level of theory.
Modeling noncovalent interactions in acyclic systems supported templating and preorganization. To conclude these studies, we looked to recapitulate the preorganization experimentally in acyclic systems. The model design was guided by crystal structures and intuition. If preorganization was important, we hypothesized that the rotamer equilibrium would shift to the rotamer that facilitated the formation of the hydrogen bond. The NMR spectra of the ethylamide model systems substantiated this belief (Figure 2). In all cases, the models favored one rotamer and the spectra were well resolved, a stark contrast to the spectra of the precursors (available in the supporting information).
The spectra also revealed some global trends. First, and indicated with a green asterisk, the model compounds showed markedly similar chemical shifts across a range of resonances. This behavior was distinct from that of the macrocycles that subdivided into two populations, indicated with the yellow- and blue-shaded regions of the spectra. Second, as indicated by the arrows, similar chemical shift differences were observed between the models and the corresponding macrocycles. These changes were pronounced for C-NH and α-NH, which engaged in hydrogen bonding interactions either indirectly or directly. The negligible changes in the chemical shifts of H+ and NNH suggested that those environments were very similar in the model and macrocycle.

3. Experimental

Computational methods. The calculations used the Gaussian 16 electronic structure package [29]. The calculations used density functional theory (DFT) [30] at the M062X/6-311++(2d,2p) [31,32,33] level using the solvation model based on density (SMD) to model the water solvent [34]. The calculations used ωB97X-D/6-31+G(d,p) DFT. All geometries were optimized in SMD solvent. The pKa calculations were validated against a set of experimental melamine derivatives pKa [27]. Additional details are available in the Supplementary Materials.
Specific experimental methods and synthesis. The details of the synthesis of these compounds and the characterization data can be found in the Supplementary Materials.
NMR Spectroscopy. Room temperature 1H NMR spectra were recorded on a 400 MHz Bruker Avance spectrometer. Chemical shifts for the 1H NMR spectra (in parts per million) were referenced to a corresponding solvent resonance (e.g., DMSO-d6, δ = 2.52 ppm). The 13C{1H} NMR spectra were recorded on the same 400 MHz Bruker spectrometer referenced to the corresponding solvent resonance. All 2D spectra were taken on the 400 MHz Bruker Avance relative to corresponding solvent resonances. Low-temperature spectra were acquired on a 500 MHz Varian NMR spectrometer at the University of North Texas in Denton. The identification of NMR signals were as follows: s = singlet, d = doublet, t = triplet, dd = doublet of doublet, m = multiplet, etc. NMR solvents were deuterated and purchased as a bottle or ampule. Structural assignments were made with and supported by two-dimensional NMR experiments (COSY, HSQC, and rOesy).
General chemistry. Flash chromatography experiments were carried out on silica gel with a porosity of 60Å, particle size of 50–63 m, surface area of 500–600 m2/g, a bulk density of 0.4 g/mL, and a pH range of 6.5–7.5. Dichloromethane/methanol was used as the eluent for chromatographic purification. Thin-layer chromatography experiments were carried out in sealed chambers and visualized with UV or submersion in ninhydrin (1.5 g ninhydrin in 100 mL of n-butanol and 3.0 mL acetic acid) followed by heating. Excess solvents were removed via rotary evaporation on a Buchi Rotavapor RII with a Welch Self-Cleaning Dry Vacuum System. All workup and purification procedures were carried out with reagent-grade solvents under ambient atmosphere.

4. Conclusions

Experimental and computational evidence support a model for templated cyclization. Computation suggests that hydrogen bonding contributes 14.2 kcal/mol of energy to the preorganized intermediate when compared with the extended model. This network of hydrogen bonds is preserved in the structures of the products. The magnitude of the contribution is consistent with the observed quantitative yields. Furthering the argument for preorganization, protonation reduces the number of rotamers sampled in the acyclic models. This preference is unsurprising given that the cost for adopting the desired rotamer in the protonated state was negligible in comparison.
The model that emerges to anchor our intuition, as advanced in Scheme 1, is validated. That is, upon condensation to form the first hydrazone, a hydrogen bonding network forms within the molecule. The entropic costs predicted to be associated with this conformational change are hypothesized to be offset by strong hydrogen bonds. These quantitative yields appear to be independent of the choice of the amino acid. Further studies will evaluate the generality of this prediction.

Supplementary Materials

Supporting information including spectra and synthetic details can be downloaded at https://www.mdpi.com/article/10.3390/molecules28031144/s1. Experimental Details; Chart S1: Compounds described in the supporting material; Figure S1: DFT-computed relative energies of protonated triazine models; Figure S2: Coding the computational data onto the rotamer diagrams; Figure S3: Relative energies of protonated triazine models; Figure S4: Folding energy of the neutral (unprotonated) intermediate; Figure S5: Folding energy of the monoprotonated (terminal) intermediate; Figure S6: Folding energy of the monoprotonated (interior) intermediate; Figure S7: Folding energy of the diprotonated intermediate; Figure S8: Folding energy vs. number of protons; Figure S9: The 400 MHz 1H NMR spectrum of 1 in DMSO-d6; Figure S10: The 100 MHz 13C{1H} NMR spectrum of 1 in DMSO-d6; Figure S11: The 400 MHz 1H NMR spectrum of 2 in DMSO-d6; Figure S12: The 100 MHz 13C{1H} NMR spectrum of 2 in DMSO-d6; Figure S13: The 400 MHz 1H NMR spectrum of 3 in DMSO-d6; Figure S14: The 100 MHz 13C{1H} NMR spectrum of 3 in DMSO-d6; Figure S15: The 400 MHz 1H NMR spectrum of GEA in DMSO-d6; Figure S16: The 100 MHz 13C{1H} NMR spectrum of GEA in DMSO-d6; Figure S17: The 400 MHz COSY of GEA in DMSO-d6; Figure S18: The 400 MHz 1H NMR spectrum of VEA in DMSO-d6; Figure S19: The 100 MHz 13C{1H} NMR spectrum of VEA in DMSO-d6; Figure S20: The 400 MHz COSY of VEA in DMSO-d6; Figure S21: The 400 MHz 1H NMR spectrum of IEA in DMSO-d6; Figure S22: The 100 MHz 13C{1H} NMR spectrum of IEA in DMSO-d6; Figure S23: The 400 MHz COSY of IEA in DMSO-d6; Figure S24: The 400 MHz 1H NMR spectrum of 4 in DMSO-d6; Figure S25: The 100 MHz 13C{1H} NMR spectrum of 4 in DMSO-d6; Figure S26: The 400 MHz 1H NMR of G-G in DMSO-d6; Figure S27: The 400 MHz 1H NMR of V-V in DMSO-d6; Figure S28. The 400 MHz 1H NMR of I-I in DMSO-d6.

Author Contributions

Conceptualization, E.E.S.; investigation, A.J.M., N.C.H., L.C.K., A.N.E. and B.G.J.; writing–original draft preparation, E.E.S.; writing–review and editing, A.J.M., B.G.J. and E.E.S.; supervision, A.J.M.; project administration, E.E.S.; funding acquisition, B.G.J. and E.E.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Robert A. Welch Foundation (EES; P-0008) and the National Institutes of Health (EES; R15GM135900).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. DeGoey, D.A.; Chen, H.J.; Cox, P.B.; Wendt, M.D. Beyond the Rule of 5: Lessons Learned from AbbVie’s Drugs and Compound Collection. J. Med. Chem. 2018, 61, 2636–2651. [Google Scholar] [CrossRef] [PubMed]
  2. Sebastiano, M.R.; Doak, B.C.; Backlund, M.; Poongavanam, V.; Over, B.; Ermondi, G.; Caron, G.; Matsson, P.; Kihlberg, J. Impact of Dynamically Exposed Polarity on Permeability and Solubility of Chameleonic Drugs Beyond the Rule of 5. J. Med. Chem. 2018, 61, 4189–4202. [Google Scholar] [CrossRef] [PubMed]
  3. Vinogradov, A.A.; Yin, Y.Z.; Suga, H. Macrocyclic Peptides as Drug Candidates: Recent Progress and Remaining Challenges. J. Am. Chem. Soc. 2019, 141, 4167–4181. [Google Scholar] [CrossRef] [PubMed]
  4. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
  5. Danelius, E.; Poongavanam, V.; Peintner, S.; Wieske, L.H.E.; Erdelyi, M.; Kihlberg, J. Solution Conformations Explain the Chameleonic Behaviour of Macrocyclic Drugs. Chem. Eur. J. 2020, 26, 5231–5244. [Google Scholar] [CrossRef] [PubMed]
  6. Sedmera, P.; Jegorov, A.; Buchta, M.; Cvak, L. 11-Demethylcyclosporins exhibit a single conformation in methanol and dimethylsulfoxide. J. Pept. Res. 2001, 58, 229–236. [Google Scholar] [CrossRef]
  7. Naylor, M.R.; Ly, A.M.; Handford, M.J.; Ramos, D.P.; Pye, C.R.; Furukawa, A.; Klein, V.G.; Noland, R.P.; Edmondson, Q.; Turmon, A.C.; et al. Lipophilic Permeability Efficiency Reconciles the Opposing Roles of Lipophilicity in Membrane Permeability and Aqueous Solubility. J. Med. Chem. 2018, 61, 11169–11182. [Google Scholar] [CrossRef]
  8. Avdeef, A.; Kansy, M. “Flexible-Acceptor” General Solubility Equation for beyond Rule of 5 Drugs. Mol. Pharm. 2020, 17, 3930–3940. [Google Scholar] [CrossRef]
  9. Ono, S.; Naylor, M.R.; Townsend, C.E.; Okumura, C.; Okada, O.; Lee, H.W.; Lokey, R.S. Cyclosporin A: Conformational Complexity and Chameleonicity. J. Chem. Inf. Model. 2021, 61, 5601–5613. [Google Scholar] [CrossRef]
  10. Wenger, R.M. Synthesis of Cyclosporine. Total Syntheses of Cyclosporin-A and Cyclosporin-H, 2 Fungal Metabolites Isolaterd from the species Tolypocladium-inflatum gams. Helv. Chim. Acta 1984, 67, 502–525. [Google Scholar] [CrossRef]
  11. Marti-Centelles, V.; Pandey, M.D.; Burguete, M.I.; Luis, S.V. Macrocyclization Reactions: The Importance of Conformational, Configurational, and Template-Induced Preorganization. Chem. Rev. 2015, 115, 8736–8834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Blankenstein, J.; Zhu, J.P. Conformation-directed macrocyclization reactions. Eur. J. Org. Chem. 2005, 10, 1949–1964. [Google Scholar] [CrossRef]
  13. Owston, P.G.; Peters, R.; Ramsammy, E.; Tasker, P.A.; Trotter, J. Non-template synthesis of N4 macrocyclic imine ligands with variable sizes—The importance of intramolecular hydrogen-bonding—x-ray crystal-structures of 3 macrocyclic and 2 open-chain ligands. J. Chem. Soc. Chem. Commun. 1980, 24, 1218–1220. [Google Scholar] [CrossRef]
  14. Clark, T.D.; Kobayashi, K.; Ghadiri, M.R. Covalent capture and stabilization of cylindrical beta-sheet peptide assemblies. Chem. Eur. J. 1999, 5, 782–792. [Google Scholar] [CrossRef]
  15. Lisowski, J. Imine- and Amine-Type Macrocycles Derived from Chiral Diamines and Aromatic Dialdehydes. Molecules 2022, 27, 4097. [Google Scholar] [CrossRef] [PubMed]
  16. Beeren, S.R.; Sanders, J.K.M. Ferrocene-amino acid macrocycles as hydrazone-based receptors for anions. Chem. Sci. 2011, 2, 1560–1567. [Google Scholar] [CrossRef]
  17. Taylor, C.M.; Kilah, N.L. Synthesis of 2 + 2 Schiff base macrocycles by a solvent templating strategy and halogen bonding directed assembly. J. Inclus. Phenom. Macro. Chem. Ournal Incl. Phenom. Macrocycl. Chem. 2022, 102, 543–555. [Google Scholar] [CrossRef]
  18. Jin, Y.H.; Wang, Q.; Taynton, P.; Zhang, W. Dynamic Covalent Chemistry Approaches toward Macrocycles, Molecular Cages, and Polymers. Acc. Chem. Res. 2014, 47, 1575–1586. [Google Scholar] [CrossRef]
  19. Jiao, T.Y.; Wu, G.C.; Zhang, Y.; Shen, L.B.; Lei, Y.; Wang, C.Y.; Fahrenbach, A.C.; Li, H. Self-Assembly in Water with N-Substituted Imines. Angew. Chem. Int. Ed. 2020, 59, 18350–18367. [Google Scholar] [CrossRef]
  20. Yepremyan, A.; Mehmood, A.; Asgari, P.; Janesko, B.G.; Simanek, E.E. Synthesis of Macrocycles Derived from Substituted Triazines. Chembiochem 2019, 20, 241–246. [Google Scholar] [CrossRef]
  21. Sharma, V.R.; Mehmood, A.; Janesko, B.G.; Simanek, E.E. Efficient syntheses of macrocycles ranging from 22–28 atoms through spontaneous dimerization to yield bis-hydrazones. RSC Adv. 2020, 10, 3217–3220. [Google Scholar] [CrossRef] [Green Version]
  22. Menke, A.J.; Gloor, C.J.; Claton, L.E.; Mekhail, M.A.; Pan, H.; Stewart, M.D.; Green, K.N.; Pavan, G.M.; Capelli, R.; Simanek, E.E. A Model for the Rapid Assessment of Solution-Structures for 24-Atom Macrocycles: The Impact of-Branched Amino Acids on Conformation. J. Org. Chem. 2022. under review. [Google Scholar]
  23. Katritzky, A.R.; Ghiviriga, I.; Steel, P.J.; Oniciu, D.C. Restricted rotations in 4,6-bis- and 2,4,6-tris-(N,N-dialkylamino)-s-triazines. JCS Perkins Trans. 2 1996, 3, 443–447. [Google Scholar] [CrossRef]
  24. Amm, M.; Platzer, N.; Guilhem, J.; Bouchet, J.P.; Volland, J.P. Structural and conformational study of substituted triazines by N-15 NMR and X-ray analysis. Magn. Res. Chem. 1998, 36, 587–596. [Google Scholar] [CrossRef]
  25. Birkett, H.E.; Harris, R.K.; Hodgkinson, P.; Carr, K.; Charlton, M.H.; Cherryman, J.C.; Chippendale, A.M.; Glover, R.P. NMR studies of exchange between triazine rotamers. Magn. Res. Chem. 2000, 38, 504–511. [Google Scholar] [CrossRef]
  26. Capelli, R.; Menke, A.J.; Pan, H.; Janesko, B.G.; Simanek, E.E.; Pavan, G.M. Well-Tempered Metadynamics Simulations Predict the Structural and Dynamic Properties of a Chiral 24-Atom Macrocycle in Solution. ACS Omega 2022, 7, 30291–30296. [Google Scholar] [CrossRef]
  27. List, M.; Puchinger, H.; Gabriel, H.; Monkowius, U.; Schwarzinger, C. N-Methylmelamines: Synthesis, Characterization, and Physical Properties. J. Org. Chem. 2016, 81, 4066–4075. [Google Scholar] [CrossRef]
  28. Jang, Y.H.; Hwang, S.; Chang, S.B.; Ku, J.; Chung, D.S. Acid Dissociation Constants of Melamine Derivatives from Density Functional Theory Calculations. J. Phys. Chem. A 2009, 113, 13036–13040. [Google Scholar] [CrossRef]
  29. Gaussian 16; Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2019.
  30. Parr, R.G.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford Univsrsity Press: Oxford, UK, 1989. [Google Scholar]
  31. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar]
  32. McLean, A.D.; Chandler, G.S. Contracted Gaussian-basis sets for molecular calculations. 1. 2nd row atoms, Z = 11–18. J. Chem. Phys. 1980, 72, 5639–5648. [Google Scholar] [CrossRef]
  33. Raghavachari, K.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-Consistent Molecular Orbital Methods. 20. Basis set for correlated wave-functions. J. Chem. Phys. 1980, 72, 650–654. [Google Scholar]
  34. Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal solvation model based on solute electron density and a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef] [PubMed]
Scheme 1. The dimerization reaction catalyzed by the addition of acid was hypothesized to proceed through a preorganized intermediate.
Scheme 1. The dimerization reaction catalyzed by the addition of acid was hypothesized to proceed through a preorganized intermediate.
Molecules 28 01144 sch001
Chart 1. Molecules described in this study. NMR labels are indicated. The X-ray structure used a morpholine derivative of the molecule shown, wherein the dimethylamine group was replaced to promote crystallization.
Chart 1. Molecules described in this study. NMR labels are indicated. The X-ray structure used a morpholine derivative of the molecule shown, wherein the dimethylamine group was replaced to promote crystallization.
Molecules 28 01144 ch001
Figure 1. The solid-state structure of the valine macrocycle. Arrows indicate hydrogen bonding interactions wherein the carbonyl from one subunit engages in hydrogen bonding with the protonated triazine of the other.
Figure 1. The solid-state structure of the valine macrocycle. Arrows indicate hydrogen bonding interactions wherein the carbonyl from one subunit engages in hydrogen bonding with the protonated triazine of the other.
Molecules 28 01144 g001
Figure 2. The 400 MHz 1H NMR spectra of macrocycles and acyclic models. Two subpopulations are indicated with yellow (V-V and both isomers of I-I) and blue boxes (G-G). Arrows indicate a common change in chemical shift between the macrocycle and its corresponding acyclic model. Elements shared across the acyclic intermediates are indicated with a green asterisk. See text for additional detail.
Figure 2. The 400 MHz 1H NMR spectra of macrocycles and acyclic models. Two subpopulations are indicated with yellow (V-V and both isomers of I-I) and blue boxes (G-G). Arrows indicate a common change in chemical shift between the macrocycle and its corresponding acyclic model. Elements shared across the acyclic intermediates are indicated with a green asterisk. See text for additional detail.
Molecules 28 01144 g002
Chart 2. Rotamers of a substituted triazine. The arrows indicate the orientation of hydrogen atoms, intended to guide the eye in distinguishing rotamers.
Chart 2. Rotamers of a substituted triazine. The arrows indicate the orientation of hydrogen atoms, intended to guide the eye in distinguishing rotamers.
Molecules 28 01144 ch002
Figure 3. Calculations of the ground-state energies of different rotamers.
Figure 3. Calculations of the ground-state energies of different rotamers.
Molecules 28 01144 g003
Scheme 2. Calculations of the ground-state energies of different rotamers.
Scheme 2. Calculations of the ground-state energies of different rotamers.
Molecules 28 01144 sch002
Figure 4. The relative energies of the protonated triazines in kcal/mol for the four rotamers arising from hindered rotation about the triazine-N bond. Color reflects the relative energetic differences with values reported in kcal/mol. See text for details. Dotted circles correspond to potential sites for protonation.
Figure 4. The relative energies of the protonated triazines in kcal/mol for the four rotamers arising from hindered rotation about the triazine-N bond. Color reflects the relative energetic differences with values reported in kcal/mol. See text for details. Dotted circles correspond to potential sites for protonation.
Molecules 28 01144 g004
Figure 5. The relative energies in the templated intermediates relative to (a) the linear extended structure including (b) the neutral intermediate, (c) one monoprotonated intermediate, (d) the other monoprotonated intermediate, and (e) the proposed intermediate. Color reflects the relative energetic differences with values reported in kcal/mol.
Figure 5. The relative energies in the templated intermediates relative to (a) the linear extended structure including (b) the neutral intermediate, (c) one monoprotonated intermediate, (d) the other monoprotonated intermediate, and (e) the proposed intermediate. Color reflects the relative energetic differences with values reported in kcal/mol.
Molecules 28 01144 g005
Figure 6. Computation revealed that intramolecular hydrogen bonding in the deprotonated macrocycle was stabilizing by 8.7 kcal/mol. Color reflects the relative energetic differences with values reported in kcal/mol.
Figure 6. Computation revealed that intramolecular hydrogen bonding in the deprotonated macrocycle was stabilizing by 8.7 kcal/mol. Color reflects the relative energetic differences with values reported in kcal/mol.
Molecules 28 01144 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Menke, A.J.; Henderson, N.C.; Kouretas, L.C.; Estenson, A.N.; Janesko, B.G.; Simanek, E.E. Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles. Molecules 2023, 28, 1144. https://doi.org/10.3390/molecules28031144

AMA Style

Menke AJ, Henderson NC, Kouretas LC, Estenson AN, Janesko BG, Simanek EE. Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles. Molecules. 2023; 28(3):1144. https://doi.org/10.3390/molecules28031144

Chicago/Turabian Style

Menke, Alexander J., Nicholas C. Henderson, Lola C. Kouretas, Anne N. Estenson, Benjamin G. Janesko, and Eric E. Simanek. 2023. "Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles" Molecules 28, no. 3: 1144. https://doi.org/10.3390/molecules28031144

APA Style

Menke, A. J., Henderson, N. C., Kouretas, L. C., Estenson, A. N., Janesko, B. G., & Simanek, E. E. (2023). Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles. Molecules, 28(3), 1144. https://doi.org/10.3390/molecules28031144

Article Metrics

Back to TopTop