Azobenzene as a photoregulator covalently attached to RNA: a quantum mechanics/molecular mechanics-surface hopping dynamics study

Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment

Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.

The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Theoretical Investigation on the Reversible Photoswitch Mechanism of Benzylidene-Oxazolone System

The design and application of molecular photoswitches have attracted much attention. Herein, we performed a detailed computational study on the photoswitch benzylidene-oxazolone system based on static electronic structure calculations and on-the-fly excited-state dynamic simulations. For the Z and E isomer, we located six and four minimum energy conical intersections (MECIs) between the first excited state (S1) and the ground state (S0), respectively. Among them, the relaxation pathway driven by ring-puckering motion is the most competitive channel with the photoisomeization process, leading to the low photoisomerization quantum yield. In the dynamic simulations, about 88 % and 66 % trajectories decay from S1 to S0 for Z and E isomer, respectively, within the total simulation time of ~2 ps. The photoisomeization quantum yields obtained in our study (0.20 for Z→E and 0.12 for E→Z) agree well with the experimental measured values (0.25 and 0.11), even though the number of trajectories is limited to 50. Our study sheds light on the complexity of the benzylidene-oxazolone system 's deactivation process and the competitive mechanisms among different reaction channels, which provides theoretical guidance for further design and development of benzylidene-oxazolone based molecular photoswitches.

Controlled destabilization of caged circularized DNA oligonucleotides predicted by replica exchange molecular dynamics simulations

Spatiotemporal control is a critical issue in the design of strategies for the photoregulation of oligonucleotide activity. Efficient uncaging, i.e., activation by removal of photolabile protecting groups (PPGs), often necessitates multiple PPGs. An alternative approach is based on circularization strategies, exemplified by intrasequential circularization, also denoted photo-tethering, as introduced in [Seyfried et al., Angew. Chem., Int. Ed., 2017, 56, 359]. Here, we develop a computational protocol, relying on replica exchange molecular dynamics (REMD), in order to characterize the destabilization of a series of circularized, caged DNA oligonucleotides addressed in the aforementioned study. For these medium-sized (32 nt) oligonucleotides, melting temperatures are computed, whose trend is in good agreement with experiment, exhibiting a large destabilization and, hence, reduction of the melting temperature of the order of ΔTm ∼ 30 K as compared with the native species. The analysis of free energy landscapes confirms the destabilization pattern experienced by the circularized oligonucleotides. The present study underscores that computational protocols that capture controlled destabilization and uncaging of oligonucleotides are promising as predictive tools in the tailored photocontrol of nucleic acids.

Dual photoisomerization mechanism of azobenzene embedded in a lipid membrane

The photoisomerization of chromophores embedded in biological environments is of high importance for biomedical applications, but it is still challenging to define the photoisomerization mechanism both experimentally and computationally. We present here a computational study of the azobenzene molecule embedded in a DPPC lipid membrane, and assess the photoisomerization mechanism by means of the quantum mechanics/molecular mechanics surface hopping (QM/MM-SH) method. We observe that while the trans-to-cis isomerization is a slow process governed by a torsional mechanism due to the strong interaction with the environment, the cis-to-trans mechanism is completed in sub-ps time scale and is governed by a pedal-like mechanism in which both weaker interactions with the environment and a different geometry of the potential energy surface play a key role.

Photochemical and thermochemical pathways to S2 and polysulfur formation in the atmosphere of Venus

Polysulfur species have been proposed to be the unknown near-UV absorber in the atmosphere of Venus. Recent work argues that photolysis of one of the (SO)₂ isomers, cis-OSSO, directly yields S₂ with a branching ratio of about 10%. If correct, this pathway dominates polysulfur formation by several orders of magnitude, and by addition reactions yields significant quantities of S₃, S₄, and S₈. We report here the results of high-level ab-initio quantum-chemistry computations that demonstrate that S₂ is not a product in cis-OSSO photolysis. Instead, we establish a novel mechanism in which S₂ is formed in a two-step process. Firstly, the intermediate S₂O is produced by the coupling between the S and Cl atmospheric chemistries (in particular, SO reaction with ClS) and in a lesser extension by O-abstraction reactions from cis-OSSO. Secondly, S₂O reacts with SO. This modified chemistry yields S₂ and subsequent polysulfur abundances comparable to the photolytic cis-OSSO mechanism through a more plausible pathway. Ab initio quantification of the photodissociations at play fills a critical data void in current atmospheric models of Venus.

Polysulfur compounds have been ascribed as the unknown near-UV absorbers in Venusian atmosphere and play a key role in the sulfur chemical cycle of this planet. Here, authors establish their production from (SO)₂ on the grounds of quantifications of photochemical and thermal pathways involved in the sulfur chemical cycle of the planet.

