Chemical Shielding Constants Research Articles

New Pharmaceutical Salts of Trazodone.

New pharmaceutically acceptable salts of trazodone (trazodone hydrogen bromide and trazodone 1-hydroxy-2-naphthonic acid) for the treatment of central nervous system disorders are synthesized and described. Although trazodone salts are poorly crystalline, single-crystal X-ray diffraction data for trazodone 1-hydroxy-2-naphthonic acid were collected and analyzed as well as compared to the previously described crystal structure of commercially available trazodone hydrochloride. The powder samples of all new salts were characterized by Fourier transform infrared spectroscopy, X-ray diffraction and 13C solid-state nuclear magnetic resonance spectroscopy. Spectroscopic studies were supported by gauge including projector augmented wave (GIPAW) calculations of carbon chemical shielding constants. The main goal of our research was to find salts with better physicochemical properties and to make an attempt to associate them with both the anion structure and the most prominent interactions exhibited by the protonated trazodone cation. The dissolution profiles of trazodone from tablets prepared from various salts with lactose monohydrate were investigated. The studies revealed that salts with simple anions show a fast release of the drug while the presence of more complex anion, more strongly interacting with the cation, effects a slow-release profile of the active substance and can be used for the preparation of the tables with a delay or prolonged mode of action.

Structure Identification of Novel Compounds Using Simple IR, 1H, and 13C NMR Spectroscopy and Computational Tools

A simple method is suggested to identify the structure of novel compounds with basic IR, 1H, and 13C nuclear magnetic resonance (NMR) spectroscopy and computational tools. With the molecular formula obtained from high‐resolution mass spectrometry, all possible isomers are calculated. The absence or presence of particular functional groups which were inferred from IR, 1H, and 13C NMR and DEPT spectra are used to sort all the possible isomers using SMARTS code and RDKit. Then, the final structure(s) are identified by comparing the correlation between Quantum mechanically calculated 13C NMR chemical shielding constants and experimental 13C chemical shifts of molecules in consideration. We have applied this protocol to five natural compounds, the number of heavy atoms ranging from 11 to 15, and correctly identified the structures of all test cases. The limitations and further consideration are discussed.

Quantum chemical calculations of 31P NMR chemical shifts of P-donor ligands in platinum(II) complexes.

This work aims to find the most suitable method that is practically applicable for the calculation of 31P NMR chemical shifts of Pt(II) complexes. The influence of various all-electron and ECP basis sets, DFT functionals, and solvent effects on the optimized geometry was tested. A variety of combinations of DFT functionals BP86, B3LYP, PBE0, TPSSh, CAM-B3LYP, and ωB97XD with all-electron basis sets 6-31G, 6-31G(d), 6-31G(d,p), 6-311G(d,p), and TZVP and ECP basis sets SDD, LanL2DZ, and CEP-31G were used. Chemical shielding constants were then calculated using BP86, PBE0, and B3LYP functionals in combination with the TZ2P basis. The magnitude of spin-orbit interactions was also evaluated.

Calculations of Geometric Parameters and Physicochemical Properties of Complexes Formed of FE(II)-Reactive 1,2,4-Trioxolane Ring and Some Anti-Malaria Drugs Via Traceless Linker

There is a method to employ aberrant levels of mobile ferrous iron (Fe(II)) for selective drug delivery in vivo. This approach makes use of a 1,2,4-trioxolane moiety, which serves as an Fe(II)- sensitive “trigger,” making drug release contingent on Fe(II)-promoted trioxolane fragmentation. In this study we focused on the prototype drug conjugate of 1,2,4-trioxolane ring joined via a traceless linker to drug species, which is conjugated to the linker via a free amine or alcohol function to form a complex.We studied four complexes formed of different drugs with various linkers. The structure and physicochemical properties of complexes were studied using the Hartree – Fock method with 6-311++G** basis set. The influence of the oxalan ring and various drugs on the properties of complexes has been studied. Geometric parameter, energies, nuclear chemical shielding constants, direction of dipole moment vector, partition coefficients, electric polarizabilities, and other physicochemical properties of these compounds have been calculated.

