Structure of Molecular Tweezer Complexes in the Solid State: NMR Experiments, X-ray Investigations, and Quantum Chemical Calculations

Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning

Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.

Synthesis of Doubly Annulated m-Terphenyl-Based Molecular Tweezers and Their Charge-Transfer Complexes with DDQ as a Guest

The synthesis of a doubly-annulated m-terphenyl-based tweezer platform has been developed, which affords ready incorporation of various pincer units from monobenzenoid to polybenzenoid electron donors. The complexation study with 2,3-dichloro-5,6-dicyano-1,4-benzoquinone (DDQ) as guest has been carried out, and the crystal structure of T-Py∩DDQ reveals the sandwich-type binding mode in the solid state.

Sensitivity of ab Initio vs Empirical Methods in Computing Structural Effects on NMR Chemical Shifts for the Example of Peptides

The structural sensitivity of NMR chemical shifts as computed by quantum chemical methods is compared to a variety of empirical approaches for the example of a prototypical peptide, the 38-residue kaliotoxin KTX comprising 573 atoms. Despite the simplicity of empirical chemical shift prediction programs, the agreement with experimental results is rather good, underlining their usefulness. However, we show in our present work that they are highly insensitive to structural changes, which renders their use for validating predicted structures questionable. In contrast, quantum chemical methods show the expected high sensitivity to structural and electronic changes. This appears to be independent of the quantum chemical approach or the inclusion of solvent effects. For the latter, explicit solvent simulations with increasing number of snapshots were performed for two conformers of an eight amino acid sequence. In conclusion, the empirical approaches neither provide the expected magnitude nor the patterns of NMR chemical shifts determined by the clearly more costly ab initio methods upon structural changes. This restricts the use of empirical prediction programs in studies where peptide and protein structures are utilized for the NMR chemical shift evaluation such as in NMR refinement processes, structural model verifications, or calculations of NMR nuclear spin relaxation rates.

The use of a selective saturation pulse to suppress t1 noise in two-dimensional (1)H fast magic angle spinning solid-state NMR spectroscopy

A selective saturation pulse at fast magic angle spinning (MAS) frequencies (60+kHz) suppresses t1 noise in the indirect dimension of two-dimensional (1)H MAS NMR spectra. The method is applied to a synthetic nucleoside with an intense methyl (1)H signal due to triisopropylsilyl (TIPS) protecting groups. Enhanced performance in terms of suppressing the methyl signal while minimising the loss of signal intensity of nearby resonances of interest relies on reducing spin diffusion--this is quantified by comparing two-dimensional (1)H NOESY-like spin diffusion spectra recorded at 30-70 kHz MAS. For a saturation pulse centred at the methyl resonance, the effect of changing the nutation frequency at different MAS frequencies as well as the effect of changing the pulse duration is investigated. By applying a pulse of duration 30 ms and nutation frequency 725 Hz at 70 kHz MAS, a good compromise of significant suppression of the methyl resonance combined with the signal intensity of resonances greater than 5 ppm away from the methyl resonance being largely unaffected is achieved. The effectiveness of using a selective saturation pulse is demonstrated for both homonuclear (1)H-(1)H double quantum (DQ)/single quantum (SQ) MAS and (14)N-(1)H heteronuclear multiple quantum coherence (HMQC) two-dimensional solid-state NMR experiments.

Dynamic nuclear polarization NMR spectroscopy of microcrystalline solids

Dynamic nuclear polarization (DNP) solid-state NMR has been applied to powdered microcrystalline solids to obtain sensitivity enhancements on the order of 100. Glucose, sulfathiazole, and paracetamol were impregnated with bis-nitroxide biradical (bis-cyclohexyl-TEMPO-bisketal, bCTbK) solutions of organic solvents. The organic solvents were carefully chosen to be nonsolvents for the compounds, so that DNP-enhanced solid-state NMR spectra of the unaltered solids could be acquired. A theoretical model is presented that illustrates that for externally doped organic solids characterized by long spin-lattice relaxation times (T(1)((1)H) > 200 s), (1)H-(1)H spin diffusion can relay enhanced polarization over micrometer length scales yielding substantial DNP enhancements (ε). ε on the order of 60 are obtained for microcrystalline glucose and sulfathiazole at 9.4 T and with temperatures of ca. 105 K. The large gain in sensitivity enables the rapid acquisition of (13)C-(13)C correlation spectra at natural isotopic abundance. It is anticipated that this will be a general method for enhancing the sensitivity of solid-state NMR experiments of organic solids.

