Next Article in Journal
Acidity Quantification and Structure Analysis of Amide-AlCl3 Liquid Coordination Complexes for C4 Alkylation Catalysis
Next Article in Special Issue
Synthesis and Structural and Optical Behavior of Dehydrohelicene-Containing Polycyclic Compounds
Previous Article in Journal
In Silico Investigation on the Molecular Behavior and Structural Stability of the Rosette Nanotubes as the Drug Vehicles for Paclitaxel, an Anti-Cancer Drug
Previous Article in Special Issue
A Chiral Relay Race: Stereoselective Synthesis of Axially Chiral Biaryl Diketones through Ring-Opening of Optical Dihydrophenan-threne-9,10-diols
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Nature of the Rotational Energy Barrier of Atropisomeric Hydrazides

1
Department of Industrial Chemistry “Toso Montanari”, Alma Mater Studiorum University of Bologna, Via P. Gobetti 85, 40129 Bologna, Italy
2
Center for the Chemical Catalysis—C3, Alma Mater Studiorum University of Bologna, Via P. Gobetti 85, 40129 Bologna, Italy
3
National Interuniversity Consortium of Materials Science and Technology (INSTM), Bologna Research Unit, Via P. Gobetti 85, 40129 Bologna, Italy
4
Department of Chemistry and Technologies of Drug, Sapienza University of Rome, Piazzale A. Moro 5, 00185 Rome, Italy
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(23), 7856; https://doi.org/10.3390/molecules28237856
Submission received: 25 October 2023 / Revised: 20 November 2023 / Accepted: 25 November 2023 / Published: 29 November 2023

Abstract

:
N-N atropisomers represent a useful class of compounds that has recently received important attention from many research groups. This article presents an in-depth analysis of the energy barrier needed for the racemization process of atropoisomeric hydrazides, combining an experimental and computational approach. The focus is on examining how electronic and steric factors impact the racemization process. The results obtained indicate that the barrier observed during the racemization process mainly arises from an increase in the p-orbital character of the nitrogen atoms.

Graphical Abstract

1. Introduction

Atropisomers are conformationally stable rotamers that have been recently rediscovered for their important role in medicinal chemistry and chemical biology. Historically, biaryls have been considered representative examples of atropisomerism since they are largely employed as catalysts and ligands for asymmetric synthesis. In recent years, novel types of atropisomers with stereogenic axes different from the C-C single bond have received massive attention from the scientific community. N-N atropisomers signify the latest frontier in atroposelective synthesis; they are valuable for their useful biological activity and their significance as noteworthy “smart” materials [1]. In the past, this class of compounds has been rarely studied by a few research groups that have analyzed the conformational properties of different types of substituted hydrazines, quinazoline-diones, and aminocamphorimides featuring a hindered N-N single bond [2,3,4,5,6,7]. Nevertheless, it is only in the last two years that important examples of the enantioselective preparation of heteroaromatic and hydrazido N-N atropisomers have been reported using consolidated catalytic strategies [8,9,10,11,12,13,14]. In most cases, the high stability of the rotation along the N-N single bond has been determined and ascribed to the steric hindrance of the substituents surrounding the stereogenic axis. When chemists face the synthesis of novel atropisomeric architectures, the experimental determination of the rotational energy barrier gives fundamental information on their stereochemical stability. Most of the time, this value is unknown, and the synthesis can be attempted without a clear indication other than that deriving from the structural comparison with similar compounds previously reported. Also, computational methods, such as Density Functional Theory (DFT), can be complementarily used to study the rotational energy barrier, and in recent years, the use of the aforementioned methods has been employed multiple times, becoming a fundamental topic of scientific reports on atroposelective synthesis. In some cases, they represent the only way to estimate the rotational energy barrier, for example, when the experimental measurement requires high temperatures, which can cause the product to decompose. The robustness of these DFT methods, which are based on the agreement between the computed rotational barrier and the experimental one, makes them a solid tool for the prediction of the rotational energy barrier of plausible atropisomeric compounds. In this way, given a target molecule, it is possible to estimate the energy barrier with a high degree of fidelity and decide whether to attempt the synthesis or not. In this article, the key role that computational calculations have in the prediction of the rotational energy barrier for the racemization of tetrasubstituted hydrazides is reported and subsequently analyzed, clarifying the effects that electronic and steric factors may have in the process. In particular, the value obtained by varying the nature of alkyl substituents surrounding the N-N single bond clearly suggests which are those required to obtain atropisomeric hydrazides and which are not suitable for this purpose. The experimental values obtained using both classical racemization experiments and dynamic approaches were comparable with those computed using DFT methods, thus confirming their usefulness and their crucial role in the design of novel atropisomeric scaffolds.

2. Results and Discussion

At the beginning of our investigation, a series of possible atropisomeric hydrazides was submitted to the computational investigation for the determination of the torsional energy barrier (Figure 1).
To give a description of the geometry, dihedral angles (ϑ), and atomic distances (r) are defined in Figure 2.

2.1. Conformational Search

An initial conformation analysis was performed on all seven hydrazides (Figure 1), showing high structural flexibility as evidenced by (i) the high number of conformers obtained in the first conformational search performed with CREST and (ii) the wide dispersion of the Boltzmann populations of the most stable conformers (Table 1).

2.2. Ground Potential Analysis

