A Review on Combination of Ab Initio Molecular Dynamics and NMR Parameters Calculations

The EFG Rosetta Stone: translating between DFT calculations and solid state NMR experiments

We present a comprehensive study on the best practices for integrating first principles simulations in experimental quadrupolar solid-state nuclear magnetic resonance (SS-NMR), exploiting the synergies between theory and experiment for achieving the optimal interpretation of both. Most high performance materials (HPMs), such as battery electrodes, exhibit complex SS-NMR spectra due to dynamic effects or amorphous phases. NMR crystallography for such challenging materials requires reliable, accurate, efficient computational methods for calculating NMR observables from first principles for the transfer between theoretical material structure models and the interpretation of their experimental SS-NMR spectra. NMR-active nuclei within HPMs are routinely probed by their chemical shielding anisotropy (CSA). However, several nuclear isotopes of interest, e.g.7Li and 27Al, have a nuclear quadrupole and experience additional interactions with the surrounding electric field gradient (EFG). The quadrupolar interaction is a valuable source of information about atomistic structure, and in particular, local symmetry, complementing the CSA. As such, there is a range of different methods and codes to choose from for calculating EFGs, from all-electron to plane wave methods. We benchmark the accuracy of different simulation strategies for computing the EFG tensor of quadrupolar nuclei with plane wave density functional theory (DFT) and study the impact of the material structure as well as the details of the simulation strategy. Especially for small nuclei with few electrons, such as 7Li, we show that the choice of physical approximations and simulation parameters has a large effect on the transferability of the simulation results. To the best of our knowledge, we present the first comprehensive reference scale and literature survey for 7Li quadrupolar couplings. The results allow us to establish practical guidelines for developing the best simulation strategy for correlating DFT to experimental data extracting the maximum benefit and information from both, thereby advancing further research into HPMs.

Unexpected Intermolecular C-H···O Hydrogen Bonds and 1H NMR Chemical Shifts in a Key Linker for Fluorine-18 Labeling of Dimeric Drugs

The compound 2-{[(trifluoromethyl)sulfonyl]oxy}propane-1,3-diyl bis(4-methylbenzenesulfonate) (TPB) is a crucial intermediate in the synthesis of 18F-radiolabeled cromolyn derivatives. In this work, we combine 1H NMR spectroscopy, X-ray crystallography, ab initio molecular dynamics, and NMR calculations to examine the structure, interactions, and solvation dynamics of the TPB molecule. In CDCl3, the CH2 groups within its glyceryl-derived linker exhibit a single set of proton signals in the 1H NMR measurements. However, when TPB is dissolved in DMSO-d6, distinct splitting patterns emerge despite its seemingly symmetric chemical structure. Crystallographic analysis further unveils the absence of overall symmetry in its three-dimensional arrangement. To elucidate these unique NMR features, we carry out ab initio molecular dynamics simulations and characterize the solvation structures and dynamics of TPB in CHCl3 and DMSO solutions. In contrast to the predominantly nonpolar nature of the CHCl3 solvents, DMSO directly participates in C-H···O hydrogen-bonding interactions with the solute molecule, leading to the splitting of its -CH2 chemical shifts into two distinct distributions. The comprehensive understanding of the structure and solvation interactions of TPB provides essential insights into its application in the radiofluorination reactions of cromolyn derivatives and holds promise for the future development of radiolabeled dimeric drugs.

High-Resolution 17O Solid-State NMR as a Unique Probe for Investigating Oxalate Binding Modes in Materials: The Case Study of Calcium Oxalate Biominerals

Oxalate ligands are found in many classes of materials, including energy storage materials and biominerals. Determining their local environments at the atomic scale is thus paramount to establishing the structure and properties of numerous phases. Here, we show that high-resolution 17O solid-state NMR is a valuable asset for investigating the structure of crystalline oxalate systems. First, an efficient 17O-enrichment procedure of oxalate ligands is demonstrated using mechanochemistry. Then, 17O-enriched oxalates were used for the synthesis of the biologically relevant calcium oxalate monohydrate (COM) phase, enabling the analysis of its structure and heat-induced phase transitions by high-resolution 17O NMR. Studies of the low-temperature COM form (LT-COM), using magnetic fields from 9.4 to 35.2 T, as well as 13C-17O MQ/D-RINEPT and 17O{1H} MQ/REDOR experiments, enabled the 8 inequivalent oxygen sites of the oxalates to be resolved, and tentatively assigned. The structural changes upon heat treatment of COM were also followed by high-resolution 17O NMR, providing new insight into the structures of the high-temperature form (HT-COM) and anhydrous calcium oxalate α-phase (α-COA), including the presence of structural disorder in the latter case. Overall, this work highlights the ease associated with 17O-enrichment of oxalate oxygens, and how it enables high-resolution solid-state NMR, for "NMR crystallography" investigations.

NMR analysis of structural geometry and molecular dynamics in perovskite-type N(CH3)4CdBr3 crystal near high-temperature phase transition

