Proton Potential in Acetylacetone

Computational and experimental examinations of new antitumor palladium(II) complex: CT-DNA-/BSA-binding, in-silico prediction, DFT perspective, docking, molecular dynamics simulation and ONIOM

Since the design of metal complexes with better biological activities is important, herein a new palladium(II) complex bearing en and acac (en and acac stand for ethylenediamine and acetylacetonato, respectively) as its ligands, [Pd(en)(acac)]NO3 complex, was synthesized and fully characterized. Quantum chemical computations of the palladium(II) complex were done via DFT/B3LYP method. Cytotoxicity activity of the new compound on leukemia cell line (K562) was assessed via MTT method. The findings indicated that the metal complex has remarkable cytotoxic effect than cisplatin. OSIRIS DataWarrior software was employed to calculate in-silico physicochemical and toxicity parameters of the synthesized complex which rendered significant results. To comprehend the interaction type of new metal compound with macromolecules, the in depth investigation of interaction of mentioned complex with CT-DNA and BSA was accomplished by fluorescence, UV-Visible absorption spectroscopy, viscosity measurement, gel electrophoresis, FRET analysis and circular dichroism (CD) spectroscopy. On the other hand, computational molecular docking was carried out and the obtained data demonstrated that H-bond and van der Waals forces are the dominant forces for the binding of the compound to the mentioned biomolecules. Molecular dynamics simulation was also done and confirmed the stability of best docked pose of palladium(II) complex inside DNA or BSA over the time and in presence of water solvent. Also, Our own N-layered Integrated molecular Orbital and molecular Mechanics (ONIOM) methodology based on the hybridization of quantum mechanics and molecular mechanics (QM/MM) methodology was accomplished to inquire about binding of Pd(II) complex with DNA or BSA.HIGHLIGHTSNew biologically active Pd(II) complex was synthesized and characterized.The in silico studies of the designed complex and its ligands were accomplished by OSIRIS DataWarrior softwareInteraction with CT-DNA and BSA was assessed by various spectroscopic methods.Molecular docking simulation supported the interaction with both macromolecules.Based on ONIOM analysis, the structures of the complex and biomolecules are altered after binding. Communicated by Ramaswamy H. Sarma.

Investigating Intramolecular H Atom Transfer Dynamics in β-Diketones with Ultrafast Infrared Spectroscopies and Theoretical Modeling

The vibrational signatures and ultrafast dynamics of the intramolecular H-bond in a series of β-diketones are investigated with 2D IR spectroscopy and computational modeling. The chosen β-diketones exhibit a range of H atom donor-acceptor distances and asymmetry along the H atom transfer coordinate that tunes the intramolecular H-bond strength. The species with the strongest H-bonds are calculated to have very soft H atom potentials, resulting in highly red-shifted OH stretch fundamental frequencies and dislocation of the H atom upon vibrational excitation. These soft potentials lead to significant coupling to the other normal mode coordinates and give rise to the very broad vibrational signatures observed experimentally. The 2D IR spectra in both the OH and OD stretch regions of the light and deuterated isotopologues reveal broadened and long-lived ground-state bleach signatures of the vibrationally hot molecules. Polarization-sensitive transient absorption measurements in the OH and OD stretch regions reveal notable isotopic differences in orientational dynamics. Orientational relaxation was measured to occur on ∼600 fs and ∼2 ps time scales for the light and deuterated isotopologues, respectively. The orientational dynamics are interpreted in terms of activated H/D atom transfer events driven by collective intramolecular structural rearrangements.

Vibrational Dynamics of the Intramolecular H-Bond in Acetylacetone Investigated with Transient and 2D IR Spectroscopy

