Bulk Chemical Shifts in Hydrogen-Bonded Systems from First-Principles Calculations and Solid-State-NMR

Cyanophycin modifications for applications in tissue scaffolding

Cyanophycin (CGP) is a polypeptide consisting of amino acids-aspartic acid in the backbone and arginine in the side chain. Owing to its resemblance to cell adhesive motifs in the body, it can be considered suitable for use in biomedical applications as a novel component to facilitate cell attachment and tissue regeneration. Although it has vast potential applications, starting with nutrition, through drug delivery and tissue engineering to the production of value-added chemicals and biomaterials, CGP has not been brought to the industry yet. To develop scaffolds using CGP powder produced by bacteria, its properties (e.g., biocompatibility, morphology, biodegradability, and mechanical strength) should be tailored in terms of the requirements of the targeted tissue. Crosslinking commonly stands for a primary modification method for renovating biomaterial features to these extents. Herein, we aimed to crosslink CGP for the first time and present a comparative study of different methods of CGP crosslinking including chemical, physical, and enzymatic methods by utilizing glutaraldehyde (GTA), UV exposure, genipin, 1-ethyl-3-[3-dimethylaminopropyl] carbodiimide hydrochloride/N-hydroxysuccinimide (EDC/NHS), and monoamine oxidase (MAO). Crosslinking efficacy varied among the samples crosslinked via the different crosslinking methods. All crosslinked CGP were non-cytotoxic to L929 cells, except for the groups with higher GTA concentrations. We conclude that CGP is a promising candidate for scaffolding purposes to be used as part of a composite with other biomaterials to maintain the integrity of scaffolds. The initiative study demonstrated the unknown characteristics of crosslinked CGP, even though its feasibility for biomedical applications should be confirmed by further examinations. KEY POINTS: • Cyanophycin was crosslinked by 5 different methods • Crosslinked cyanophycin is non-cytotoxic to L929 cells • Crosslinked cyanophycin is a promising new material for scaffolding purposes.

Distinguishing between COOH, COO- , and hydrogen disordered COOH sites in solids with 13 C chemical shift anisotropy and T1 measurements

Since 1993, it has been known that 13 C chemical shift tensor (i.e., δ11 , δ22 , and δ33 ) provides information sufficient to distinguish between COOH and COO- sites. Herein, four previously unreported metrics are proposed for differentiating COOH/COO- moieties. A new relationship is also introduced that correlates the asymmetry (i.e., δ11 -δ22 ) of COOH sites to the proximity of hydrogen bond donating partners within 2.6 Å with high accuracy (±0.05 Å). Conversely, a limitation to all proposed metrics is that they fail to distinguish between COO- and hydrogen disordered COOH sites. To reconcile this omission, a new approach is proposed based on T1 measurements of both 1 H and 13 C. The 13 C T1 values are particularly sensitive with the T1 for hydrogen disordered COOH moieties found to be nearly six times smaller than T1 's from COO- sites.

Opportunities and Challenges in Applying Solid-State NMR Spectroscopy in Organic Mechanochemistry

In recent years it is shown that mechanochemical strategies can be beneficial in directed conversions of organic compounds. Finding new reactions proved difficult, and due to the lack of mechanistic understanding of mechanochemical reaction events, respective efforts have mostly remained empirical. Spectroscopic techniques are crucial in shedding light on these questions. In this overview, the opportunities and challenges of solid-state nuclear magnetic resonance (NMR) spectroscopy in the field of organic mechanochemistry are discussed. After a brief discussion of the basics of high-resolution solid-state NMR under magic-angle spinning (MAS) conditions, seven opportunities for solid-state NMR in the field of organic mechanochemistry are presented, ranging from ex situ approaches to structurally elucidated reaction products obtained by milling to the potential and limitations of in situ solid-state NMR approaches. Particular strengths of solid-state NMR, for instance in differentiating polymorphs, in NMR-crystallographic structure-determination protocols, or in detecting weak noncovalent interactions in molecular-recognition events employing proton-detected solid-state NMR experiments at fast MAS frequencies, are discussed.

van der Wel P. Selective observation of semi-rigid non-core residues in dynamically complex mutant huntingtin protein fibrils

•

Mutant huntingtin exon 1 fibrils feature a broad range of molecular dynamics.

•

Molecular motion is coupled to water dynamics outside the fiber core.

•

