Uracil-water interaction revisited - in search of single H-bonded secondary minima.

Monohydrated uracil (UW) complexes are stabilized by both O⋯HO and NH⋯O hydrogen bonds (H-bonds), simultaneously participating in forming three stable cyclic structures. The role and contribution of these individual H-bonds (O⋯HO and NH⋯O) to the stability of the three UW complexes are still not understood, because of the technical problems in obtaining their optimized structures by standard geometry optimization. The present study explores a non-standard approach to identify three single H-bonded local minima structures without imaginary frequency using DFT (M06-2X, B3LYP and B3LYP-D3), MP2 and CCSD(T) theories and Dunning's correlation-consistent aug-cc-pVTZ basis set, in both vacuum and aqueous media (CPCM method). The results reveal that these new structures are very shallow secondary minima between two deep wells or next to a deep well of primary minima (double H-bonded structures) in the potential energy surface. The H-bond energy of these single H-bonded complexes is found to be less sensitive to a wide range (about 15-20 degrees) of O⋯HO and NH⋯O angles, and the linearity is preferred in the stable three single H-bonded structures. The technical method used to locate such a shallow minimum is described in detail and may be useful for identifying local minima in other cases where consecutive multiple H-bonded structures are global minima. Energy decomposition (using symmetry adapted perturbation theory (SAPT)) of interaction energy, electron redistribution, and relevant vibrational modes are discussed.

Electron Correlation or Basis Set Quality: How to Obtain Converged and Accurate NMR Shieldings for the Third-Row Elements?

The quality of theoretical NMR shieldings calculated at the quantum-chemical level depends on various theoretical aspects, of which the basis set type and size are among the most important factors. Nevertheless, not much information is available on the basis set effect on theoretical shieldings of the NMR-active nuclei of the third row. Here, we report on the importance of proper basis set selection to obtain accurate and reliable NMR shielding parameters for nuclei from the third row of the periodic table. All calculations were performed on a set of eleven compounds containing the elements Na, Mg, Al, Si, P, S, or Cl. NMR shielding tensors were calculated using the SCF-HF, DFT-B3LYP, and CCSD(T) methods, combined with the Dunning valence aug-cc-pVXZ, core-valence aug-cc-pCVXZ, Jensen polarized-convergent aug-pcSseg-n and Karlsruhe x2c-Def2 basis set families. We also estimated the complete basis set limit (CBS) values of the NMR parameters. Widely scattered nuclear shieldings were observed for the Dunning polarized-valence basis set, which provides irregular convergence. We show that the use of Dunning core-valence or Jensen basis sets effectively reduces the scatter of theoretical NMR results and leads to their exponential-like convergence to CBS. We also assessed the effect of vibrational, temperature, and relativistic corrections on the predicted shieldings. For systems with single bonds, all corrections are relatively small, amounting to less than 4% of the CCSD(T)/CBS value. Vibrational and temperature corrections were less reliable for H3PO and HSiCH due to the high anharmonicity of the molecules. An abnormally high relativistic correction was observed for phosphorus in PN, reaching ~20% of the CCSD(T)/CBS value, while the correction was less than 7% for other tested molecules.

On the impact of side methyl groups on the structure and vibrational properties of ß-carotenoids. The case of butadiene and isoprene.

Theoretical consideration about the impact of methyl groups on the structure and vibrational properties of ß-carotenoids, using medium size molecules of trans-butadiene and trans-isoprene, are reported. Density functional theory (DFT) calculations with correlation-consistent and polarization-consistent basis sets were applied to trans-1,3-butadiene and trans-isoprene as the smallest building bricks of ß-carotenoids. Their structure and harmonic vibrations were estimated in the complete basis set limit (CBS) using the non-linear least square fit. Optimized geometries and harmonic frequencies, obtained with B3LYP and BLYP density functionals and large basis sets, were favorably reproduced by a significantly faster approach, using a recently modified STO(1M)-3G Slater-type basis set. Selected density functionals with STO(1M)-3G and 6-311++G** basis sets were also successful in predicting ß-carotene structures and harmonic vibrations. This work demonstrates the potential applicability of the proposed level of theory for larger molecules, including ß-carotenoids, present in numerous natural food sources. The proposed scheme of molecular modeling, applied to biologically active compounds in food, could provide a deeper insight into their function in vivo, which is directly related to their structure and spectroscopic properties. It could also support the experimental qualitative analysis, based on peak assignment of ß-carotenoids in various food sources.