Structural Characterization and Pharmaceutical Properties of Three Novel Cocrystals of Ethenzamide With Aliphatic Dicarboxylic Acids

Ethenzamide (ET) was screened in cocrystallization experiments with pharmaceutically acceptable coformer molecules to discover materials of improved physicochemical properties, that is, higher solubility and better stability. Three novel cocrystals of ET with glutaric, malonic, and maleic acids were obtained by neat grinding and slow evaporation from solution. The purpose of the study was to notice the changes in the geometry and interactions of ET molecule in crystalline phase introduced by different acid and relate them to physicochemical properties of pure ET. Therefore, the crystal structure of the cocrystals was determined by single crystal X-ray diffraction analysis. The powder samples were characterized by differential scanning calorimetry, Fourier-transform infrared spectroscopy, and 13C and 15N solid-state nuclear magnetic resonance spectroscopy. Spectroscopic studies were supported by gauge including projector augmented wave calculations of chemical shielding constants. The high stability of cocrystals during direct compression was proved. The solubility in simulated gastric fluids for studied cocrystals appeared to be approximately 1.6 times-fold higher than ET. The dissolution rates of all ET cocrystals were not faster than the pure drug, but after 240 min, more drugs were released.

New look at the Badger-Bauer rule: Correlations of spectroscopic IR and NMR parameters with hydrogen bond energy and geometry. FHF complexes

In this work correlation dependencies between hydrogen bond energy ΔE for complexes with FH⋯F hydrogen bond and their spectroscopic characteristics of the IR and NMR spectra are presented. We considered 26 complexes in a wide hydrogen bond energy range 0.2–47 kcal/mol. For each complex we calculated complexation energy (MP2/6-311++G(d,p)), IR spectroscopic parameters (FH stretching frequency ν, FH stretching frequency in local mode approximation νLM at MP2/6-311++G(d,p) level) and NMR parameters (chemical shift of hydrogen δH and fluorine nuclei δF, Nuclear Independent Chemical Shielding and spin-spin coupling constants 1JFH, 1hJH...F, 2hJFF at B3LYP/pcSseg-2 level). It was shown that changes of parameters upon complexation, i.e. changes of the stretching frequency in local mode approximation ΔνLM, change of the proton chemical shift ΔδH and change of the absolute value of spin-spin coupling constant 1JFH could be used for estimation of corresponding hydrogen bond strength. Furthermore, we build correlation dependencies between abovementioned spectroscopic characteristics and geometric ones, such as the asymmetry of bridging proton position q1 = 0.5·(rFH – rH…F).

Self-Consistent Field Calculation of Nuclear Magnetic Resonance Chemical Shielding Constants Using Gauge-Including Atomic Orbitals and Approximate Two-Electron Integrals.

The chain-of-spheres method (COS) for approximating two-electron integrals is applied to Hartree-Fock and density functional theory calculations of nuclear magnetic resonance chemical shielding tensors, based on gauge-including atomic orbitals. The accuracy of the approximation is compared to that of the resolution of the identity (RI) approach, using a benchmark test set of 15 small molecules. Reasonable auxiliary basis sets and grid sizes are selected on the basis of a careful investigation of how approximating each of the two-electron terms in the self-consistent field (SCF) and coupled perturbed SCF equations affects the calculated shielding constants. It is found that the errors are linearly additive but can have either sign. The mean absolute relative error due to applying the RI/COS approximations with the chosen settings to all two-electron terms is on the order of 0.01% and therefore negligible compared to the errors due to basis set incompleteness (∼1%) and the method used (10-50%). Several larger organic systems are used to assess the efficiency of the RI approximation for both Coulomb- and exchange-type integrals (RIJK) as well as a combination of RI for Coulomb and COS for exchange contributions (RIJCOSX). The RIJK approximation is more efficient for small molecules, while for systems of over 100 electrons and 1000 basis functions, the RIJCOSX approximation is superior.

Calculation of Magnetic Shielding Constants with meta-GGA Functionals Employing the Multipole-Accelerated Resolution of the Identity: Implementation and Assessment of Accuracy and Efficiency.