Applications of high-resolution 1H solid-state NMR

This article reviews the large increase in applications of high-resolution (1)H magic-angle spinning (MAS) solid-state NMR, in particular two-dimensional heteronuclear and homonuclear (double-quantum and spin-diffusion NOESY-like exchange) experiments, in the last five years. These applications benefit from faster MAS frequencies (up to 80 kHz), higher magnetic fields (up to 1 GHz) and pulse sequence developments (e.g., homonuclear decoupling sequences applicable under moderate and fast MAS). (1)H solid-state NMR techniques are shown to provide unique structural insight for a diverse range of systems including pharmaceuticals, self-assembled supramolecular structures and silica-based inorganic-organic materials, such as microporous and mesoporous materials and heterogeneous organometallic catalysts, for which single-crystal diffraction structures cannot be obtained. The power of NMR crystallography approaches that combine experiment with first-principles calculations of NMR parameters (notably using the GIPAW approach) are demonstrated, e.g., to yield quantitative insight into hydrogen-bonding and aromatic CH-π interactions, as well as to generate trial three-dimensional packing arrangements. It is shown how temperature-dependent changes in the (1)H chemical shift, linewidth and DQ-filtered signal intensity can be analysed to determine the thermodynamics and kinetics of molecular level processes, such as the making and breaking of hydrogen bonds, with particular application to proton-conducting materials. Other applications to polymers and biopolymers, inorganic compounds and bioinorganic systems, paramagnetic compounds and proteins are presented. The potential of new technological advances such as DNP methods and new microcoil designs is described.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