For each molecule, the lower-energy conformer was then optimized at the DFT level. All the optimized structures show nearly perfectly planar sp2-hybridized nitrogen atoms, as can be seen from dihedral angles ϑ1 and ϑ2 in Table 2. In addition, the best conformer for all molecules except E5 and G7 has almost perfect C2 symmetry. Three of these structures (B2, D4, and G7) have a ϑ3 dihedral angle around the N-N bond of about 60°, while for the other four structures, this angle is closer to the orthogonal arrangement, varying from 77 to about 88°.
These favored conformations appear to be preferred due to two main non-covalent intramolecular interactions (Figure 3a): (a) an interaction between one or two methyl protons of one t-butyl group with the π-system of an aromatic substituent on the opposite nitrogen atom, and (b) the interaction between protons of opposite BOC groups. In addition, it is observed that the relative orientation of the carbonyl groups of the BOCs generates an s-trans configuration (regarding the NNCO fragment) apart from E5, which presents both carbonyls in an s-cis configuration.
Furthermore, an analysis of the total electric dipole moment (Table 3) has been deployed to better understand why conformers prefer the s-trans (outside) configuration against the s-cis (inside) configuration (Table 3). The opposite conformation had been manually generated by the manual rotation of 180° of the ϑ5 and ϑ6 dihedral angles, followed by an optimization at the same level. It is then possible to observe that they all tend to minimize the dipole moment except for E5, probably thanks to the gain in NCIs between methyl groups and carbonyl oxygen atoms.
To obtain an initial guess of the geometry of the torsional TS, like common practices [15,16], a relaxed scan of the ϑ3 dihedral angle around the N-N bond was performed (Figure 3b) both in clockwise and counterclockwise modes. The nomenclature adopted refers to GP1 as the previous optimized structure and GP2 as a conformer of its enantiomer obtained from the scan; cis refers to TS geometries with the BOC groups on the same side of the σ N-N bond, as opposed to trans that have the BOC groups on opposite sides.
Figure 3. (a) B2 best conformer geometry at ωB97x-D/6-31g(d) level and main NCI-index surfaces [17]. (b) Scans of the dihedral θ3 angle of the B2 molecule.
Figure 3. (a) B2 best conformer geometry at ωB97x-D/6-31g(d) level and main NCI-index surfaces [17]. (b) Scans of the dihedral θ3 angle of the B2 molecule.
Molecules 28 07856 g003

2.3. Transition State Optimizations

The guess geometries obtained from the scan did not result in geometries particularly close to the actual TSs, greatly increasing the difficulty of the search. Therefore, it became necessary to tweak the maximum displacement parameter for each cycle of optimization in Gaussian. The efficacy of this strategy may be due to a particularly flat PES of these structures, which, in the search for a transition state, can easily lead to loss of the correct frequency at an excessive motion. This workflow resulted in finding most of the correct TS structures except those for E5 and F6. For these elusive TSs, new predicted geometries were obtained using the Growing String Method [14] (GSM) coupled with xTB. This approach allowed us to obtain TS geometries at the DFT level for all remaining structures. Geometrical parameters for the so-found TSs are reported in Table 4.
A common feature was noted in the spacing of the carbonyl oxygen atoms of 4–5 Å, except for some cis TS geometries. Similarly, there are N-N bond distances around 1.4 Å, slightly longer than the bonds in the minimal structures. Both nitrogen atoms also show strong opposite pyramidalization, as displayed by dihedral angles ϑ1 and ϑ2 (Table 4 and orange circle in Figure 4). Compared with equilibrium structures, TS geometries exhibit a mixed s-cis and s-trans configuration for BOC fragments, as shown by ϑ5 and ϑ6.
From NCI surface analysis (Figure 4 and Supplementary Materials), some important interactions can be detected: (i) π-π stacking between aromatic systems; (ii) CH-π interaction of the methyl groups pointing towards an aromatic ring; (iii) a carbonyl interaction of the oxygen atoms with the hydrogen of the opposite CH2; (iv) a BOC interaction between the carbonyl oxygen of the BOCs that can also give hydrogen bonds on its own t-butyl.
The HOMO orbital of all the TSs is mainly localized in the π-system (when present), as was in the GPs analyzed. In contrast, it is possible to observe a decrease in carbamate fragment conjugation (see Supplementary Materials for more information).

2.4. Energy Barriers from Experimental Data

The atropisomeric interconversion of structurally symmetric samples (A1–F6) was investigated by liquid chromatography. (In the case of compound G7, decomposition was observed during the experimental determination of the rotational barrier).
Two approaches were used: (i) the off-column racemization process monitored by enantioselective HPLC and (ii) the on-line enantiomerization process on chiral stationary phase at variable temperature. After initial screening, aiming for the optimization of analytical conditions for each racemate, racemization and/or enantiomerization experiments were achieved. As reported in Table S1, samples F6 and E5 were baseline separated at low temperatures (10 °C and 0 °C, respectively), attesting to the fast-interconverting speed of each atropisomeric pair. In fact, when analyzed at room temperature, they showed a typical plateau between two peaks. Chromatographic UV and CD profiles were reported in Figures S1 and S2.
Energy values of enantiomerization processes, temperature, and speed constants are listed in Table 5. Notably, the ΔG of racemization (ΔGrac) comes from the off-column experiment, while the ΔG of enantiomerization (ΔGenant) comes from the on-column one. The relationship between two values is: k racemization = 2k enantiomerization. To make the comparison easier with computational data, all racemization-free energies were converted into the corresponding enantiomerization values.
For A1 to D4 samples, an off-column approach was employed at 70 degrees, and the values of the higher free activation energy were expressed by sample C3 (ΔGenant 26.01 kcal/mol). For sample A1, the value of ΔGenant < 24.3 kcal/mol is an estimated value due to the fast racemization (complete in less than 30 min) and to the low column performance in terms of enantioselectivity (α = 1.12). A more flexible N-N bond is the one in sample F6. The plateau is clearly visible at 40 °C (see chromatogram in Figure S3), and by Auto DHPLC y2k, the simulated chromatogram provided a ΔGenant of 23.10 kcal/mol. Unfortunately, it was not possible to acquire chromatograms at higher temperatures to preserve the integrity of the polysaccharide-based stationary phase. A lower energy value has been recorded from sample E5. The ΔGenant of 20.61 kcal/mol at 25 °C has been calculated by simulation of the corresponding dynamic chromatogram (chromatographic trace on the top of Figure 5 right).