We present a highly efficient implementation for density functional calculations of chemical shielding constants. It employs the multipole-accelerated resolution of the identity for the calculation of the Coulomb part, which complements the usage of low order scaling routines for the evaluation of the exchange-correlation part and the Hartree-Fock exchange part. Introduced errors for shifts of chemical shielding constants of H, C, F, and P are evaluated for respective test sets of molecules and are related to the accuracy of shifts obtained with hybrid and nonhybrid functionals of the generalized gradient approximation type as well as for meta-GGA functionals themselves. Efficiency is demonstrated for α-d-glucose chains with more than 2500 atoms on a single CPU as well as with an OpenMP parallelized version.

Exploring the pnicogen bond non-covalent interactions in 4-XPhNH2:PFnH3-n complexes (n = 1–3, X = H, F, CN, CHO, NH2, CH3, NO2 and OCH3)

Quantum chemical calculations were carried out to explore the effects of substituents in para position of aniline as well as phosphine on the strength of interactions in pnicogen-bonded (PB) complexes of the 4-XPhNH2:PFnH3-n (n=1–3 and X=H, F, CN, CHO, NH2, CH3, NO2 and OCH3). The P⋯N pnicogen bond strength in these complexes was examined at M06-2X/6–311++G(d,p), MP2/6–311++G(d,p), MP2/cc-pVDZ, CCSD/cc-pVDZ and MP2/aug-cc-pVTZ levels of theory. The dispersion corrected M06-2X-GD3, B2PLYP-GD3 and mPW2PLYP-GD2 functionals were also used to calculate the binding energies. The changes in the pnicogen bond strength due to the change of the substituents were well monitored by the changes in the interaction energy, structural parameters, natural charge, charge transfer, electron density topology, NMR chemical shielding and spin-spin coupling constants. The BSSE-corrected BEs change from 14.5 to 21.4kJmol−1 for the 4-XPhNH2:PF3, 21.5 to 30.5kJmol−1 for the 4-XPhNH2:PF2H and 24.2 to 33.1kJmol−1 for the XPhNH2:PFH2 complex at MP2/aug-cc-pVTZ level and 12.5 to 18.9kJmol−1 for the 4-XPhNH2:PF3, 17.4 to 26.5kJmol−1 for the 4-XPhNH2:PF2H and 18.9 to 26.5kJmol−1 for the XPhNH2:PFH2 complexes at MP2/6–311++G(d,p) level. For each PFnH3-n monomer, results demonstrated that the strength of the pnicogen bonds increases by introducimg electron-donating substituents in 4-XPhNH2 molecule, whereas a reverse situation was found upon introducing electron-accepting substituents. Natural bond orbital analysis confirmed that the charge transfer takes place from 4-XPhNH2 to PFnH3-n. Electron density properties based on atoms in molecules theory were also utilized to characterize the pnicogen bonds.

Bifurcatriol, a New Antiprotozoal Acyclic Diterpene from the Brown Alga Bifurcaria bifurcata.

Linear diterpenes that are commonly found in brown algae are of high chemotaxonomic and ecological importance. This study reports bifurcatriol (1), a new linear diterpene featuring two stereogenic centers isolated from the Irish brown alga Bifurcaria bifurcata. The gross structure of this new natural product was elucidated based on its spectroscopic data (IR, 1D and 2D-NMR, HRMS). Its absolute configuration was identified by experimental and computational vibrational circular dichroism (VCD) spectroscopy, combined with the calculation of 13C-NMR chemical shielding constants. Bifurcatriol (1) was tested for in vitro antiprotozoal activity towards a small panel of parasites (Plasmodium falciparum, Trypanosoma brucei rhodesiense, T. cruzi, and Leishmania donovani) and cytotoxicity against mammalian primary cells. The highest activity was exerted against the malaria parasite P. falciparum (IC50 value 0.65 μg/mL) with low cytotoxicity (IC50 value 56.6 μg/mL). To our knowledge, this is the first successful application of VCD and DP4 probability analysis of the calculated 13C-NMR chemical shifts for the simultaneous assignment of the absolute configuration of multiple stereogenic centers in a long-chain acyclic natural product.

