Małgorzata A. Broda

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

AbstractThe 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 CH2 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. FigureSensitivity of water B3LYP calculated harmonic and anharmonic vs(OH) frequencies on selected Pople and polarization consistent basis sets size. The results for pc-n basis sets were fitted with two parameter formula and the CBS(2,3,4) estimated

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

AbstractFormamide 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 CCl4 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, NH2 scissoring, and NH2 in plane bend) was observed, probably due to the presence of severe anharmonicity and Fermi resonance. FigureFormamide 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 analysed for all fundamental formamide vibrational modes in the gas phase and solution.

Structural and Molecular Characterization of meso-Substituted Zinc Porphyrins: A DFT Supported Study

Structural parameters of a range of over 100 meso-substituted zinc porphyrins were reviewed and compared to show how far the nature of the functional group may affect the interatomic distances and bond angles within the porphyrin core. It was proved that even despite evident deformations of the molecular structure, involving twisting of the porphyrins central plane, the coupled π-bonding system remains flexible and stable. DFT calculations were applied to a number of selected porphyrins representative for the reviewed compounds to emphasize the relevance of theoretical methods in structural investigations of complex macrocyclic molecular systems. Experimental and DFT-simulated IR spectral data were reported and analyzed in context of the individual molecular features introduced by the meso substituents into the porphyrin moiety base. Raw experimental spectral data, including 1H- and 13C-NMR, UV-Vis, FTIR, XRD, and other relevant physicochemical details have been provided for a specially chosen reference zinc porphyrin functionalized by tert-butylphenyl groups.

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‐NMe2) 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. Copyright

Effect of tin and lead chlorotriphenyl analogues on selected living cells

Three kinds of living cells, human embryonic kidney cells, Saccharomyces cerevisiae, and Escherichia coli, were tested for their sensitivity to chlorotriphenyltin and chlorotriphenyllead. The tin compound proved definitely more toxic than the lead derivative, particularly in the case of the human embryonic kidney cells devoid of any protective cell wall. Electron paramagnetic resonance (EPR) comparative studies carried out by using a natural model liposome system (egg yolk lecithin) confirmed considerable changes within the lipid bilayer upon doping by the aforementioned additives, which may be crucial to the mechanism of the observed cell cleavage. The individual dopants revealed diverse impact upon the membranes condition, chlorotriphenyltin distinctly fluidized the lipid system, whereas chlorotriphenyllead stiffened the medium within the membrane. A theoretical approach concerning such different behaviors of studied tin and lead analogues because of their high toxicity in living cells has been presented.

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.

The conformational properties of dehydrobutyrine and dehydrovaline: theoretical and solid-state conformational studies.

Dehydrobutyrine is the most naturally occurring dehydroamino acid. It is also the simplest dehydroamino acid having the geometrical isomers E/Z. To investigate its conformational properties, a theoretical analysis was performed on N‐acetyl‐α,β‐dehydrobutyrine N′‐methylamides, Ac‐(E)‐ΔAbu‐NHMe and Ac‐(Z)‐ΔAbu‐NHMe, as well as the dehydrovaline derivative Ac‐ΔVal‐NHMe. The ϕ, ψ potential energy surfaces and the localised conformers were calculated at the B3LYP/6‐311 + + G(d,p) level of theory both in vacuo and with inclusion of the solvent (chloroform, water) effect (SCRF method). The X‐ray crystal structures of Ac‐(Z)‐ΔAbu‐NHMe and Ac‐ΔVal‐NHMe were determined at 85 and 100 K, respectively. The solid‐state conformational preferences for the studied residues have been analysed and compared with the other related structures. Despite the limitations imposed by the Cα = Cβ double bond on the topography of the side chains, the main chains of the studied dehydroamino acids are more flexible than in standard alanine. The studied dehydroamino acids differ in their conformational preferences, which depend on the polarity of the environment. This might be a reason why the nature quite precisely differentiates between ΔVal and each of the ΔAbu isomers, and why, particularly so with the latter, they are used as a conformational tool to influence the biological action of usually small, cyclic dehydropeptides. Copyright

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

AbstractIn 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 AbstractIn 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

β‐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‐NMe2 (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‐NMe2), both in vacuo and in the polar environment. On the contrary, N‐methylated E isomer (Ac‐Gly‐(E)‐ΔPhe‐NMe2) can easier adopt β‐turn conformation, but the backbone torsion angles (ϕ1, ψ1, ϕ2, ψ2) are off the limits for common β‐turn types.

The conformational properties of α,β-dehydroamino acids with a C-terminal ester group

α,β‐Dehydroamino acid esters occur in nature. To investigate their conformational properties, a systematic theoretical analysis was performed on the model molecules Ac‐ΔXaa‐OMe [ΔXaa = ΔAla, (E)‐ΔAbu, (Z)‐ΔAbu, ΔVal] at the B3LYP/6‐311+ + G(d,p) level in the gas phase as well as in chloroform and water solutions with the self‐consistent reaction field‐polarisable continuum model method. The Fourier transform IR spectra in CCl4 and CHCl3 have been analysed as well as the analogous solid state conformations drawn from The Cambridge Structural Database. The ΔAla residue has a considerable tendency to adopt planar conformations C5 (ϕ, ψ ≈ − 180°, 180°) and β2 (ϕ, ψ ≈ − 180°, 0°), regardless of the environment. The ΔVal residue prefers the conformation β2 (ϕ, ψ ≈ − 120°, 0°) in a low polar environment, but the conformations α (ϕ, ψ ≈ − 55°, 35°) and β (ϕ, ψ ≈ − 55°, 145°) when the polarity increases. The ΔAbu residues reveal intermediate properties, but their conformational dispositions depend on configuration of the side chain of residue: (E)‐ΔAbu is similar to ΔAla, whereas (Z)‐ΔAbu to ΔVal. Results indicate that the low‐energy conformation β2 is the characteristic feature of dehydroamino acid esters. The studied molecules constitute conformational patterns for dehydroamino acid esters with various side chain substituents in either or both Z and E positions. Copyright