Dynamics-based spectral editing ssNMR reveals mobile non-core residues.

•

Intermediate-motion selection via dipolar dephasing of rigid sites.

•

Semi-mobile glutamines outside the fiber core observed and identified.

Many amyloid-forming proteins, which are normally intrinsically disordered, undergo a disorder-to-order transition to form fibrils with a rigid β-sheet core flanked by disordered domains. Solid-state NMR (ssNMR) and cryogenic electron microscopy (cryoEM) excel at resolving the rigid structures within amyloid cores but studying the dynamically disordered domains remains challenging. This challenge is exemplified by mutant huntingtin exon 1 (HttEx1), which self-assembles into pathogenic neuronal inclusions in Huntington disease (HD). The mutant protein’s expanded polyglutamine (polyQ) segment forms a fibril core that is rigid and sequestered from the solvent. Beyond the core, solvent-exposed surface residues mediate biological interactions and other properties of fibril polymorphs. Here we deploy magic angle spinning ssNMR experiments to probe for semi-rigid residues proximal to the fibril core and examine how solvent dynamics impact the fibrils’ segmental dynamics. Dynamic spectral editing (DYSE) 2D ssNMR based on a combination of cross-polarization (CP) ssNMR with selective dipolar dephasing reveals the weak signals of solvent-mobilized glutamine residues, while suppressing the normally strong background of rigid core signals. This type of ‘intermediate motion selection’ (IMS) experiment based on cross-polarization (CP) ssNMR, is complementary to INEPT- and CP-based measurements that highlight highly flexible or highly rigid protein segments, respectively. Integration of the IMS-DYSE element in standard CP-based ssNMR experiments permits the observation of semi-rigid residues in a variety of contexts, including in membrane proteins and protein complexes. We discuss the relevance of semi-rigid solvent-facing residues outside the fibril core to the latter’s detection with specific dyes and positron emission tomography tracers.

A toolbox for improving the workflow of NMR crystallography

NMR crystallography is a powerful tool with applications in structural characterization and crystal structure verification, to name two. However, applying this tool presents several challenges, especially for industrial users, in terms of consistency, workflow, time consumption, and the requirement for a high level of understanding of experimental solid-state NMR and GIPAW-DFT calculations. Here, we have developed a series of fully parameterized scripts for use in Materials Studio and TopSpin, based on the .magres file format, with a focus on organic molecules (e.g. pharmaceuticals), improving efficiency, robustness, and workflow. We separate these tools into three major categories: performing the DFT calculations, extracting & visualizing the results, and crystallographic modelling. These scripts will rapidly submit fully parameterized CASTEP jobs, extract data from the calculations, assist in visualizing the results, and expedite the process of structural modelling. Accompanied with these tools is a description on their functionality, documentation on how to get started and use the scripts, and links to video tutorials for guiding new users. Through the use of these tools, we hope to facilitate NMR crystallography and to harmonize the process across users.

5-amino-2-methylpyridinium hydrogen fumarate: An XRD and NMR crystallography analysis

Single-crystal X-ray diffraction structures of the 5-amino-2-methylpyridinium hydrogen fumarate salt have been solved at 150 and 300 K (CCDC 1952142 and 1952143). A base-acid-base-acid ring is formed through pyridinium-carboxylate and amine-carboxylate hydrogen bonds that hold together chains formed from hydrogen-bonded hydrogen fumarate ions. ¹ H and ¹³ C chemical shifts as well as ¹⁴ N shifts that additionally depend on the quadrupolar interaction are determined by experimental magic angle spinning (MAS) solid-state nuclear magnetic resonance (NMR) and gauge-including projector-augmented wave (GIPAW) calculation. Two-dimensional homonuclear ¹ H-¹ H double-quantum (DQ) MAS and heteronuclear ¹ H-¹³ C and ¹⁴ N-¹ H spectra are presented. Only small differences of up to 0.1 and 0.6 ppm for ¹ H and ¹³ C are observed between GIPAW calculations starting with the two structures solved at 150 and 300 K (after geometry optimisation of atomic positions, but not unit cell parameters). A comparison of GIPAW-calculated ¹ H chemical shifts for isolated molecules and the full crystal structures is indicative of hydrogen bonding strength.

NMR crystallography of molecular organics

Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.

Hydrogen bonding networks in gabapentin protic pharmaceutical salts: NMR and in silico studies

