Role of the Solvent and Intramolecular Hydrogen Bonds in the Antioxidative Mechanism of Prenylisoflavone from Leaves of Vatairea guianensis.

This work discusses the electron structure, antioxidative properties, and solvent contribution of two new antioxidant molecules discovered, named S10 and S11, extracted from a medicinal plant called Vatairea guianensis, found in the Amazon rain-forest. To gain a better understanding, a study using density functional theory coupled with the polarizable-continuum model and the standard 6-311++G(d,p) basis set was conducted. The results indicate that S10 has a higher antioxidant potential than S11, confirming the experimental expectations. In the gas phase, the hydrogen atom transfer route dominates the hydrogen scavenging procedure. However, in the water solvents, the antioxidant mechanism prefers the sequential proton loss electron transfer mechanism. Furthermore, the solvent plays a fundamental role in the antioxidant mechanism. The formation of an intramolecular OH···OCH3 hydrogen bond is crucial for accurately describing the hydrogen scavenging phenomenon, better aligning with the experimental data. The results suggest that the two isoflavones investigated are promising for the pharmacologic and food industries.

Biotransformation of 1-nitro-2-phenylethane [Formula: see text] 2-phenylethanol from fungi species of the Amazon biome: an experimental and theoretical analysis.

CONTEXT: Natural products and their biotransformation procedures are a powerful source of new chromophores with potential applications in fields like biology, pharmacology and materials science. Thus, this work discusses about the extraction procedure of 1-nitro-2-phenylethane (1N2PE) from Aniba canelilla, its biotransformation setup into 2-phenylethanol (2PE) using four fungi, Lasiodiplodia caatinguensis (phytopathogenic fungus from Citrus sinensis), Colletotrichum sp. (phytopathogenic fungus from Euterpe oleracea), Aspergillus flavus and Rigidoporus lineatus isolated from copper mining waste located in the interior of the Brazilian Amazon. A detailed experimental and theoretical vibrational analysis (IR and Raman) have allowed us to perform some charge transfer effects on the title compounds (push-pull effect) by monitoring specific vibrational modes of their electrophilic and nucleophilic molecular sites. The solvent interactions promote molecular conformations that affect the vibrational spectra of the donor and acceptor groups, as can be seen comparatively in the gas and aqueous solution spectra, an effect possibly related to the bathochromic shift in the calculated optical spectrum of the compounds. The nonlinear optical behavior shows that while the solvent reduces the response of 1N2PE, the response of 2PE increases the optical parameters, which presents low refractive index (n) and first hyperpolarizability. ([Formula: see text]) is almost eight times that reported for urea (42.79 a.u.), a common nonlinear optical material. Furthermore, the bioconversion goes from an electrophilic to a nucleophilic compound, affecting its molecular reactivity. METHODS: 1N2PE was obtained from Aniba canelilla, whose essential oil is constituted of [Formula: see text] of 2PE. The A. canelilla essential oil was extracted under hydrodistillation. The biotransformation reactions were performed in autoclaved liquid media (100 mL) composed of malt extract (2%) in 250 mL Erlenmeyer flask. Each culture was incubated in an orbital shaker (130 rpm) at [Formula: see text]C during 7 days and after that, 50 mg of 1N2PE (80%) were diluted in 100 [Formula: see text]L of dimethylsulfoxide (DMSO) and added to the reactions flasks. Aliquots (2 mL) were removed using ethyl acetate (2 mL) and analyzed by GC-MS (fused silica capillary col1umn, Rtx -5MS 30 m [Formula: see text] 0.25 mm [Formula: see text] 0.25 [Formula: see text]m) in order to determine the amount of 1N2PE biotransformation. FTIR 1N2PE and 2PE spectra were obtained by attenuated total reflectance (ATR), using a Agilent CARY 630 spectrometer, in the spectral region 4000-650 cm[Formula: see text]. The quantum chemical calculations were carried out in the Gaussian 09 program while the DICE code was used to perform the classical Monte Carlo simulations and generate the liquid environment using the classical All-Atom Optimized parameters for Liquid Simulations (AA-OPLS). All nonlinear optical properties, reactive parameters, and electronic excitations were calculated using the Density Functional Theory framework coupled to the standard 6-311++G(d,p) basis set.

The Solute Polarization and Structural Effects on the Nonlinear Optical Response of Based Chromone Molecules.