The NMR chemical shifts, linewidths, spin-lattice relaxation times in the rotating system T1ρ, and spin-lattice relaxation times in the laboratory system T1 were evaluated for the perovskite-type N(CH3)4CdBr3 crystal, aiming to understand the changes in the structural geometry and molecular dynamics from phase I to phase II. From the temperature-dependence of the 1H, 13C, 14N, and 113Cd NMR chemical shifts, the structural geometry underwent a continuous change, without anomalous changes around (TC = 390 K). However, the linewidths in phase I were narrower than those in phase II, indicating that the motional averaging effects were caused by the rapid rotation of the N(CH3)4 group. Sudden changes in T1 and T1ρ were observed near TC, for which the activation energy Ea in phase I was approximately 12 times larger than that in phase II; the small Ea values in phase II indicate a large degree of freedom for the methyl group and CdBr6 octahedra, whereas the large Ea in phase I was primarily attributed to the overall N(CH3)4 and the 113Cd in the CdBr6 groups. Consequently, the phase transition mechanisms of N(CH3)4CdBr3 are related to reorientation of the N(CH3)4 group and the arrangement of the CdBr6 groups.

Structures and Dynamics of Complex Guest Molecules in Confinement, Revealed by Solid-State NMR, Molecular Dynamics, and Calorimetry

This review gives an overview of current trends in the investigation of confined molecules such as water, small and higher alcohols, carbonic acids, ethylene glycol, and non-ionic surfactants, such as polyethylene glycol or Triton-X, as guest molecules in neat and functionalized mesoporous silica materials employing solid-state NMR spectroscopy, supported by calorimetry and molecular dynamics simulations. The combination of steric interactions, hydrogen bonds, and hydrophobic and hydrophilic interactions results in a fascinating phase behavior in the confinement. Combining solid-state NMR and relaxometry, DNP hyperpolarization, molecular dynamics simulations, and general physicochemical techniques, it is possible to monitor these confined molecules and gain deep insights into this phase behavior and the underlying molecular arrangements. In many cases, the competition between hydrogen bonding and electrostatic interactions between polar and non-polar moieties of the guests and the host leads to the formation of ordered structures, despite the cramped surroundings inside the pores.

Molecular dynamics and NMR reveal the coexistence of H-bond-assisted and through-space JFH coupling in fluorinated amino alcohols

The JFH coupling constants in fluorinated amino alcohols were investigated through experimental and theoretical approaches. The experimental JFH couplings were only reproduced theoretically when explicit solvation through molecular dynamics (MD) simulations was conducted in DMSO as the solvent. The combination of MD conformation sampling and DFT NMR spin-spin coupling calculations for these compounds reveals the simultaneous presence of through-space (TS) and hydrogen bond (H-bond) assisted JFH coupling between fluorine and hydrogen of the NH group. Furthermore, MD simulations indicate that the hydrogen in the amino group participates in both an intermolecular bifurcated H-bond with DMSO and in transmitting the observed JFH coupling. The contribution of TS to the JFH coupling is due to the spatial proximity of the fluorine and the NH group, aided by a combination of the non-bonding transmission pathway and the hydrogen bonding pathway. The experimental JFH coupling observed for the molecules studied should be represented as 4TS/1hJFH coupling.

Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment

NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.

Dynamic and relativistic effects on Pt-Pt indirect spin-spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent-solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin-orbit Kohn-Sham (2c-ZKS) and four-component Dirac-Kohn-Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1-Pt2(NH3)4(Am)2-L2]n+ (Am = α-pyrrolidonate and pivalamidate; L = H2O, Cl-, and Br-) are studied. The computed Pt-Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N-C-O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X-) and aquahalo (L1 = H2O; L2 = X-) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

The accuracy limit of chemical shift predictions for species in aqueous solution

Interpreting NMR experiments benefits from first-principles predictions of chemical shifts. Reaching the accuracy limit of theory is relevant for unambiguous structural analysis and dissecting theoretical approximations. Since accurate chemical shift measurements are based on using internal reference compounds such as trimethylsilylpropanesulfonate (DSS), a detailed comparison of experimental with theoretical data requires simultaneous consideration of both target and reference species ensembles in the same solvent environment. Here we show that ab initio molecular dynamics simulations to generate liquid-state ensembles of target and reference compounds, including explicitly their short-range solvation environments and combined with quantum-mechanical solvation models, allows for predicting highly accurate 1H (∼0.1-0.5 ppm) and aliphatic 13C (∼1.5 ppm) chemical shifts for aqueous solutions of the model compounds trimethylamine N-oxide (TMAO) and N-methylacetamide (NMA), referenced to DSS without any system-specific adjustments. This encompasses the two peptide bond conformations of NMA identified by NMR. The results are used to derive a general-purpose guideline set for predictive NMR chemical shift calculations of NMA in the liquid state and to identify artifacts of force field models. Accurate predictions are only obtained if a sufficient number of explicit water molecules is included in the quantum-mechanical calculations, disproving a purely electrostatic model of the solvent effect on chemical shifts.