Acetylacetone (AcAc) has proven to be a fruitful but highly challenging model system for the experimental and computational interrogation of strong intramolecular hydrogen bonds. Key questions remain, however, regarding the identity of the minimum-energy structure of AcAc and the dynamics of intramolecular proton transfer. Here, we investigate the OH/OD stretch and bend regions of the enol tautomer of AcAc and its deuterated isotopologue with transient absorption and 2D IR spectroscopy. The OH bend region reveals a single dominant diagonal transition near 1625 cm^-1 with intense cross peaks to lower-frequency modes, demonstrating highly mixed fingerprint transitions that contain OH bend character. The anharmonic coupling of the OH bend results in a highly elongated OH bend excited-state absorption transition that indicates a large manifold of OH bend overtone/combination bands in the OH stretch region that leads to strong bend-stretch Fermi resonance interactions. The OH and OD stretch regions consist of broad ground-state bleach signals, but there is no clear evidence of ω₂₁ excited-state absorptions due to rapid population relaxation arising from strong intramolecular coupling to bending, fingerprint, and low-frequency H-bond modes. Orientational relaxation dynamics persist for timescales longer than the vibrational lifetimes, with polarization anisotropy components decaying within approximately 2 and 10 periods of the O-O oscillation for the OH and OD stretch, respectively. The significant isotopic dependence of the orientational dynamics is discussed in the context of intramolecular mode coupling, diffusional processes, and contributions from proton/deuteron transfer dynamics.

Full-dimensional potential energy surface for acetylacetone and tunneling splittings

New, full-dimensional potential energy surface for acetylacetone allows for description of H-tunneling dynamics and characterization of stationary points.

NMR and IR Investigations of Strong Intramolecular Hydrogen Bonds

For the purpose of this review, strong hydrogen bonds have been defined on the basis of experimental data, such as OH stretching wavenumbers, νOH, and OH chemical shifts, δOH (in the latter case, after correction for ring current effects). Limits for O-H···Y systems are taken as 2800 > νOH > 1800 cm-1, and 19 ppm > δOH > 15 ppm. Recent results as well as an account of theoretical advances are presented for a series of important classes of compounds such as β-diketone enols, β-thioxoketone enols, Mannich bases, proton sponges, quinoline N-oxides and diacid anions. The O···O distance has long been used as a parameter for hydrogen bond strength in O-H···O systems. On a broad scale, a correlation between OH stretching wavenumbers and O···O distances is observed, as demonstrated experimentally as well as theoretically, but for substituted β-diketone enols this correlation is relatively weak.

Proton transfer in acetylacetone and its α-halo derivatives

A two-dimensional potential energy function has been applied to study the bent intramolecular H-bonds within acetylacetone and its α-halo derivatives. The theoretically predicted proton transfer barrier heights correlate very well with geometrical parameters and electronic properties related to the H-bond strength.

Role of topological charge stabilization in protomeric tautomerism

Protomeric tautomerism is analyzed in view of the topological charge stabilization rules. Based on Hückel molecular orbital considerations and modern DFT calculations, it was found that the branching of amino or hydroxyl groups significantly contributes to the stability of major species through the first- and second-order perturbations with respect to the isoelectronic hydrocarbon. While amino-imino tautomerism is almost completely dominated by topological charge stabilization, hydroxyl-oxo tautomerism is affected by changes in the resonance integral of C-O/C═O bonds. Nevertheless, apart from side effects such as hydrogen bonds or solvent effects, a quantitative preference rule for the prediction of the tautomeric stability can be developed using topological π-electron energetics. As well as the analyses of simple bases, applications to complex or extended systems are exemplified analyzing purine bases, polyguanide, and polyuret. The present approach can be useful in conjunction with chemical intuition that comes from conventional valence bond theory.

Intramolecular OHO bonding in dibenzoylmethane: symmetry and spectral manifestations

Calculations reveal that both symmetric and asymmetric structures of dibenzoylmethane are present in the gas phase at room temperature.

Pharmacophore-based database mining for probing fragmental drug-likeness of diketo acid analogues

A number of the structurally diverse chemical compounds with functional diketo acid (DKA) subunit(s) have been revealed by combined online and MoStBiodat 3D pharmacophore-guided ZINC and PubChem database screening. We used the structural data available from such screening to analyse the similarities of the compounds containing the DKA fragment. Generally, the analysis by principal component analysis and self-organizing neural network approaches reveals four families of compounds complying with the chemical constitution (aromatic, aliphatic) of the compounds. From a practical point of view, similar studies may reveal potential bioisosteres of known drugs, e.g. raltegravir/elvitegravir. In this context, it seems that mono-halogenated aryl substructures with para group show the closest similarity to these compounds, in contrast to structures where the aromatic ring is halogenated in both ortho- and para-locations.

Car-Parrinello simulation of hydrogen bond dynamics in sodium hydrogen bissulfate

