A QM/MM Study of the Absorption Spectrum of Harmane in Water Solution and Interacting with DNA: The Crucial Role of Dynamic Effects

Thermodynamics and Intermolecular Interactions during the Insertion of Anionic Naproxen into Model Cell Membranes

The insertion process of Naproxen into model dimyristoylphosphatidylcholine (DMPC) membranes is studied by resorting to state-of-the-art classical and quantum mechanical atomistic computational approaches. Molecular dynamics simulations indicate that anionic Naproxen finds an equilibrium position right at the polar/nonpolar interphase when the process takes place in aqueous environments. With respect to the reference aqueous phase, the insertion process faces a small energy barrier of ≈5 kJ mol^-1 and yields a net stabilization of also ≈5 kJ mol^-1. Entropy changes along the insertion path, mainly due to a growing number of realizable microstates because of structural reorganization, are the main factors driving the insertion. An attractive fluxional wall of noncovalent interactions is characterized by all-quantum descriptors of chemical bonding (natural bond orbitals, quantum theory of atoms in molecules, noncovalent interaction, density differences, and natural charges). This attractive wall originates in the accumulation of tiny transfers of electron densities to the interstitial region between the fragments from a multitude of individual intermolecular contacts stabilizing the tertiary drug/water/membrane system.

Optical properties of photodynamic therapy drugs in different environments: the paradigmatic case of temoporfin

Computational tools have been used to study the photophysical and photochemical features of photosensitizers in photodynamic therapy (PDT) - a minimally invasive, less aggressive alternative for cancer treatment. PDT is mainly based on the activation of molecular oxygen through the action of a photoexcited sensitizer (photosensitizer). Temoporfin, widely known as mTHPC, is a second-generation photosensitizer, which produces the cytotoxic singlet oxygen when irradiated with visible light and hence destroys tumor cells. However, the bioavailability of the mostly hydrophobic photosensitizer, and hence its incorporation into cells, is fundamental to achieve the desired effect on malignant tissues via PDT. In this study, we focus on the optical properties of the temoporfin chromophore in different environments -in vacuo, in solution, encapsulated in drug delivery agents, namely cyclodextrin, and interacting with a lipid bilayer.

2,3-(Dibenzimidazol-2-yl)quinoxalines: Unexpected Dynamical Effect on Steady-State Electronic Absorption Spectra

We report on the electronic absorption spectra, conformational behavior, and intra- and intermolecular hydrogen bonds of 2,3-(dibenzimidazol-2-yl)-quinoxaline (DBIQ). The experimentally found strong solvent dependence of the absorption spectra of DBIQ solutions cannot be assigned to electronic excitations of the equilibrium ground-state DBIQ structure. Extended consideration including the nonequilibrium structures within the framework of ab initio molecular dynamics (MD) revealed the importance of torsion molecular motions not covered by the static case. The strong impact of solute-solvent hydrogen bonding on stabilization of these nonequilibrium structures and on conformational composition of DBIQ was demonstrated. A presence of twisted nonplanar geometries along the whole MD trajectory was shown to drastically influence not only energies but also characters of electronic excitations, resulting in a change of local π-π* character in a solution of 1,2-dichloroethane to charge-transfer character in polar dimethylsulfoxide.

Interplay between conformational and solvent effects in UV-visible absorption spectra: curcumin tautomers as a case study

We present a combined theoretical and experimental study on the UV-vis spectra of enol-keto (EK) and keto-keto (KK) tautomeric forms of curcumin dissolved in aqueous solution. Solvent effects have been investigated by resorting to the implicit polarizable continuum model (QM/PCM) and non-polarizable and fully polarizable QM/MM approaches, the latter based on the fluctuating charges (FQ) force-field. In particular, all methods are challenged to rationalize the contribution of conformational, electrostatic and polarization effects in the calculation of the vertical excitation spectra of curcumin tautomers. The obtained results highlight that for both tautomers specific solute-solvent hydrogen-bond interactions play a minor role with respect to conformational and electrostatic effects.

Novel Molecular-Dynamics-Based Protocols for Phase Space Sampling in Complex Systems

