Fluorescence of PRODAN in water: A computational QM/MM MD study

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Reconstruction of Nuclear Ensemble Approach Electronic Spectra Using Probabilistic Machine Learning

The theoretical prediction of molecular electronic spectra by means of quantum mechanical (QM) computations is fundamental to gain a deep insight into many photophysical and photochemical processes. A computational strategy that is attracting significant attention is the so-called Nuclear Ensemble Approach (NEA), that relies on generating a representative ensemble of nuclear geometries around the equilibrium structure and computing the vertical excitation energies (ΔE) and oscillator strengths (f) and phenomenologically broadening each transition with a line-shaped function with empirical full-width δ. Frequently, the choice of δ is carried out by visually finding the trade-off between artificial vibronic features (small δ) and over-smoothing of electronic signatures (large δ). Nevertheless, this approach is not satisfactory, as it relies on a subjective perception and may lead to spectral inaccuracies overall when the number of sampled configurations is limited due to an excessive computational burden (high-level QM methods, complex systems, solvent effects, etc.). In this work, we have developed and tested a new approach to reconstruct NEA spectra, dubbed GMM-NEA, based on the use of Gaussian Mixture Models (GMMs), a probabilistic machine learning algorithm, that circumvents the phenomenological broadening assumption and, in turn, the use of δ altogether. We show that GMM-NEA systematically outperforms other data-driven models to automatically select δ overall for small datasets. In addition, we report the use of an algorithm to detect anomalous QM computations (outliers) that can affect the overall shape and uncertainty of the NEA spectra. Finally, we apply GMM-NEA to predict the photolysis rate for HgBrOOH, a compound involved in Earth's atmospheric chemistry.

Solvent-Dependent Excited-State Evolution of Prodan Dyes

Excited-state character and dynamics of two 6-(dimethylamino)-2-acylnaphthalene dyes (Prodan and Badan-SCH₂CH₂OH) were studied by picosecond time-resolved IR spectroscopy (TRIR) in solvents of different polarity and relaxation times: hexane, CD₃OD, and glycerol-d₈. In all these solvents, near-UV excitation initially produced the same S₁(ππ*) excited state characterized by a broad TRIR signal. A very fast decay (3, ∼100 ps) followed in hexane, whereas conversion to a distinct IR spectrum with a ν(C═O) band downshifted by 76 cm^-1 occurred in polar/H-bonding solvents, slowing down on going from CD₃OD (1, 23 ps) to glycerol-d₈ (5.5, 51, 330 ps). The final relaxed excited state was assigned as planar Me₂N → C═O intramolecular charge transfer S₁(ICT) by comparing experimental and TDDFT-calculated spectra. TRIR conversion kinetics are comparable to those of early stages of multiexponential fluorescence decay and dynamic fluorescence red-shift. This work presents a strong evidence that Prodan-type dyes undergo solvation-driven charge separation in their S₁ state, which is responsible for the dynamic fluorescence Stokes shift observed in polar/H-bonding solvents. The time evolution of the optically prepared S₁(ππ*) state to the S₁(ICT) final state reflects environment relaxation and solvation dynamics. This finding rationalizes the widespread use of Prodan-type dyes as probes of environment dynamics and polarity.

Photophysical Properties of BADAN Revealed in the Study of GGBP Structural Transitions

The fluorescent dye BADAN (6-bromoacetyl-2-dimetylaminonaphtalene) is widely used in various fields of life sciences, however, the photophysical properties of BADAN are not fully understood. The study of the spectral properties of BADAN attached to a number of mutant forms of GGBP, as well as changes in its spectral characteristics during structural changes in proteins, allowed to shed light on the photophysical properties of BADAN. It was shown that spectral properties of BADAN are determined by at least one non-fluorescent and two fluorescent isomers with overlapping absorbing bands. It was found that BADAN fluorescence is determined by the unsolvated "PICT" (planar intramolecular charge transfer state) and solvated "TICT" (twisted intramolecular charge transfer state) excited states. While "TICT" state can be formed both as a result of the "PICT" state solvation and as a result of light absorption by the solvated ground state of the dye. BADAN fluorescence linked to GGBP/H152C apoform is quenched by Trp 183, but this effect is inhibited by glucose intercalation. New details of the changes in the spectral characteristics of BADAN during the unfolding of the protein apo and holoforms have been obtained.

