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.

Thiazole-amino acids: influence of thiazole ring on conformational properties of amino acid residues.

Post-translational modified thiazole-amino acid (Xaa-Tzl) residues have been found in macrocyclic peptides (e.g., thiopeptides and cyanobactins), which mostly inhibit protein synthesis in Gram + bacteria. Conformational study of the series of model compounds containing this structural motif with alanine, dehydroalanine, dehydrobutyrine and dehydrophenylalanine were performed using DFT method in various environments. The solid-state crystal structure conformations of thiazole-amino acid residues retrieved from the Cambridge Structural Database were also analysed. The studied structural units tend to adopt the unique semi-extended ß2 conformation; which is stabilised mainly by N-H⋯NTzl hydrogen bond, and for dehydroamino acids also by π-electron conjugation. The conformational preferences of amino acids with a thiazole ring were compared with oxazole analogues and the role of the sulfur atom in stabilising the conformations of studied peptides was discussed.

Oxidation of diclofenac in the presence of iron(II) octacarboxyphthalocyanine.

This paper presents the results of the research on the influence of catalytic activity of iron(II) octacarboxyphthalocyanines (FePcOC) on the transformation of diclofenac (DCF) which is the most popular anti-inflammatory analgesic. Diclofenac poses a serious threat to the natural environment. The paper demonstrates that diclofenac, in the presence a monomeric form of iron octacarboxyphthalocyanine and hydroxyl radicals (HO•) (from H2O2), undergoes a transformation into diclofenac-2,5-iminoquinone (DCF-2,5-IQ), causing distinct changes in the UV-Vis absorption spectrum. In the presence of iron octacarboxyphthalocyanine and H2O2, the previously colourless diclofenac solution becomes intense orange. As a result, a new band at approx. 450 nm appears in the absorption spectrum. HPLC analysis has shown that the concentration of diclofenac decreases with time. TD-DFT calculations using the CAM-B3LYP/6-31+G (d, p) method have been conducted to confirm experimental data concerning the formation of a new band at λmax = 450 nm.

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.

Hydrogen Bond-Mediated Conjugates Involving Lanthanide Diphthalocyanines and Trifluoroacetic Acid (Lnpc2@TFA): Structure, Photoactivity, and Stability.

The interaction between lanthanide diphthalocyanine complexes, LnPc2 (Ln = Nd, Sm, Eu, Gd, Yb, Lu; Pc = C32H16N8, phthalocyanine ligand) and trifluoroacetic acid (TFA) was investigated in benzene, and the stability of the resulting molecular system was assessed based on spectral (UV-Vis) and kinetic measurements. Structural Density Functional Theory (DFT) calculations provided interesting data regarding the nature of the bonding and allowed estimating the interaction energy between the LnPc2 and TFA species. Conjugates are created between the LnPc2 and TFA molecules via hydrogen bonds of moderate strength (>N∙∙H··) at the meso- -bridges of the Pc moieties, which renders the sandwich system to flatten. Attachment of TFA is followed by rearrangement of electronic density within the chromophore system of the macrocycles manifested in considerable changes in their UV-Vis spectra and consequently the color of the studied solutions (from green to orange). The LnPc2@TFA conjugates including Nd, Sm, Eu, and Gd appeared evidently less photostable when exposed to UV radiation than the related mother compounds, whereas in the case of Yb and Lu derivatives some TFA-prompted stabilizing effect was noticed. The conjugates displayed the capacity for singlet oxygen generation in contrast to the LnPc2s itself. Photon upconversion through sensitized triplet-triplet annihilation was demonstrated by the TFA conjugates of Nd, Sm, Eu, and Gd.

What is the form of muscimol from fly agaric mushroom (Amanita muscaria) in water? An insight from NMR experiment supported by molecular modeling.