The adequate exploration of the phase space of a chromophore is a fundamental necessity for the simulation of their optical and photophysical properties, taking into account the effects of vibrational motion and, most importantly, the coupling with a (non-homogeneous) molecular environment. A representative set of conformational snapshots around the Franck-Condon region is also required to perform non-adiabatic molecular dynamics, for instance in the framework of surface hopping. Indeed, in the latter case one needs to prepare a set of initial conditions providing a meaningful and complete statistical base for the subsequent trajectory propagation. In this contribution, we propose two new protocols for molecular dynamics-based phase space sampling, called "local temperature adjustment" and "individual QM/MM-based relaxation." These protocols are intended for situations in which the popular Wigner distribution sampling procedure is not applicable-as it is the case when anharmonic or nonlinear vibrations are present-and where regular molecular dynamics sampling might suffer from an inaccurate distribution of internal energy or from inaccurate force fields. The new protocols are applied to the case of phase space sampling of [Re(CO)3(Im)(Phen)]+ (im, imidazole; phen, phenanthroline) in aqueous solution, showing the advantages and limitations of regular Wigner and molecular dynamics sampling as well as the strengths of the new protocols.

Vibrational Sampling and Solvent Effects on the Electronic Structure of the Absorption Spectrum of 2-Nitronaphthalene

The influence of vibrational motion on electronic excited state properties is investigated for the organic chromophore 2-nitronaphtalene in methanol. Specifically, the performance of two vibrational sampling techniques - Wigner sampling and sampling from an ab initio molecular dynamics trajectory- is assessed, in combination with implicit and explicit solvent models. The effects of the different sampling/solvent combinations on the energy and electronic character of the absorption bands are analyzed in terms of charge transfer and exciton size, computed from the electronic transition density. The absorption spectra obtained using sampling techniques and its underlying properties are compared to those of the electronic excited states calculated at the Franck-Condon equilibrium geometry. It is found that the absorption bands of the vibrational ensembles are red-shifted compared to the Franck-Condon bright states, and this red-shift scales with the displacement from the equilibrium geometry. Such displacements are found larger and better described when using ensembles from the harmonic Wigner distribution than snapshots from the molecular dynamics trajectory. Particularly relevant is the torsional motion of the nitro group that quenches the charge transfer character of some of the absorption bands. This motion, however, is better described in the molecular dynamics trajectory. Thus, none of the vibrational sampling approaches can satisfactorily capture all important aspects of the nuclear motion. The inclusion of solvent also red-shifts the absorption bands with respect to the gas phase. This red-shift scales with the charge-transfer character of the bands and is found larger for the implicit than for the explicit solvent model. The advantages and drawbacks of the different sampling and solvent models are discussed to guide future research on the calculation of UV-vis spectra of nitroaromatic compounds.

Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules

Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.

DNA damage photo-induced by chloroharmine isomers: hydrolysis versus oxidation of nucleobases

Photodynamic therapy (PDT) is an emerging clinical treatment currently being used against a wide range of both cancerous and noncancerous diseases. The search for new active photosensitizers as well as the development of novel selective delivery systems are the major challenges faced in the application of PDT. We investigated herein three chloroharmine derivatives (6-, 8- and 6,8-dichloroharmines) with quite promising intrinsic photochemical tunable properties and their ability to photoinduce DNA damage in order to elucidate the underlying photochemical mechanisms. Data revealed that the three compounds are quite efficient photosensitizers. The overall extent of photo-oxidative DNA damage induced by both 8-chloro-substituted β-carbolines is higher than that induced by 6-chloro-harmine. The predominant type of lesion generated also depends on the position of the chlorine atom in the β-carboline ring. Both 8-chloro-substituted β-carbolines mostly oxidize purines via type I mechanism, whereas 6-chloro-harmine mainly behaves as a "clean" artificial photonuclease inducing single-strand breaks and site of base loss via proton transfer and concerted (HO^--mediated) hydrolytic attack. The latter finding represents an exception to the general photosensitizing reactions and, to the best of our knowledge, this is the first time that this process is well documented. The controlled and selective production of different oxygen-independent lesions could be fine-tuned by simply changing the substituent groups in the β-carboline ring. This could be a promising tool for the design and development of novel photo-therapeutic agents aimed to tackle hypoxic conditions shown in certain types of tumours.