Effects of Enzyme-Ligand Interactions on the Photoisomerization of a Light-Regulated Chemotherapeutic Drug

Molecular photoswitches permit using light to control protein activity with high spatiotemporal resolutions, thereby alleviating the side effects of conventional chemotherapy. However, due to the challenges in probing ultrafast photoisomerization reactions in biological environments, it remains elusive how the protein influences the photochemistry of the photoswitches, which hampers the rational design of light-regulated therapeutics. To overcome this challenge, we employed first-principles nonadiabatic dynamics simulations to characterize the photodynamics of the phototrexate (PTX), a recently developed photoswitchable anticancer chemotherapeutic that reversibly inhibits its target enzyme dihydrofolate reductase (DHFR). Our simulations show that the protein environment impedes the trans to cis photoisomerization of the PTX. The confinement in the ligand-binding cavity slows down the isomerization kinetics and quantum yield of the photoswitch by reshaping its conical intersection, increasing its excited-state free-energy barrier and quenching its local density fluctuations. Also, the protein environment results in a suboptimal binding mode of the photoproduct that needs to undergo large structural rearrangement to effectively inhibit the enzyme. Therefore, we predict that the PTX's trans → cis photoisomerization in solution precedes its binding with the protein, despite the favorable binding energy of the trans isomer. Our findings highlight the importance of the protein environment on the photochemical reactions of the molecular photoswitches. As such, our work represents an important step toward the rational design of light-regulated drugs in photopharmacology.

Photoinduced azobenzene-modified DNA dehybridization: insights into local and cooperativity effects from a molecular dynamics study

Photoresponsive azobenzene-modified DNA (RNA) has become a very fruitful material for nanotechnology due to the capability of switching on and off hybridization (i.e., duplex formation) in smart nanostructures. This nanomaterial exploits the well-known azobenzene trans/cis photo-isomerization. In fact, it has been found that DNA tethered with trans-azobenzene shows normal nucleic acid recognition and hybridization, while the cis form destabilizes the duplex configuration, eventually leading to DNA unzipping. However, while the working principle of the light-triggered DNA dehybridization is apparent, specific details of this mechanism still remain elusive to experiments. Previous in silico studies successfully addressed some aspects (e.g., local structural effects, thermal stability, and early events of azobenzene photoisomerization) of this challenging molecular process characterized by timescales spanning several orders of magnitude, from picoseconds (i.e., azobenzene photoisomerization) to micro- and milli-seconds (i.e., complete strand detachment). In this work, inspired by a recent report by Asanuma and coworkers, we focus on the local and cooperativity effects played by multiple azobenzene units on a 10-mer azobenzene-modified DNA duplex. Using molecular dynamics (MD) simulations, we investigated nine systems equipped with a variable number (from 1 to 7) of photoswitch units and different configurations, focusing our analysis on the initial events (from few ps to hundreds of ns) characterizing DNA destabilization upon trans-to-cis isomerization, such as hydrogen bonding breakage and base pair misalignment. Results highlight, on one hand, the local effects of single azobenzene units on DNA duplex structure and, on the other hand, the cooperative role that multiple photoswitches show in enhancing and accelerating DNA dehybridization following trans-to-cis conversion, in agreement with previously reported data and observations.

Photochemistry and Non-adiabatic Photodynamics of the HOSO Radical

Hydroxysulfinyl radical (HOSO) is important due to its involvement in climate geoengineering upon SO₂ injection and generation of the highly hygroscopic H₂SO₄. Its photochemical behavior in the upper atmosphere is, however, uncertain. Here we present the ultraviolet-visible photochemistry and photodynamics of this species by simulating the atmospheric conditions with high-level quantum chemistry methods. Photocleavage to HO + SO arises as the major solar-induced channel, with a minor contribution of H + SO₂ photoproducts. The efficient generation of SO is relevant due to its reactivity with O₃ and the consequent depletion of ozone in the stratosphere.

Theoretical insights into the effect of size and substitution patterns of azobenzene derivatives on the DNA G-quadruplex