2.5. Energy Barrier Values at Different Temperatures for Samples 3 and 6: Eyring Plots

Exploring the dependence of free energy versus temperature, two Eyring plots were built from samples B2 and E5 (Table S2 and Figure S4). Data were obtained by off-column racemization experiments for sample B2 and by on-column enantiomerization experiments for sample E5. In both cases, a slight increment of ΔG was recorded as a temperature increase. In addition, good linearity was observed when plotting 1/T vs. ΔG/T (R2 = 0.9936 and R2 = 0.9979 for samples B2 and E5, respectively), and the corresponding equations have shown a not negligible entropic effect in the N-N interconversion process. The ΔH values are 15.45 kcal/mol and 17.01 kcal/mol, while the ΔS values are −27.9 u.e. and −12.0 u.e., respectively, for sample B2 and sample E5. Figure 5 reports as an example the off-column and on-column approaches employed in this study.

2.6. Comparison Analysis

After these analyses, the obtained energy barrier values have been compared with the experimental ones (Table 6). To better compare the two data, the mRRHO approximation [18] were computed at the temperature of the experiment. It has been noticed that trans-TS is the one determining the epimerization barrier due to its lower energy for all compounds analyzed.
To obtain as much information as possible on the nature of the rotation barrier, further analyses were conducted, such as the Distortion/Interaction Model [19,20] (DIM) and Natural Bond Order [21,22] (NBO). The employment of the DIM allowed an energy estimate to be attributed to the contribution of the various fragments into which the molecule was divided (Figure 2) to the energy of the TS. The plot of Figure 6a, each fragment distortion energy (Equation (1)) and interaction energy value (Equation (2)) have been stacked (on the left) to compare them with the TS energy (on the right). In many of the TSs analyzed, the distortion energy of the fragments is destabilizing; only rarely is this value mildly stabilizing. This phenomenon can be rooted back to internal fragment tension release, thanks to the more room the fragment is given in such an open TS as these ones.
The largest contribution to the TS energy is the deformation of the N-N fragment, while the interaction energies can be stabilizing or destabilizing (Figure 6a and others in Supplementary Materials).
To have a better understanding of the actual nature of the barrier [23], a correlation between the experimental ∆G and the calculated distortion energy that seems to characterize the ∆G (Figure 6b) has been reported. A fair correlation came out: the not-so-high R2 correlation parameter is, for the most part, due to points F6 and E5, for which their ∆Gexp were obtained at different temperatures than the others.
The correlation found indicates that the barrier observed during the racemization process is not only due to the torsion about the σ N-N bond and the clash of the substituents generating the entire rotational barrier, but also to an increase in the p-orbital character of the nitrogen atoms. This hypothesis is also supported by the loss of conjugation of the nitrogen atoms with the carbonyl system of the BOC and their pyramidalization, as it is notable from the increase of the r5, r6, r7, and r8 distances and from the surfaces of the HOMO orbitals.
To further confirm this assumption, a NBO population analysis was carried out. This study allowed us to define the percentage of p orbitals present in the bonds in which nitrogen atoms are involved. In Table 7, the differences in the percentage of the p orbitals of the nitrogen atom between TSs and GPs are correlated to the difference in the bond distance. As previously supposed, slightly more p characters are present in TSs.
In addition to these parameters, bond angles were measured for each nitrogen in the molecules. This, along with ϑ1 and ϑ2, facilitated the assessment of nitrogen pyramidalization and deviation from perfect sp2 hybridization, indicating proximity to sp3 characteristics.
Theoretical expectations dictate that the sum of bond angles around a perfectly sp2 nitrogen should be 360°, while for a perfectly sp3 nitrogen (taking ammonia as a reference), it is 321° (three times the H–N–H bond angle of 107°).
To quantify this pyramidalization of the nitrogen atoms in the TSs geometries, for each one of the two nitrogen atoms, the difference between the sum of bond angles in the GP and the sum of bond angles in the TS geometry was calculated. These values are reported in Table 8 in absolute terms as ∆A and as the %ratio (∆A%) between ∆A and the difference between perfectly sp2 and sp3 observed for ammonia (39°, as shown in Table 8-top) taken as a reference.
In summary, the consideration of the nitrogen atom pyramidalization involves two key parameters: the sum of the bond angles of the substituents on each nitrogen atom (∆A%) and the overall increase in the p character in the bonds around the nitrogen atoms. It is noteworthy that all geometries exhibit a high ∆A%, and each bond demonstrates a higher p percentile in the TS rather than in the GP. This observation is supported by the good correlation between these two factors: the more the sum of the p orbital into the bond, the more the sum of the bond angles approaches the ideal sp3 geometry (Figure 7).

3. Methods and Materials

3.1. Computational Methods