The solute polarization due to solvent is a an electrostatic quantum effect that impacts diverse molecular properties, including the nonlinear optical response of a material. An iterative procedure that allows updating the solute charge distribution in the presence of the solvent is combined with a sequential Monte Carlo/Quantum Mechanics methodology and Density Functional Theory methods to evaluate the nonlinear optical (NLO) response using the hyper Rayleigh scattering (HRS) of a series of chromones recently identified in Chamaecrista diphylla, an herbaceous plant abundant throughout the Americas and used in folk medicine. From this study, it is determined that from gas to solvent environment, the systems acquire low refractive index (n) and an improvement of the first hyperpolarizability (ßHRS ), signaling potential NLO uses. It is shown that the octupolar contributions (ßJ=3 ) superate the dipolar ones (ßJ=1 ) and dominate the second-order optical response in both gas and liquid phases, which indicate nontrivial optical materials. Moreover, the solvent environment and structural changes in the periphery can tune significantly the dipolar/octupolar balance, showing a key to control the decoupling between these contributions.

The Importance of the Density Functional Theory Exchange-Correlation Hartree-Fock Term in Magnetic Resonance: Application to an Aqueous Environment.

Within the framework of Density Functional Theory (DFT), the relevance of the term Hartree-Fock exchange (HFE) for a variety of molecular properties is a critical point. For this reason, we spend efforts to understand these relationships in nuclear magnetic resonance (NMR) parameters in a water solvent. This work takes advantage of the appropriate aug-cc-pVTZ-J basis set and the Minnesota family of DFT methods, which consider different portions of HFE contributions. With regard to solvent participation, the results are based on a sequential Monte Carlo/Quantum Mechanics procedure, which builds the structures of the liquid under realistic thermodynamic conditions. Compared to the accurate results of second-order polarization propagator approximation (SOPPA) and experimental data, all NMR parameters show a huge dependence on the size of the HFE contribution. For instance, the inclusion of this term in 1JOH and 2JHH indirect spin-spin couplings does vary with 49.661 and 25.459 Hz, respectively. The M06-HF method accounts for 100% of HFE and better matches the σO and σH shielding constants. On the other hand, 1JOH and 2JHH demand a medium contribution (54% of HFE), the best description being associated with the M06-2X method. Thus, the dependence varies regarding the phenomenology of the property in focus and the order for independent treatments. For elements that participate in hydrogen bonds simultaneously as donor and acceptor actors, the results indicate that explicit solvent molecules must be considered in the quantum mechanical calculations for better modeling of paramagnetic shielding constants.

Physical achirality in geometrically chiral rotamers of hydrazine and boranylborane molecules.

The diagonal components and the trace of tensors which account for chiroptical response of the hydrazine molecule N2 H4 , that is, static anapole magnetizability and frequency-dependent electric dipole-magnetic dipole polarisability, are a function of the ϕ ≡ ∠ H─N─N─H dihedral angle. They vanish for symmetry reasons at ϕ = 0° and ϕ = 180°, corresponding respectively to C2v and C2h point group symmetries, that is, cis and trans conformers characterized by the presence of molecular symmetry planes. Nonetheless, vanishing diagonal components have been observed also in the proximity of ∠ H─N─N─H = 90°, in which the point group symmetry is C2 and hydrazine is unquestionably chiral. In the boranylborane molecule B2 H4 , assuming the B─B bond in the y direction, the ayy component of the anapole magnetizability tensor approximately vanishes for dihedral angles ∠ H─B─B─H corresponding to chiral rotamers which belong to D2 symmetry. Such anomalous effects have been ascribed to physical achirality of these conformers, that is, to their inability to sustain electronic current densities inducing either anapole moments, or electric and magnetic dipole moments, about the chiral axis connecting heavier atoms, as well as perpendicular directions. In other terms, the structure of certain geometrically chiral rotamers may be such that neither toroidal nor helical flow, which determine chiroptical phenomenology, can take place in the presence of perturbing fields parallel or orthogonal to the chiral axis.

A Density Functional Theory Study of Optical Rotation in Some Aziridine and Oxirane Derivatives.

