A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution

Computational Spectroscopy of Aqueous Solutions: The Underlying Role of Conformational Sampling

Fully atomistic multiscale polarizable quantum mechanics (QM)/molecular mechanics (MM) approaches, combined with techniques to sample the solute-solvent phase space, constitute the most accurate method to compute spectral signals in aqueous solution. Conventional sampling strategies, such as classical molecular dynamics (MD), may encounter drawbacks when the conformational space is particularly complex, and transition barriers between conformers are high. This can lead to inaccurate sampling, which can potentially impact the accuracy of spectral calculations. For this reason, in this work, we compare classical MD with enhanced sampling techniques, i.e., replica exchange MD and metadynamics. In particular, we show how the different sampling techniques affect computed UV, electronic circular dichroism, nuclear magnetic resonance shielding, and optical rotatory dispersion of N-acetylproline-amide in aqueous solution. Such a system is a model peptide characterized by complex conformational variability. Calculated values suggest that spectral properties are influenced by solute conformers, relative population, and solvent effects; therefore, particular care needs to be paid for when choosing the sampling technique.

Substituent Effect versus Aromaticity─A Curious Case of Fulvene Derivatives

A computational study on amino- and nitro-substituted penta- and heptafulvenes reveals the interplay between the aromaticity and the substituent effect (SE). Ring substitution alone has little influence on the aromaticity, but in combination with an exo substituent of opposite properties, it substantially enhances the cyclic π-electron delocalization. Despite the SE being stronger for β substitution, only γ substitution leads to higher aromaticity. An explanation is provided by the electron density of delocalized bonds (EDDB) method, which proves to be a valuable tool in analyzing both cyclic delocalization and the SE.

Restriction of photoinduced electron transfer as a mechanism for the aggregation-induced emission of a trityl-functionalised maleimide fluorophore

The restriction of intramolecular motion (RIM) and restricted access to a conical intersection (RACI) have been accepted as general working mechanisms for aggregation-induced emission (AIE) phenomena. However, as the family of AIE molecules grows, the RIM and RACl mechanisms cannot be used to fully understand some AIE phenomena. Herein, the restriction of the photoinduced electron transfer (RPET) state is proposed to rationalize the AIE phenomena of trityl-functionalised maleimide molecule based on density functional theory calculations. The "state-crossing from a locally excited to an electron transfer state" (SLEET) model was employed to predict the ON/OFF molecular PET in solution and solid states. According to the SLEET model, we showed that a non-emissive electron transfer excited state leads to the fluorescence quenching of trityl-functionalised maleimide in solution. However, due to the reduced polarity of the environment in aggregates, the electron transfer state is thermodynamically inaccessible, and a low-lying locally excited state exhibits intense emission. These findings provide a theoretical foundation to understand the working mechanisms of AIE molecules and the design of new AIEgens. We expect that the RPET mechanism can be used to screen potential AIEgens using the SLEET model.

Unraveling the acid-base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand

Understanding the interactions and the solvent effects on the distribution of several species in equilibrium and how it can influence the ¹H-NMR properties, spectroscopy (UV-vis absorption), and the acid-base equilibria can be especially challenging. This is the case of a bis-bidentate bridging ligand bis(2-pyridyl)-benzo-bis(imidazole), where the two pyridyl and four imidazolyl nitrogen atoms can be protonated in different ways, depending on the solvent, generating many isomeric/tautomeric species. Herein, we report a combined theoretical-experimental approach based on a sequential quantum mechanics/molecular mechanics procedure that was successfully applied to describe in detail the acid-base characterization and its effects on the electronic properties of such a molecule in solution. The calculated free-energies allowed the identification of the main species present in solution as a function of the solvent polarity, and its effects on the magnetic shielding of protons (¹H-NMR chemical shifts), the UV-vis absorption spectra, and the acid-base equilibrium constants (pK_as) in aqueous solution. Three acid-base equilibrium constants were experimentally/theoretically determined (pK_a1 = 1.3/1.2, pK_a2 = 2.1/2.2 and pK_a5 = 10.1/11.3) involving mono-deprotonated and mono-protonated cis and trans species. Interestingly, other processes with pK_a3 = 3.7 and pK_a4 = 6.0 were also experimentally determined and assigned to the protonation and deprotonation of dimeric species. The dimerization of the most stable neutral species was investigated by Monte Carlo simulations and its electronic effects were considered for the elucidation of the UV-vis absorption bands, revealing transitions mainly with the charge-transfer characteristic and involving both the monomeric species and the dimeric species. The good matching of the theoretical and experimental results provides an atomistic insight into the solvent effects on the electronic properties of this bis-bidentate bridging ligand.

Spectrophotometric studies of charge-transfer complexes formed with ions N,N'-alkyldiyl-bis(pyridinium) derivatives and iodide

The charge-transfer complexes (CTC) between electron donor iodide and pyridinium dimers and monomers as acceptors have been studied spectrophotometrically in acetonitrile, DMSO and water. The dimers were: N,N'-alkyldiyl-bis(4-cyanopyridinium) and N,N'-alkyldiyl-bis(2-bromopyridinium), and the monomers were: N-alkyl-4-cyanopyridinium and N-alkyl-2-bromopyridinium, bridged by n methylene spacers. The formation constant (K_CTC), molar absorptivity (ε_CTC), fluorescence-quenching constant (K_SV) were determined. The results indicate that the stoichiometry of the CTC for dimers and monomers is 1:1 (equimolar ratio), and its formation is favored for dimers more than for monomers, especially dimers with short methylene spacers forming a "sandwiched type-complex" in which the iodide is close to the two pyridinium rings. Solvents with low polarity favored the complex, which was destroyed by the presence of water. The experimental studies were complemented with theoretical studies with quantum mechanics (QM) calculations using density functional theory (DFT) and molecular mechanics (MM) simulations with Molecular Dynamics for identify the most stable conformers in acetonitrile solution. The electronic excitations were calculated with sequential QM/MM hybrid method, showing a charge-transfer wavelength in agreement with the UV-Vis absorption spectra obtained experimentally. It confirms the "sandwiched type-complex" conformers favoring the interaction of the iodide with one pyridinium rings and simultaneously forming one, or two, hydrogen bonds with the alkyl chain and additionally, for the case of dimers, with the other pyridinium ring.

Free Energy Computation for an Isomerizing Chromophore in a Molecular Cavity via the Average Solvent Electrostatic Configuration Model: Applications in Rhodopsin and Rhodopsin-Mimicking Systems

We present a novel technique for computing the free energy differences between two chromophore "isomers" hosted in a molecular environment (a generalized solvent). Such an environment may range from a relatively rigid protein cavity to a flexible solvent environment. The technique is characterized by the application of the previously reported "average electrostatic solvent configuration" method, and it is based on the idea of using the free energy perturbation theory along with a chromophore annihilation procedure in thermodynamic cycle calculations. The method is benchmarked by computing the ground-state room-temperature relative stabilities between (i) the cis and trans isomers of prototypal animal and microbial rhodopsins and (ii) the analogue isomers of a rhodopsin-like light-driven molecular switch in methanol. Furthermore, we show that the same technology can be used to estimate the activation free energy for the thermal isomerization of systems i-ii by replacing one isomer with a transition state. The results show that the computed relative stability and isomerization barrier magnitudes for the selected systems are in line with the available experimental observation in spite of their widely diverse complexity.