NMR shielding and spin–rotation constants in XCO (X = Ni, Pd, Pt) molecules

Ab initio nonrelativistic and four-component relativistic DFT (density functional theory) methods are combined to study the spin–rotation and absolute nuclear magnetic resonance (NMR) shielding constants of group 10 transition metal monocarbonyls. Good agreement is obtained between the calculated and available experimental data for the spin–rotation constants and shielding spans for PdCO and PtCO. These data allow us to determine accurate absolute chemical shielding constants for all the nuclei, as well as for the unknown spin–rotation constants. We compare the four-component shielding constants with those obtained from the spin–orbit zeroth-order regular approximation, together with an assessment of the performance of different basis sets. For the first time, relativistically optimised basis sets for the heavy atoms used in the four-component calculations are shown to give converged results for both magnetic properties studied.

Probing theoretical level effect on fluorine chemical shielding calculations

The extensive quantum‐chemical calculations are assessed for prediction of the 19F nuclear magnetic shielding constants. The methods within domain of density functional theory (DFT) and ab initio, Hartree–Fock (HF), and second‐order Møller–Plesset perturbation theory (MP2) are employed for calculation of the 19F chemical shielding constants for a series of 30 small fluorine‐containing molecules. The importance of inclusion of four factors, namely, electron correlation treatment, triple‐ξ valence shell, diffuse function, and polarization function on calculating the fluorine NMR chemical shieldings for a variety of Pople's standard basis sets at gauge invariant atomic orbital (GIAO) and continuous set of gauge transformation (CSGT) conditions are discussed by employing the linear regression analysis. The systematic investigation performed here on the influence of the method, level of theory, and basis set size on the accuracy of 19F chemical shielding demonstrates that the presence of a high level of electron correlation factor in level as well as the polarization function in basis set leads to an accurate wave function to compute the reliable fluorine chemical shielding. The results analysis illustrates that the polarization function has a significant effect on the basis sets to minimize the difference between theoretical and experimental values in several fluorine molecules. An additional comparison of the GIAO and CSGT results is used as a supplementary study to confirm the validity of the results. The basis set dependence of fluorine chemical shielding constants verifies that GIAO provides results that are often more accurate than those that are calculated with CSGT approach at the same basis set size, as expected. © 2013 Wiley Periodicals, Inc. Concepts Magn Reson Part A 42A: 140‐153, 2013.

Errors in the calculation of (27)Al nuclear magnetic resonance chemical shifts.

Computational chemistry is an important tool for signal assignment of 27Al nuclear magnetic resonance spectra in order to elucidate the species of aluminum(III) in aqueous solutions. The accuracy of the popular theoretical models for computing the 27Al chemical shifts was evaluated by comparing the calculated and experimental chemical shifts in more than one hundred aluminum(III) complexes. In order to differentiate the error due to the chemical shielding tensor calculation from that due to the inadequacy of the molecular geometry prediction, single-crystal X-ray diffraction determined structures were used to build the isolated molecule models for calculating the chemical shifts. The results were compared with those obtained using the calculated geometries at the B3LYP/6-31G(d) level. The isotropic chemical shielding constants computed at different levels have strong linear correlations even though the absolute values differ in tens of ppm. The root-mean-square difference between the experimental chemical shifts and the calculated values is approximately 5 ppm for the calculations based on the X-ray structures, but more than 10 ppm for the calculations based on the computed geometries. The result indicates that the popular theoretical models are adequate in calculating the chemical shifts while an accurate molecular geometry is more critical.

High-field multinuclear solid-state nuclear magnetic resonance (NMR) and first principle calculations in MgSO4 polymorphs