For each substrate, an initial conformational search was performed using the CREST [24] software utility (version 2.11.1), which is based on the xTB [25] engine (version 6.4.1), and it executes an interactive meta-dynamics with a cross-genetic grow algorithm (called iMTD-GC). This analysis was conducted at 298.15 K at the GFN2-xTB [26] level, discarding all conformers with an energy higher than the predefined cutoff value of 6 kcal/mol relative to the lowest-energy conformer.
The seven conformer ensembles obtained were then refined at the DFT-theory level using the CENSO [27] (version 1.2.0) framework interfaced with the ORCA quantum chemistry program package (version 5.0.1) [28] according to the following protocol:
  • part0 cheap prescreening: b97-d3 [29]/def2-SV(P)//GFN2-xTB (Input geometry)
  • part1 prescreening: r2scan-3c [30] + C-PCM[hexane] + GmRRHO [18] (GFN2[ALPB [31]]-bhess [32])//GFN2-xTB (Input geometry)
  • part2 optimization: r2scan-3c + C-PCM[hexane] + GmRRHO(GFN2[ALPB]-bhess)//r2scan-3c[SMD]
This framework enabled faster processing of the vast conformational ensemble of a large number of conformers since each part analyzed the conformers and sorted out those with a relative energy greater than the default (for part0 4.0 kcal/mol, part1 3.5 kcal/mol, and part2 2.5 kcal/mol). For the last two parts, the free energy was calculated via the mRRHO approximation [18].
All following calculations have been performed using the Gaussian16 (rev. A.03) [33] program at the DFT level with the ωB97x-D/6-31g(d) model chemistry and C-PCM implicit solvation in n-hexane. These parameters were chosen on the basis of previous works on similar compounds [34]. Free energy has been calculated with a vertical excitation energy correction, modifying the functional and the basis set to M06-2X/def2-TZVP, and using the previously calculated thermochemical correction [35].
To accurately describe non-covalent interactions, the NCI index was calculated with the NCIPLOT [36] software (version 4.0) using the FINE grid of integration.
To study the energetic contributions to the transition states (TS), the distortion-interaction model (DIM) was employed [19,20]. This analysis was carried out by fragmenting each molecule into five groups (Figure 2): each of four substituents and the N-N core. For each fragment, the distortion energy ( E d i ) was calculated as the difference between the energy of the segment at its TS geometry ( E T S i * ) and at its ground potential (GP1) geometry ( E G P i ). Thus, the interaction energy ( E i ) was obtained by the difference between the single point energy of the whole TS and the sum of all distortion energies.
E d i = E T S i * E G P i
E i = E T S fragment E d i

3.2. Synthesis and Experimental Part

Hydrazides have been synthesized with an achiral methodology, as reported in Figure 8. The synthetic strategy can be divided into three categories: the first two employed the di-tert-butyl hydrazine-1,2-dicarboxylate and the corresponding bromide derivates or anhydrides to give A1-F6 compounds. In alternative, di-tert-butyl azodicarboxylate has been employed in a sequential pathway to obtain compound G7. Only the first one has already been discussed in the literature [37]; the others were not optimized procedures.

4. Conclusions

During this study, a systematic computational and experimental analysis of seven atropoisomeric hydrazides was conducted. For the herein studied molecules, both the GPs and the transition state geometries have been thoroughly investigated, in good accordance with the experimental racemization studies. Furthermore, with both the DIM and the NBO analyses, it has been possible to assign as the cause of the rotational impeding the predominant distortion energy of the N–N bonds and the partial rehybridization of the two nitrogen atoms.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28237856/s1; File S1 Supporting information, refs [38,39,40,41,42,43,44,45] are cited in Supplementary Materials.

Author Contributions