We studied proton dynamics of a short hydrogen bond of the crystalline sodium hydrogen bissulfate, a hydrogen-bonded ferroelectric system. Our approach was based on the established Car-Parrinello molecular dynamics (CPMD) methodology, followed by an a posteriori quantization of the OH stretching motion. The latter approach is based on snapshot structures taken from CPMD trajectory, calculation of proton potentials, and solving of the vibrational Schrodinger equation for each of the snapshot potentials. The so obtained contour of the OH stretching band has the center of gravity at about 1540 cm(-1) and a half width of about 700 cm(-1), which is in qualitative agreement with the experimental infrared spectrum. The corresponding values for the deuterated form are 1092 and 600 cm(-1), respectively. The hydrogen probability densities obtained by solving the vibrational Schrodinger equation allow for the evaluation of potential of mean force along the proton transfer coordinate. We demonstrate that for the present system the free energy profile is of the single-well type and features a broad and shallow minimum near the center of the hydrogen bond, allowing for frequent and barrierless proton (or deuteron) jumps. All the calculated time-averaged geometric parameters were in reasonable agreement with the experimental neutron diffraction data. As the present methodology for quantization of proton motion is applicable to a variety of hydrogen-bonded systems, it is promising for potential use in computational enzymology.

Car-Parrinello simulation of an O-H stretching envelope and potential of mean force of an intramolecular hydrogen bonded system: application to a Mannich base in solid state and in vacuum

Car-Parrinello molecular dynamics (CPMD) study was performed for an anharmonic system-an intramolecularly hydrogen bonded Mannich-base-type compound, 4,5-dimethyl-2(N,N-dimethylaminemethyl)phenol, to investigate the vibrational spectrum associated with the O-H stretching. Calculations were carried out for the solid state and for an isolated molecule. The classical CPMD simulation was performed and then the proton potential snapshots were extracted from the trajectory. The vibrational Schrodinger equation for the snapshots was solved numerically, and the (O-H) envelope was calculated as a superposition of the 0-->1 transitions. The potential of mean force for the proton stretching mode was calculated from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling, nuclear quantum effects, and effects of the environment. Perspectives for application of the presented methodology in the computational support of biocatalysis are given in the study.

Luminescence and Raman Spectra of Acetylacetone at Low Temperatures

Raman spectra of acetylacetone were recorded for molecules isolated in an argon matrix at 10 K and for a polycrystalline sample. In the solid sample, broad bands appear superimposed on a much weaker Raman spectrum corresponding mainly to the stable enol form. The position of these bands depends on the excitation wavelength (514.5 and 488.8 nm argon ion laser lines were used), sample temperature, and cooling history. They are attributed to transitions from an excited electronic state to various isomer states in the ground electronic state. Laser photons have energies comparable to energies of a number of excited triplet states predicted for a free acetylacetone molecule (Chen, X.-B.; Fang, W.-H.; Phillips, D. L. J. Phys. Chem. A 2006, 110, 4434). Since singlet-to-triplet photon absorption transitions are forbidden, states existing in the solid have mixed singlet/triplet character. Their decay results in population of different isomer states, which except for the lowest isomers SYN enol, TS2 enol (described in Matanović I.; Doslić, N. J. Phys. Chem. A 2005, 109, 4185), and the keto form, which can be detected in the Raman spectra of the solid, are not vibrationally resolved. Differential scanning calorimetry detected two signals upon cooling of acetylacetone, one at 229 K and one at 217 K, while upon heating, they appear at 254 and 225 K. The phase change at higher temperature is attributed to a freezing/melting transition, while the one at lower temperature seems to correspond to freezing/melting of keto domains, as suggested by Johnson et al. (Johnson, M. R.; Jones, N. H.; Geis, A; Horsewill. A. J.; Trommsdorff, H. P. J. Chem. Phys. 2002, 116, 5694). Using matrix isolation in argon, the vibrational spectrum of acetylacetone at 10 K was recorded. Strong bands at 1602 and 1629 cm(-1) are assigned as the SYN enol bands, while a weaker underlying band at 1687 cm(-1) and a medium shoulder at 1617 cm(-1) are assigned as TS2 enol bands.

Experimental and Theoretical UV Characterizations of Acetylacetone and Its Isomers