We present time-dependent density functional theory (TDDFT) calculations of the electronic optical rotation (ORP) for seven oxirane and two aziridine derivatives in the gas phase and in solution and compare the results with the available experimental values. For seven of the studied molecules it is the first time that their optical rotation was studied theoretically and we have therefore investigated the influence of several settings in the TDDFT calculations on the results. This includes the choice of the one-electron basis set, the exchange-correlation functional or the particular polarizable continuum model (PCM). We can confirm that polarized quadruple zeta basis sets augmented with diffuse functions are necessary for converged results and find that the aug-pc-3 basis set is a viable alternative to the frequently employed aug-cc-pVQZ basis set. Based on our study, we cannot recommend the generalized gradient functional KT3 for calculations of the ORP in these compounds, whereas the hybrid functional PBE0 gives results quite similar to the long-range correct CAM-B3LYP functional. Finally, we observe large differences in the solvent effects predicted by the integral equation formalism of PCM and the SMD variant of PCM. For the majority of solute/solvent combinations in this study, we find that the SMD model in combination with the PBE0 functional and the aug-pc-3 basis set gives the best agreement with the experimental values.

Noniterative Doubles Corrections to the Random Phase and Higher Random Phase Approximations: Singlet and Triplet Excitation Energies.

The second-order noniterative doubles-corrected random phase approximation (RPA) method has been extended to triplet excitation energies and the doubles-corrected higher RPA method as well as a shifted version for calculating singlet and triplet excitation energies are presented here for the first time. A benchmark set consisting of 20 molecules with a total of 117 singlet and 71 triplet excited states has been used to test the performance of the new methods by comparison with previous results obtained with the second-order polarization propagator approximation (SOPPA) and the third order approximate coupled cluster singles, doubles and triples model CC3. In general, the second-order doubles corrections to RPA and HRPA significantly reduce both the mean deviation as well as the standard deviation of the errors compared to the CC3 results. The accuracy of the new methods approaches the accuracy of the SOPPA method while using only 10-60% of the calculation time. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

RPA(D) and HRPA(D): Two new models for calculations of NMR indirect nuclear spin-spin coupling constants.

In this article, the RPA(D) and HRPA(D) models for the calculation of linear response functions are presented. The performance of the new RPA(D) and HRPA(D) models is compared to the performance of the established RPA, HRPA, and SOPPA models in calculations of indirect nuclear spin-spin coupling constants using the CCSD model as a reference. The doubles correction offers a significant improvement on both the RPA and HRPA models; however, the improvement is more dramatic in the case of the RPA model. For all coupling types investigated in this study, the results obtained using the HRPA(D) model are comparable in accuracy to those given by the SOPPA model, while requiring between 30% and 90% of the calculation time needed for SOPPA. The RPA(D) model, while of slightly lower accuracy compared to the CCSD model than HRPA(D), offered calculation times of only approximately 25% of those required for SOPPA for all the investigated molecules. © 2018 Wiley Periodicals, Inc.

Computational study of basis set and electron correlation effects on anapole magnetizabilities of chiral molecules.

In the presence of a static, nonhomogeneous magnetic field, represented by the axial vector B at the origin of the coordinate system and by the polar vector C=∇×B, assumed to be spatially uniform, the chiral molecules investigated in this paper carry an orbital electronic anapole, described by the polar vector A. The electronic interaction energy of these molecules in nonordered media is a cross term, coupling B and C via a¯, one third of the trace of the anapole magnetizability aαß tensor, that is, WBC=-a¯B·C. Both A and W(BC) have opposite sign in the two enantiomeric forms, a fact quite remarkable from the conceptual point of view. The magnitude of a¯ predicted in the present computational investigation for five chiral molecules is very small and significantly biased by electron correlation contributions, estimated at the density functional level via three different functionals. © 2016 Wiley Periodicals, Inc.

On the truncation of the number of excited states in density functional theory sum-over-states calculations of indirect spin spin coupling constants.

It is investigated, whether the number of excited (pseudo)states can be truncated in the sum-over-states expression for indirect spin-spin coupling constants (SSCCs), which is used in the Contributions from Localized Orbitals within the Polarization Propagator Approach and Inner Projections of the Polarization Propagator (IPPP-CLOPPA) approach to analyzing SSCCs in terms of localized orbitals. As a test set we have studied the nine simple compounds, CH4, NH3, H2O, SiH4, PH3, SH2, C2H2, C2H4, and C2H6. The excited (pseudo)states were obtained from time-dependent density functional theory (TD-DFT) calculations with the B3LYP exchange-correlation functional and the specialized core-property basis set, aug-cc-pVTZ-J. We investigated both how the calculated coupling constants depend on the number of (pseudo)states included in the summation and whether the summation can be truncated in a systematic way at a smaller number of states and extrapolated to the total number of (pseudo)states for the given one-electron basis set. We find that this is possible and that for some of the couplings it is sufficient to include only about 30% of the excited (pseudo)states.

