Investigation of structure, vibrational and NMR spectra of oxycodone and naltrexone: a combined experimental and theoretical study

In this work, two important opioid antagonists, naltrexone and oxycodone, were prepared from thebaine and were characterized by IR, (1)H NMR and (13)C NMR spectroscopy. Moreover, computational NMR and IR parameters were obtained using density functional theory (DFT) at B3LYP/6-311++G** level of theory. Complete NMR and vibrational assignment were carried out using the observed and calculated spectra. The IR frequencies and NMR chemical shifts, determined experimentally, were compared with those obtained theoretically from DFT calculations, showed good agreements. The RMS errors observed between experimental and calculated data for the IR absorptions are 85 and 105 cm(-1), for the (1)H NMR peaks are 0.87 and 0.17 ppm and for those of (13)C NMR are 5.6 and 5.3 ppm, respectively for naltrexone and oxycodone.     

Correlation between 19F environment and isotropic chemical shift in barium and calcium fluoroaluminates

High-speed MAS (19)F NMR spectra are recorded and reconstructed for 10 compounds from BaF(2)-AlF(3) and CaF(2)-AlF(3) binary systems which leads to the determination of 77 isotropic (19)F chemical shifts in various environments. A first attribution of NMR lines is performed for 8 compounds using a superposition model as initially proposed by B. Bureau et al. The phenomenological parameters of this model are then refined to improve the NMR line assignment. A satisfactory reliability is reached with a root-mean-square (RMS) deviation between calculated and measured values equal to 6 ppm. The refined parameters are then successfully tested on alpha-BaCaAlF(7) whose structure was recently determined. Finally, the isotropic chemical shift ranges are defined for shared, unshared, and "free" fluorine atoms encountered in the investigated binary systems. So, the fluorine surroundings can be deduced from the NMR line positions in compounds whose structure is unknown. Such an approach can also be applied to fluoride glasses.     

Cluster models and ab initio calculations of (19)F NMR isotropic chemical shifts for inorganic fluorides

(19)F NMR isotropic chemical shift (delta(iso)) calculations are performed in crystallized compounds using the GIAO method with the B3LYP hybrid functional at DFT level. Clusters centered on the studied fluorine atoms mimic the crystalline structures. The 6-311+G(d) basis set is chosen for the central fluorine atom, and the LanL2DZ basis set for the others. The metal atoms are described by the 3-21G(2d) basis set or, when not available, by the CRENBL basis set with the corresponding ECP, and augmented with 2d polarization functions when existing. First, for high-symmetry systems (MF, MF(2), and MF(3) compounds), a systematization of the cluster building up from coordination spheres is proposed, generalized to fluoroperovskites and fluoroaluminates KAlF(4) and RbAlF(4). When applied to rather low symmetry systems such as barium fluorometalates BaMgF(4), BaZnF(4), and Ba(2)ZnF(6), the definition of the coordination spheres is far from easy. Then, for structures built up from a MF(6) octahedron network, we may define different "starting clusters": [FM(2)F(8)] for the shared fluorine atoms, [FMF(4)] for the unshared ones, and [FBa(4)](7+) for the "free" ones. Analogous "starting clusters" are then tested on compounds from the NaF-AlF(3), BaF(2)-AlF(3), and CaF(2)-AlF(3) binary systems and for alpha-BaCaAlF(7) that are also built up from a MF(6) octahedron network. For each of these corresponding fluorine sites, delta(iso) values are calculated with the "starting clusters" and several larger clusters and compared to the experimental delta(iso) values. For the barium-containing clusters, the RMS deviation is equal to 51 ppm. It is suggested that this result may be related to the poor quality of the barium basis sets for which no polarization functions are available for the moment. In total, chemical shifts were calculated for 122 fluorine sites, in a various range of compounds. For the clusters without barium, the ab initio method leads to a RMS equal to 22 ppm, which is a quite nice result keeping in mind that the (19)F chemical shift range is larger than 200 ppm.     

On the accuracy of the GIAO‐DFT calculation of 15N NMR chemical shifts of the nitrogen‐containing heterocycles – a gateway to better agreement with experiment at lower computational cost