NMR crystallography of campho[2,3-c]pyrazole (Z' = 6): combining high-resolution 1H-13C solid-state MAS NMR spectroscopy and GIPAW chemical-shift calculations

(1)H-(13)C two-dimensional magic-angle spinning (MAS) solid-state NMR correlation spectra, recorded with the MAS-J-HMQC experiment, are presented for campho[2,3-c]pyrazole. For each (13)C moiety, there are six resonances associated with the six distinct molecules in the asymmetric unit cell (Z' = 6). The one-bond C-H correlations observed in the 2D (1)H-(13)C MAS-J-HMQC spectra allow the experimental determination of the (1)H and (13)C chemical shifts associated with the separate CH, CH(2), and CH(3) groups. (1)H and (13)C chemical shifts calculated by using the GIPAW (Gauge Including Projector Augmented Waves) plane-wave pseudopotential approach are presented. Calculations for the whole unit cell (12 × 29 = 348 atoms, with geometry optimization of all atoms) allow the assignment of the experimental (1)H and (13)C chemical shifts to the six distinct molecules. The calculated chemical shifts for the full crystal structure are compared with those for isolated molecules as extracted from the geometry-optimized crystal structure. In this way, the effect of intermolecular interactions on the observed chemical shifts is quantified. In particular, the calculations are sufficiently precise to differentiate the small (<1 ppm) differences between the (1)H chemical shifts of the six resonances associated with each distinct CH or CH(2) moiety.

Complete (1)H resonance assignment of beta-maltose from (1)H-(1)H DQ-SQ CRAMPS and (1)H (DQ-DUMBO)-(13)C SQ refocused INEPT 2D solid-state NMR spectra and first principles GIPAW calculations

A disaccharide is a challenging case for high-resolution (1)H solid-state NMR because of the 24 distinct protons (14 aliphatic and 10 OH) having (1)H chemical shifts that all fall within a narrow range of approximately 3 to 7 ppm. High-resolution (1)H (500 MHz) double-quantum (DQ) combined rotation and multiple pulse sequence (CRAMPS) solid-state NMR spectra of beta-maltose monohydrate are presented. (1)H-(1)H DQ-SQ CRAMPS spectra are presented together with (1)H (DQ)-(13)C correlation spectra obtained with a new pulse sequence that correlates a high-resolution (1)H DQ dimension with a (13)C single quantum (SQ) dimension using the refocused INEPT pulse-sequence element to transfer magnetization via one-bond (13)C-(1)H J couplings. Compared to the observation of only a single broad peak in a (1)H DQ spectrum recorded at 30 kHz magic-angle spinning (MAS), the use of DUMBO (1)H homonuclear decoupling in the (1)H DQ CRAMPS experiment allows the resolution of distinct DQ correlation peaks which, in combination with first-principles chemical shift calculations based on the GIPAW (Gauge Including Projector Augmented Waves) plane-wave pseudopotential approach, enables the assignment of the (1)H resonances to the 24 distinct protons. We believe this to be the first experimental solid-state NMR determination of the hydroxyl OH (1)H chemical shifts for a simple sugar. Variable-temperature (1)H-(1)H DQ CRAMPS spectra reveal small increases in the (1)H chemical shifts of the OH resonances upon decreasing the temperature from 348 K to 248 K.

Preparation, structures and preliminary host-guest studies of fluorinated syn-bis-quinoxaline molecular tweezers

A series of polycyclic frameworks with fluorinated syn-facial quinoxaline sidewalls has been prepared as potential molecular tweezers for electron-rich guest compounds. Our synthetic route to the cyclooctadiene-derived scaffolds 16a-d takes advantage of the facile isolation of a novel spirocyclic precursor 9b with the crucial syn-orientation of its two alkene moieties. The crystal structure of 16c displays two features typical of a molecular tweezer: inclusion of a solvent molecule in the molecular cleft and self-association of the self-complementary scaffolds. Furthermore, host-guest NMR studies of compound 16c in solution show chemical exchange between the unbound and bound electron-rich guest, N,N,N',N'-tetramethyl-p-phenylenediamine.

Determining relative proton-proton proximities from the build-up of two-dimensional correlation peaks in 1H double-quantum MAS NMR: insight from multi-spin density-matrix simulations

The build-up of intensity-as a function of the number, n(rcpl), of POST-C7 elements used for the excitation and reconversion of double-quantum (DQ) coherence (DQC)-is analysed for the fifteen distinct DQ correlation peaks that are observed experimentally for the eight separate (1)H resonances in a (1)H (500 MHz) DQ CRAMPS solid-state (12.5 kHz MAS) NMR spectrum of the dipeptide beta-AspAla (S. P. Brown, A. Lesage, B. Elena, and L. Emsley, J. Am. Chem. Soc., 2004, 126, 13230). The simulation in SPINEVOLUTION (M. Veshtort and R. G. Griffin, J. Magn. Reson., 2006, 178, 248) of t(1) ((1)H DQ evolution) FIDs for clusters of eight dipolar-coupled protons gives separate simulated (1)H DQ build-up curves for the CH(2)(a), CH(2)(b), CH(Asp), CH(Ala), NH and OH (1)H single-quantum (SQ) (1)H resonances. An analysis of both the simulated and experimental (1)H DQ build-up leads to the following general observations: (i) considering the build-up of (1)H DQ peaks at a particular SQ frequency, maximum intensity is observed for the DQC corresponding to the shortest H-H distance; (ii) for the maximum intensity (1)H DQ peak at a particular SQ frequency, the recoupling time for the observed maximum intensity depends on the corresponding H-H distance, e.g., maximum intensity for the CH(2)(a)-CH(2)(b) (H-H distance = 1.55 A) and OH-CH(Asp) (H-H distance = 2.49 A) DQ peaks is observed at n(rcpl) = 2 and 3, respectively; (iii) for DQ peaks involving a CH(2) proton at a non-CH(2) SQ frequency, there is much reduced intensity and a maximum intensity at a short recoupling time; (iv) for the other lower intensity (1)H DQ peaks at a particular SQ frequency, maximum intensity is observed for the same (or close to the same) recoupling time, but the relative intensity of the DQ peaks is a reliable indicator of the relative H-H distance-the ratio of the maximum intensities for the peaks at the CH(Ala) SQ frequency due to the two DQCs with the NH and OH protons are found to be approximately in the ratio of the squares of the corresponding dipolar coupling constants. While the simulated (1)H DQ build-up curves reproduce most of the features of the experimental curves, maximum intensity is often observed at a longer recoupling time in simulations. In this respect, simulations for two to eight spins show a trend towards a faster decay for an increasing number of considered spins. Finally, simulations show that increasing either the Larmor frequency (to 1 GHz) or the MAS frequency (to 125 kHz) does not lead to changes in the marked differences between the (1)H DQ build-up curves at the CH(Asp) SQ frequency for DQCs to the CH(2)(a) and OH protons that correspond to similar H-H distances (2.39 A and 2.49 A, respectively).

Computational Studies of a Protein-based Nanoactuator for Nanogripping Applications

The design hypothesis, architectures, and computational modeling of a novel peptide-based nanoactuator are presented in this paper. We engineered the α-helical coiled-coil portion of the yeast transcriptional activator peptide called GCN4 to obtain an environmentally responsive nanoactuator. The dimeric coiled-coil peptide consists of two identical approximately 4.5 nm long and approximately 3 nm wide polypeptide chains. The actuation mechanism depends on the modification of electrostatic charges along the peptide by varying the pH of the solution resulting in the reversible movement of helices and, therefore, creating the motion of an actuator. Using molecular dynamics simulations we showed that pH changes led to a reversible opening of up to 1.5 nm which is approximately 150% of the initial separation of the nanoactuator. We also investigated the forces generated by the nanoactuator upon pH actuation, using a new method based on a modified steered molecular dynamics technique. Owing to its open and close motion resembling that of tweezers, the new nanoactuator can potentially be used as a nanogripper in various nanomanipulation tasks such as detection and removal of heavy metal ions during nanofabrication processes or as a molecular switch.

Linear-scaling method for calculating nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals within Hartree-Fock and density-functional theory

Details of a new density matrix-based formulation for calculating nuclear magnetic resonance chemical shifts at both Hartree-Fock and density functional theory levels are presented. For systems with a nonvanishing highest occupied molecular orbital-lowest unoccupied molecular orbital gap, the method allows us to reduce the asymptotic scaling order of the computational effort from cubic to linear, so that molecular systems with 1000 and more atoms can be tackled with today's computers. The key feature is a reformulation of the coupled-perturbed self-consistent field (CPSCF) theory in terms of the one-particle density matrix (D-CPSCF), which avoids entirely the use of canonical MOs. By means of a direct solution for the required perturbed density matrices and the adaptation of linear-scaling integral contraction schemes, the overall scaling of the computational effort is reduced to linear. A particular focus of our formulation is to ensure numerical stability when sparse-algebra routines are used to obtain an overall linear-scaling behavior.

Selectivity of guest–host interactions in self-assembled hydrogen-bonded nanostructures observed by NMR

We studied the incorporation of various small guest molecules into calix[4]hydroquinone nanotubes and nanoclusters using solid-state proton NMR spectroscopy in combination with quantum chemical calculations. While the molecules exhibit different types of hydrogen bonding and van der Waals interactions, they show different affinities to the nanotube host structures. As the guest molecules are located inside the complexes, they experience a shift in the NMR resonance line caused by screening effects from the aromatic electrons of the host superstructure. The abilities to fill the otherwise empty space within the tubes can hence be measured indirectly by the displacement of the NMR lines relative to the free molecules. In this way, we can probe which guest molecules are recognized by the calix[4]hydroquinones as suitable for filling their nanoporous superstructures. Selective guest-host interactions have been successfully achieved for the three component mixture of water and acetone with either 2-methyl-2-propanol or 2-propanol. In both cases, the alcohols were superior to acetone in filling the CHQ tubes.

Molecular recognition in molecular tweezers systems: quantum-chemical calculation of NMR chemical shifts

Quantum-chemical calculations for molecular tweezers systems are presented, where the focus is not only on the recognition process in the host-guest systems, but on the self aggregation of the tweezers host as well. Such intermolecular interactions influence the corresponding NMR spectra strongly by up to 6 ppm for proton chemical shifts, since ring-current effects are particularly important. The quantum-chemical results allow one to reliably assign the spectra and to gain information both on the structure and on the importance of intra- and intermolecular interactions. In addition, we study the accuracy of a variety of density functionals for describing the present host-guest systems, where we observe a considerable underestimation of ring-current effects on (1)H NMR chemical shifts at the density functional theory (DFT) level using smaller basis sets such as 6-31G**, so that larger bases like TZP are required. This stands in contrast to the behavior of the Hartree-Fock scheme, where small basis sets, such as 6-31G**, provide reliable (1)H NMR shieldings for molecular tweezers systems.