On Complex Formation between 5-Fluorouracil and ß-Cyclodextrin in Solution and in the Solid State: IR Markers and Detection of Short-Lived Complexes by Diffusion NMR.

In this work, the nuclear magnetic resonance (NMR) and IR spectroscopic markers of the complexation between 5-fluorouracil (5-FU) and ß-cyclodextrin (ß-CD) in solid state and in aqueous solution are investigated. In the attenuated total reflectance(ATR) spectra of 5-FU/ß-CD products obtained by physical mixing, kneading and co-precipitation, we have identified the two most promising marker bands that could be used to detect complex formations: the C=O and C-F stretching bands of 5-FU that experience a blue shift by ca. 8 and 2 cm-1 upon complexation. The aqueous solutions were studied by NMR spectroscopy. As routine NMR spectra did not show any signs of complexation, we have analyzed the diffusion attenuation of spin-echo signals and the dependence of the population factor of slowly diffusing components on the diffusion time (diffusion NMR of pulsed-field gradient (PFG) NMR). The analysis has revealed that, at each moment, ~60% of 5-FU molecules form a complex with ß-CD and its lifetime is ca. 13.5 ms. It is likely to be an inclusion complex, judging from the independence of the diffusion coefficient of ß-CD on complexation. The obtained results could be important for future attempts of finding better methods of targeted anticancer drug delivery.

Factors Governing the Chemical Stability and NMR Parameters of Uracil Tautomers and Its 5-Halogen Derivatives.

We report on the density functional theory (DFT) modelling of structural, energetic and NMR parameters of uracil and its derivatives (5-halogenouracil (5XU), X = F, Cl, Br and I) in vacuum and in water using the polarizable continuum model (PCM) and the solvent model density (SMD) approach. On the basis of the obtained results, we conclude that the intramolecular electrostatic interactions are the main factors governing the stability of the six tautomeric forms of uracil and 5XU. Two indices of aromaticity, the harmonic oscillator model of aromaticity (HOMA), satisfying the geometric criterion, and the nuclear independent chemical shift (NICS), were applied to evaluate the aromaticity of uracil and its derivatives in the gas phase and water. The values of these parameters showed that the most stable tautomer is the least aromatic. A good performance of newly designed xOPBE density functional in combination with both large aug-cc-pVQZ and small STO(1M)-3G basis sets for predicting chemical shifts of uracil and 5-fluorouracil in vacuum and water was observed. As a practical alternative for calculating the chemical shifts of challenging heterocyclic compounds, we also propose B3LYP calculations with small STO(1M)-3G basis set. The indirect spin-spin coupling constants predicted by B3LYP/aug-cc-pVQZ(mixed) method reproduce the experimental data for uracil and 5-fluorouracil well.

Performance of polarization-consistent vs. correlation-consistent basis sets for CCSD(T) prediction of water dimer interaction energy.

Detailed study of Jensen's polarization-consistent vs. Dunning's correlation-consistent basis set families performance on the extrapolation of raw and counterpoise-corrected interaction energies of water dimer using coupled cluster with single, double, and perturbative correction for connected triple excitations (CCSD(T)) in the complete basis set (CBS) limit are reported. Both 3-parameter exponential and 2-parameter inverse-power fits vs. the cardinal number of basis set, as well as the number of basis functions were analyzed and compared with one of the most extensive CCSD(T) results reported recently. The obtained results for both Jensen- and Dunning-type basis sets underestimate raw interaction energy by less than 0.136 kcal/mol with respect to the reference value of - 4.98065 kcal/mol. The use of counterpoise correction further improves (closer to the reference value) interaction energy. Asymptotic convergence of 3-parameter fitted interaction energy with respect to both cardinal number of basis set and the number of basis functions are closer to the reference value at the CBS limit than other fitting approaches considered here. Separate fits of Hartree-Fock and correlation interaction energy with 3-parameter formula additionally improved the results, and the smallest CBS deviation from the reference value is about 0.001 kcal/mol (underestimated) for CCSD(T)/aug-cc-pVXZ calculations. However, Jensen's basis set underestimates such value to 0.012 kcal/mol. No improvement was observed for using the number of basis functions instead of cardinal number for fitting. Graphical Abstract.

