QM and QM/MM Studies on Excited-State Relaxation Mechanisms of Unnatural Bases in Vacuo and Base Pairs in DNA

Photocyclization and Photoisomerization Mechanisms of an Indolylfulgide Derivative in Acetonitrile Solution by Using the QM (MS-CASPT2)/MM Method

The indolylfulgide systems have been extensively investigated due to their potential applications as photochromic materials. In this work, the photoinduced ring-closure/opening and isomerization reactions of a photochromic indolylfulgide in vacuum and acetonitrile solvent have been investigated by means of MS-CASPT2//CASSCF and QM(MS-CASPT2)//CASSCF/MM. The deactivation mechanisms of indolylfulgide have been proposed based on the optimized structures in the S0 and S1 states, S1/S0 conical intersections, and the calculated minimum-energy paths. After excitation into the first singlet excited-state, which is spectroscopically bright in the Franck-Condon point of the E, the photoprocesses proceed toward a nearby S1 minimum. Then, two possible nonadiabatic relaxation paths exist to repopulate the ground state. In the ring closure reaction, the S1 E isomer evolves directly into one S1/S0 conical intersection and decays to the ground state with bifurcation toward C or E. In the E → Z tautomerization pathway, the excited system can deactivate to the S0 state via a distinct conical intersection. The minimum-energy paths of the indolylfulgide revealed that the ring closure reaction in the solvent is more facile to take place than the E → Z isomerization after irradiation of the same E. Furthermore, for the ring opening reaction from the C side, there exists an energy barrier (11.1 kcal/mol) in the S1 state before arriving at the conical intersection. The computational results showed that the solvent has some influence on the system compared with that in the gas phase. The present work could contribute to comprehending the photoreactions of indolylfulgide and its derivatives in solution.

"On-the-Fly" Nonadiabatic Dynamics Simulation on the Ultrafast Photoisomerization of a Molecular Photoswitch Iminothioindoxyl: An RMS-CASPT2 Investigation

Iminothioindoxyl (ITI) is a new class of photoswitch that exhibits many excellent properties including well-separated absorption bands in the visible region for both conformers, ultrafast Z to E photoisomerization as well as the millisecond reisomerization at room temperature for the E isomer, and switchable ability in both solids and various solvents. However, the underlying ultrafast photoisomerization mechanism at the atomic level remains unclear. In this work, we have employed a combination of high-level RMS-CASPT2-based static electronic structure calculations and nonadiabatic dynamics simulations to investigate the ultrafast photoisomerization dynamics of ITI. Based on the minimum-energy structures, minimum-energy conical intersections, linear interpolation internal coordinate paths, and nonadiabatic dynamics simulations, the overall photoisomerization scenario of ITI upon excitation is established. Upon excitation around 416 nm, the molecule will be excited to the S2 state considering its close energy to the experimentally measured absorption maximum and larger oscillator strength, from which ultrafast decay of S2 to S1 state can take place efficiently with a time constant of 62 fs. However, the photoisomerization is not likely to complete in the S2 state since the dihedral associated with the Z to E isomerization changes little during the relaxation. Upon relaxing to the S1 state, the molecule will decay to the S0 state ultrafast with a time constant of 232 fs. In contrast, the decay of the S1 state is important for the isomerization considering that the dihedral related to the isomerization of the hopping structures is close to 90°. Therefore, the S1/S0 intersection region should be important for the isomerization of ITI. Arriving at the S0 state, the molecule can either go back to the original Z reactant or isomerize to the E products. At the end of the 500 fs simulation time, the E configuration accounts for nearly 37% of the final structures. Moreover, the photoisomerization mechanism is different from the isomerization mechanism in the ground state; i.e., instead of the inversion mechanism in the ground state, the photoisomerization prefers the rotation mechanism. Our results not only agree well with previous experimental studies but also provide some novel insights that could be helpful for future improvements in the performance of the ITI photoswitches.