Conceptualization, A.P. and G.B.; Methodology, A.P., L.M., P.R., S.M. and A.C.; Software, A.P. and P.R.; Validation, P.R.; Formal analysis, G.C., C.P., M.C., S.B.O., S.M. and A.C.; Investigation, A.P., L.M., G.C., C.P., M.C., S.B.O. and A.C.; Data curation, S.M. and A.C.; Writing—original draft, A.P. and G.B.; Writing—review & editing, P.R., A.C. and G.B.; Supervision, G.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union, Next Generation EU, DM 737/2021, University of Bologna Project ALMArieCurie 2022, SUpER, and the “Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali” (INSTM).

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Centonze, G.; Portolani, C.; Righi, P.; Bencivenni, G. Enantioselective Strategies for The Synthesis of N-N Atropisomers. Angew. Chem. 2023, 135, e202303966. [Google Scholar] [CrossRef]
  2. Dewar, M.J.S.; Jennings, B. Rotation-Inversion Barriers in Hydrazines. J. Am. Chem. Soc. 1969, 91, 3655–3656. [Google Scholar] [CrossRef]
  3. Verma, S.M.; Prasad, R. Conformational Analysis by Nuclear Magnetic Resonance Spectroscopy. N′-Derivatives of N-Aminocamphorimides. J. Org. Chem. 1973, 38, 1004–1010. [Google Scholar] [CrossRef]
  4. Atkinson, R.S.; Barker, E.; Price, C.J.; Russell, D.R. The N–N Bond as a Chiral Axis: 3-(Diacylamino)Quinazolin-4(3H)Ones as Chiral Acylating Agents. J. Chem. Soc. Chem. Commun. 1994, 1159–1160. [Google Scholar] [CrossRef]
  5. Atkinson, R.S.; Barker, E.; Edwards, P.J.; Thomson, G.A. The N–N Bond as a Chiral Axis: 3-Diacylaminoquinazolinones as Chiral Acylating Agents. J. Chem. Soc. Perkin Trans. 1996, 1, 1047–1055. [Google Scholar] [CrossRef]
  6. Coogan, M.P.; Smart, E.; Hibbs, D.E. Asymmetric Transformation (Deracemisation) of an Atropisomeric Bisheterocyclic Amine. Chem. Commun. 1999, 1991–1992. [Google Scholar] [CrossRef]
  7. Coogan, M.P.; Passey, S.C. Tetraacyl Hydrazines and 3,3′-Biquinazoline-4,4′-Diones; Synthesis, Studies of Rotational Barriers and Deracemisation. J. Chem. Soc. Perkin Trans. 2000, 2, 2060–2066. [Google Scholar] [CrossRef]
  8. Wang, X.-M.; Zhang, P.; Xu, Q.; Guo, C.-Q.; Zhang, D.-B.; Lu, C.-J.; Liu, R.-R. Enantioselective Synthesis of Nitrogen–Nitrogen Biaryl Atropisomers via Copper-Catalyzed Friedel–Crafts Alkylation Reaction. J. Am. Chem. Soc. 2021, 143, 15005–15010. [Google Scholar] [CrossRef]
  9. Chen, K.-W.; Chen, Z.-H.; Yang, S.; Wu, S.-F.; Zhang, Y.-C.; Shi, F. Organocatalytic Atroposelective Synthesis of N−N Axially Chiral Indoles and Pyrroles by De Novo Ring Formation. Angew. Chem. Int. Ed. 2022, 61, e202116829. [Google Scholar] [CrossRef]
  10. Gao, Y.; Wang, L.-Y.; Zhang, T.; Yang, B.-M.; Zhao, Y. Atroposelective Synthesis of 1,1′-Bipyrroles Bearing a Chiral N−N Axis: Chiral Phosphoric Acid Catalysis with Lewis Acid Induced Enantiodivergence. Angew. Chem. Int. Ed. 2022, 61, e202200371. [Google Scholar] [CrossRef]
  11. Lin, W.; Zhao, Q.; Li, Y.; Pan, M.; Yang, C.; Yang, G.-H.; Li, X. Asymmetric Synthesis of N–N Axially Chiral Compounds via Organocatalytic Atroposelective N-Acylation. Chem. Sci. 2021, 13, 141–148. [Google Scholar] [CrossRef] [PubMed]
  12. Mei, G.-J.; Wong, J.J.; Zheng, W.; Nangia, A.A.; Houk, K.N.; Lu, Y. Rational Design and Atroposelective Synthesis of N–N Axially Chiral Compounds. Chem 2021, 7, 2743–2757. [Google Scholar] [CrossRef]
  13. Chen, Z.; Li, T.; Wang, N.; Ma, X.; Ni, S.; Zhang, Y.; Shi, F. Organocatalytic Enantioselective Synthesis of Axially Chiral N,N′-Bisindoles. Angew. Chem. Int. Ed. 2023, 62, e202300419. [Google Scholar] [CrossRef]
  14. Wang, L.-Y.; Miao, J.; Zhao, Y.; Yang, B.-M. Chiral Acid-Catalyzed Atroposelective Indolization Enables Access to 1,1′-Indole-Pyrroles and Bisindoles Bearing a Chiral N–N Axis. Org. Lett. 2023, 25, 1553–1557. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, Q.; Mándi, A.; Li, S.; Chen, Y.; Zhang, W.; Tian, X.; Zhang, H.; Li, H.; Zhang, W.; Zhang, S.; et al. N-N-Coupled Indolo-Sesquiterpene Atropo-Diastereomers from a Marine-Derived Actinomycete. Eur. J. Org. Chem. 2012, 2012, 5256–5262. [Google Scholar] [CrossRef]
  16. LaPlante, S.R.; Edwards, P.J.; Fader, L.D.; Jakalian, A.; Hucke, O. Revealing Atropisomer Axial Chirality in Drug Discovery. ChemMedChem 2011, 6, 505–513. [Google Scholar] [CrossRef]
  17. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  18. Grimme, S. Supramolecular Binding Thermodynamics by Dispersion-Corrected Density Functional Theory. Chem.–Eur. J. 2012, 18, 9955–9964. [Google Scholar] [CrossRef]
  19. Houk, K.N.; Bickelhaupt, M. Analyzing Reaction Rates with the Distortion Interaction-Activation Strain Model. Angew. Chem. Int. Ed. 2017, 56, 10070–10086. [Google Scholar] [CrossRef]
  20. Velasco-Juárez, E.; Arpa, E.M. A Novel Partitioning Scheme for the Application of the Distortion/Interaction–Activation Strain Model to Intramolecular Reactions. Theor. Chem. Acc. 2021, 140, 107. [Google Scholar] [CrossRef]
  21. Gloaguen, E.; Brenner, V.; Alauddin, M.; Tardivel, B.; Mons, M.; Zehnacker-Rentien, A.; Declerck, V.; Aitken, D.J. Direct Spectroscopic Evidence of Hyperconjugation Unveils the Conformational Landscape of Hydrazides. Angew. Chem. Int. Ed. 2014, 53, 13756–13759. [Google Scholar] [CrossRef] [PubMed]
  22. Weinhold, F.; Landis, C.R.; Glendening, E.D. What Is NBO Analysis and How Is It Useful? Int. Rev. Phys. Chem. 2016, 35, 399–440. [Google Scholar] [CrossRef]
  23. He, C.Q.; Simon, A.; Lam, Y.; Brunskill, A.P.J.; Yasuda, N.; Tan, J.; Hyde, A.M.; Sherer, E.C.; Houk, K.N. Model for the Enantioselectivity of Asymmetric Intramolecular Alkylations by Bis-Quaternized Cinchona Alkaloid-Derived Catalysts. J. Org. Chem. 2017, 82, 8645–8650. [Google Scholar] [CrossRef] [PubMed]
  24. Pracht, P.; Bohle, F.; Grimme, S. Automated Exploration of the Low-Energy Chemical Space with Fast Quantum Chemical Methods. Phys. Chem. Chem. Phys. 2020, 22, 7169–7192. [Google Scholar] [CrossRef] [PubMed]
  25. Bannwarth, C.; Caldeweyher, E.; Ehlert, S.; Hansen, A.; Pracht, P.; Seibert, J.; Spicher, S.; Grimme, S. Extended TIGHT‐BINDING Quantum Chemistry Methods. WIREs Comput. Mol. Sci. 2021, 11, e1493. [Google Scholar] [CrossRef]
  26. Bannwarth, C.; Ehlert, S.; Grimme, S. GFN2-xTB—An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions. J. Chem. Theory Comput. 2019, 15, 1652–1671. [Google Scholar] [CrossRef] [PubMed]
  27. Grimme, S.; Bohle, F.; Hansen, A.; Pracht, P.; Spicher, S.; Stahn, M. Efficient Quantum Chemical Calculation of Structure Ensembles and Free Energies for Nonrigid Molecules. J. Phys. Chem. A 2021, 125, 4039–4054. [Google Scholar] [CrossRef]
  28. Neese, F. Software Update: The ORCA Program System—Version 5.0. WIREs Comput. Mol. Sci. 2022, 12, e1606. [Google Scholar] [CrossRef]
  29. Brandenburg, J.G.; Bannwarth, C.; Hansen, A.; Grimme, S. B97-3c: A Revised Low-Cost Variant of the B97-D Density Functional Method. J. Chem. Phys. 2018, 148, 064104. [Google Scholar] [CrossRef]
  30. Grimme, S.; Hansen, A.; Ehlert, S.; Mewes, J.-M. r2SCAN-3c: A “Swiss Army Knife” Composite Electronic-Structure Method. J. Chem. Phys. 2021, 154, 064103. [Google Scholar] [CrossRef]
  31. Ehlert, S.; Stahn, M.; Spicher, S.; Grimme, S. Robust and Efficient Implicit Solvation Model for Fast Semiempirical Methods. J. Chem. Theory Comput. 2021, 17, 4250–4261. [Google Scholar] [CrossRef] [PubMed]
  32. Spicher, S.; Grimme, S. Single-Point Hessian Calculations for Improved Vibrational Frequencies and Rigid-Rotor-Harmonic-Oscillator Thermodynamics. J. Chem. Theory Comput. 2021, 17, 1701–1714. [Google Scholar] [CrossRef] [PubMed]
  33. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision A.03; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  34. Portolani, C.; Centonze, G.; Luciani, S.; Pellegrini, A.; Righi, P.; Mazzanti, A.; Ciogli, A.; Sorato, A.; Bencivenni, G. Synthesis of Atropisomeric Hydrazides by One-Pot Sequential Enantio- and Diastereoselective Catalysis. Angew. Chem. Int. Ed. 2022, 61, e202209895. [Google Scholar] [CrossRef] [PubMed]
  35. Simón, L.; Goodman, J.M. How Reliable Are DFT Transition Structures? Comparison of GGA, Hybrid-Meta-GGA and Meta-GGA Functionals. Org. Biomol. Chem. 2011, 9, 689–700. [Google Scholar] [CrossRef]
  36. Contreras-García, J.; Johnson, E.R.; Keinan, S.; Chaudret, R.; Piquemal, J.-P.; Beratan, D.N.; Yang, W. NCIPLOT: A Program for Plotting Noncovalent Interaction Regions. J. Chem. Theory Comput. 2011, 7, 625–632. [Google Scholar] [CrossRef]
  37. Rasmussen, L.K. Facile Synthesis of Mono-, Di-, and Trisubstituted Alpha-Unbranched Hydrazines. J. Org. Chem. 2006, 71, 3627–3629. [Google Scholar] [CrossRef] [PubMed]
  38. Bredihhin, A.; Mäeorg, U. Effective Strategy for the Systematic Synthesis of Hydrazine Derivatives. Tetrahedron 2008, 64, 6788–6793. [Google Scholar] [CrossRef]
  39. Gasparrini, F.; Lunazzi, L.; Mazzanti, A.; Pierini, M.; Pietrusiewicz, K.M.; Villani, C. Comparison of Dynamic HPLC and Dynamic NMR in the Study of Conformational Stereodynamics: Case of the Enantiomers of a Hindered Secondary Phosphine Oxide1. J. Am. Chem. Society 2000, 122, 4776–4780. [Google Scholar] [CrossRef]
  40. Gasparrini, F.; D’Acquarica, I.; Pierini, M.; Villani, C. Chromatographic Resolution and Enantiomerization Barriers of Axially Chiral 1-Naphthamides. J. Sep. Sci. 2001, 24, 941–946. [Google Scholar] [CrossRef]
  41. Sabia, R.; Ciogli, A.; Pierini, M.; Franzini, R.; Iazzetti, A.; Villani, C. Chromatographic Separation of the Interconverting Enantiomers of Imidazo- and Triazole-Fused Benzodiazepines. J. Chromatogr. A 2021, 1647, 462148. [Google Scholar] [CrossRef]
  42. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera? A Visualization System for Exploratory Research and Analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed]
  43. Zimmerman, P.M. Growing String Method with Interpolation and Optimization in Internal Coordinates: Method and Examples. J. Chem. Phys. 2013, 138, 184102. [Google Scholar] [CrossRef] [PubMed]
  44. Luchini, G.; Alegre-Requena, J.V.; Funes-Ardoiz, I.; Paton, R.S. GoodVibes: Automated Thermochemistry for Heterogeneous Computational Chemistry Data. F1000Research 2020, 9, 291. [Google Scholar] [CrossRef]
  45. Zhao, Q.; Hsu, H.-H.; Savoie, B. Conformational Sampling for Transition State Searches on a Computational Budget. J. Chem. Theory Comput. 2021, 18, 3006–3016. [Google Scholar] [CrossRef]