Investigation of the resonance-assisted hydrogen bond in model ß-diketones through localized molecular orbital analysis of the spin-spin coupling constants related to the O-H···O hydrogen bond.

The resonance-assisted hydrogen bond (HB) phenomenon has been studied theoretically by a localized molecular orbital (LMO) decomposition of the spin-spin coupling constants between atoms either involved or close to the O-H···O system of some ß-diketones and their saturated counterparts. The analysis, carried out at the level of the second-order polarization propagator approximation, shows that the contributions in terms of LMO to the paramagnetic spin orbital and the spin dipolar Ramsey terms proof the importance of the delocalized π-electron structure supporting the idea of the existence of the resonance-assisted HB phenomenon phenomenon. The LMO contributions to the Fermi contact term indicate mainly the presence of the HB that may or not be linked to the π-electrons.

Communication: Localized molecular orbital analysis of the effect of electron correlation on the anomalous isotope effect in the NMR spin-spin coupling constant in methane.

We discuss the effect of electron correlation on the unexpected differential sensitivity (UDS) in the (1)J(C-H) coupling constant of CH4 using a decomposition into contributions from localized molecular orbitals and compare with the (1)J(N-H) coupling constant in NH3. In particular, we discuss the well known fact that uncorrelated coupled Hartree-Fock (CHF) calculations are not able to reproduce the UDS in methane. For this purpose we have implemented for the first time a localized molecular orbital analysis for the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes--SOPPA(CCSD) in the DALTON program. Comparing the changes in the localized orbital contributions at the correlated SOPPA and SOPPA(CCSD) levels and at the uncorrelated CHF level, we find that the latter overestimates the effect of stretching the bond between the coupled atoms on the contribution to the coupling from the localized bonding orbital between these atoms. This disturbs the subtle balance between the molecular orbital contributions, which lead to the UDS in methane.

Monte Carlo-quantum mechanics study of magnetic properties of hydrogen peroxide in liquid water.

A theoretical study of magnetic properties of hydrogen peroxide in water has been carried out by means of Monte Carlo simulation and quantum mechanics calculations. The solvent effects were evaluated in supermolecular structures generated by simulations in the NPT ensemble. The solute-solvent structure was analyzed in terms of radial distribution functions, and the solute-solvent hydrogen bonds were identified with geometric and energetic criteria. Approximately three water molecules are hydrogen bonded to H2O2 (0.6 and 0.8 in each hydrogen and oxygen atom, respectively, of the H2O2). Although, on average, both hydroxyls of the peroxide are equivalent, the distribution of hydrogen-bonded water molecules is highly asymmetric. Analyzing the statistics of the hydrogen bonds, we identify that only 34% of the configurations give symmetric distributions around the two hydroxyls of the H2O2 simultaneously. The magnetic shieldings and the indirect spin-spin coupling constants were calculated at the B3LYP/aug-cc-pVTZ and aug-cc-pVTZ-J computational level. We find that the solvent shields the oxygen and unshields the hydrogen atoms of the peroxide (+5.5 and -2.9 ppm, respectively), with large fluctuation from configuration to configuration in the oxygen case, an effect largely accounted for in terms of a single hydrogen bond with H2O2 as the proton donor. The most sensitive coupling in the presence of the solvent is observed to be the one-bond J(O,H).

Electric field effects on nuclear magnetic shielding of the 1:1 and 2:1 (homo and heterochiral) complexes of XOOX' (X, X' = H, CH3) with lithium cation and their chiral discrimination.

The set of 1:1 and 2:1 complexes of XOOX' (X, X' = H, CH(3)) with lithium cation has been studied to determine if they are suitable candidates for chiral discrimination in an isotropic medium via nuclear magnetic resonance spectroscopy. Conventional nuclear magnetic resonance is unable to distinguish between enantiomers in the absence of a chiral solvent. The criterion for experimental detection is valuated by the isotropic part of nuclear shielding polarisability tensors, related to a pseudoscalar of opposite sign for two enantiomers. The study includes calculations at coupled Hartree-Fock and density functional theory schemes for (17)O nucleus in each compound. Additional calculations for (1)H are also included for some compounds. A huge static homogeneous electric field, perpendicular to the magnetic field of the spectromer, as big as ≈1.7 × 10(8) V m(-1) should be applied to observe a shift of ≈1 ppm for (17)O magnetic shielding in the proposed set of complexes.

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.

