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

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.