TD-DFT analysis of excitation of the closed-form spironaphthopyran in methanol solution: The contribution of vibronic transitions and intermolecular hydrogen bonds

A set of forty hybrid functionals, combined with a 6-31++G(d,p) basis set and an IEFPCM solvent description, was applied to calculate the S0→S1 excitation of the spironaphthopyran (SNP) photochromic compound. None of the hybrid functionals resulted in Cspiro-O bond cleavage upon optimization of the S1 excited state of the SNP, in contrast to the spiropyran BIPS which we studied earlier. At the same time, the deformations of the SNP during excitation turned out to be very small (in contrast to those for BIPS) which made it possible to calculate its vibronic absorption spectrum. The best agreement between the calculated and experimental absorption bands was demonstrated by the SOGGA11X functional. The main result of this work is confirmation of the vibration participation in the breaking of the Cspiro-O bond in the S1 excited state of the SNP. This assumption based on experimental data was previously put forward by other authors. Strong H-bonds between the SNP and two methanol molecules enhance the photoinduced distortions of the latter and are slightly weakened upon excitation.

Intramolecular excimers of open forms of 2H-benzopyran, 2H- and 3H-naphthopyrans in solution: TD-DFT/DFT analysis

We hypothesized that in the open forms of diphenyl-substituted photochromic compounds, immediately after the photoinduced cleavage of the C-O bond, one of the phenyl rings forms a stack with an aromatic system at the other end of the alkyl linker. The formation of these intramolecular stacking excimers is made possible by a known increase in the rotational mobility of the linker double bonds in the excited state. The flexibilities of linkers are analyzed in terms of changes in their dihedral angles. After internal conversion, stacking is partially broken but it still prevents the recovery of the C-O bond. This explains the significantly greater stability of the open forms of diphenyl-substituted benzopyran and naphthopyran compared to dimethyl-substituted ones. Open structures with residual stacking are intermediate and by thermal motion, they transform into stable trans-trans and cis-trans forms. However, in our opinion, they significantly affect the equilibrium of closed and open isomers in solution. In particular, the calculated stacking energy of 2H-naphthopyran turned out to be higher than that of 3H-naphthopyran, which partly explains the significantly greater stability of the open form of the former. For the DFT/TD-DFT calculations performed in this work, the following functionals were selected to provide reliable stacking for the open forms of all three studied compounds: APFD, M052X, M06HF, as well as B3LYP with D3(BJ) Grimme dispersion correction. The 6-311++G(d,p) basis set and cyclohexane solvent modeled using IEFPCM were also used.

4a,4b-Dihydrophenanthrene → cis-stilbene photoconversion: TD-DFT/DFT study

CONTEXT

DHP → CS photoconversion was analyzed in terms of electron density redistribution for the first time. The following explanation for the non-recovery of the C4a-C4b bond upon CS relaxation is proposed: during this process, the Coulomb repulsion energy between these pairs of atoms increases by almost one and a half times, and their bonding by an electron at LUMO is insufficient to recover the C4a-C4b bond. According to calculations, upon CS relaxation, the linker connecting the benzene rings undergoes significant structural changes. In this case, the distance between the C4a and C4b atoms increases from 3.00 Å to 3.28 Å. Calculations showed that the C4a-C4b vibration of the DHP bond has a very low intensity. Therefore, thermal motion does not contribute to the rupture of this bond.

METHODS

All calculations were performed using the Gaussian16 software package at the B3LYP/6-311 +  + G(d,p)/IEFPCM theory level. B3LYP was the only hybrid functional supported by Gaussian16, which ensured the cleavage of the C4a-C4b bond of DHP while optimizing its S1 excited state. A quantitative description of the redistribution of electron density in the studied conformers was carried out using the analysis of the NPA of atomic charges. Cyclohexane was used as an implicitly specified non-polar solvent. Visualization of molecular orbitals, and electron densities, as well as plotting of calculated IR spectra, were performed using the Gaussview6 software package.