Optimized basis sets for the calculation of indirect nuclear spin-spin coupling constants involving the atoms B, Al, Si, P, and Cl.

The aug-cc-pVTZ-J series of basis sets for indirect nuclear spin-spin coupling constants has been extended to the atoms B, Al, Si, P, and Cl. The basis sets were obtained according to the scheme previously described by Provasi et al. [J. Chem. Phys. 115, 1324 (2001)]. First, the completely uncontracted correlation consistent aug-cc-pVTZ basis sets were extended with four tight s and three tight d functions. Second, the s and p basis functions were contracted with the molecular orbital coefficients of self-consistent-field calculations performed with the uncontracted basis sets on the simplest hydrides of each atom. As a first illustration, we have calculated the one-bond indirect spin-spin coupling constants in BH(4)(-), BF, AlH, AlF, SiH(4), SiF(4), PH(3), PF(3), H(2)S, SF(6), HCl, and ClF at the level of density functional theory using the Becke three parameter Lee-Yang-Parr and the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes.

NMR chemical shielding and spin-spin coupling constants of liquid NH3: a systematic investigation using the sequential QM/MM method.

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the (1)J(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, sigma((15)N) calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Analysis of isotope effects in NMR one-bond indirect nuclear spin-spin coupling constants in terms of localized molecular orbitals.

We recently showed, by analyzing contributions from localized molecular orbitals, that the anomalous deuterium isotope effect in the one-bond indirect nuclear spin-spin coupling constant of methane, also called the unexpected differential sensitivity, can be explained by the transfer of s-orbital character from the stretched bond to the other unchanged bonds [ChemPhysChem, 2008, 9, 1259]. We now extend this analysis of isotope effects to the molecules BH(4)(-), NH(4)(+), SiH(4), H(2)O and NH(3) in order to test our conclusions on a wider rage of XH(4) compounds and to investigate whether the lone-pair orbitals are really responsible for the absence of a similar effect in water and ammonia as proposed earlier [J. Chem. Phys., 2000, 113, 3121].

Stereochemical dependence of NMR geminal spin-spin coupling constants.

In this work it was sought to explore the versatility of geminal spin-spin coupling constants, (2)J(XY) SSCCs, as probes for stereochemical studies. A set of compounds, where their experimental (2)J(XY) SSCCs through the X-C-Y molecular fragment are predicted to be sensitive to hyperconjugative interactions involving either bonding or antibonding orbitals containing the C carbon atom ('coupling pathway'), were analyzed. SSCC calculations were performed for some selected examples using the second order polarization propagator approximation (SOPPA) method or within the DFT-B3LYP framework. Hyperconjugative interactions were calculated within the Natural Bond Orbital (NBO) approach. Results are condensed in two qualitative rules: Rule I(M)-hyperconjugative interactions transferring charge into the coupling pathway yield a positive increase to the Fermi contact (FC), contribution to (2)K(XY) reduced spin-spin coupling constants (RSSCC), and Rule II(M)-hyperconjugative interactions transferring charge from the coupling pathway yield a negative increase to the FC contribution to (2)K(XY) RSSCC.

Resolving an apparent discrepancy between theory and experiment: spin-spin coupling constants for FCCF.

Ab initio equation of motion coupled cluster singles and doubles (EOM-CCSD) and second-order polarization propagator approximation (SOPPA) calculations have been performed to evaluate spin-spin coupling constants for FCCF (difluoroethyne). The computed EOM-CCSD value of (3)J(F-F) obtained at the experimental geometry of this molecule supports the previously reported experimental value of 2.1 Hz, thereby resolving an apparent discrepancy between theory and experiment. This coupling constant exhibits a strong dependence on the C-C and C-F distances, and its small positive value results from a sensitive balance of paramagnetic spin-orbit (PSO) and spin-dipole (SD) terms. The three other unique FCCF coupling constants (1)J(C-C), (1)J(C-F), and (2)J(C-F) have also been reported and compared with experimental data. While (1)J(C-F) is in agreement with experiment, the computed value of (2)J(C-F) is larger than our estimate of the experimental coupling constant.