Theoretical insights into the photophysics of an unnatural base Z: A MS-CASPT2 investigation

We have employed the highly accurate multistate complete active space second-order perturbation theory (MS-CASPT2) method to investigate the photoinduced excited state relaxation properties of one unnatural base, namely Z. Upon excitation to the S2 state of Z, the internal conversion to the S1 state would be dominant. From the S1 state, two intersystem crossing paths leading to the T2 and T1 states and one internal conversion path to the S0 state are possible. However, considering the large barrier to access the S1 /S0 conical intersection and the strong spin-orbit coupling between S1 and T2 states (>40 cm-1 ), the intersystem crossing to the triplet manifolds is predicted to be more preferred. Arriving at the T2 state, the internal conversion to the T1 state and the intersystem crossing back to the S1 state are both possible considering the S1 /T2 /T1 three-state intersection near the T2 minimum. Upon arrival at the T1 state, the deactivation to S0 can be efficient after overcoming a small barrier to access T1 /S0 crossing point, where the spin-orbit coupling (SOC) is as large as 39.7 cm-1 . Our present work not only provides in-depth insights into the photoinduced process of unnatural base Z, but can also help the future design of novel unnatural bases with better photostability.

Investigation of dynamical flexibility of D5SIC-DNAM inside DNA duplex in aqueous solution: a systematic classical MD approach

Incorporation of artificial 3rd base pairs (unnatural base pairs, UBPs) has emerged as a fundamental technique in pursuit of expanding the genetic alphabet. 2,6-Dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN), a potential unnatural base pair (UBP) developed by Romesberg and colleagues, has been shown to have remarkable capability for replication within DNA. Crystal structures of a Taq polymerase/double-stranded DNA (ds-DNA) complex containing a DS-DN pair in the 3' terminus showed a parallelly stacked geometry for the pre-insertion, and an intercalated geometry for the post-insertion structure. Unconventional orientations of DS-DN inside a DNA duplex have inspired scientists to investigate the conformational orientations and structural properties of UBP-incorporated DNA. In recent years, computational simulations have been used to investigate the geometry of DS-DN within the DNA duplex; nevertheless, unresolved questions persist owing to inconclusive findings. In this work, we investigate the structural and dynamical properties of DS and DN inside a ds-DNA strand in aqueous solution considering both short and long DNA templates using polarizable, and non-polarizable classical MD simulations. Flexible conformational change of UBP with major populations of Watson-Crick-Franklin (WCF) and three distinct non-Watson-Crick-Franklin (nWCFP1, nWCFP2, nWCFO) conformations through intra and inter-strand flipping have been observed. Our results suggest that a dynamical conformational change leads to the production of diffierent conformational distribution for the systems. Simulations with a short ds-DNA duplex suggest nWCF (P1 and O) as the predominant structures, whereas long ds-DNA duplex simulations indicate almost equal populations of WCF, nWCFP1, nWCFO. DS-DN in the terminal position is found to be more flexible with occasional mispairing and fraying. Overall, these results suggest flexibility and dynamical conformational change of the UBP as well as indicate varied conformational distribution irrespective of starting orientation of the UBP and length og DNA strand.

Combined QM (MS-CASPT2)/MM studies on photocyclization and photoisomerization of a fulgide derivative in toluene solution

Photocyclization and photoisomerization of fulgides have been extensively studied experimentally and computationally due to their significant potential applications for example as photoswitches in memory devices. However, the reported excited-state decay mechanisms of fulgides do not include the effects of solvation explicitly to date. Herein, calculations using the high-level MS-CASPT2//CASSCF method were conducted to explore the photoinduced excited-state decay processes of the E_α conformer of a fulgide derivative in toluene with solvent effects treated by implicit PCM and explicit QM/MM models, respectively. Several minima and conical intersections were optimized successfully in and between the S₀ and S₁ states; then, two nonadiabatic excited-state decay channels that could efficiently drive the system to the ground state were proposed based on the excited-state ring-closure and isomerization paths. In addition, we also found that in the ring-closure path, the potential energy surface is essentially barrierless before approaching the conical intersection, while it needs to overcome a small energy barrier along the E → Z photoisomerization path for the nonadiabatic S₁ → S₀ internal conversion process. The present computational results could provide useful mechanistic insights into the photoinduced cyclization and isomerization reactions of fulgide and its derivatives.

Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations

Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.