Local aromaticity mapping in the vicinity of planar and nonplanar molecules.

We report on nucleus-independent magnetic shielding (NICS) scans over the centers of six- and five-membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H2 Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, 3 He, and 7 Li+ , approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and 3 He chemical shifts in the middle of neutral C60 and C606- anion agreed well with the reported experimental NMR values for He@C60 and He@C606- .

Convergence of nuclear magnetic shieldings and one-bond 1 J(11 B1 H) indirect spin-spin coupling constants in small boron molecules.

Self-consistent field Hartree-Fock, density functional theory, and coupled-cluster calculations of the nuclear magnetic shielding constants of BH and BH3 molecules have been conducted to characterize the convergence of individual results obtained with correlation-consistent and polarization-consistent basis sets. The individual 11 B and 1 H NMR parameters were estimated in the complete basis set limit and compared with benchmark results. Only the KT3 density functional accurately reproduced 11 B shielding in BH molecule.

DFT studies on the structural and vibrational properties of polyenes.

Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C = C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. BLYP/6-311++G** calculated ν(C = C) and ν(C-C) frequencies for all-trans and all-cis polyenes containing up to 14 carbon-carbon double bonds were used to estimate these values for very long chains. Correction of raw, unscaled vibrational data was performed by comparing theoretical and experimental wavenumbers for polyenes chains containing 3 to 12 conjugated C = C units with both ends substituted by tert-butyl groups. The corrected ν(C = C) and ν(C-C) wavenumbers for all-trans molecules were used to estimate the presence of 9 - 12 C = C units in all-trans polyene pigment in red coral.

Predicting the structure and vibrational frequencies of ethylene using harmonic and anharmonic approaches at the Kohn-Sham complete basis set limit.

In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn-Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested.The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies. Graphical Abstract In this study, the performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, and 6) was tested. For the first time, we demonstrated regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn-Sham complete basis set (KS CBS) limit.

Toward engineering efficient peptidomimetics. Screening conformational landscape of two modified dehydroaminoacids.

Effective peptidomimetics should posses structural rigidity and appropriate interaction pattern leading to potential spatial and electronic matching to the target receptor site. Rational design of such small bioactive molecules could push chemical synthesis and molecular modeling toward faster progress in medicinal chemistry. Conformational properties of N-t-butoxycarbonyl-glycine-(E/Z)-dehydrophenylalanine N',N'-dimethylamides (Boc-Gly-(E/Z)-ΔPhe-NMe2 ) in chloroform were studied by NMR and IR spectroscopy. The experimental findings were supported by extensive calculations at DFT(B3LYP, M06-2X) and MP2 levels of theory and the ß-turn tendency for both isomers of the studied dipeptide were determined in vacuum and in solution. The theoretical data and experimental IR results were used as an additional information for the NMR-based determination of the detailed solution conformations of the peptides. The obtained results reveal that N-methylation of C-terminal amide group changes dramatically the conformational properties of studied dehydropeptides. Theoretical conformational analysis reveals that the tendency to adopt ß-turn conformations is much weaker for the N-methylated Z isomer (Boc-Gly-(Z)-ΔPhe-NMe2 ), both in vacuum and in polar environment. On the contrary, N-methylated E isomer (Boc-Gly-(E)-ΔPhe-NMe2 ) can easily adopt ß-turn conformation, but the backbone torsion angles (φ1, ψ1, φ2, ψ2) are off the limits for common ß-turn types.

ß-turn tendency in N-methylated peptides with dehydrophenylalanine residue: DFT study.

