1
|
Valverde D, Mai S, Canuto S, Borin AC, González L. Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water. JACS AU 2022; 2:1699-1711. [PMID: 35911449 PMCID: PMC9327080 DOI: 10.1021/jacsau.2c00250] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S2 state that converts internally to the dark S1 state, from which the population is transferred to the triplet T2 state via intersystem crossing and finally to the lowest T1 state. This S2 → S1 → T2 → T1 deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T1 lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T1/S0 crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ1 = 131 fs and τ2 = 191 fs) that are in excellent agreement with the experimentally reported value (τexp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.
Collapse
Affiliation(s)
- Danillo Valverde
- Department
of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Avenida Professor Lineu Prestes, 748, São Paulo, São Paulo CEP 05508-000, Brazil
- Institute
of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, São Paulo CEP 05508-090, Brazil
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, Vienna 1090, Austria
| | - Sylvio Canuto
- Institute
of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, São Paulo CEP 05508-090, Brazil
| | - Antonio Carlos Borin
- Department
of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Avenida Professor Lineu Prestes, 748, São Paulo, São Paulo CEP 05508-000, Brazil
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, Vienna 1090, Austria
| |
Collapse
|
2
|
Xie BB, Tang XF, Liu XY, Chang XP, Cui G. Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations. Phys Chem Chem Phys 2021; 23:27124-27149. [PMID: 34849517 DOI: 10.1039/d1cp03994f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China.
| | - Xiu-Fang Tang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China.
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, Sichuan, 610068, China
| | - Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China.
| |
Collapse
|
3
|
Kadaoluwa Pathirannahalage SP, Meftahi N, Elbourne A, Weiss ACG, McConville CF, Padua A, Winkler DA, Costa Gomes M, Greaves TL, Le TC, Besford QA, Christofferson AJ. Systematic Comparison of the Structural and Dynamic Properties of Commonly Used Water Models for Molecular Dynamics Simulations. J Chem Inf Model 2021; 61:4521-4536. [PMID: 34406000 DOI: 10.1021/acs.jcim.1c00794] [Citation(s) in RCA: 64] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Water is a unique solvent that is ubiquitous in biology and present in a variety of solutions, mixtures, and materials settings. It therefore forms the basis for all molecular dynamics simulations of biological phenomena, as well as for many chemical, industrial, and materials investigations. Over the years, many water models have been developed, and it remains a challenge to find a single water model that accurately reproduces all experimental properties of water simultaneously. Here, we report a comprehensive comparison of structural and dynamic properties of 30 commonly used 3-point, 4-point, 5-point, and polarizable water models simulated using consistent settings and analysis methods. For the properties of density, coordination number, surface tension, dielectric constant, self-diffusion coefficient, and solvation free energy of methane, models published within the past two decades consistently show better agreement with experimental values compared to models published earlier, albeit with some notable exceptions. However, no single model reproduced all experimental values exactly, highlighting the need to carefully choose a water model for a particular study, depending on the phenomena of interest. Finally, machine learning algorithms quantified the relationship between the water model force field parameters and the resulting bulk properties, providing insight into the parameter-property relationship and illustrating the challenges of developing a water model that can accurately reproduce all properties of water simultaneously.
Collapse
Affiliation(s)
- Sachini P Kadaoluwa Pathirannahalage
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia.,Laboratoire de Chimie, Ecole Normale Supérieure de Lyon, CNRS, Lyon 69342, France
| | - Nastaran Meftahi
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Aaron Elbourne
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Alessia C G Weiss
- Leibniz-Institut für Polymerforschung e.V., Hohe Straße 6, 01069 Dresden, Germany
| | - Chris F McConville
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia.,Institute for Frontier Materials, Deakin University, Geelong, Victoria 3220, Australia
| | - Agilio Padua
- Laboratoire de Chimie, Ecole Normale Supérieure de Lyon, CNRS, Lyon 69342, France
| | - David A Winkler
- School of Biochemistry and Genetics, La Trobe Institute for Molecular Science, La Trobe University, Bundoora, Victoria 3086, Australia.,Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria 3052, Australia.,School of Pharmacy, University of Nottingham, Nottingham NG7 2QL, U.K
| | | | - Tamar L Greaves
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Tu C Le
- School of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
| | - Quinn A Besford
- Leibniz-Institut für Polymerforschung e.V., Hohe Straße 6, 01069 Dresden, Germany
| | - Andrew J Christofferson
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia.,ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| |
Collapse
|
4
|
Fang YG, Valverde D, Mai S, Canuto S, Borin AC, Cui G, González L. Excited-State Properties and Relaxation Pathways of Selenium-Substituted Guanine Nucleobase in Aqueous Solution and DNA Duplex. J Phys Chem B 2021; 125:1778-1789. [PMID: 33570942 PMCID: PMC8023715 DOI: 10.1021/acs.jpcb.0c10855] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
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 S1 1(nSeπ5*) and S2 1(πSeπ5*) states are predicted to be the spectroscopically dark and bright states, respectively. Two triplet states, T1 3(πSeπ5*) and T2 3(nSeπ5*), are found energetically below the S2 state. Extending the QM region to include the 6SeG-Cyt base pair slightly stabilizes the S2 state and destabilizes the S1, 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 S2 state, 6SeG evolves to the S2 minimum, near to a S2/S1 conical intersection that allows for internal conversion to the S1 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 S2/S1 conical intersection. From the S1 state, 6SeG can further decay to the T1 3(πSeπ5*) state via intersystem crossing, where it will be trapped due to the existence of a sizable energy barrier between the T1 minimum and the T1/S0 crossing point. Although this general S2 → T1 mechanism takes place in both media, the presence of DNA induces a steeper S2 potential energy surface, that it is expected to accelerate the S2 → S1 internal conversion.
Collapse
Affiliation(s)
- Ye-Guang Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Danillo Valverde
- Institute of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, SP 05508-090, Brazil
| | - Sebastian Mai
- Photonics Institute, Vienna University of Technology, Gußhausstraße 27-29, 1040 Vienna, Austria.,Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Sylvio Canuto
- Institute of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, SP 05508-090, Brazil
| | - Antonio Carlos Borin
- Department of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Av. Prof. Lineu Prestes 748, 05508-000. São Paulo, SP Brazil
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| |
Collapse
|
5
|
Peng Q, Zhu YH, Zhang TS, Liu XY, Fang WH, Cui G. Selenium substitution effects on excited-state properties and photophysics of uracil: a MS-CASPT2 study. Phys Chem Chem Phys 2020; 22:12120-12128. [PMID: 32440669 DOI: 10.1039/d0cp01369b] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The photophysics of selenium-substituted nucleobases has attracted recent experimental attention because they could serve as potential photosensitizers in photodynamic therapy. Herein, we present a comprehensive MS-CASPT2 study on the spectroscopic and excited-state properties, and photophysics of 2-selenouracil (2SeU), 4-selenouracil (4SeU), and 2,4-selenouracil (24SeU). Relevant minima, conical intersections, crossing points, and excited-state relaxation paths in the lowest five electronic states (i.e., S0, S1, S2, T2, and T1) are explored. On the basis of these results, their photophysical mechanisms are proposed. Upon photoirradiation to the bright S2 state, 2SeU quickly relaxes to its S2 minimum and then moves in an essentially barrierless way to a nearby S2/S1 conical intersection near which the S1 state is populated. Next, the S1 system arrives at an S1/T2/T1 intersection where a large S1/T1 spin-orbit coupling of 430.8 cm-1 makes the T1 state populated. In this state, a barrier of 6.8 kcal mol-1 will trap 2SeU for a while. In parallel, for 4SeU or 24SeU, the system first relaxes to the S2 minimum and then overcomes a small barrier to approach an S2/S1 conical intersection. Once hopping to the S1 state, there exists an extended region with very close S1, T2, and T1 energies. Similarly, a large S1/T1 spin-orbit coupling of 426.8 cm-1 drives the S1→ T1 intersystem crossing process thereby making the T1 state populated. Similarly, an energy barrier heavily suppresses electronic transition to the S0 state. The present work manifests that different selenium substitutions on uracil can lead to a certain extent of different vertical and adiabatic excitation energies, excited-state properties, and relaxation pathways. These insights could help understand the photophysics of selenium-substituted nucleobases.
Collapse
Affiliation(s)
- Qin Peng
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | | | | | | | | | | |
Collapse
|
6
|
Sun H, Zhao L, Peng S, Huang N. Incorporating replacement free energy of binding-site waters in molecular docking. Proteins 2014; 82:1765-76. [DOI: 10.1002/prot.24530] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Revised: 01/17/2014] [Accepted: 01/28/2014] [Indexed: 12/24/2022]
Affiliation(s)
- Hanzi Sun
- College of Life Sciences; Beijing Normal University; Beijing 100875 China
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| | - Lifeng Zhao
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| | - Shiming Peng
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| | - Niu Huang
- College of Life Sciences; Beijing Normal University; Beijing 100875 China
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| |
Collapse
|
7
|
Vijay Solomon R, Angeline Vedha S, Venuvanalingam P. A new turn in codon–anticodon selection through halogen bonds. Phys Chem Chem Phys 2014; 16:7430-40. [DOI: 10.1039/c3cp54442g] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
8
|
Sheng J, Gan J, Soares AS, Salon J, Huang Z. Structural insights of non-canonical U*U pair and Hoogsteen interaction probed with Se atom. Nucleic Acids Res 2013; 41:10476-87. [PMID: 24013566 PMCID: PMC3905866 DOI: 10.1093/nar/gkt799] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Unlike DNA, in addition to the 2′-OH group, uracil nucleobase and its modifications play essential roles in structure and function diversities of non-coding RNAs. Non-canonical U•U base pair is ubiquitous in non-coding RNAs, which are highly diversified. However, it is not completely clear how uracil plays the diversifing roles. To investigate and compare the uracil in U-A and U•U base pairs, we have decided to probe them with a selenium atom by synthesizing the novel 4-Se-uridine (SeU) phosphoramidite and Se-nucleobase-modified RNAs (SeU-RNAs), where the exo-4-oxygen of uracil is replaced by selenium. Our crystal structure studies of U-A and U•U pairs reveal that the native and Se-derivatized structures are virtually identical, and both U-A and U•U pairs can accommodate large Se atoms. Our thermostability and crystal structure studies indicate that the weakened H-bonding in U-A pair may be compensated by the base stacking, and that the stacking of the trans-Hoogsteen U•U pairs may stabilize RNA duplex and its junction. Our result confirms that the hydrogen bond (O4…H-C5) of the Hoogsteen pair is weak. Using the Se atom probe, our Se-functionalization studies reveal more insights into the U•U interaction and U-participation in structure and function diversification of nucleic acids.
Collapse
Affiliation(s)
- Jia Sheng
- Department of Chemistry, Georgia State University, Atlanta, GA, 30303, USA and Department of Biology, Brookhaven National Laboratory, Upton, NY, 11973, USA
| | | | | | | | | |
Collapse
|