A combination of solid-state nuclear magnetic resonance (NMR) and first principles calculations was applied to obtain 17O, 25Mg, and 33S NMR parameters for two polymorphs of anhydrous magnesium sulfate. Working at the very high magnetic field of 21.14 T results in a dramatic improvement of resolution through a reduction of the effects of quadrupolar interactions and significant improvement in sensitivity. Experimental 25Mg and 33S spectra are dominated by quadrupolar interactions with quadrupolar parameters unique for each polymorph. In the case of 17O, there is a substantial contribution of the chemical shift anisotropy. The use of multiple-quantum magic-angle spinning (MQMAS) experiments allows the resolution of distinct oxygen species and assignment of signals in the experimental 17O spectrum. Chemical shielding constants and quadrupolar parameters for all three nuclei were calculated using plane wave pseudopotential density functional theory as implemented in the CASTEP computational package. The calculated NMR parameters are in very good agreement with the experimental results and help in signal assignment of the 17O spectrum. The results suggest applicability of such a combined computational and experimental solid-state NMR approach for the refinement of crystallographic data.

Theoretical Study of the 1:1 and 2:1 (Homo- and Heterochiral) Complexes of XOOX′ (X, X′ = H, CH3) with Lithium Cation

A theoretical study of the 1:1 and 2:1 complexes of XOOX' (X, X' = H, CH(3)) with the lithium cation has been carried out by means of ab initio computational methods up to the MP2/aug-cc-pVTZ level. The optical rotatory power and NMR parameters (absolute chemical shielding and indirect coupling constants) have been calculated. In addition, the racemization barriers within the complexes formed have been evaluated. Special attention has been paid concerning the differences between the 2:1 homo- and heterochiral complexes.

Theoretical study on the electronic, structural, properties and reactivity of a series of mono-, di-, tri- and tetrachlorothiophenes as well as corresponding radical cation forms as monomers for conducting polymers.

In this paper, electrical and structural properties of mono-, di-, tri- and tetrachlorothiophenes and their radical cations have been studied using the density functional theory and B3LYP method with 6-311++G** basis set. The effects of the number and position of the substituent of chlorine atoms on the properties of the thiophene ring for all chlorothiophenes and their radical cations have been studied. Vibrational frequencies, nuclear chemical shielding constants, spin-density distribution, size and direction of dipole moment vector, ionization potential, electric polarizabilities and NICS values of these compounds have been calculated as well. The analysis of these data showed that double bonds in 3-chlorothiophene are more delocalized and it is the best possible candidate monomer among all chlorothiophenes for the synthesis of corresponding conducting polymers with modified characteristics.

Quantitative prediction of gas-phase N15 and P31 nuclear magnetic shielding constants

High-level ab initio benchmark calculations of the (15)N and (31)P NMR chemical shielding constants for a representative set of molecules are presented. The computations have been carried out at the Hartree-Fock self-consistent field (HF-SCF), density functional theory (DFT) (B-P86 and B3-LYP), second-order Moller-Plesset perturbation theory (MP2), coupled cluster singles and doubles (CCSD), and CCSD augmented by a perturbative treatment of triple excitations [CCSD(T)] level of theory using basis sets of triple zeta quality or better. The influence of the geometry, the treatment of electron correlation, as well as basis set and zero-point vibrational effects on the shielding constants are discussed and the results are compared to gas-phase experimental shifts. As for the first time a study using high-level post-HF methods is carried out for a second-row element, we also propose a family of basis sets suitable for the computation of (31)P shielding constants. The mean deviations observed for (15)N and (31)P are 0.9 [CCSD(T)/13s9p4d3f] and -3.3 ppm [CCSD(T)/15s12p4d3f2g], respectively, when corrected for zero-point vibrational effects. Results obtained at the DFT level of theory are of comparable accuracy to MP2 for (15)N and of comparable accuracy to HF-SCF for (31)P. However, they are not improved by inclusion of zero-point vibrational effects. The PN molecule is an especially interesting case with exceptionally large electron correlation effects on shielding constants beyond MP2 which, therefore, represents an excellent example for further benchmark studies.

Solvent Dependence of (14)N Nuclear Magnetic Resonance Chemical Shielding Constants as a Test of the Accuracy of the Computed Polarization of Solute Electron Densities by the Solvent.