Photophysics of chlorin e6: from one- and two-photon absorption to fluorescence and phosphorescence

Linear and non-linear optical properties of a known photosensitizer producing singlet oxygen, chlorin e6, have been studied, including dynamics effects.

Deciphering the photosensitization mechanisms of hypericin towards biological membranes

Molecular modeling and simulation allow unravelling the molecular basis of the photosensitization of biological membranes by the hypericin drug used in photodynamic therapy.

Two-photon-absorption DNA sensitization via solvated electron production: unraveling photochemical pathways by molecular modeling and simulation

Infrared optical properties and photochemistry of a potential therapeutic dye interacting with DNA: a multiscale computational approach.

Transition matrices and orbitals from reduced density matrix theory

In this contribution, we report two different methodologies for characterizing the electronic structure reorganization occurring when a chromophore undergoes an electronic transition. For the first method, we start by setting the theoretical background necessary to the reinterpretation through simple tensor analysis of (i) the transition density matrix and (ii) the natural transition orbitals in the scope of reduced density matrix theory. This novel interpretation is made more clear thanks to a short compendium of the one-particle reduced density matrix theory in a Fock space. The formalism is further applied to two different classes of excited states calculation methods, both requiring a single-determinant reference, that express an excited state as a hole-particle mono-excited configurations expansion, to which particle-hole correlation is coupled (time-dependent Hartree-Fock/time-dependent density functional theory) or not (configuration interaction single/Tamm-Dancoff approximation). For the second methodology presented in this paper, we introduce a novel and complementary concept related to electronic transitions with the canonical transition density matrix and the canonical transition orbitals. Their expression actually reflects the electronic cloud polarisation in the orbital space with a decomposition based on the actual contribution of one-particle excitations from occupied canonical orbitals to virtual ones. This approach validates our novel interpretation of the transition density matrix elements in terms of the Euclidean norm of elementary transition vectors in a linear tensor space. A proper use of these new concepts leads to the conclusion that despite the different principles underlying their construction, they provide two equivalent excited states topological analyses. This connexion is evidenced through simple illustrations of (in)organic dyes electronic transitions analysis.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

DNA Electronic Circular Dichroism on the Inter-Base Pair Scale: An Experimental-Theoretical Case Study of the AT Homo-Oligonucleotide

A successful elucidation of the near-ultraviolet electronic circular dichroism spectrum of a short double-stranded DNA is reported. Time-dependent density functional theory methods are shown to accurately predict spectra and assign bands on the microscopic base-pair scale, a finding that opens the field for using circular dichroism spectroscopy as a sensitive nanoscale probe of DNA to reveal its complex interactions with the environment.

Photophysics of acetophenone interacting with DNA: why the road to photosensitization is open

Deoxyribonucleic acid photosensitization, i.e. the photoinduced electron- or energy-transfer of chromophores interacting with DNA, is a crucial phenomenon that triggers important DNA lesions such as pyrimidine dimerization, even upon absorption of relatively low-energy radiation. Oxidative lesions may also be produced via the photoinduced production of reactive oxygen species. Aromatic ketones, and acetophenone in particular, are well known for their sensitization effects. In this contribution we model the structural and dynamical properties of the acetophenone/DNA aggregates as well as their spectroscopic and photophysical properties using high-level hybrid quantum mechanics/molecular mechanics methods. We show that the key steps of the photochemistry of acetophenone in gas phase are conserved in the macromolecular environment and thus an ultrafast singlet-triplet conversion of acetophenone is expected prior to the transfer to DNA.

Spectroscopic behavior of loratadine and desloratadine in different aqueous media conditions studied by means of TD-DFT calculations

In this study, we explained the influence of the stepwise protonation of two antihistaminic drugs on their experimental absorption spectra. We demonstrated the capability of the TD-CAM-B3LYP method, combined with a polarizable continuum model, to produce good performance for the calculated spectra. The lowest energy transitions and the molecular orbital plots were analyzed in detail. The calculated UV spectra are proposed as potential alternatives to initialize the well-known MCR-ALS algorithm, especially when the spectra of the pure analytes are not available. Moreover, it can be a useful strategy for planning an experimental methodology oriented to multiway analysis when the drug species exhibit acid-base properties.