The tendency to adopt ß-turn conformation by model dipeptides with α,ß-dehydrophenylalanine (ΔPhe) residue in the gas phase and in solution is investigated by theoretical methods. We pay special attention to a dependence of conformational properties on the side-chain configuration of dehydro residue and the influence of N-methylation on ß-turn stability. An extensive computational study of the conformational preferences of Z and E isomers of dipeptides Ac-Gly-(E/Z)-ΔPhe-NHMe (1a / 1b) and Ac-Gly-(E/Z)-ΔPhe-NMe(2) (2a/2b) by B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) methods is reported. It is shown that, in agreement with experimental data, Ac-Gly-(Z)-ΔPhe-NHMe has a great tendency to adopt ß-turn conformation. In the gas phase the type II ß-turn is preferred, whereas in the polar environment, the type I. On the other hand, dehydro residue in Ac-Gly-(E)-ΔPhe-NHMe has a preference to adopt extended conformations in all environments. N-methylation of C-terminal amide group, which prevents the formation of 1←4 intramolecular hydrogen bond, change dramatically the conformational properties of studied dehydropeptides. Especially, the tendency to adopt ß-turn conformations is much weaker for the N-methylated Z isomer (Ac-Gly-(Z)-ΔPhe-NMe(2) ), both in vacuo and in the polar environment. On the contrary, N-methylated E isomer (Ac-Gly-(E)-ΔPhe-NMe(2) ) can easier adopt ß-turn conformation, but the backbone torsion angles (ϕ(1) , ψ(1) , ϕ(2) , ψ(2) ) are off the limits for common ß-turn types.

Estimating the carbonyl anharmonic vibrational frequency from affordable harmonic frequency calculations.

A linear correlation between harmonic and anharmonic frequencies of water calculated at B3LYP level of theory was observed with a number of basis sets. Similar relationships were found in both the gas phase and solution for several small molecules. The best correlation was found for C = O stretch mode in formaldehyde, formamide and N-methylacetamide. The average difference between B3LYP harmonic and anharmonic ν(C = O) frequencies calculated with several basis sets in these molecules was 30 cm(-1). The ad hoc correction of -30 cm(-1), added to harmonic frequencies of two different carbonyl groups present in a structure of a larger molecule was tested as a fast way of predicting anharmonic frequencies without elaborated calculations. The proposed approach was tested successfully on a larger molecule of E and Z isomers of N-acetyl-α,ß-dehydrophenylalanine N',N'-dimethylamide [Ac-(E/Z)-ΔPhe-NMe(2)] and the estimated anharmonic ν(C = O) frequencies were close to directly calculated results.

Experimental and theoretical NMR and IR studies of the side-chain orientation effects on the backbone conformation of dehydrophenylalanine residue.

Conformation of N-acetyl-(E)-dehydrophenylalanine N', N'-dimethylamide (Ac-(E)-ΔPhe-NMe(2)) in solution, a member of (E)-α, ß-dehydroamino acids, was studied by NMR and infrared spectroscopy and the results were compared with those obtained for (Z) isomer. To support the spectroscopic interpretation, the Φ, Ψ potential energy surfaces were calculated at the MP2/6-31 + G(d,p) level of theory in chloroform solution modeled by the self-consistent reaction field-polarizable continuum model method. All minima were fully optimized by the MP2 method and their relative stabilities were analyzed in terms of π-conjugation, internal H-bonds and dipole interactions between carbonyl groups. The obtained NMR spectral features were compared with theoretical nuclear magnetic shieldings, calculated using Gauge Independent Atomic Orbitals (GIAO) approach and rescaled to theoretical chemical shifts using benzene as reference. The calculated indirect nuclear spin-spin coupling constants were compared with available experimental parameters.

Effects of side-chain orientation on the backbone conformation of the dehydrophenylalanine residue. Theoretical and X-ray study.

Two E isomers of α,ß-dehydro-phenylalanine, Ac-(E)-ΔPhe-NHMe (1a) and Ac-(E)-ΔPhe-NMe(2) (2a), have been synthesized and their low temperature structures determined by single-crystal X-ray diffraction. A systematic theoretical analysis was performed on these molecules and their Z isomers (1b and 2b). The ϕ,ψ potential energy surfaces were calculated at the MP2/6-31+G(d,p) and B3LYP/6-31+G(d,p) levels in the gas phase and at the B3LYP/6-31+G(d,p) level in the chloroform and water solutions with the SCRF-PCM method. All minima were fully optimized by the MP2 and DFT methods, and their relative stabilities were analyzed in terms of π-conjugation, internal H-bonds, and dipole-dipole interactions between carbonyl groups. The results indicate that all the studied compounds can adopt the conformation H (ϕ, ψ ≈ ±40°, ∓120°) which is atypical for standard amino acids residues. A different arrangement of the side chain in the E and Z isomers causes them to have different conformational preferences. In the presence of a polar solvent both Z isomers of ΔPhe (1b and 2b) are found to adopt the 3(10)-helical conformation (left- and right-handed are equally likely). On the other hand, this conformation is not accessible or highly energetic for E isomers of ΔPhe (1a and 2a). Those isomers have an intrinsic inclination to have an extended conformation. The conformational space of the Z isomers is much more restricted than that of the E derivative both in the gas phase and in solution. In the gas phase the E isomers of ΔPhe have lower energies than the Z ones, but in the aqueous solution the energy order is reversed.