Multiscale simulation of fluids: coupling molecular and continuum

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.

Rapid Measurement of Heteronuclear Coupling Constants in Complex NMR Spectra

The NMR analysis of fluorine-containing molecules, increasingly widespread due to their importance in pharmaceuticals and biochemistry, poses significant challenges. Severe peak overlap in the proton spectrum often hinders the extraction of critical structural information in the form of heteronuclear scalar coupling constants, which are crucial for determining pharmaceutical properties and biological activity. Here, a new method, IPAP-FESTA, is reported that drastically simplifies measurements of the signs and magnitudes of proton-fluorine couplings. Its usefulness is demonstrated for the structural study of the steroidal drug fluticasone propionate extracted from a commercial formulation and for assessing solvent effects on the conformational equilibrium in a physically inseparable fluorohydrin mixture.

The need for operando modelling of 27Al NMR in zeolites: the effect of temperature, topology and water

Solid state (ss-) 27Al NMR is one of the most valuable tools for the experimental characterization of zeolites, owing to its high sensitivity and the detailed structural information which can be extracted from the spectra. Unfortunately, the interpretation of ss-NMR is complex and the determination of aluminum distributions remains generally unfeasible. As a result, computational modelling of 27Al ss-NMR spectra has grown increasingly popular as a means to support experimental characterization. However, a number of simplifying assumptions are commonly made in NMR modelling, several of which are not fully justified. In this work, we systematically evaluate the effects of various common models on the prediction of 27Al NMR chemical shifts in zeolites CHA and MOR. We demonstrate the necessity of operando modelling; in particular, taking into account the effects of water loading, temperature and the character of the charge-compensating cation. We observe that conclusions drawn from simple, high symmetry model systems such as CHA do not transfer well to more complex zeolites and can lead to qualitatively wrong interpretations of peak positions, Al assignment and even the number of signals. We use machine learning regression to develop a simple yet robust relationship between chemical shift and local structural parameters in Al-zeolites. This work highlights the need for sophisticated models and high-quality sampling in the field of NMR modelling and provides correlations which allow for the accurate prediction of chemical shifts from dynamical simulations.

DFT approach to predict 13C NMR chemical shifts of hydrocarbon species adsorbed on Zn-modified zeolite

13C MAS NMR spectroscopy is a powerful technique to study the mechanisms of hydrocarbon transformations on heterogeneous catalysts. It can reliably identify the surface intermediates and the adsorbed products based...

The Physical Chemistry and Chemical Physics (PCCP) Section of the International Journal of Molecular Sciences in Its Publications: The First 300 Thematic Articles in the First 3 Years

The Physical Chemistry and Chemical Physics Section (PCCP Section) is one of the youngest among the sections of the International Journal of Molecular Sciences (IJMS)-the year 2021 will only mark three years since its inception [...].

Application of Molecular Dynamics Simulations in the Analysis of Cyclodextrin Complexes

Cyclodextrins (CDs) are highly respected for their ability to form inclusion complexes via host-guest noncovalent interactions and, thus, ensofance other molecular properties. Various molecular modeling methods have found their applications in the analysis of those complexes. However, as showed in this review, molecular dynamics (MD) simulations could provide the information unobtainable by any other means. It is therefore not surprising that published works on MD simulations used in this field have rapidly increased since the early 2010s. This review provides an overview of the successful applications of MD simulations in the studies on CD complexes. Information that is crucial for MD simulations, such as application of force fields, the length of the simulation, or solvent treatment method, are thoroughly discussed. Therefore, this work can serve as a guide to properly set up such calculations and analyze their results.

Structure Solution of Nano-Crystalline Small Molecules Using MicroED and Solid-State NMR Dipolar-Based Experiments

Three-dimensional electron diffraction crystallography (microED) can solve structures of sub-micrometer crystals, which are too small for single crystal X-ray crystallography. However, R factors for the microED-based structures are generally high because of dynamic scattering. That means R factor may not be reliable provided that kinetic analysis is used. Consequently, there remains ambiguity to locate hydrogens and to assign nuclei with close atomic numbers, like carbon, nitrogen, and oxygen. Herein, we employed microED and ssNMR dipolar-based experiments together with spin dynamics numerical simulations. The NMR dipolar-based experiments were 1H-14N phase-modulated rotational-echo saturation-pulse double-resonance (PM-S-RESPDOR) and 1H-1H selective recoupling of proton (SERP) experiments. The former examined the dephasing effect of a specific 1H resonance under multiple 1H-14N dipolar couplings. The latter examined the selective polarization transfer between a 1H-1H pair. The structure was solved by microED and then validated by evaluating the agreement between experimental and calculated dipolar-based NMR results. As the measurements were performed on 1H and 14N, the method can be employed for natural abundance samples. Furthermore, the whole validation procedure was conducted at 293 K unlike widely used chemical shift calculation at 0 K using the GIPAW method. This combined method was demonstrated on monoclinic l-histidine.