Theoretical Study of Light-Induced Crosslinking Reaction Between Pyrimidine DNA Bases and Aromatic Amino Acids

Low-lying electronic excited states and their relaxation pathways as well as energetics of the crosslinking reaction between uracil as a model system for pyrimidine-type building blocks of DNA and RNA and benzene as a model system for aromatic groups of tyrosine (Tyr) and phenylalanine (Phe) amino acids have been studied in the framework of density functional theory. The equilibrium geometries of the ground and electronic excited states as well as the crossing points between the potential energy surfaces of the uracil–benzene complex were computed. Based on these results, different relaxation pathways of the electronic excited states that lead to either back to the initial geometry configuration or the dimerization between the six-membered rings of the uracil–benzene complex have been identified, and the energetic conditions for their occurrence are discussed. It can be concluded that the DNA–protein crosslinking reaction can be induced by the external electromagnetic field via the dimerization reaction between the six-membered rings of the uracil–benzene pair at the electronic excited-state level of the complex. In the case of the uracil–phenol complex, the configuration of the cyclic adduct (dimerized) conformation is less likely to be formed.

All-electron real-time and imaginary-time time-dependent density functional theory within a numeric atom-centered basis function framework

Real-time time-dependent density functional theory (RT-TDDFT) is an attractive tool to model quantum dynamics by real-time propagation without the linear response approximation. Sharing the same technical framework of RT-TDDFT, imaginary-time time-dependent density functional theory (it-TDDFT) is a recently developed robust-convergence ground state method. Presented here are high-precision all-electron RT-TDDFT and it-TDDFT implementations within a numerical atom-centered orbital (NAO) basis function framework in the FHI-aims code. We discuss the theoretical background and technical choices in our implementation. First, RT-TDDFT results are validated against linear-response TDDFT results. Specifically, we analyze the NAO basis sets' convergence for Thiel's test set of small molecules and confirm the importance of the augmentation basis functions for adequate convergence. Adopting a velocity-gauge formalism, we next demonstrate applications for systems with periodic boundary conditions. Taking advantage of the all-electron full-potential implementation, we present applications for core level spectra. For it-TDDFT, we confirm that within the all-electron NAO formalism, it-TDDFT can successfully converge systems that are difficult to converge in the standard self-consistent field method. We finally benchmark our implementation for systems up to ∼500 atoms. The implementation exhibits almost linear weak and strong scaling behavior.

Position of the Benzene Ring Substituent Regulates the Excited-State Deactivation Process of the Benzyluracil Systems

A systematic theoretical study of the regulating effect of the substituent position on the photoinduced deactivation process of the benzyluracil systems has been performed based on the high-level static electronic structure calculations and on-the-fly full-dimensional excited-state dynamics simulations. Similarities and differences coexist for the two systems by comparative studies on the photoinduced deactivation process of the 5-benzyluracil (5-BU) and 6-benzyluracil (6-BU) systems. They both obey an S₂ → S₁ → S₀ two-step decay pattern, and the decay coordinates of the S₂ → S₁ and S₁ → S₀ processes are mainly driven by the elongation of the bridging bond and the out-of-plane ring deformation motion, respectively. However, the puckering motion occurring at the C2 atom in the uracil fragment dominates the decay pathway of the 5-BU system. On the contrary, the puckering motion at the C5 atom in the benzene fragment mainly drives the decay coordinate of the 6-BU system. Therefore, the substituent position could play significant roles in the deactivation process of the benzyluracil systems. Moreover, the S₁ → S₀ decay process of the 6-BU system consists of five pathways, possessing a more complex deactivation picture than the 5-BU system. The fitted time scale of the puckering motion is compatible with the experimentally observed lifetimes. This work provides a fundamental understanding of the photophysical and photochemical properties of the benzyluracil systems and can give rational suggestions to further design or regulate the bionic molecular systems.