Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

An Expedited Route to Optical and Electronic Properties at Finite Temperature via Unsupervised Learning

Electronic properties and absorption spectra are the grounds to investigate molecular electronic states and their interactions with the environment. Modeling and computations are required for the molecular understanding and design strategies of photo-active materials and sensors. However, the interpretation of such properties demands expensive computations and dealing with the interplay of electronic excited states with the conformational freedom of the chromophores in complex matrices (i.e., solvents, biomolecules, crystals) at finite temperature. Computational protocols combining time dependent density functional theory and ab initio molecular dynamics (MD) have become very powerful in this field, although they require still a large number of computations for a detailed reproduction of electronic properties, such as band shapes. Besides the ongoing research in more traditional computational chemistry fields, data analysis and machine learning methods have been increasingly employed as complementary approaches for efficient data exploration, prediction and model development, starting from the data resulting from MD simulations and electronic structure calculations. In this work, dataset reduction capabilities by unsupervised clustering techniques applied to MD trajectories are proposed and tested for the ab initio modeling of electronic absorption spectra of two challenging case studies: a non-covalent charge-transfer dimer and a ruthenium complex in solution at room temperature. The K-medoids clustering technique is applied and is proven to be able to reduce by ∼100 times the total cost of excited state calculations on an MD sampling with no loss in the accuracy and it also provides an easier understanding of the representative structures (medoids) to be analyzed on the molecular scale.

Nature of the Ultrafast Interligands Electron Transfers in Dye-Sensitized Solar Cells

Charge-transfer dynamics and interligand electron transfer (ILET) phenomena play a pivotal role in dye-sensitizers, mostly represented by the Ru-based polypyridyl complexes, for TiO2 and ZnO-based solar cells. Starting from metal-to-ligand charge-transfer (MLCT) excited states, charge dynamics and ILET can influence the overall device efficiency. In this letter, we focus on N34- dye ( [Ru(dcbpy)2(NCS)2]4-, dcbpy = 4,4'-dicarboxy-2,2'-bipyridine) to provide a first direct observation with high time resolution (<20 fs) of the ultrafast electron exchange between bpy-like ligands. ILET is observed in water solution after photoexcitation in the ∼400 nm MLCT band, and assessment of its ultrafast time-scale is here given through a real-time electronic dynamics simulation on the basis of state-of-the-art electronic structure methods. Indirect effects of water at finite temperature are also disentangled by investigating the system in a symmetric gas-phase structure. As main result, remarkably, the ILET mechanism appears to be based upon a purely electronic evolution among the dense, experimentally accessible, MLCT excited states manifold at ∼400 nm, which rules out nuclear-electronic couplings and proves further the importance of the dense electronic manifold in improving the efficiency of dye sensitizers in solar cell devices.

Electronic and Vibrational Manifold of Tetracyanoethylene-Chloronaphthalene Charge Transfer Complex in Solution: Insights from TD-DFT and Ab Initio Molecular Dynamics

The interplay between light absorption and the molecular environment has a central role in the observed photophysics of a wide range of photoinduced chemical and biological phenomena. The understanding of the interplay between vibrational and electronic transitions is the focus of this work, since it can provide a rationale to tune the optical properties of charge transfer (CT) materials used for technological applications. A clear description of these processes poses a nontrivial challenge from both the theoretical and experimental points of view, where the main issue is how to accurately describe and probe drastic changes in the electronic structure and the ultrafast molecular relaxation and dynamics. In this work we focused on the intermolecular CT reaction that occurs upon photon absorption in a π-stacked model system in dichloromethane solution, in which the 1-chloronaphthalene (1ClN) acts as the electron donor and tetracyanoethylene (TCNE) is the electron acceptor. Density functional theory calculations have been carried out to characterize both the ground-state properties and more importantly the low-lying CT electronic transition, and excellent agreement with recently available experimental results [Mathies, R. A.; et al. J. Phys. Chem. A 2018, 122 (14), 3594] was obtained. The minima of the ground state and first singlet excited state have been accurately characterized in terms of spatial arrangements and vibrational Raman frequencies, and the CT natures of the first two low-lying electronic transitions in the absorption spectra have been addressed and clarified too. Finally, by modeling the possible coordination sites of the TCNE electron acceptor with respect to monovalent ions (Na+, K+) in an implicit solution of acetonitrile, we find that TCNE can accommodate a counterion in two different arrangements, parallel and orthogonal to the C═C axis, leading to the formation of a contact ion pair. The nature of the counterion and its relative position entail structural modifications of the TCNE radical anion, mainly the central C═C and C≡N bonds, compared to the isolated case. An important red shift of the C═C stretching frequency was observed when the counterion is orthogonal to the double bond, to a greater extent for Na+. On the contrary, in the second case, where the counterion ion lies along the internuclear C═C axis, we find that K+ polarizes the electron density of the double bond more, resulting in a greater red shift than with Na+.