Hydrogen bonds (HBs) play a key role in the supramolecular arrangement of crystalline solids and, although they have been extensively studied, the influence of their strength and geometry on crystal packing remains poorly understood. Here we describe the crystal structures of two novel protic gabapentin (GBP) pharmaceutical salts prepared with the coformers methanesulfonic acid (GBP:METHA) and ethanesulfonic acid (GBP:ETHA). This study encompasses experimental and computational electronic structure analyses of ¹ H NMR chemical shifts (CSs), upon in silico HB cleavage. GBP:METHA and GBP:ETHA crystal packing comprise two main structural domains: an ionic layer (characterized by the presence of charge-assisted ⁺ NH_GBP ⋯O^-_METHA/ETHA HB interactions) and a neutral layer generated in a different way for each salt, mainly due to the presence of bifurcated HB interactions. A comprehensive study of HB networks is presented for GBP:METHA, by isolating molecular fragments involved in distinct HB types (NH⋯O, OH⋯O, and CH⋯O) obtained from in silico disassembling of an optimized three-dimensional packing structure. Formation of HB leads to calculated ¹ H NMR CS changes from 0.4 to ~5.8 ppm. This study further attempts to assess how ¹ H NMR CS of protons engaged in certain HB are affected when other nearby HB, involving bifurcated or geminal/vicinal hydrogen atoms, are removed.

Investigation of structure and dynamics in a photochromic molecular crystal by NMR crystallography

A photochromic anil, N-(3,5-di-t-butylsalicylidene)-4-amino-pyridine, has been studied by single-crystal X-ray diffraction, multinuclear magic-angle spinning NMR, and first-principles density functional theory (DFT) calculations. Interpretation of the solid-state NMR data on the basis of calculated chemical shifts confirms the structure is primarily composed of molecules in the ground-state enol tautomer, whereas thermally activated cis-keto and photoisomerised trans-keto states exist as low-level defects with populations that are too low to detect experimentally. Variable temperature ¹³ C NMR data reveal evidence for solid-state dynamics, which is found to be associated with fast rotational motion of t-butyl groups and 180° flips of the pyridine ring, contrasting the time-averaged structure obtained by X-ray diffraction. Comparison of calculated chemical shifts for the full crystal structure and an isolated molecule also reveals evidence for an intermolecular hydrogen bond involving the pyridine ring and an adjacent imine carbon, which facilitates the flipping motion. The DFT calculations also reveal that the molecular conformation in the crystal structure is very close to the energetic minimum for an isolated molecule, indicating that the ring dynamics arise as a result of considerable steric freedom of the pyridine ring and which also allows the molecule to adopt a favourable conformation for photochromism.

An NMR crystallography investigation of furosemide

This paper presents an NMR crystallography study of three polymorphs of furosemide. Experimental magic-angle spinning (MAS) solid-state NMR spectra are reported for form I of furosemide, and these are assigned using density-functional theory (DFT)-based gauge-including projector augmented wave (GIPAW) calculations. Focusing on the three known polymorphs, we examine the changes to the NMR parameters due to crystal packing effects. We use a recently developed formalism to visualise which regions are responsible for the chemical shielding of particular sites and hence understand the variation in NMR parameters between the three polymorphs.

An XRD and NMR crystallographic investigation of the structure of 2,6-lutidinium hydrogen fumarate

A crystallographic study highlighting the benefits of a combined XRD and NMR approach in investigating both stability and variation within an organic multicomponent crystal.

Predicting anisotropic thermal displacements for hydrogens from solid-state NMR: a study on hydrogen bonding in polymorphs of palmitic acid

The hydrogen-bonding environments at the COOH moiety in eight polycrystalline polymorphs of palmitic acid are explored using solid-state NMR. Although most phases have no previously reported crystal structure, measured 13C chemical shift tensors for COOH moieties, combined with DFT modeling establish that all phases crystallize with a cyclic dimer (R22(8)) hydrogen bonding arrangement. Phases A2, Bm and Em have localized OH hydrogens while phase C has a dynamically disordered OH hydrogen. The phase designated As is a mix of five forms, including 27.4% of Bm and four novel phases not fully characterized here due to insufficient sample mass. For phases A2, Bm, Em, and C the anisotropic uncertainties in the COOH hydrogen atom positions are established using a Monte Carlo sampling scheme. Sampled points are retained or rejected at the ±1σ level based upon agreement of DFT computed 13COOH tensors with experimental values. The collection of retained hydrogen positions bear a remarkable resemblance to the anisotropic displacement parameters (i.e. thermal ellipsoids) from diffraction studies. We posit that this similarity is no mere coincidence and that the two are fundamentally related. The volumes of NMR-derived anisotropic displacement ellipsoids for phases with localized OH hydrogens are 4.1 times smaller than those derived from single crystal X-ray diffraction and 1.8 times smaller than the volume of benchmark single crystal neutron diffraction values.