Enzymatic construction of metal-mediated nucleic acid base pairs

Artificial metal base pairs have become increasingly important in nucleic acids chemistry due to their high thermal stability, water solubility, orthogonality to natural base pairs, and low cost of production. These interesting properties combined with ease of chemical and enzymatic synthesis have prompted their use in several practical applications, including the construction of nanomolecular devices, ions sensors, and metal nanowires. Chemical synthesis of metal base pairs is highly efficient and enables the rapid screening of novel metal base pair candidates. However, chemical synthesis is limited to rather short oligonucleotides and requires rather important synthetic efforts. Herein, we discuss recent progress made for the enzymatic construction of metal base pairs that can alleviate some of these limitations. First, we highlight the possibility of generating metal base pairs using canonical nucleotides and then describe how modified nucleotides can be used in this context. We also provide a description of the main analytical techniques used for the analysis of the nature and the formation of metal base pairs together with relevant examples of their applications.

Quantum Mechanics/Molecular Mechanics Studies on the Photophysical Mechanism of Methyl Salicylate

Methyl salicylate (MS) as a subunit of larger salicylates found in commercial sunscreens has been shown to exhibit keto-enol tautomerization and dual fluorescence emission via excited-state intramolecular proton transfer (ESIPT) after the absorption of ultraviolet (UV) radiation. However, its excited-state relaxation mechanism is unclear. Herein, we have employed the quantum mechanics(CASPT2//CASSCF)/molecular mechanics method to explore the ESIPT and excited-state relaxation mechanism of MS in the lowest three electronic states, that is, S₀, S₁, and T₁ states, in a methanol solution. Based on the optimized geometric and electronic structures, conical intersections and crossing points, and minimum-energy paths combined with the computed linearly interpolated Cartesian coordinate paths, the photophysical mechanism of MS has been proposed. The S₁ state is a spectroscopically bright ¹ππ* state in the Franck-Condon region. From the initially populated S₁ state, there exist three nonradiative relaxation paths to repopulate the S₀ state. In the first one, the S₁ system (i.e., ketoB form) first undergoes an ESIPT path to generate an S₁ tautomer (i.e., enol form) that exhibits a large Stokes shift in experiments. The generated S₁ enol tautomer further evolves toward the nearby S₁/S₀ conical intersection and then hops to the S₀ state, followed by the backward ground-state intramolecular proton transfer (GSIPT) to the initial ketoB form S₀ state. In the second one, the S₁ system first hops through the S₁ → T₁ intersystem crossing (ISC) to the T₁ state, which then further decays to the S₀ state via T₁ → S₀ ISC at the T₁/S₀ crossing point. In the third path, the T₁ system that stems from the S₁ → T₁ ISC process via the S₁/T₁ crossing point first takes place a T₁ ESIPT to generate a T₁ enol tautomer, which can further decay to the S₀ state via T₁-to-S₀ ISC. Finally, the GSIPT occurs to back the system to the initial ketoB form S₀ state. Our present work could contribute to understanding the photophysics of MS and its derivatives.

Excited-State Properties and Relaxation Pathways of Selenium-Substituted Guanine Nucleobase in Aqueous Solution and DNA Duplex

The excited-state properties and relaxation mechanisms after light irradiation of 6-selenoguanine (6SeG) in water and in DNA have been investigated using a quantum mechanics/molecular mechanics (QM/MM) approach with the multistate complete active space second-order perturbation theory (MS-CASPT2) method. In both environments, the S₁ ¹(n_Seπ₅^*) and S₂ ¹(π_Seπ₅^*) states are predicted to be the spectroscopically dark and bright states, respectively. Two triplet states, T₁ ³(π_Seπ₅^*) and T₂ ³(n_Seπ₅^*), are found energetically below the S₂ state. Extending the QM region to include the 6SeG-Cyt base pair slightly stabilizes the S₂ state and destabilizes the S₁, due to hydrogen-bonding interactions, but it does not affect the order of the states. The optimized minima, conical intersections, and singlet-triplet crossings are very similar in water and in DNA, so that the same general mechanism is found. Additionally, for each excited state geometry optimization in DNA, three kind of structures ("up", "down", and "central") are optimized which differ from each other by the orientation of the C═Se group with respect to the surrounding guanine and thymine nucleobases. After irradiation to the S₂ state, 6SeG evolves to the S₂ minimum, near to a S₂/S₁ conical intersection that allows for internal conversion to the S₁ state. Linear interpolation in internal coordinates indicate that the "central" orientation is less favorable since extra energy is needed to surmount the high barrier in order to reach the S₂/S₁ conical intersection. From the S₁ state, 6SeG can further decay to the T₁ ³(π_Seπ₅^*) state via intersystem crossing, where it will be trapped due to the existence of a sizable energy barrier between the T₁ minimum and the T₁/S₀ crossing point. Although this general S₂ → T₁ mechanism takes place in both media, the presence of DNA induces a steeper S₂ potential energy surface, that it is expected to accelerate the S₂ → S₁ internal conversion.

Enzymatic Formation of an Artificial Base Pair Using a Modified Purine Nucleoside Triphosphate

The expansion of the genetic alphabet with additional, unnatural base pairs (UBPs) is an important and long-standing goal in synthetic biology. Nucleotides acting as ligands for the coordination of metal cations have advanced as promising candidates for such an expansion of the genetic alphabet. However, the inclusion of artificial metal base pairs in nucleic acids mainly relies on solid-phase synthesis approaches, and very little is known about polymerase-mediated synthesis. Herein, we report the selective and high yielding enzymatic construction of a silver-mediated base pair (dIm^C-Ag^I-dPur^P) as well as a two-step protocol for the synthesis of DNA duplexes containing such an artificial metal base pair. Guided by DFT calculations, we also shed light into the mechanism of formation of this artificial base pair as well as into the structural and energetic preferences. The enzymatic synthesis of the dIm^C-Ag^I-dPur^P artificial metal base pair provides valuable insights for the design of future, more potent systems aiming at expanding the genetic alphabet.

Enzymatic Construction of Artificial Base Pairs: The Effect of Metal Shielding

Th formation of metal base pairs is a versatile method for the introduction of metal cations into nucleic acids that has been used in numerous applications including the construction of metal nanowires, development of energy, charge-transfer devices and expansion of the genetic alphabet. As an alternative, enzymatic construction of metal base pairs is an alluring strategy that grants access to longer sequences and offers the possibility of using such unnatural base pairs (UBPs) in SELEX experiments for the identification of functional nucleic acids. This method remains rather underexplored, and a better understanding of the key parameters in the design of efficient nucleotides is required. We have investigated the effect of methylation of the imidazole nucleoside (dIm^nMe TP) on the efficiency of the enzymatic construction of metal base pairs. The presence of methyl substituents on dImTP facilitates the polymerase-driven formation of dIm^4Me -Ag^I -dIm and dIm^2Me TP-Cr^III -dIm base pairs. Steric factors rather than the basicity of the imidazole nucleobase appear to govern the enzymatic formation of such metal base pairs. We also demonstrate the compatibility of other metal cations rarely considered in the construction of artificial metal bases by enzymatic DNA synthesis under both primer extension reaction and PCR conditions. These findings open up new directions for the design of nucleotide analogues for the development of metal base pairs.

QM/MM Studies on the Photophysical Mechanism of a Truncated Octocrylene Model

Methyl 2-cyano-3,3-diphenylacrylate (MCDPA) shares the same molecular skeleton with octocrylene (OCR) that is one of the most common molecules used in commercially available sunscreens. However, its excited-state relaxation mechanism is unclear. Herein, we have used the QM(CASPT2//CASSCF)/MM method to explore spectroscopic properties, geometric and electronic structures, relevant conical intersections and crossing points, and excited-state relaxation paths of MCDPA in methanol solution. We found that in the Franck-Condon (FC) region, the V(¹ππ*) state is energetically lower than the V'(¹ππ*) state only by 2.8 kcal/mol and is assigned to experimentally observed maximum absorption band. From these two initially populated singlet states, there exist three nonradiative relaxation paths to repopulate the S₀ state. In the first one, when the V(¹ππ*) state is populated in the FC region, the system diabatically evolves along the V(¹ππ*) state into its minimum where the internal conversion to S₀ occurs. In the second one, the V'(¹ππ*) state is populated in the FC region and the system adiabatically overcomes a barrier of ca. 3.0 kcal/mol to approach the V(¹ππ*) minimum eventually leading to a V(¹ππ*)-to-S₀ internal conversion. In the third one, the V'(¹ππ*) state first hops via the intersystem crossing to the T₂ state, which then decays through the internal conversion to the T₁ state. The T₁ state is finally converted to the S₀ state via the T₁/S₀ crossing point. Our present work contributes to understanding the photophysics of OCR and its variants.

Optimization of Replication, Transcription, and Translation in a Semi-Synthetic Organism

Previously, we reported the creation of a semi-synthetic organism (SSO) that stores and retrieves increased information by virtue of stably maintaining an unnatural base pair (UBP) in its DNA, transcribing the corresponding unnatural nucleotides into the codons and anticodons of mRNAs and tRNAs, and then using them to produce proteins containing noncanonical amino acids (ncAAs). Here we report a systematic extension of the effort to optimize the SSO by exploring a variety of deoxy- and ribonucleotide analogues. Importantly, this includes the first in vivo structure-activity relationship (SAR) analysis of unnatural ribonucleoside triphosphates. Similarities and differences between how DNA and RNA polymerases recognize the unnatural nucleotides were observed, and remarkably, we found that a wide variety of unnatural ribonucleotides can be efficiently transcribed into RNA and then productively and selectively paired at the ribosome to mediate the synthesis of proteins with ncAAs. The results extend previous studies, demonstrating that nucleotides bearing no significant structural or functional homology to the natural nucleotides can be efficiently and selectively paired during replication, to include each step of the entire process of information storage and retrieval. From a practical perspective, the results identify the most optimal UBP for replication and transcription, as well as the most optimal unnatural ribonucleoside triphosphates for transcription and translation. The optimized SSO is now, for the first time, able to efficiently produce proteins containing multiple, proximal ncAAs.

Towards the enzymatic formation of artificial metal base pairs with a carboxy-imidazole-modified nucleotide

The identification of synthetic nucleotides that sustain the formation of orthogonal, unnatural base pairs is an important goal in synthetic biology. Such artificial synthons have been used for the generation of semi-synthetic organisms as well as functional nucleic acids with enhanced binding properties. The enzymatic formation of artificial metal-base pairs is a vastly underexplored and alluring alternative to existing systems. Here, we report the synthesis and biochemical characterization of 1‑(2-deoxy‑β‑d‑ribofuranosyl) imidazole‑4‑carboxylate nucleoside triphosphate (dIm^CTP) which is equipped with a carboxylic acid moiety on the imidazole moiety in order to increase the coordination environment to [2 + 2] and [2 + 1]. A clear metal dependence was observed for the single incorporation of the modified nucleotide into DNA by the DNA polymerase from Thermus aquaticus (Taq). The presence of Ag^I in primer extension reactions conducted with combinations of 1‑(2‑deoxy‑β‑d‑ribofuranosyl) imidazole nucleoside triphosphate (dImTP) and dIm^CTP supported the unusual [2 + 1] coordination pattern. The efficiency of the tailing reactions mediated by the terminal deoxynucleotidyl transferase (TdT) was markedly improved when using dIm^CTP instead of dImTP. Even though products with multiple modified nucleotides were not observed, the appendage of additional metal binding ligands on the imidazole nucleobase appears to be a valid approach to improve the biochemical properties of modified triphosphates in the context of an expansion of the genetic alphabet with metal base pairs.

Progress Toward a Semi-Synthetic Organism with an Unrestricted Expanded Genetic Alphabet

We have developed a family of unnatural base pairs (UBPs), exemplified by the pair formed between dNaM and dTPT3, for which pairing is mediated not by complementary hydrogen bonding but by hydrophobic and packing forces. These UBPs enabled the creation of the first semisynthetic organisms (SSOs) that store increased genetic information and use it to produce proteins containing noncanonical amino acids. However, retention of the UBPs was poor in some sequence contexts. Here, to optimize the SSO, we synthesize two novel benzothiophene-based dNaM analogs, dPTMO and dMTMO, and characterize the corresponding UBPs, dPTMO-dTPT3 and dMTMO-dTPT3. We demonstrate that these UBPs perform similarly to, or slightly worse than, dNaM-dTPT3 in vitro. However, in the in vivo environment of an SSO, retention of dMTMO-dTPT3, and especially dPTMO-dTPT3, is significantly higher than that of dNaM-dTPT3. This more optimal in vivo retention results from better replication, as opposed to more efficient import of the requisite unnatural nucleoside triphosphates. Modeling studies suggest that the more optimal replication results from specific internucleobase interactions mediated by the thiophene sulfur atoms. Finally, we show that dMTMO and dPTMO efficiently template the transcription of RNA containing TPT3 and that their improved retention in DNA results in more efficient production of proteins with noncanonical amino acids. This is the first instance of using performance within the SSO as part of the UBP evaluation and optimization process. From a general perspective, the results demonstrate the importance of evaluating synthetic biology "parts" in their in vivo context and further demonstrate the ability of hydrophobic and packing interactions to replace the complementary hydrogen bonding that underlies the replication of natural base pairs. From a more practical perspective, the identification of dMTMO-dTPT3 and especially dPTMO-dTPT3 represents significant progress toward the development of SSOs with an unrestricted ability to store and retrieve increased information.

Photochemical and Photodynamical Properties of Sulfur-Substituted Nucleic Acid Bases

Sulfur-substituted nucleobases (a.k.a., thiobases) are among the world's leading prescriptions for chemotherapy and immunosuppression. Long-term treatment with azathioprine, 6-mercaptopurine and 6-thioguanine has been correlated with the photoinduced formation of carcinomas. Establishing an in-depth understanding of the photochemical properties of these prodrugs may provide a route to overcoming these carcinogenic side effects, or, alternatively, a basis for developing effective compounds for targeted phototherapy. In this review, a broad examination is undertaken, surveying the basic photochemical properties and excited-state dynamics of sulfur-substituted analogs of the canonical DNA and RNA nucleobases. A molecular-level understanding of how sulfur substitution so remarkably perturbs the photochemical properties of the nucleobases is presented by combining experimental results with quantum-chemical calculations. Structure-property relationships demonstrate the impact of site-specific sulfur substitution on the photochemical properties, particularly on the population of the reactive triplet state. The value of fundamental photochemical investigations for driving the development of ultraviolet-A chemotherapeutics is showcased. The most promising photodynamic agents identified thus far have been investigated in various carcinoma cell lines and shown to decrease cell proliferation upon exposure to ultraviolet-A radiation. Overarching principles have been elucidated for the impact that sulfur substitution of the carbonyl oxygen has on the photochemical properties of the nucleobases.