Cryogenic matrix isolation experiments have allowed the measurement of the UV absorption spectra of the high-energy non-chelated isomers of acetylacetone, these isomers being produced by UV irradiation of the stable chelated form. Their identification has been done by coupling selective UV-induced isomerization, infrared spectroscopy, and harmonic vibrational frequency calculations using density functional theory. The relative energies of the chelated and non-chelated forms of acetylacetone in the S0 state have been obtained using density functional theory and coupled-cluster methods. For each isomer of acetylacetone, we have calculated the UV transition energies and dipole oscillator strengths using the excited-state coupled-cluster methods, including EOMCCSD (equation-of-motion coupled-cluster method with singles and doubles) and CR-EOMCCSD(T) (the completely renormalized EOMCC approach with singles, doubles, and non-iterative triples). For dipole-allowed transition energies, there is a very good agreement between experiment and theory. In particular, the CR-EOMCCSD(T) approach explains the blue shift in the electronic spectrum due to the formation of the non-chelated species after the UV irradiation of the chelated form of acetylacetone. Both experiment and CR-EOMCCSD(T) theory identify two among the seven non-chelated forms to be characterized by red-shifted UV transitions relative to the remaining five non-chelated isomers.

The C2v Structure of Enolic Acetylacetone

It has often been postulated that the lowest energy enolic form of Acetylacetone (AcAc) assumes C(s) symmetry, i.e., has a double-minimum potential possibly exhibiting a low barrier to internal proton transfer and not a single minimum, C(2v). Recent theoretical calculations and experimental work support the C(s) hypothesis but the literature on this fascinating molecule is divided. Toward this objective, the high-resolution rotational spectra of enolic acetylacetone and 3 isotopologues have been obtained, revealing C(2v) symmetry. The two methyl groups exhibit a very low barrier to internal rotation, thus making AcAc internally highly dynamic.

H-bonded complexes between acetylacetone and two molecules of methanol: HF and DFT level study

Five stable H-bonded complexes (supersystems) between acetylacetone and two methanol molecules were investigated at the B3LYP and HF levels of theory using the 6-311G** and 6-11++G** basis sets. The most stable complex was found as the one with the highest relative bonding and interaction energies. All vibrational frequencies resulting from calculations with the 6-311++G** basis set were compared with the recorded IR spectrum of acetylacetone/methanol mixture in a molar ratio 1:2.

Theoretical IR spectroscopy based on QM/MM calculations provides changes in charge distribution, bond lengths, and bond angles of the GTP ligand induced by the Ras-protein

The GTPase Ras p21 is a crucial switch in cellular signal transduction. Fourier transform infrared (FTIR) spectra of the substrate guanosine triphosphate (GTP) show remarkable changes when it binds to the enzyme. The reduced band widths indicate that the flexible GTP molecule is guided by the protein into a preferred conformation. The delocalized phosphate vibrations of unbound GTP become localized. The frequency shifts show an electron movement toward beta-phosphate, which probably contributes to catalysis by reducing the free activation energy. To quantify these qualitative observations we performed QM/MM molecular dynamics simulations of Ras.GTP and GTP in water. The triphosphate part of GTP was treated quantum mechanically using density functional theory (DFT). Vibrational spectra were calculated in harmonic approximation with an average deviation of 3% from the experimental frequencies. This provides a high confidence in the computational results as vibrational spectra are highly sensitive to conformation and charge distribution. As compared to GTP in water, Ras-bound GTP shows a shift of negative charge of approximately 0.2 e toward the beta-phosphate from gamma-phosphate and from alpha-phosphate due to the positive charge of the magnesium ion, to a lesser extent of Lys-16, and surprisingly without any effect of the P-loop backbone. Magnesium and Gly-13 twist and bend the gamma-O-beta bonds such that the crucial bond is stretched before cleaving.

Infrared Spectroscopy of the Intramolecular Hydrogen Bond in Acethylacetone: A Computational Approach

The intramolecular hydrogen bond in the enol-acethylacetone (ACAC) is investigated by performing reduced-dimensional quantum calculations. To analyze the shared proton vibrations, two sets of coordinates were employed: normal mode coordinates describing the motion in the vicinity of the most stable configuration, and internal coordinates accounting for the double minimum proton motion. It is proved that the extreme broadness of the OH-stretch band in ACAC is a consequence of the coexistence of two enol-ACAC structures: the global minimum and the transition state for rotation of the distal methyl group. Further, a ground-state tunneling splitting of 116 cm(-1) is found, and it is shown that the inclusion of the kinematic coupling is mandatory when treating large-amplitude proton motion. In the OH-stretch direction a splitting of 853 cm(-1) was predicted.