The biologically active alkaloid muscimol is present in fly agaric mushroom (Amanita muscaria), and its structure and action is related to human neurotransmitter γ -aminobutyric acid (GABA). The current study reports on determination of muscimol form present in water solution using multinuclear 1 H and 13 C nuclear magnetic resonance (NMR) experiments supported by density functional theory molecular modeling. The structures of three forms of free muscimol molecule both in the gas phase and in the presence of water solvent, modeled by polarized continuous model, and nuclear magnetic isotropic shieldings, the corresponding chemical shifts, and indirect spin-spin coupling constants were calculated. Several J-couplings observed in proton and carbon NMR spectra, not available before, are reported. The obtained experimental spectra, supported by theoretical calculations, favor the zwitterion form of muscimol in water. This structure differs from NH isomer, previously determined in dimethyl sulfoxide (DMSO) solution. In addition, positions of signals C3 and C5 are reversed in both solvents.

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.

Performance of revised STO(1M)-3G basis set for prediction of 5-fluorocytosine chemical shifts.

Nuclear shieldings and chemical shifts of 5-fluorocytosine (5FC) were predicted in the gas phase and DMSO solution modeled by polarizable continuum model using B3LYP density functional and revised STO(1M)-3G basis set. For comparison, eight arbitrary selected basis sets including STO-3G and medium-size Pople-type and larger dedicated Jensen-type ones were applied. The former basis sets were significantly smaller, but the calculated structural parameters, harmonic vibrational frequencies, were very accurate and close to those obtained with larger, polarization-consistent ones. The predicted 13 C and 1 H chemical shieldings of 5FC and cytosine, selected as parent molecule, were acceptable (root mean square for 13 C chemical shifts in DMSO of about 5 ppm and less) though less accurate than those calculated with large basis sets, dedicated for prediction of nuclear magnetic resonance parameters.

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.

Method and basis set dependence of the NICS indexes of aromaticity for benzene.

The role of theory level in prediction of benzene magnetic indexes of aromaticity is analysed and compared with calculated nuclear magnetic shieldings of 3 He used as NMR probe. Three closely related nucleus-independent chemical shift (NICS) based indexes were calculated for benzene at SCF-HF, MP2, and DFT levels of theory and the impact of basis set on these quantities was studied. The changes of benzene NICS(0), NICS(1), and NICS(1)zz parameters calculated using SCF-HF, MP2 and several density functionals were within 1 to 3 ppm. Similar deviations between magnetic indexes of aromaticity were observed for values calculated with selected basis sets. Only very small effect of polar solvent on benzene aromaticity was predicted. The 3 He nuclear magnetic isotropic shielding (σ) and its zz-components (σzz ) of helium atom approaching the centre of benzene ring from above produced similar curves versus benzene-He distance to NICS parameters calculated for similarly moving Bq ghost atom. We also propose an experimental verification of NICS calculations by designing the 3 He NMR measurement for benzene saturated with helium gas or in low temperature matrices.

The influence of solvent on conformational properties of peptides with Aib residue-a DFT study.

The conformational propensities of the Aib residue on the example of two model peptides Ac-Aib-NHMe (1) and Ac-Aib-NMe2 (2), were studied by B3LYP and M06-2X functionals, in the gas phase and in the polar solvents. To verify the reliability of selected functionals, we also performed MP2 calculations for the tested molecules in vacuum. Polarizable continuum models (PCM and SMD) were used to estimate the solvent effect. Ramachandran maps were calculated to find all energy minima. Noncovalent intramolecular interactions due to hydrogen-bonds and dipole attractions between carbonyl groups are responsible for the relative stabilities of the conformers. In order to verify the theoretical results, the available conformations of similar X-ray structures from the Cambridge Crystallographic Data Center (CCDC) were analyzed. The results of the calculations show that both derivatives with the Aib residue in the gas phase prefer structures stabilized by intramolecular N-H⋯O hydrogen bonds, i.e., C5 and C7 conformations, while polar solvent promotes helical conformation with φ, ψ values equal to +/-60°, +/-40°. In addition, in the case of molecule 2, the helical conformation is the only one available in the polar environment. This result is fully consistent with the X-ray data. Graphical abstract Effect of solvent on the Ramachandran maps of the model peptides with Aib residue.

Pyrazole amino acids: hydrogen bonding directed conformations of 3-amino-1H-pyrazole-5-carboxylic acid residue.