Introducing photoswitches into the DNA G-quadruplex provides excellent opportunities to control folding and unfolding of these assemblies, demonstrating their potential in the development of novel nanodevices with medical and nanotechnology applications. Using a quantum mechanics/molecular mechanics (QM/MM) scheme, we carried out a series of simulations to identify the effect of the size and substitution patterns of three azobenzene derivatives (AZ1, AZ2 and AZ3) on the excitation energies of the two lowest excited states of the smallest photoswitchable G-quadruplex reported to date. We demonstrated that the size and the substitution pattern do not affect the ultrafast cis-trans photoiomerization mechanism of the azobenzene derivatives significantly, in agreement with the experiment. However, molecular dynamics simulations revealed that while AZ2 and AZ3 G-quadruplexes are structurally stable during the simulations, the AZ1 G-quadruplex undergoes larger structural changes and shows two ground state populations that differ in the azobenzene backbone adopting two different conformations. AZ1, with para-para substitution pattern, provides more flexibility to the whole G-quadruplex structure compared to AZ2 and AZ3, and can thus facilitate the photoisomerization reaction between a nonpolymorphic, stacked, tetramolecular G-quadruplex and an unstructured state after trans-cis isomerization occurring in a longer time dynamics, in agreement with the experimental findings. The QM/MM simulations of the absorption spectra indicated that the thermal fluctuation plays a more crucial role in the main absorption band of the azobenzene derivatives than the inclusion of the G-quadruplex, implying that the influence of the G-quadruplex environment is minimal. We propose that the latter is attributed to the position of the azobenzene linkers in the G-quadruplexes, i.e. the edgewise loops containing the azobenzene moieties that are located above the G-quartets, not being fully embedded inside or involved in the stacked structure. Our theoretical findings provide support to a recent study of the photoresponsive formation of photoswitchable G-quadruplex motifs.

Spectroscopic Properties of Two 5'-(4-Dimethylamino)Azobenzene Conjugated G-Quadruplex Forming Oligonucleotides

The synthesis of two 5'-end (4-dimethylamino)azobenzene conjugated G-quadruplex forming aptamers, the thrombin binding aptamer (TBA) and the HIV-1 integrase aptamer (T30695), was performed. Their structural behavior was investigated by means of UV, CD, fluorescence spectroscopy, and gel electrophoresis techniques in K+-containing buffers and water-ethanol blends. Particularly, we observed that the presence of the 5'-(4-dimethylamino)azobenzene moiety leads TBA to form multimers instead of the typical monomolecular chair-like G-quadruplex and almost hampers T30695 G-quadruplex monomers to dimerize. Fluorescence studies evidenced that both the conjugated G-quadruplexes possess unique fluorescence features when excited at wavelengths corresponding to the UV absorption of the conjugated moiety. Furthermore, a preliminary investigation of the trans-cis conversion of the dye incorporated at the 5'-end of TBA and T30695 showed that, unlike the free dye, in K+-containing water-ethanol-triethylamine blend the trans-to-cis conversion was almost undetectable by means of a standard UV spectrophotometer.

Finding the true pathway for reversible isomerization of a single azobenzene molecule tumbling on Au(111) surface

Switchable trans-cis isomerization of azobenzene (AB) and its derivatives on metallic surfaces have offered rich possibilities to functionalize molecular devices. However, the lack of a good understanding of the isomerization pathway has severely limited our ability for rational design. One of the long-debated issues is the cis configuration of the parental AB on the Au(111) surface, for which the experimentally inferred structure differs from the theoretically predicted global minimum. Here, we theoretically identify a new in situ metastable configuration for cis-AB on Au(111) that can reproduce all the observations reported in the scanning tunneling microscopy experiments. It reveals that the bistability of AB on the Au(111) surface is attributed to the significantly increased kinetic stability of the newly discovered cis-AB isomer. A fascinating tumbling pathway that overcomes two energy barriers stimulated by tunneling electrons for the trans-cis AB isomerization on Au(111) has been verified, suggesting a new type of molecular motion based on the AB systems.

Aminoazobenzene@Ag modified meshes with large extent photo-response: towards reversible oil/water removal from oil/water mixtures

A large-extent photo-responsive wettability transformation is realized by aminoazobenzene@Ag to selectively remove oil/water.

Photo-responsive materials with superwetting properties, especially in the azo-based class, have been used in water treatment because of their smart performance on wettability changes. However, their transformation extent in wettability has always troubled researchers. Here, we modified nano-Ag pine needles and aminoazobenzene (AABN) on polydopamine (PDA) pre-treated porous meshes, realizing a large-extent reversible photo-responsive wettability transformation from highly hydrophobic to highly hydrophilic. The contact angle is about 150.0° after being exposed to visible light, and is about 10.0° under 365 nm UV light. Accordingly, the modified mesh can achieve photo-responsive removal between oil and water from oil/water mixtures. This facile and universal approach based on trans–cis isomerization of AABN could be endowed to various commercial conductive meshes. Moreover, the modified meshes exhibit satisfactory removal efficiency, reusability and physical/chemical stability, which are more promising for practical applications such as fuel recycling, remote controlled oil/water separation and astronautical resource regeneration.