Figure 1. Substrates considered.
Figure 1. Substrates considered.
Molecules 28 07856 g001
Figure 2. Geometrical parameter definition and fragmentation are used to apply the distortion/interaction model.
Figure 2. Geometrical parameter definition and fragmentation are used to apply the distortion/interaction model.
Molecules 28 07856 g002
Figure 4. (a) B2-Cis TS at ωB97x-D/6-31g(d) level obtained from the counterclockwise relaxed scan of the dihedral angle around the N-N bond. (b) B2-Trans TS at ωB97x-D/6-31g(d) level obtained from the clockwise relaxed scan of the same dihedral [17].
Figure 4. (a) B2-Cis TS at ωB97x-D/6-31g(d) level obtained from the counterclockwise relaxed scan of the dihedral angle around the N-N bond. (b) B2-Trans TS at ωB97x-D/6-31g(d) level obtained from the clockwise relaxed scan of the same dihedral [17].
Molecules 28 07856 g004
Figure 5. Thermal racemization of B2 (left) and enantiomerization of E5 (right).
Figure 5. Thermal racemization of B2 (left) and enantiomerization of E5 (right).
Molecules 28 07856 g005
Figure 6. (a) DIM analysis of the trans TS for B2 molecule at ωB97x-D/6-31g(d) level. (b) Correlation between experimental ∆G and NN distortion energy.
Figure 6. (a) DIM analysis of the trans TS for B2 molecule at ωB97x-D/6-31g(d) level. (b) Correlation between experimental ∆G and NN distortion energy.
Molecules 28 07856 g006
Figure 7. Correlation between the %ratio of the bond angles and the overall sum of the p character of the nitrogen atoms.
Figure 7. Correlation between the %ratio of the bond angles and the overall sum of the p character of the nitrogen atoms.
Molecules 28 07856 g007
Figure 8. Synthetic procedures for investigated molecules.
Figure 8. Synthetic procedures for investigated molecules.
Molecules 28 07856 g008
Table 1. Number of conformers generated from the CREST-CENSO protocol. Values in parentheses represent the Boltzmann population of the best conformer at that stage.
Table 1. Number of conformers generated from the CREST-CENSO protocol. Values in parentheses represent the Boltzmann population of the best conformer at that stage.
ParameterCRESTPart0Part1Part2
A193 (67.0%)897848 (31.5%)
B2241 (22.2%)17712796 (4.5%)
C3125 (46.9%)723731 (29.4%)
D471 (24.7%)473730 (17.0%)
E586 (30.1%)605951 (11.2%)
F666 (9.8%)404035 (7.5%)
G7115 (38.1%)836250 (13.1%)
Table 2. Geometrical parameters of the best conformer for each compound at ωB97x-D/6-31g(d) level.
Table 2. Geometrical parameters of the best conformer for each compound at ωB97x-D/6-31g(d) level.
ParameterA1B2C3D4E5F6G7
ϑ1 [°]−0.759.26−1.26−10.414.351.48−8.31
ϑ2 [°]−0.758.84−1.25−10.41−5.821.48−11.52
ϑ3 [°]77.89−61.5577.5862.7387.70−84.1262.11
ϑ5 [°]179.87171.67179.23−171.3819.01179.24−170.04
ϑ6 [°]179.87171.35179.23−171.38−20.29179.24−172.92
r4 [Å]1.371.381.371.381.371.381.38
r5 [Å]1.381.381.381.381.411.411.38
r6 [Å]1.381.381.381.381.411.411.38
Table 3. Dipole momentum of the optimized structures at ωB97x-D/6-31g(d) level.
Table 3. Dipole momentum of the optimized structures at ωB97x-D/6-31g(d) level.
Molecules 28 07856 i001
InsideOutsideΔG kcal/mol
(In-Out)
A11.9848 [b]0.9121 [a]+5.85
B21.5692 [b]0.8129 [a]+8.21
C31.9374 [b]0.8278 [a]+8.47
D41.5456 [b]0.8058 [a]+3.32
E52.1177 [a]5.2156 [b]−1.75
F62.7641 [b]5.7090 [a]+3.92
G71.7223 [b]0.6977 [a]+5.43
[a] structures generated from the CREST-CENSO protocol. [b] manually generated by manually rotating by 180° the ϑ4 and ϑ5 dihedral angles of the structures.
Table 4. Geometrical parameters of the transition state determine the torsional barrier at the ωB97x-D/6-31g(d) level.
Table 4. Geometrical parameters of the transition state determine the torsional barrier at the ωB97x-D/6-31g(d) level.
ParameterA1B2C3D4E5F6G7
ϑ1 [°]20.71−21.3335.4923.46−25.12−31.56−31.60
ϑ2 [°]−29.5328.85−20.30−34.3728.2524.5319.85
ϑ3 [°]−173.62−161.50170.61−177.37169.33−166.89−176.93
ϑ5 [°]162.93−166.00−179.6011.77−129.29−179.69−2.07
ϑ6 [°]−1.7113.9739.47−172.14−40.38−8.82169.36
r4 [Å]1.441.441.441.441.421.431.44
r5 [Å]1.391.381.441.401.461.411.41
r6 [Å]1.411.421.391.401.471.401.39
Table 5. Energy values of enantiomerization processes, temperature, speed constants, and half-life for each sample were listed.
Table 5. Energy values of enantiomerization processes, temperature, speed constants, and half-life for each sample were listed.
SampleT [°C]krac [min−1]kenant [min−1]ΔGenant [kcal/mol]t1/2 [min]
D4700.0111 (R2 = 0.9974)-25.5731
C3710.0064 (R2 = 0.9989)-26.0154
B2710.0250 (R2 = 0.9991)-25.0914
A1700.0645-<24.3 a5.4
F640-0.029623.10 b12
E525-0.3560 b20.61 b2
a estimated value: the racemization is completed in 30 min. b averaged value between the free energy of the direct (ΔG1-2) and reversed (ΔG2-1) enantiomerization processes.
Table 6. Single-point refined rotational ΔG at M06-2x/def2-TZVP//ωB97x-D/6-31g(d) level.
Table 6. Single-point refined rotational ΔG at M06-2x/def2-TZVP//ωB97x-D/6-31g(d) level.
Molecule∆G Exp
[kcal/mol]
T exp [°C]∆G Computational [kcal/mol]Deviation from
Experimental Value [%]
A1<24.37029.2320.29 [a]
B225.097126.314.87
C326.017136.2339.29
D425.577024.81−2.96
E520.612522.579.52
F623.104028.4623.20
G7N.A.26.87 [b]N.A.
N.A.: not available. [a]: considering experimental ΔG = 24.3 kcal/mol; [b]: calculated at 298.15 K.
Table 7. NBO analysis at ωB97x-D/6-31g(d) level.
Table 7. NBO analysis at ωB97x-D/6-31g(d) level.
Bond [a]A1B2C3D4E5F6G7
∆bd [Å]p [%]∆bd [Å]p [%]∆bd [Å]p [%]∆bd [Å]p [%]∆bd [Å]p [%]∆bd [Å]p [%]∆bd [Å]p [%]
1–20.072.560.062.100.074.880.063.550.053.020.052.960.063.98
2–10.073.900.064.340.072.530.063.390.052.840.052.410.062.26
1–30.011.560.012.020.035.810.033.16−0.012.010.077.230.045.18
1–50.012.950.001.100.067.530.023.090.056.600.004.900.044.44
2–40.045.200.033.640.011.330.022.85−0.013.400.067.030.020.85
2–60.044.790.045.650.012.640.023.240.065.730.003.090.001.75
∆bd: difference in bond distance between TS and GP. ∆p [%]: difference in %p orbital involved in the bond (%pTS-%pGS). [a]: The numbers defining the bond refer to the indexing encoded in Figure 2.
Table 8. Difference of angles (GP-TS). The index of the nitrogen is according to Figure 2 numeration.
Table 8. Difference of angles (GP-TS). The index of the nitrogen is according to Figure 2 numeration.
Molecules 28 07856 i002
N1N2
Molecule∆A∆A [%]∆A∆A [%]
A110.2626.3110.2626.31
B27.3318.7917.3844.56
C329.3475.239.0023.08
D413.8035.3813.7035.13
E516.9343.4117.6645.28
F622.4457.5416.1941.51
G717.8545.776.6817.13
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