The main factors affecting the accuracy and computational cost of the gauge‐independent atomic orbital density functional theory (GIAO‐DFT) calculation of ¹⁵N NMR chemical shifts in the representative series of key nitrogen‐containing heterocycles – azoles and azines – have been systematically analyzed. In the calculation of ¹⁵N NMR chemical shifts, the best result has been achieved with the KT3 functional used in combination with Jensen's pcS‐3 basis set (GIAO‐DFT‐KT3/pcS‐3) resulting in the value of mean absolute error as small as 5 ppm for a range exceeding 270 ppm in a benchmark series of 23 compounds with an overall number of 41 different ¹⁵N NMR chemical shifts. Another essential finding is that basically, the application of the locally dense basis set approach is justified in the calculation of ¹⁵N NMR chemical shifts within the 3–4 ppm error that results in a dramatic decrease in computational cost. Based on the present data, we recommend GIAO‐DFT‐KT3/pcS‐3//pc‐2 as one of the most effective locally dense basis set schemes for the calculation of ¹⁵N NMR chemical shifts. Copyright © 2014 John Wiley & Sons, Ltd.     

A combined semiempirical-DFT study of oligomers within the finite-chain approximation, evolution from oligomers to polymers

A computationally cheap approach combining time-independent density functional theory (TIDFT) and semiempirical methods with an appropriate extrapolation procedure is proposed to accurately estimate geometrical and electronic properties of conjugated polymers using just a small set of oligomers. The highest occupied molecular orbital-lowest unoccupied molecular orbital gap (HLG) obtained at a TIDFT level (B3PW91) for two polymers, trans-polyacetylene--the simplest conjugated polymer, and a much larger poly(2-methoxy-5-(2,9-ethyl-hexyloxy)-1,4-phenylenevinylene (MEH-PPV) polymer converge to virtually the same asymptotic value than the excitation energy obtained with time-dependent DFT (TDDFT) calculations using the same functional. For TIDFT geometries, the HLG is found to converge to a value within the experimentally accepted range for the band gap of these polymers, when an exponential extrapolation is used; however if semiempirical geometries are used, a linear fit of the HLG versus 1/n is found to produce the best results. Geometrical parameters are observed to reach a saturation value in good agreement with experimental information, within the length of oligomers calculated here and no extrapolation was considered necessary. Finally, the performance of three different semiempirical methods (AM1, PM3, and MNDO) and for the TIDFT calculations, the performance of 7 different full electron basis sets (6-311+G**, 6-31+ +G**, 6-311+ +G**, 6-31+G**, 6-31G**, 6-31+G*, and 6-31G) is compared and it is determined that the choice of semiempirical method or the basis set does not significantly affect the results.     

Ab initio and semiempirical computational studies on 1‐{(2E)‐2‐[(aminocarbonothioyl)hydrazono]‐2‐(3‐mesityl‐3‐methylcyclobutyl)ethyl}‐pyrrolidine‐2,5‐dione

The molecular geometry, vibrational frequencies, and gauge including atomic orbital (GIAO) ¹H‐ and ¹³C NMR chemical shift values of the title compound in the ground state have been calculated using the Hartree‐Fock (HF) and density functional theory (DFT) methods with 6‐31G(d) basis sets, and compared with the experimental data. The calculated results show that the optimized geometries can well reproduce the crystal structural parameters and the theoretical vibrational frequencies, and ¹H‐ and ¹³C NMR chemical shift values show good agreement with experimental data. To determine conformational flexibility, the molecular energy profile of the title compound was obtained by semiempirical (AM1) calculations with respect to the selected torsion angle, which was varied from ?180° to +180° in steps of 10°. The energetic behavior of the title compound in solvent media was examined using the B3LYP method with the 6‐31G(d) basis set by applying the Onsager and the polarizable continuum model (PCM). The results obtained with these methods reveal that the PCM method provided more stable structure than Qnsager's method. By using TD‐DFT method, electronic absorption spectra of the title compound have been predicted and a good agreement with the TD‐DFT method and the experimental one is determined. The predicted nonlinear optical properties of the title compound are much greater than ones of urea. In addition, the molecular electrostatic potential (MEP), frontier molecular orbitals (FMO) analysis, NBO analysis and thermodynamic properties of the title compound were investigated using theoretical calculations. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

(E)-4-((4-chlorobenzylidene)amino)-N-(thiazole-2yl) benzenesulfonamide: Synthesis,characterization and electronic structure theory and docking studies

Spectroscopic methods (NMR, UV–vis) were utilized to get an understanding of the molecular structure of the (E)-4-((4-chlorobenzylidene)amino)-N-(thiazole-2yl) benzenesulfonamide (4CLBTH). The DFT/B3LYP/cc-pVDZ, were used in the theoretical study of the compound. Calculations using DFT can produce geometrical parameters such as bond lengths and bond angles. The GIAO method was utilized in DMSO to calculate the ¹H and ¹³C NMR for TMS. These results were then compared to the observed data. The UV–Vis spectrum of 4CLBTH was computed by using the B3LYP with the cc-pVDZ basis set. Using theoretical calculations, an investigation into the FMO analysis and the MEP were carried out. The theoretical findings and the experimental findings are in reasonable agreement with one another. The molecular docking was investigated by the autodock suite and visualised using the discovery studio visualizer. The multiwfn software was used to calculate the wavefunction studies.     

Towards the versatile DFT and MP2 computational schemes for 31P NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the calculation of ³¹P NMR chemical shifts in the representative series of organophosphorous compounds are examined at the density functional theory (DFT) and second‐order Møller–Plesset perturbation theory (MP2) levels. At the DFT level, the best functionals for the calculation of ³¹P NMR chemical shifts are those of Keal and Tozer, KT2 and KT3. Both at the DFT and MP2 levels, the most reliable basis sets are those of Jensen, pcS‐2 or larger, and those of Pople, 6‐311G(d,p) or larger. The reliable basis sets of Dunning's family are those of at least penta‐zeta quality that precludes their practical consideration. An encouraging finding is that basically, the locally dense basis set approach resulting in a dramatic decrease in computational cost is justified in the calculation of ³¹P NMR chemical shifts within the 1–2‐ppm error. Relativistic corrections to ³¹P NMR absolute shielding constants are of major importance reaching about 20–30 ppm (ca 7%) improving (not worsening!) the agreement of calculation with experiment. Further better agreement with the experiment by 1–2 ppm can be obtained by taking into account solvent effects within the integral equation formalism polarizable continuum model solvation scheme. We recommend the GIAO‐DFT‐KT2/pcS‐3//pcS‐2 scheme with relativistic corrections and solvent effects taken into account as the most versatile computational scheme for the calculation of ³¹P NMR chemical shifts characterized by a mean absolute error of ca 9 ppm in the range of 550 ppm. Copyright © 2014 John Wiley & Sons, Ltd.     

Phosphorus mononitride: A difficult case for theory

Phosphorus nitride (PN) is the simplest molecule formed solely by phosphorus and nitrogen. It represents an interesting model for materials, where phosphorus is directly attached to nitrogen. Nevertheless, both theoretical and experimental studies often provide an incomplete picture on the structural, electronic, and spectral properties of PN. Theoretical predictions often suffer from insufficient level of theory, incomplete basis set, or from neglecting several effects, for example, zero-point vibrational correction (ZPVC). Therefore, we performed an extensive benchmark study on structural, electronic, and spectral properties of PN at the Hartree-Fock, density functional theory (DFT), or even the coupled-cluster levels. We paid special attention to the basis set effect. We tested three variants of Dunning's aug-cc-pVXZ basis sets with the size from double-ζ to sextuple-ζ, as well as Jensen's aug-pc-n, aug-pcJ-n, and aug-pcSseg-n basis sets, where n = 1-4. Obtained energetics, PN distance, dipole moment, vibrational frequencies, and nuclear magnetic resonance (NMR) parameters were extrapolated to the complete basis set limit (CBS) using three- or two-parameter formulas. The ³¹P NMR shieldings estimated with the aug-cc-pVXZ and aug-cc-pV(X + d)Z basis sets strongly depend on the basis set size providing scattered convergence patterns toward CBS. The Hartree-Fock self-consistent field (HF-SCF) NMR parameters evinced similar behavior as the coupled-cluster data. The only smooth convergence was achieved using the aug-cc-pCVXZ basis sets that include core-valence effects. The KT3 functional underestimated the phosphorus CBS shieldings by about 12 ppm compared to coupled cluster with singles and doubles (CCSD) (T). Nevertheless, KT3 unambiguously surpasses the HF-SCF and CCSD levels that provide ³¹P shieldings that are lower by about 150 ppm and 24 ppm compared to CCSD(T). The convergence of nitrogen shieldings was regular for all basis set hierarchies and all theoretical methods. Relativistic and vibrational effects on selected properties were also discussed.     

Calculating the response of NMR shielding tensor σ(31P) and 2J(31P,13C) coupling constants in nucleic acid phosphate to coordination of the Mg2+ cation

Dependence of NMR (31)P shielding tensor and (2)J(P,C) coupling constants on solvation of nucleic acid phosphate by Mg(2+) and water was studied using methods of bioinformatic structural analyses of crystallographic data and DFT B3LYP calculations of NMR parameters. The effect of solvent dynamics on NMR parameters was calculated using molecular dynamic. The NMR calculations for representative solvation patterns determined in crystals of B-DNA and A-RNA molecules pointed out the crucial importance of local Mg(2+) coordination geometry, including hydration by explicit water molecules and necessity of dynamical averaging over the solvent reorientation. The dynamically averaged (31)P chemical shift decreased by 2-9.5 ppm upon Mg(2+) coordination, the chemical shielding anisotropy increased by 0-20 ppm, and the (2)J(P,C5') coupling magnitude decreased by 0.2-1.8 Hz upon Mg(2+) coordination. The calculated decrease of the (31)P chemical shift is in excellent agreement with the 1.5-10 ppm decrease of the phosphorothioate (31)P chemical shift upon Cd(2+) coordination probed experimentally in hammerhead ribozyme (Suzumura; et al. J. Am. Chem. Soc. 2002, 124, 8230-8236; Osborne; et al., Biochemistry 2009, 48, 10654-10664). None of the dynamically averaged NMR parameters unequivocally distinguishes the site-specific Mg(2+) coordination to one of the two nonesterified phosphate oxygen atoms of the phosphate determined by bioinformatic analyses. By comparing the limit cases of static and dynamically averaged solvation, we propose that mobility of the solvent has a dramatic impact on NMR parameters of nucleic acid phosphate and must be taken into account for their accurate modeling.     

1H and 13C NMR scaling factors for the calculation of chemical shifts in commonly used solvents using density functional theory

Calculation of NMR chemical shifts and coupling constants using quantum mechanical calculations [density functional theory (DFT)], has become a very popular tool for the determination of conformation and the assignment of stereochemistry within a molecule. We present the scaling factors (linear regression parameters) from 10 DFT methods for 10 commonly used NMR solvents using the same set of reference compounds. The results were compared with the corresponding gas‐phase calculations to assess the inclusion of the polarizable continuum model for solvent effects. © 2014 Wiley Periodicals, Inc.     

Computational structural determination and energy landscape analysis of the hepatic carcinogen 2-(acetylamino)fluorene

 2-(Acetylamino)fluorene (AAF), a potent mutagen and a prototypical example of the mutagenic aromatic amines, forms covalent adducts to DNA after metabolic activation in the liver. A benchmark study of AAF is presented using a number of the most widely used molecular mechanics and semiempirical computational methods and models. The results are compared to higher-level quantum calculations and to experimentally obtained crystal structures. Hydrogen bonding between AAF molecules in the crystal phase complicates the direct comparison of gas-phase theoretical calculations with experiment, so Hartree–Fock (HF) and Becke–Perdew (BP) density functional theory (DFT) calculations are used as benchmarks for the semiempirical and molecular mechanics results. Systematic conformer searches and dihedral energy landscapes were carried out for AAF using the SYBYL and MMFF94 molecular mechanics force fields; the AM1, PM3 and MNDO semiempirical quantum mechanics methods; HF using the 3-21G*and 6-31G* basis sets; and DFT using the nonlocal BP functional and double numerical polarization basis sets. MMFF94, AM1, HF and DFT calculations all predict the same planar structures, whereas SYBYL, MNDO and PM3 all predict various nonplanar geometries. The AM1 energy landscape is in substantial agreement with HF and DFT predictions; MMFF94 is qualitatively similar to HF and DFT; and the MNDO, PM3 and SYBYL results are qualitatively different from the HF and DFT results and from each other. These results are attributed to deficiencies in MNDO, PM3 and SYBYL. The MNDO, PM3 and SYBYL models may be unreliable for compounds in which an amide group is immediately adjacent to an aromatic ring. Received: 26 May 2002 / Accepted: 12 December 2002 / Published online: 14 February 2003     

A comparative DFT study of tetracyanoquinodimethane and its difluoro and tetrafluoro analogs

The optimized molecular structures, harmonic vibrational wavenumbers and corresponding vibrational assignments, frontier molecular orbitals and UV data of 7,7,8,8-tetracyanoquinodimethane, 2,5-difluoro-7,7,8,8-tetracyanoquinodimethane and 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane were computed using quantum mechanical code. Calculations were carried out at Becke-3-Lee-Yang-Parr (B3LYP) functional with density functional theory (DFT) and time-dependent density functional theory (TD-DFT) using the 6–311++G (d,p) basis set. The theoretical results were successfully compared with some of the available experimental data. The influence of fluorine on the structural, vibrational and electronic properties were investigated. The addition of fluorine reduces both the electrical and optical band gaps. The findings of this research can be useful for analogs of the molecules studied which have potential applications in the design of organic semiconductors for electronic devise applications.     

Spectroscopic investigation (IR and NMR) and HOMO-LUMO analysis of aromatic imines using theoretical approach

In this work, the theoretical studies on the structure, FT-IR, NMR, and UV–Vis spectroscopy of (E)-N-benzylidenebenzenamine (A1) and (E)-N-(2, 4′-dichlorobenzylidene) propan-1-amine (A2) are presented. The optimized structure of the molecules, NMR and UV–Vis spectra analysis were determined by the Density Functional Theory (DFT) method using B3LYP/6-311G (d, p) basis set. For FT-IR analysis, both the HF and DFT methods were used in order to determine their accuracy and reliability in theoretical calculations. The computed result of DFT calculations in comparison with the experimental results showed that the DFT method gives a more accurate prediction. The infrared (IR) spectra for the imine molecules have been recorded in the region of 500–4000 cm^?1. The gauge-independent atomic orbital (GIAO) method has been used to evaluate the ¹³C and ¹H nuclear magnetic resonance (NMR) chemical shifts of the molecules. The computed results of NMR spectra of the molecules was found to be in good agreement with the experimental data. The UV–Vis spectra of the molecules were computed to determine the HOMO-LUMO energies in order to gain insight into their electronic properties. Mulliken population analysis on atomic charges of the molecules was also calculated using the HF (Hartree-Fock) and B3LYP method. All the computed results indicated that the B3LYP method provides satisfactory results and, therefore, can be employed to support experimental data. It also demonstrated a reliable approach towards characterization of molecules in chemical science.     

Conformational study of (8alpha,8'beta)-bis(substituted phenyl)-lignano-9,9'-lactones by means of combined computational, database mining, NMR, and chemometric approaches

Beta-(3,4-Methylenedioxybenzyl)-gamma-butyrolactone (MDBL) and (-)-hinokinin (HK) were obtained by partial synthesis and characterized by 1H NMR and computational methods (conformational analysis, molecular modeling, structural data mining and chemometrics). Three conformers were detected for MDBL and nine were found for HK. The energy differences are around 1 and 2 kcal mol(-1) and rotation barriers are less than 3 and 5 kcal mol(-1) for MDBL and HK conformers, respectively. The geometries of these conformers, obtained from semiempirical PM3 and density functional theory (DFT) B3LYP 6-31G** calculations agree satisfactorily with 1H NMR data (vicinal proton-proton coupling constants) and structures retrieved from the Cambridge Structural Database (torsion angles). DFT combined with some variants of the Haasnoot-de Leeuuw-Altona equations gives the best predictions for the coupling constants. The molecular conformation of MDBL, of HK, and of related systems depends not only on intramolecular interactions but also on crystal packing forces and solvent-solute interactions, in particular hydrogen bonds and polar interactions. Hydration favors more stable HK conformers, which can be important for their behavior in chemical and biological systems.     

51V solid-state magic angle spinning NMR spectroscopy and DFT studies of oxovanadium(V) complexes mimicking the active site of vanadium haloperoxidases

A series of 11 oxovanadium(V) complexes mimicking the active site of vanadium haloperoxidases have been investigated by (51)V magic angle spinning NMR spectroscopy and density functional theory (DFT). The MAS spectra are dominated by the anisotropic quadrupolar and chemical shielding interactions; for these compounds, C(Q) ranges from 3 to 8 MHz, and delta(sigma) is in the range 340-730 ppm. The quadrupolar coupling and chemical shielding tensors as well as their relative orientations have been determined by numerical simulations of the spectra. The spectroscopic NMR observables appear to be very sensitive to the details of the electronic and geometric environment of the vanadium center in these complexes. For the four crystallographically characterized compounds from the series, the quadrupolar and chemical shielding anisotropies were computed at the DFT level using two different basis sets, and the calculated tensors were in general agreement with the experimental solid-state NMR data. A combination of (51)V solid-state NMR and computational methods is thus beneficial for investigation of the electrostatic and geometric environment in diamagnetic vanadium systems with moderate quadrupolar anisotropies.