TD-DFT study of BIPS spiropyran: Effects of functionals and high-polarity solvent on Cspiro O bond dissociation and recovery

We performed a theoretical analysis of the BIPS photochemical cycle using an extensive set of forty hybrid functionals and taking into account a highly polar solvent (methanol). The functionals with a small fraction of the exact Hartree-Fock exchange (%HF) showed the predominant S₀  → S₂ transition with the strengthening of the C_spiro O bond. At the same time, functionals with medium and high %HF (including those with long-range correction) gave a dominant S₀  → S₁ transition with weakening or breaking of the C_spiro O bond, which corresponds to the experimental results. The influence of a highly polar solvent on the photochemical electrocyclic transformations of BIPS turned out to be significant. The number of functionals causing dissociation of the C_spiro O bond decreased from 10 to 7 compared to the gas phase. The magnitude of the oscillator strength has increased by approximately one and a half times. Structural distortions of the BIPS molecule during excitation (both with and without C_spiro O bond cleavage) significantly decreased in methanol compared to the gas phase. The two strong hydrogen bonds of methanol molecules with the oxygen and nitrogen atoms of spiropyran also have a significant effect on its excitation. They lead to a change in the predominant transition from S₀  → S₂ to S₀  → S₁ for five functionals. The number of functionals giving dissociation of the C_spiro O bond decreased from seven to four (M08HX, M052X, CAM-B3LYP, and M11). After the opening of the excited BIPS molecule, both of its strong H-bonds with methanol are preserved. Of this set of four functionals, only M052X and CAM-B3LYP exhibited the dominant HOMO-1 → LUMO configuration observed in high-level computations by other authors. Therefore, both of these functionals are recommended for modeling the photochemical cycle of this spiropyran. The photochemical cycle of BIPS was theoretically analyzed. The redistribution of the electron density in this cycle was quantitatively described using the differences in NPA of the atomic charges. The most important result of this analysis was the electrostatic mechanism of the approach of C_spiro and oxygen atoms at the fourth stage, which causes further reduction of the C_spiro O bond.

Theoretical study of the excitation of proflavine H-dimer in an aqueous solution: the effect of functionals and dispersion corrections

Analysis of the ground and excited states of the H-dimer of the proflavine dye cation (PF) in an aqueous solution was performed by the DFT/TD-DFT method with an implicit specification...

Coumarin 102 excitation in aqueous media: contributions of vibronic coupling and hydration

For the first time, vibronic coupling was considered when analyzing the excitation of coumarin C102.

Photoexcitation of oxazine 170 dye in aqueous solution: TD-DFT study

The vibronic absorption spectra of OX170 dye in an aqueous solution using 40 hybrid functionals, the 6-31 +  + G(d,p) basis set, and the SMD solvent model were calculated. It turned out that the long-range corrected ωB97XD functional provided the best agreement with the experiment in the positions of the main maximum and the short-wavelength shoulder. Calculations showed that this shoulder is vibronic and is not caused by a separate electronic transition. At the same time, the shoulder intensity in the calculated spectrum turned out to be lower than in the experimental one. Various parameters of the OX170 cation in the ground and excited states (IR spectra, atomic charges, dipole moments, and transition moment) were calculated. Maps of the distribution of electron density and electrostatic potential have been built. The influence of four strong hydrogen bonds of the dye with water molecules on the absorption spectrum was analyzed. It was shown that these bonds are strengthened upon OX170 excitation. It was found that explicit assignment of water molecules strongly bound to the dye leads to a redshift of the calculated spectrum by ≈15 nm as a whole, and worsened its shape. Photoexcitation of the dye leads to a noticeable polarization of only one of the four considered water molecules (associated with the endocyclic nitrogen atom in the central ring of the chromophore, the electron density on which increases the most).

The vibronic absorption spectrum and electronic properties of Azure B in aqueous solution: TD-DFT/DFT study