Pellegrini, A.; Marcon, L.; Righi, P.; Centonze, G.; Portolani, C.; Capodiferro, M.; Oljira, S.B.; Manetto, S.; Ciogli, A.; Bencivenni, G. On the Nature of the Rotational Energy Barrier of Atropisomeric Hydrazides. Molecules 2023, 28, 7856. https://doi.org/10.3390/molecules28237856

AMA Style

Pellegrini A, Marcon L, Righi P, Centonze G, Portolani C, Capodiferro M, Oljira SB, Manetto S, Ciogli A, Bencivenni G. On the Nature of the Rotational Energy Barrier of Atropisomeric Hydrazides. Molecules. 2023; 28(23):7856. https://doi.org/10.3390/molecules28237856

Chicago/Turabian Style

Pellegrini, Andrea, Laura Marcon, Paolo Righi, Giovanni Centonze, Chiara Portolani, Marco Capodiferro, Shilashi Badasa Oljira, Simone Manetto, Alessia Ciogli, and Giorgio Bencivenni. 2023. "On the Nature of the Rotational Energy Barrier of Atropisomeric Hydrazides" Molecules 28, no. 23: 7856. https://doi.org/10.3390/molecules28237856

APA Style

Pellegrini, A., Marcon, L., Righi, P., Centonze, G., Portolani, C., Capodiferro, M., Oljira, S. B., Manetto, S., Ciogli, A., & Bencivenni, G. (2023). On the Nature of the Rotational Energy Barrier of Atropisomeric Hydrazides. Molecules, 28(23), 7856. https://doi.org/10.3390/molecules28237856

Article Metrics

Back to TopTop