Mechanistic analysis of light-driven overcrowded alkene-based molecular motors by multiscale molecular simulations

Computational Study on the Dynamics of a Bis(benzoxazole)-Based Overcrowded Alkene

Understanding and controlling molecular motions is of pivotal importance for designing molecular machinery and functional molecular systems, capable of performing complex tasks. Herein, we report a comprehensive theoretical study to elucidate the dynamic behavior of a bis(benzoxazole)-based overcrowded alkene displaying several coupled and uncoupled molecular motions. The benzoxazole moieties give rise to 4 different stable conformers that interconvert through single-bond rotations. By performing excited- and ground-state molecular dynamics simulations, DFT calculations, and NMR studies, we found that the photochemical E-Z isomerization of the central double bond of each stable conformer is directional and leads to a mixture of metastable isomers. This transformation is analogous to the classical Feringa-type molecular motors, with the notable difference that, during the photochemical isomerization and the subsequent thermal helix inversion (THI) steps, multiple possible pathways take place that involve single-bond rotations that can be both coupled and uncoupled to the rotation of the naphthyl half of the molecule.

On the supramolecular interactions into a pH- and Metal-Actuated Molecular Shuttle: some insights from QTAIM modeling

Supramolecular contacts responsible for chemical interaction of cucurbit[7]uril (CB[7]) macrocycle on a Tolyl-Viologen-Phenylene-Imidazole (T-VPI) molecular thread, at acid pH (T-VPI-H+) or after Ag+ cation addition (T-VPI-Ag+), are analytically addressed in a computational framework combining Quantum Theory of Atoms in Molecules (QTAIM) with Density Functional Theory (DFT). In this respect, the crystallographic structure (CCDC number 2217466) is taken as reference condition for addressing the nature of the chemical interactions driving the shuttling of the CB[7] between T and P stations recently observed in dilute water solutions. Beside the host(CB[7]) vs guest(T-VPI-H+ or T-VPI-Ag+) complexation, the coordination sphere of the Ag+ cation is also investigated by means of local electronic energy density - H(r) - descriptors. The derived non-covalent interaction patterns are found to support diagnostic 1H NMR signals used for detecting the mutual position of the CB[7] along the axle. This work highlights the potentialities of a QTAIM based approach in the characterization of supramolecular and metal-complexation effects in molecular aggregates such as not-interlocked synthetic molecular shuttles.

Excited State Dynamics in Unidirectional Photochemical Molecular Motors

Unidirectional photochemically driven molecular motors (PMMs) convert the energy of absorbed light into continuous rotational motion. As such they are key components in the design of molecular machines. The prototypical and most widely employed class of PMMs is the overcrowded alkenes, where rotational motion is driven by successive photoisomerization and thermal helix inversion steps. The efficiency of such PMMs depends upon the speed of rotation, determined by the rate of ground state thermal helix inversion, and the quantum yield of photoisomerization, which is dependent on the excited state energy landscape. The former has been optimized by synthetic modification across three generations of overcrowded alkene PMMs. These improvements have often been at the expense of photoisomerization yield, where there remains room for improvement. In this perspective we review the application of ultrafast spectroscopy to characterize the excited state dynamics in PMMs. These measurements lead to a general mechanism for all generations of PMMs, involving subpicosecond decay of a Franck-Condon excited state to populate a dark excited state which decays within picoseconds via conical intersections with the electronic ground state. The model is discussed in the context of excited state dynamics calculations. Studies of PMM photochemical dynamics as a function of solvent suggest exploitation of intramolecular charge transfer and solvent polarity as a route to controlling photoisomerization yield. A test of these ideas for a first generation motor reveals a high degree of solvent control over isomerization yield. These results suggest a pathway to fine control over the performance of future PMMs.

Excited-State Dynamics Simulations of a Light-Driven Molecular Motor in Solution

Molecular motors, where light can be transformed into motion, are promising in the design of nanomechanical devices. For applications, however, finding relationships between molecular motion and the environment is important. Here, we report the study of excited-state dynamics of an overcrowded alkene in solution using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach combined with excited-state molecular dynamics simulations. Using QM/MM surface-hopping trajectories, we calculated time-resolved emission and transient absorption spectra. These show the rise of a short-lived Franck-Condon state, followed by the formation of a dark state in the first 150 fs before the molecular motor relaxes to the ground state in about 1 ps. From the analysis of radial distribution functions, we infer that the orientation of the solvent with respect to the molecular motor in the electronic excited state is similar to that in the ground state during the photoisomerization.

Absolute binding free energy calculations improve enrichment of actives in virtual compound screening

We determined the effectiveness of absolute binding free energy (ABFE) calculations to refine the selection of active compounds in virtual compound screening, a setting where the more commonly used relative binding free energy approach is not readily applicable. To do this, we conducted baseline docking calculations of structurally diverse compounds in the DUD-E database for three targets, BACE1, CDK2 and thrombin, followed by ABFE calculations for compounds with high docking scores. The docking calculations alone achieved solid enrichment of active compounds over decoys. Encouragingly, the ABFE calculations then improved on this baseline. Analysis of the results emphasizes the importance of establishing high quality ligand poses as starting points for ABFE calculations, a nontrivial goal when processing a library of diverse compounds without informative co-crystal structures. Overall, our results suggest that ABFE calculations can play a valuable role in the drug discovery process.

Solvation Effects on the Thermal Helix Inversion of Molecular Motors from QM/MM Calculations

Molecular motors convert light and thermal energies into mechanical work, offering good opportunities to design novel molecular devices. Among them, molecular motors alternate a photoisomerization and a thermal helix inversion to achieve unidirectional rotation. The rotational speed is limited by the helix inversion step, which in turn is governed by a barrier in the electronic ground state. In this work, we systematically study the solvation effect on the thermal process of selected molecular motors, comparing reaction barriers obtained from both density functional theory (DFT) in the isolated system and umbrella sampling within a hybrid quantum mechanics/molecular mechanics (QM/MM) model in solution. We find more prominent solvation effects on those molecular motors with larger dipole moments. The results could provide insight into how to functionalize molecular motors to speed up their rotation.

Modeling Molecular Emitters in Organic Light-Emitting Diodes with the Quantum Mechanical Bespoke Force Field

Combined molecular dynamics (MD) and quantum mechanics (QM) simulation procedures have gained popularity in modeling the spectral properties of functional organic molecules. However, the potential energy surfaces used to propagate long-time scale dynamics in these simulations are typically described using general, transferable force fields designed for organic molecules in their electronic ground states. These force fields do not typically include spectroscopic data in their training, and importantly, there is no general protocol for including changes in geometry or intermolecular interactions with the environment that may occur upon electronic excitation. In this work, we show that parameters tailored for thermally activated delayed fluorescence (TADF) emitters used in organic light-emitting diodes (OLEDs), in both their ground and electronically excited states, can be readily derived from a small number of QM calculations using the QUBEKit (QUantum mechanical BEspoke toolKit) software and improve the overall accuracy of these simulations.