The vibronic absorption spectrum of Azure B (AB) in an aqueous solution is calculated using the time-dependent density functional theory (TD-DFT). The results of calculations are analyzed using all hybrid functionals supported by Gaussian16, the 6-31++G(d,p) basis set, and the IEFPCM and SMD solvent models. The solvent model IEFPCM gave significantly underestimated values of λ_max in comparison with the experiment. This is a manifestation of the TD-DFT "cyanine failure". However, the SMD model made it possible to obtain good agreement between the calculation results and experimental data. The best fit was achieved using the X3LYP functional. According to our calculations, the shoulder in the visible absorption spectrum of AB has a vibronic origin. However, the calculated shoulder is weaker than the experimental one. Explicit assignment of two water molecules, which form strong hydrogen bonds with a dye molecule, leads to a shift of the calculated absorption spectrum to longer wavelengths by approximately 17 nm but does not lead to an improvement in its shape. Comparative analysis of the calculated vibronic absorption spectra of Azure B with those obtained earlier for Azure A and methylene blue showed that the presence and intensity of the short-wavelength shoulder are determined by the location of the bands of higher vibronic transitions relative to the band of the 0₀ → 0⁰ main transitions. Photoexcitation leads to an increase in the dipole moment of the dye molecule. An insignificant photoinduced electron transfer was found in the central ring of the chromophore of the dye molecule.

The vibronic absorption spectra and electronic states of acridine orange in aqueous solution

The time-dependent density functional theory (TD-DFT) was used to obtain vibronic absorption spectra of acridine orange dye (AO) in an aqueous solution that were in good agreement with the experiment. The protonated and neutral forms of the dye have been investigated. The results of calculations using various functionals and basis sets have been analyzed. The best agreement with experiment was given by the level of theory X3LYP/6-31G(d,p). AO molecular orbitals involved in electronic transitions due light absorption in the visible region of the spectrum have been obtained. The dipole moments and atomic charges of the ground and excited states of the AO molecule have been calculated. Maps of the electrostatic potential have been drawn. An insignificant photoinduced electron transfer was found in the central ring of the chromophore of the dye molecule. According to our calculations, the vibronic coupling and the Boltzmann distribution play a significant role in the absorption spectra of the AO.

Antitumor effects and hematotoxicity of C60-Cis-Pt nanocomplex in mice with Lewis lung carcinoma

BACKGROUND

Cisplatin (Cis-Pt) is a widely used anticancer drug but its therapeutic efficiency is limited by hemato-, cardio-, hepato-, nephro- and neurotoxicity. Complexation of Cis-Pt with C₆₀ fullerene nanoparticle will allow to enhance the antitumor activity of the drug and to reduce its side toxic effects.

AIM

To estimate the antitumor effects of С60-Cis-Pt nanocomplex in Lewis lung carcinoma (LLC) and analyze hematological toxicity in tumor-bearing mice.

MATERIALS AND METHODS

Complexation of C₆₀ fullerene and Cis-Pt molecule was studied by computer simulation. С60-Cis-Pt nanocomplex was i.p. injected to LLC-bearing mice in a total dose of 7.5 mg/kg (C₆₀:Cis-Pt as 3.75:3.75 mg/kg). The survival of tumor-bearing mice and the relative reduction of tumor weight was recorded. Blood indices were determined using the Particle Counter PCE-210 automatic hematology analyzer.

RESULTS

Computer simulation demonstrated the formation of С60-Cis-Pt nanocomplex in physiological medium and its stability due to the hydrophobic interactions. Treatment with C₆₀-Cis-Pt nanocomplex increased survival time of LLC-bearing mice by 32%, normalized hemoglobin content (up to 100 g/l), erythrocyte and platelet count as compared to the untreated LLC-bearing mice. Tumor weight decreased by 35.5%; the mitotic index of tumor cells decreased by 78%, and apoptotic index increased by 75%. The revealed effects of the C₆₀-Cis-Pt nanocomplex were more pronounced than the effects of Cis-Pt or C₆₀ fullerene alone in equivalent dose.