Direct observation of the solvent organization and nuclear vibrations of [Ru(dcbpy)2(NCS)2]4-, [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], via ab initio molecular dynamics

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(II) complex, [Ru(dcbpy)₂(NCS)₂]^4- [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C₂ symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies.

Structural Origin and Vibrational Fingerprints of the Ultrafast Excited State Proton Transfer of the Pyranine-Acetate Complex in Aqueous Solution

The excited state proton transfer (ESPT) reaction from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid (HPTS or pyranine) to an acetate molecule has been investigated in explicit aqueous solution via excited state ab initio molecular dynamics simulations based on hybrid quantum/molecular mechanics (QM/MM) potentials. In all the trajectories, the direct proton transfer has been observed in the excited state within 1 ps. We find that the initial structural configuration extracted from the ground state distribution strongly affects the ESPT kinetics. Indeed, the relative orientation of the proton donor-acceptor pair and the presence of a water molecule hydrogen bonded to the phenolic acid group of the pyranine are the key factors to facilitate the ESPT. Furthermore, we analyze the vibrational fingerprints of the ESPT reaction, reproducing the blue shift of the acetate CO stretching (CO_ac), from 1666 to 1763 cm^-1 testifying the transformation of acetate to acetic acid. Finally, our findings suggest that the acetate CC stretching (CC_ac) is also sensitive to the progress of the ESPT reaction. The CC_ac stretching is indeed ruled by the two vibrational modes (928 and 1426 cm^-1), that in the excited state are alternately activated when the proton is shared or bound to the donor/acceptor, respectively.

Water-Mediated Excited State Proton Transfer of Pyranine-Acetate in Aqueous Solution: Vibrational Fingerprints from Ab Initio Molecular Dynamics

In this work, we simulate the excited state proton transfer (ESPT) reaction involving the pyranine photoacid and an acetate molecule as proton acceptor, connected by a bridge water molecule. We employ ab initio molecular dynamics combined with an hybrid quantum/molecular mechanics (QM/MM) framework. Furthermore, a time-resolved vibrational analysis based on the wavelet-transform allows one to identify two low frequency vibrational modes that are fingerprints of the ESPT event: a ring wagging and ring breathing. Their composition suggests their key role in optimizing the structure of the proton donor–acceptor couple and promoting the ESPT event. We find that the choice of the QM/MM partition dramatically affects the photoinduced reactivity of the system. The QM subspace was gradually extended including the water molecules directly interacting with the pyranine–water–acetate system. Indeed, the ESPT reaction takes place when the hydrogen bond network around the reactive system is taken into account at full QM level.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

Use of an Artificial Miniaturized Enzyme in Hydrogen Peroxide Detection by Chemiluminescence

Advanced oxidation processes represent a viable alternative in water reclamation for potable reuse. Sensing methods of hydrogen peroxide are, therefore, needed to test both process progress and final quality of the produced water. Several bio-based assays have been developed so far, mainly relying on peroxidase enzymes, which have the advantage of being fast, efficient, reusable, and environmentally safe. However, their production/purification and, most of all, batch-to-batch consistency may inherently prevent their standardization. Here, we provide evidence that a synthetic de novo miniaturized designed heme-enzyme, namely Mimochrome VI*a, can be proficiently used in hydrogen peroxide assays. Furthermore, a fast and automated assay has been developed by using a lab-bench microplate reader. Under the best working conditions, the assay showed a linear response in the 10.0-120 μM range, together with a second linearity range between 120 and 500 μM for higher hydrogen peroxide concentrations. The detection limit was 4.6 μM and quantitation limits for the two datasets were 15.5 and 186 μM, respectively. In perspective, Mimochrome VI*a could be used as an active biological sensing unit in different sensor configurations.

Unveiling anharmonic coupling by means of excited state ab initio dynamics: application to diarylethene photoreactivity

In this work, excited state ab initio molecular dynamics together with a time resolved vibrational analysis is employed to shed light on the vibrational photoinduced dynamics of a well-known diarylethene molecule experiencing a ring opening reaction upon electronic excitation. The photoreactivity of diarylethenes is recognized to be controlled by a non-adiabatic intersection point between the ground and the first excited state surfaces. The computation of an energy scan, along a suitable reaction coordinate, allows us to identify the region of potential energy surfaces in which the ground (S0) and the first excited (S1) state are well separated. The adiabatic sampling of that region in S1 shows that in the first 3 picoseconds, the central CC bond, which is subject to break, oscillates in an antiphase with respect to the energy gap ΔE(S1 - S0). A multiresolution analysis based on the wavelet transform was then applied to the structural parameters extracted from the excited state dynamics. The wavelet maps show characteristic oscillations of the frequencies, mainly CC stretching and CCC bending localized on the central 4-ring moiety. Moreover, we have identified the main frequency (methyl wagging motion) involved in the modulation of these oscillations. The anharmonic coupling within a group of vibrational modes was therefore highlighted, in good agreement with experimental evidence. For the first time, a quantitative analysis of time resolved signals from a wavelet transform/ab initio molecular dynamics approach was performed.

Iron(III) Complexes for Highly Efficient and Sustainable Ketalization of Glycerol: A Combined Experimental and Theoretical Study

The growing production of biodiesel as a promising alternative and renewable fuel led as the main problem the dramatic increase of its by-product: glycerol. Different strategies for glycerol derivatization have been reported so far, some more efficient or sustainable than others. Herein, we report a very promising and eco-friendly transformation of glycerol in nontoxic solvents and chemicals (i.e., solketal, ketals), proposing three new families of Fe(III) compounds capable of catalysing glycerol acetalization with unpublished turn over frequencies (TOFs), and adhering most of the principles of green chemistry. The comparison between the activity of complexes of formula [FeCl₃(1-R)] (1-R = substituted pyridinimine), [FeCl(2-R,R')] (2-R,R' = substituted O,O'-deprotonated salens) and their corresponding simple salts reveals that the former are extremely convenient because they are able to promote solketal formation with excellent TOFs, up to 10⁵ h^-1. Satisfactory performances were shown with respect to the entire range of substrates, with results being competitive to those reported in the literature so far. Moreover, the experimental activity was supported by an accurate and complete ab initio study, which disclosed the fundamental role of iron(III) as Lewis acid in promoting the catalytic activity. The unprecedented high activity and the low loading of the catalyst, combined with the great availability and the good eco-toxicological profile of iron, foster future applications of this catalytic process for the sustainable transformation of an abundant by-product in a variety of chemicals.