Ab-Initio Molecular Dynamics Simulations and Calculations of Spectroscopic Parameters in Hydrogen-Bonding Liquids in Confinement (Project 8)

We investigate the effect of several nanoscale confinements on structural and dynamical properties of liquid water and binary aqueous mixtures. By means of molecular dynamics simulations based on density functional theory and atomistic force fields. Our main focus is on the dependence on the structure and the hydrogen-bonding-network of the liquids near the confinement interface at atomistic resolution. As a complementary aspect, spatially resolved profiles of the proton NMR chemical shift values are used to quantify the local strength of the hydrogen-bond-network.

Role of Anomalous Water Constraints in the Efficacy of Pharmaceuticals Probed by 1H Solid-State NMR

Water plays a complex and central role in determining the structural and reactive properties in numerous chemical systems. In crystalline materials with structural water, the primary focus is often to relate hydrogen bonding motifs to functional properties such as solubility, which is highly relevant in pharmaceutical applications. Nevertheless, understanding the full electrostatic landscape is necessary for a complete structure-function picture. Herein, a combination of tools including 1H magic angle spinning NMR and X-ray crystallography are employed to evaluate the local landscape of water in crystalline hydrates. Two hydrates of an anti-leukemia drug mercaptopurine, which exhibit dramatically different dehydration temperatures (by 90°C) and a three-fold difference in the in vivo bioavailability, are compared. The results identify an electrosteric caging mechanism for a kinetically trapped water in the hemihydrate form, which is responsible for the dramatic differences in properties.

Single-crystal X-ray diffraction and NMR crystallography of a 1:1 cocrystal of dithianon and pyrimethanil

A single-crystal X-ray diffraction structure of a 1:1 cocrystal of two fungicides, namely dithianon (DI) and pyrimethanil (PM), is reported [systematic name: 5,10-dioxo-5H,10H-naphtho[2,3-b][1,4]dithiine-2,3-dicarbonitrile-4,6-dimethyl-N-phenylpyrimidin-2-amine (1/1), C14H4N2O2S2·C12H13N2]. Following an NMR crystallography approach, experimental solid-state magic angle spinning (MAS) NMR spectra are presented together with GIPAW (gauge-including projector augmented wave) calculations of NMR chemical shieldings. Specifically, experimental 1H and 13C chemical shifts are determined from two-dimensional 1H-13C MAS NMR correlation spectra recorded with short and longer contact times so as to probe one-bond C-H connectivities and longer-range C...H proximities, whereas H...H proximities are identified in a 1H double-quantum (DQ) MAS NMR spectrum. The performing of separate GIPAW calculations for the full periodic crystal structure and for isolated molecules allows the determination of the change in chemical shift upon going from an isolated molecule to the full crystal structure. For the 1H NMR chemical shifts, changes of 3.6 and 2.0 ppm correspond to intermolecular N-H...O and C-H...O hydrogen bonding, while changes of -2.7 and -1.5 ppm are due to ring current effects associated with C-H...π interactions. Even though there is a close intermolecular S...O distance of 3.10 Å, it is of note that the molecule-to-crystal chemical shifts for the involved sulfur or oxygen nuclei are small.

Phycocyanobilin in solution – a solvent triggered molecular switch

The chromophore phycocyanobilin changes its spectroscopic behaviour upon solvent change. Our calculations trace this effect back to conformational switching, induced by changes in the hydrogen bonding pattern.

Selective measurements of long-range homonuclear J-couplings in solid-state NMR

We demonstrate here that the principle of frequency-selective spin-echoes can be extended to the measurements of long-range homonuclear scalar J-couplings in the solid-state. Singly or doubly frequency-selective pulses were used to generate either a J-modulated experiment (S) or a reference experiment (S0). The combination of these two distinct experiments provides experimental data that, in favorable cases, are insensitive to incoherent relaxation effects, and which can be used to estimate long-range homonuclear J-couplings in multiple spin-systems. The concept is illustrated in the case of a uniformly (13)C and (15)N labeled sample of L-histidine, where the absolute value of homonuclear J-couplings between two spins separated by one, two or three covalent bonds are measured. Moreover, we show that a (2)J((15)N-C-(15)N) coupling as small as 0.9 Hz can be precisely measured with the method presented here.