A series of model compounds containing 3-amino-1H-pyrazole-5-carboxylic acid residue with N-terminal amide/urethane and C-terminal amide/hydrazide/ester groups were investigated by using NMR, Fourier transform infrared, and single-crystal X-ray diffraction methods, additionally supported by theoretical calculations. The studies demonstrate that the most preferred is the extended conformation with torsion angles ϕ and ψ close to ±180°. The studied 1H-pyrazole with N-terminal amide/urethane and C-terminal amide/hydrazide groups solely adopts this energetically favored conformation confirming rigidity of that structural motif. However, when the C-terminal ester group is present, the second conformation with torsion angles ϕ and ψ close to ±180° and 0°, respectively, is accessible. The conformational equilibrium is observed in NMR and Fourier transform infrared studies in solution in polar environment as well as in the crystal structures of other related compounds. The observed conformational preferences are clearly related to the presence of intramolecular interactions formed within the studied residue. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Interactions of amino acids with aluminum octacarboxyphthalocyanine hydroxide. Experimental and DFT studies.

The influence of albumin and amino acids (L-serine, glycine, L-histidine, L-tryptophan, L-cysteine) on the properties of aluminum octacarboxyphthalocyanine hydroxide (Al(OH)PcOC) was investigated in a phosphate buffer (pH 8.0). Particular attention was paid to the spectroscopic properties and photostability of Al(OH)PcOC. The effect of albumin or amino acids on the photodegradation of Al(OH)PcOC was examined in water using red light: 685 nm and daylight irradiation. Analysis of kinetic curves indicated that interaction with those molecules increases the photostability of Al(OH)PcOC. The molecular structure of Al(OH)PcOC complexes (in vacuum and in water) with axially or equatorially coordinated amino acids was studied by the B3LYP/6-31G* method, and the effects on molecular structure and electronic absorption spectrum were investigated on the basis of the density functional theory. The calculation results revealed that axial coordination significantly reduces the non-planarity of the phthalocyanine ring, and, thus, alters the electronic structure. On the other hand, hydrogen bonding of phthalocyanine side COOH groups with amino acids, in equatorial complexes, does not change the structure within the center of the phthalocyanine, and causes only a slight increase in UV-vis bands intensity, which is in perfect agreement with experimental data. Graphical abstract Structure of equatorial complex of Al(OH)PcOC with L-histidine calculated byB3LYP/6-31G(d) method. Dotted lines H-bonds.

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.

Conformational preferences and synthesis of isomers Z and E of oxazole-dehydrophenylalanine.

Dehydrophenylalanine, ΔPhe, is the most commonly studied α,ß-dehydroamino acid. In nature, further modifications of the α,ß-dehydroamino acids were found, for example, replacement of the C-terminal amide group by oxazole ring. The conformational properties of oxazole-dehydrophenylalanine residue (ΔPhe-Ozl), both isomers Z and E, were investigated. To determine all possible conformations, theoretical calculations were performed using Ac-(Z/E)-ΔPhe-Ozl(4-Me) model compounds at M06-2X/6-31++G(d,p) level of theory. Ac-(Z/E)-ΔPhe-Ozl-4-COOEt compounds were synthesized and the conformational preferences of each isomer, Z and E, were investigated using FTIR and NMR-NOE in solutions of increasing polarity (CHCl3 , DMSO-d6). The solid-state low-temperature structures of Ac-(Z)-ΔPhe-Ozl-4-COOEt and its intermediate analog Ac-(Z)-ΔPhe-Ozn(4-OH)-4-COOEt were also determined. In a weakly polar environment, the ΔPhe-Ozl residue has a tendency to adopt the conformation ß2 with the calculated φ and ψ angles of -127° and 0° for the isomer Z and -170° and 26° for the isomer E. The increase of environment polarity favors the helical conformation α and the beta-turn like conformation ß, but the conformation ß2 seems to be still accessible. The (E)-ΔPhe-Ozl residue can be obtained from the isomer Z in photoisomerization reaction. However, hydroxyl-oxazoline-dehydrophenylalanine ΔPhe-Ozn(4-OH) decomposes in such conditions. Alternatively, (E)-ΔPhe-NH2 can be applied as a substrate in the Hantzsch reaction. © 2016 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 283-294, 2016.

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.