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 the stability of D5SIC-DNAM-incorporated DNA duplex in Taq polymerase binary system: a systematic classical MD approach

DNA polymerases are fundamental enzymes that play a crucial role in processing DNA with high fidelity and accuracy ensuring the faithful transmission of genetic information. The recognition of unnatural base pairs (UBPs) by polymerases, enabling their replication, represents a significant and groundbreaking discovery with profound implications for genetic expansion. Romesberg et al. examined the impact of DNA containing 2,6-dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN) UBPs bound to T. aquaticus DNA polymerase (Taq) through crystal structure analysis. Here, we have used polarizable and nonpolarizable classical molecular dynamics (MD) simulations to investigate the structural aspects and stability of Taq in complex with a DNA duplex including a DS-DN pair in the terminal 3' and 5' positions. Our results suggest that the flexibility of UBP-incorporated DNA in the terminal position is arrested by the polymerase, thus preventing fraying and mispairing. Our investigation also reveals that the UBP remains in an intercalated conformation inside the active site, exhibiting two distinct orientations in agreement with experimental findings. Our analysis pinpoints particular residues responsible for favorable interactions with the UBP, with some relying on van der Waals interactions while other on Coulombic forces.

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.

N-Methoxy-1,3-oxazinane nucleic acids (MOANAs) - a configurationally flexible backbone modification allows post-synthetic incorporation of base moieties

(2R,3S)-4-(Methoxyamino)butane-1,2,3-triol was converted into a protected phosphoramidite building block and incorporated into the middle of a short DNA oligonucleotide. O1 and O3 of the (2R,3S)-4-(methoxyamino)butane-1,2,3-triol were engaged in phosphodiester linkages, leaving O2 and the methoxyamino function available to form an N-methoxy-1,3-oxazinane ring through reaction with an aldehyde. In modified oligonucleotides thus obtained, the oxazinane ring formally replaces the furanose ring and the aldehyde, the base moiety of natural nucleosides. The feasibility of synthesizing base-modified oligonucleotides by this approach was demonstrated with several aromatic and aliphatic aldehydes featuring various functional groups.