Probing the relation between charge transport and supramolecular organization down to ångström resolution in a benzothiadiazole-cyclopentadithiophene copolymer

Molecular modeling shows that longitudinal displacement of the backbones by a couple of ångströms has a profound impact on the electronic coupling mediating charge transport in a conjugated copolymer. These changes can be probed by monitoring the calculated X-ray scattering patterns and NMR chemical shifts as a function of sliding of the polymer chains and comparing them to experiment.

Interpretation of crystal effects on NMR chemical shift tensors: electron and shielding deformation densities

The relationship between the NMR observables and the supramolecular structure of any system is not straightforward. In this work we examine the influence of the crystal packing for three purine derivatives (hypoxanthine, theobromine, and 6-(2-methoxy)benzylaminopurine) on the principal components of the NMR chemical shift tensors (CSTs). We employ density functional calculations to obtain various molecular properties (the ground-state electron density, the magnitudes and orientations of the components of NMR chemical shift tensor, and the spatial distribution of the isotropic magnetic shielding) for the isolated molecules and for the molecules embedded in supramolecular clusters modeling the crystal environment and evaluate their differences. The concept has enabled us to rationalize the effect of the crystal packing on the NMR CSTs in terms of the redistribution of the ground-state electron density induced by intermolecular interactions in the solid state.

First Principles Calculations of NMR Chemical Shifts of Liquid Water at an Amorphous Silica Interface

We investigate the anomalous structure and hydrogen bond network of water molecules confined inside a silica nanopore (MCM-41 type). In addition to geometric data, we use proton NMR chemical shifts as a measure for the strength of the H-bonding network. We compute the 1H NMR shifts of confined water based on a first principle approach in the framework of density functional perturbation theory under periodic boundary conditions. The hydrophilic character of the silica is well manifested in the water density profile. Our calculations illustrate both the modifications of the 1H NMR chemical shifts of the water with respect to bulk water and a considerable slowing down of water diffusion. In the vicinity of silanols, weakly hydrogen bonded liquid water is observed, while at the center region of the pore, the hydrogen bonding network is enhanced with respect to bulk water.

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.

Identifying guanosine self assembly at natural isotopic abundance by high-resolution 1H and 13C solid-state NMR spectroscopy

By means of the (1)H chemical shifts and the proton-proton proximities as identified in (1)H double-quantum (DQ) combined rotation and multiple-pulse spectroscopy (CRAMPS) solid-state NMR correlation spectra, ribbon-like and quartet-like self-assembly can be identified for guanosine derivatives without isotopic labeling for which it was not possible to obtain single crystals suitable for diffraction. Specifically, characteristic spectral fingerprints are observed for dG(C10)(2) and dG(C3)(2) derivatives, for which quartet-like and ribbon-like self-assembly has been unambiguously identified by (15)N refocused INADEQUATE spectra in a previous study of (15)N-labeled derivatives (Pham, T. N.; et al. J. Am. Chem. Soc.2005, 127, 16018). The NH (1)H chemical shift is observed to be higher (13-15 ppm) for ribbon-like self-assembly as compared to 10-11 ppm for a quartet-like arrangement, corresponding to a change from NH···N to NH···O intermolecular hydrogen bonding. The order of the two NH(2)(1)H chemical shifts is also inverted, with the NH(2) proton closest in space to the NH proton having a higher or lower (1)H chemical shift than that of the other NH(2) proton for ribbon-like as opposed to quartet-like self-assembly. For the dG(C3)(2) derivative for which a single-crystal diffraction structure is available, the distinct resonances and DQ peaks are assigned by means of gauge-including projector-augmented wave (GIPAW) chemical shift calculations. In addition, (14)N-(1)H correlation spectra obtained at 850 MHz under fast (60 kHz) magic-angle spinning (MAS) confirm the assignment of the NH and NH(2) chemical shifts for the dG(C3)(2) derivative and allow longer range through-space N···H proximities to be identified, notably to the N7 nitrogens on the opposite hydrogen-bonding face.

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.