Ultrafast Formation of the Charge Transfer State of Prodan Reveals Unique Aspects of the Chromophore Environment

Lipophilic dyes such as laurdan and prodan are widely used in membrane biology due to a strong bathochromic shift in emission that reports the structural parameters of the membrane such as area per molecule. Disentangling of the factors which control the spectral shift is complicated by the stabilization of a charge-transfer-like excitation of the dye in polar environments. Predicting the emission therefore requires modeling both the relaxation of the environment and the corresponding evolution of the excited state. Here, an approach is presented in which (i) the local environment is sampled by a classical molecular dynamics (MD) simulation of the dye and solvent, (ii) the electronically excited state of prodan upon light absorption is predicted by numerical quantum mechanics (QM), (iii) the iterative relaxation of the environment around the excited dye by MD coupled with the evolution of the excited state is performed, and (iv) the emission properties are predicted by QM. The QM steps are computed using the many-body Green's function in the GW approximation and the Bethe-Salpeter equation with the environment modeled as fixed point charges, sampled in the MD simulation steps. The comparison to ultrafast time-resolved transient absorption measurements demonstrates that the iterative molecular mechanics (MM)/QM approach agrees quantitatively with both the polarity-dependent shift in emission and the time scale over which the charge transfer state is stabilized. Together the simulations and experimental measurements suggest that the evolution into the charge transfer state is slower in amphiphilic solvents.

Quantitative Design of Bright Fluorophores and AIEgens by the Accurate Prediction of Twisted Intramolecular Charge Transfer (TICT)

Inhibition of TICT can significantly increase the brightness of fluorescent materials. Accurate prediction of TICT is thus critical for the quantitative design of high-performance fluorophores and AIEgens. TICT of 14 types of popular organic fluorophores were modeled with time-dependent density functional theory (TD-DFT). A reliable and generalizable computational approach for modeling TICT formations was established. To demonstrate the prediction power of our approach, we quantitatively designed a boron dipyrromethene (BODIPY)-based AIEgen which exhibits (almost) barrierless TICT rotations in monomers. Subsequent experiments validated our molecular design and showed that the aggregation of this compound turns on bright emissions with ca. 27-fold fluorescence enhancement, as TICT formation is inhibited in molecular aggregates.

A solvent–solute cooperative mechanism for symmetry-breaking charge transfer

Random rotational fluctuation of solvents generates a large transient electric field that breaks the symmetry required for symmetry breaking charge transfer.

Orientation of Laurdan in Phospholipid Bilayers Influences Its Fluorescence: Quantum Mechanics and Classical Molecular Dynamics Study

Fluidity of lipid membranes is known to play an important role in the functioning of living organisms. The fluorescent probe Laurdan embedded in a lipid membrane is typically used to assess the fluidity state of lipid bilayers by utilizing the sensitivity of Laurdan emission to the properties of its lipid environment. In particular, Laurdan fluorescence is sensitive to gel vs liquid–crystalline phases of lipids, which is demonstrated in different emission of the dye in these two phases. Still, the exact mechanism of the environment effects on Laurdan emission is not understood. Herein, we utilize dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC) lipid bilayers, which at room temperature represent gel and liquid–crystalline phases, respectively. We simulate absorption and emission spectra of Laurdan in both DOPC and DPPC bilayers with quantum chemical and classical molecular dynamics methods. We demonstrate that Laurdan is incorporated in heterogeneous fashion in both DOPC and DPPC bilayers, and that its fluorescence depends on the details of this embedding.