CONCLUSION

Treatment with C₆₀-Cis-Pt nanocomplex prolonged the survival of LLC-bearing mice and reduced anemia in LLC-bearing mice.

Interaction of pseudoephedrine and azithromycin with losartan: Spectroscopic, dissolution and permeation studies

This study aims at investigating the potential effect of selected cationic drugs (azithromycin (AZN) and pseudoephedrine sulfate (PSD) on the dissolution profile and intestinal permeation of losartan potassium (LOS) that might occur due to ion pair salt formation. DSC, FT-IR and ¹H NMR indicated the formation of ion pair salts between LOS and each of AZN and PSD. Based on NMR chemical shifts calculations, utilizing specialized software, the most likely structures of the salt were proposed and revealed interesting structural features. The obtained ion pair products were shown to have lower aqueous solubilities (water and phosphate buffer pH 6.8) and higher apparent partition coefficient values compared to the parent compound. Neither of the cations affected the dissolution of LOS tablet (Cozaar® 100 mg) in the studied media (HCl pH 1.2 and phosphate buffer pH 6.8). Interestingly, AZN significantly increased the dissolution of LOS in phosphate buffer pH 4.5 (f₂ = 33), and an explanation based on distinguished association pattern between AZN and LOS (CH/π) was offered. Employing permeation test across Caco-2 cells monolayer, the apparent permeability coefficient (Papp) of LOS increased significantly (from 0.9 × 10^-5 cm/s to 1.8 × 10^-5 cm/s) in the presence of the selected cations. Therefore, while the employed cationic drugs were not shown to form ion pair salts under the in-vitro dissolution conditions, they may still participate in significant in-vivo interaction with LOS.

Comparative Analysis of the Antineoplastic Activity of C60 Fullerene with 5-Fluorouracil and Pyrrole Derivative In Vivo

The antitumor activity of pristine C60 fullerene aqueous solution (C60FAS) compared to 5-fluorouracil (5-FU) and pyrrole derivative 1-(4-Cl-benzyl)-3-Cl-4-(CF3-fenylamino)-1H-pyrrol-2.5-dione (MI-1) cytostatic drugs was investigated and analyzed in detail using the model of colorectal cancer induced by 1.2-dimethylhydrazine (DMH) in rats. The number, size, and location of the tumors were measured, and the pathology was examined. It was found that the number of tumors and total lesion area decreased significantly under the action of C60FAS and MI-1. Because these drugs have different mechanisms of action, their simultaneous administration can potentially increase the effectiveness and significantly reduce the side effects of antitumor therapy.

Electrostatic contribution to the energy of binding of aromatic ligands with DNA

The method of solution of the nonlinear Poisson-Boltzmann equation was used to calculate electrostatic energy of binding of various aromatic ligands with DNA oligomers of different length. Analysis of the electrostatic contribution was made in terms of a two-step DNA binding process: formation of the intercalation cavity and insertion of the ligand. The total electrostatic energy was also partitioned into components: the energy of atom-atom coulombic interactions and the energy of interaction with surrounding water. The results indicate that electrostatic interactions are, as a whole, unfavorable to the intercalation process and that a correct analysis of structure-energy interrelation for Ligand-DNA interactions should only be accomplished at the level of the components rather than at the level of total electrostatic energy.

On the origin of the decrease in stability of the DNA hairpin d(GCGAAGC) on complexation with aromatic drugs

Molecular dynamics simulations of drug-DNA complexes have been carried out in order to explain the experimentally observed decrease in thermal stability of the DNA hairpin d(GCGAAGC) on binding the aromatic drug molecules, daunomycin, ethidium bromide, novantrone and proflavine. This complexation behavior is in contrast to the stabilizing effect of the same aromatic drug molecules on DNA duplexes. Analysis of the energy parameters and the hydration properties of the complexes shows that the main factor correlating with the decrease in melting temperatures of the drug-hairpin complexes is the number of water bridges, with a reduction of at least 40% on ligand binding.