Although continuum solvation models have now been shown to provide good quantitative accuracy for calculating free energies of solvation, questions remain about the accuracy of the perturbed solute electron densities and properties computed from them. Here we examine those questions by applying the SM8, SM8AD, SMD, and IEF-PCM continuum solvation models in combination with the M06-L density functional to compute the (14)N magnetic resonance nuclear shieldings of CH3CN, CH3NO2, CH3NCS, and CH3ONO2 in multiple solvents, and we analyze the dependence of the chemical shifts on solvent dielectric constant. We examine the dependence of the computed chemical shifts on the definition of the molecular cavity (both united-atom models and models based on superposed individual atomic spheres) and three kinds of treatments of the electrostatics, namely the generalized Born approximation with the Coulomb field approximation, the generalized Born model with asymmetric descreening, and models based on approximate numerical solution schemes for the nonhomogeneous Poisson equation. Our most systematic analyses are based on the computation of relative (14)N chemical shifts in a series of solvents, and we compare calculated shielding constants relative to those in CCl4 for various solvation models and density functionals. While differences in the overall results are found to be reasonably small for different solvation models and functionals, the SMx models SM8, and SM8AD, using the same cavity definitions (which for these models means the same atomic radii) as those employed for the calculation of free energies of solvation, exhibit the best agreement with experiment for every functional tested. This suggests that in addition to predicting accurate free energies of solvation, the SM8 and SM8AD generalized Born models also describe the solute polarization in a manner reasonably consistent with experimental (14)N nuclear magnetic resonance spectroscopy. Models based on the nonhomogeneous Poisson equation show slightly reduced accuracy. Scaling the intrinsic Coulomb radii to larger values (as has sometimes been suggested in the past) does not uniformly improve the results for any kind of solvent model; furthermore it uniformly degrades the results for generalized Born models. Use of a basis set that increases the outlying charge diminishes the accuracy of continuum models that solve the nonhomogeneous Poisson equation, which we ascribe to the inability of the numerical schemes for approximately solving the nonhomogeneous Poisson equation to fully account for the effects of electronic charge outside the solute cavity.

Protein NMR Chemical Shift Calculations Based on the Automated Fragmentation QM/MM Approach

An automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach has been developed to routinely calculate ab initio protein NMR chemical shielding constants. The AF-QM/MM method is linear-scaling and trivially parallel. A general fragmentation scheme is employed to generate each residue-centric region which is treated by quantum mechanics, and the environmental electrostatic field is described with molecular mechanics. The AF-QM/MM method shows good agreement with standard self-consistent field (SCF) calculations of the NMR chemical shieldings for the mini-protein Trp cage. The root-mean-square errors (RMSEs) for 1H, 13C, and 15N NMR chemical shieldings are equal to or less than 0.09, 0.32, and 0.78 ppm, respectively, for all Hartree-Fock (HF) and density functional theory (DFT) calculations reported in this work. The environmental electrostatic potential is necessary to accurately reproduce the NMR chemical shieldings using the AF-QM/MM approach. The point-charge models provided by AMBER, AM1/CM2, PM3/CM1, and PM3/CM2 all effectively model the electrostatic field. The latter three point-charge models are generated via semiempirical linear-scaling SCF calculations of the entire protein system. The correlations between experimental 1H NMR chemical shifts and theoretical predictions are >0.95 for AF-QM/MM calculations using B3LYP with the 6-31G**, 6-311G**, and 6-311++G** basis sets. Our study, not unexpectedly, finds that conformational changes within a protein structure play an important role in the accurate prediction of experimental NMR chemical shifts from theory.

Improved Description of Nuclear Magnetic Resonance Chemical Shielding Constants Using the M06-L Meta-Generalized-Gradient-Approximation Density Functional

The performance of the M06-L density functional has been tested for four databases of NMR isotropic chemical shielding constants. Comparison with the B3LYP, BLYP, HCTH, KT1, KT2, LSDA, OPBE, OLYP, PBE, TPSS, and VSXC functionals shows that M06-L has improved performance for calculating NMR chemical shielding constants, especially for highly correlated systems. We also found that VSXC and M06-L have encouraging accuracy for calculating (13)C chemical shielding constants, and both functionals perform very well for the chemical shielding constants in the o-benzyne molecule.