Estimation of formamide harmonic and anharmonic modes in the Kohn-Sham limit using the polarization consistent basis sets.

Formamide harmonic and anharmonic frequencies of fundamental vibrations in the gas phase and in several solvents were successfully estimated in the B3LYP Kohn-Sham complete basis set limit (KS CBS). CBS results were obtained by extrapolating a power function (two-parameter formula) to the results calculated with polarization-consistent basis sets. Anharmonic corrections using the second order perturbation treatment (PT2) and hybrid B3LYP functional combined with polarization consistent pc-n (n=0, 1, 2, 3, 4) and several Pople's basis sets were analyzed for all fundamental formamide vibrational modes in the gas phase and solution. Solvent effects were modeled within a PCM method. The anharmonic frequency of diagnostic amide vibration C=O in the gas phase and the CCl(4) solution calculated with the VPT2 method was significantly closer to experimental data than the corresponding harmonic frequency. Both harmonic and anharmonic frequencies of C=O stretching mode decreased linearly with solvent polarity, expressed by relative environment permittivity (ε) ratio (ε-1)/(2ε+1). However, an unphysical behavior of solvent dependence of some low frequency anharmonic amide modes of formamide (e.g., CN stretch, NH(2) scissoring, and NH(2) in plane bend) was observed, probably due to the presence of severe anharmonicity and Fermi resonance.

Extrapolation of water and formaldehyde harmonic and anharmonic frequencies to the B3LYP/CBS limit using polarization consistent basis sets.

The harmonic and anharmonic frequencies of fundamental vibrations in formaldehyde and water were successfully estimated using the B3LYP Kohn-Sham limit. The results obtained with polarization- and correlation-consistent basis sets were fitted with a two-parameter formula. Anharmonic corrections were obtained by a second order perturbation treatment (PT2). We compared the performance of the PT2 scheme on the two title molecules using SCF, MP2 and DFT (BLYP, B3LYP, PBE and B3PW91 functionals) methods combined with polarization consistent pc-n (n = 0, 1, 2, 3, 4) basis sets, Dunning's basis sets (aug)-cc-pVXZ where X = D, T, Q, 5, 6 and Pople's basis sets up to 6-311++G(3df,2pd). The influence of SCF convergence level and density grid size on the root mean square of harmonic and anharmonic frequency deviations from experimental values was tested. The wavenumber of formaldehyde CH(2) anharmonic asymmetric stretching mode is very sensitive to grid size for large basis sets; this effect is not observed for harmonic modes. BLYP-calculated anharmonic frequencies consistently underestimate observed wavenumbers. On the basis of formaldehyde anharmonic frequencies, we show that increasing the Pople basis set size does not always lead to improved agreement between anharmonic frequencies and experimental values.

The effect of beta-methylation on the conformation of alpha, beta-dehydrophenylalanine: a DFT study.

Dehydroamino acids are non-coded amino acids that offer unique conformational properties. Dehydrophenylalanine (DeltaPhe) is most commonly used to modify bioactive peptides to constrain the topography of the phenyl ring in the side chain, which commonly serves as a pharmacophore. The Ramachandran maps (in the gas phase and in CHCl(3) mimicking environments) of DeltaPhe analogues with methyl groups at the beta position of the side chain as well as at the C-terminal amide were calculated using the B3LYP/6-31 + G** method. Unexpectedly, beta-methylation alone results in an increase of conformational freedom of the affected DeltaPhe residue. However, further modification by introducing an additional methyl group at C-terminal methyl amide results in a steric crowding that fixes the torsion angle psi of all conformers to the value 123 degrees , regardless of the Z or E position of the phenyl ring. The number of conformers is reduced and the accessible conformational space of the residues is very limited. In particular, (Z)-Delta(betaMe)Phe with the tertiary C-terminal amide can be classified as the amino acid derivative that has a single conformational state as it seems to adopt only the beta conformation.