
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

38851775
48871
10.1038/s41467-024-48871-4
Article
Sub-100-fs energy transfer in coenzyme NADH is a coherent process assisted by a charge-transfer state
http://orcid.org/0000-0002-5090-7984
Jaiswal Vishal Kumar 1
http://orcid.org/0000-0003-0747-6266
Aranda Ruiz Daniel 2
Petropoulos Vasilis 3
http://orcid.org/0000-0003-4591-5100
Kabaciński Piotr 3
Montorsi Francesco 1
http://orcid.org/0000-0002-6969-4894
Uboldi Lorenzo 3
http://orcid.org/0009-0002-2750-8093
Ugolini Simone 1
http://orcid.org/0000-0002-6015-3135
Mukamel Shaul 4
http://orcid.org/0000-0002-9534-2702
Cerullo Giulio giulio.cerullo@polimi.it

3
http://orcid.org/0000-0002-0796-289X
Garavelli Marco marco.garavelli@unibo.it

1
http://orcid.org/0000-0003-4402-2685
Santoro Fabrizio 5
http://orcid.org/0000-0003-3071-5341
Nenov Artur artur.nenov@unibo.it

1
1 https://ror.org/01111rn36 grid.6292.f 0000 0004 1757 1758 Dipartimento di Chimica industriale “Toso Montanari”, Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
2 https://ror.org/043nxc105 grid.5338.d 0000 0001 2173 938X ICMol, Universidad de Valencia, Catedrático José Beltrán Martínez, 2, 46980 Paterna, Spain
3 https://ror.org/01nffqt88 grid.4643.5 0000 0004 1937 0327 Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy
4 grid.266093.8 0000 0001 0668 7243 Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, CA 92697 USA
5 Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
8 6 2024
8 6 2024
2024
15 490018 9 2023
15 5 2024
© The Author(s) 2024, corrected publication 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Excitation energy transfer (EET) is a key photoinduced process in biological chromophoric assemblies. Here we investigate the factors which can drive EET into efficient ultrafast sub-ps regimes. We demonstrate how a coherent transport of electronic population could facilitate this in water solvated NADH coenzyme and uncover the role of an intermediate dark charge-transfer state. High temporal resolution ultrafast optical spectroscopy gives a 54±11 fs time constant for the EET process. Nonadiabatic quantum dynamical simulations computed through the time-evolution of multidimensional wavepackets suggest that the population transfer is mediated by photoexcited molecular vibrations due to strong coupling between the electronic states. The polar aqueous solvent environment leads to the active participation of a dark charge transfer state, accelerating the vibronically coherent EET process in favorably stacked conformers and solvent cavities. Our work demonstrates how the interplay of structural and environmental factors leads to diverse pathways for the EET process in flexible heterodimers and provides general insights relevant for coherent EET processes in stacked multichromophoric aggregates like DNA strands.

Excitation energy transfer is important for many photoinduced biological processes in systems with multiple chromophores. Here, the authors elucidate this process for the coenzyme NADH using ultrafast spectroscopy and quantum dynamics.

Subject terms

Excited states
Energy transfer
Quantum chemistry
Optical spectroscopy
https://doi.org/10.13039/100006151 DOE | SC | Basic Energy Sciences (BES) DE-SC0022225 DE-SC0022225 DE-SC0022225 DE-SC0022225 DE-SC0022225 Jaiswal Vishal Kumar Montorsi Francesco Mukamel Shaul Garavelli Marco Nenov Artur PRIN: PROGETTI DI RICERCA DI RILEVANTE INTERESSE NAZIONALE-Bando 2022, Prot. 2022HL9PRP. PRIN: PROGETTI DI RICERCA DI RILEVANTE INTERESSE NAZIONALE-Bando 2022, Prot. P20224AWLB.PNRR MUR project ECS_00000033_ECOSISTER Si-Fi-MYSTERY project, PRIN: PROGETTI DI RICERCA DI RILEVANTE INTERESSE NAZIONALE-Bando 2022, Prot. P2022WSJKSissue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Following the excitation energy transfer (EET) between chromophoric units in multimeric aggregates is of key importance for understanding biological processes like photosynthesis and DNA photodamage, as well as for the design of efficient light-harvesting molecular assemblies.1–3 EET in these systems falls between the regimes of coherent and incoherent energy transport.4 Where a process is placed in this sliding scale is dictated by: a) the energetic separation of relevant electronic states versus the coupling strength to each other and with the environment5–9; b) the spatial separation of the monomers in the macro-structure.

In homoaggregates, coherent population oscillations of quantum origin between the individual sites coupled to each other with a fixed phase relation in so-called excitonic excited states (ES) arise due to purely electronic coupling between the monomers. Coupling the electronic dynamics to a bath of nuclear degrees of freedom facilitates an efficient and ultrafast unidirectional population transfer to the lowest excitonic state. The role of such electronic coherences in disordered biological systems, such as the reaction centers of light-harvesting complexes, is being intensely debated10–13.

In heterodimers with energetically well-separated electronic states, the motion of the nuclear wavepacket can strongly tune the energy gaps promoting crossings where electronic and vibronic couplings, i.e. wave function mixing along nuclear degrees of freedom, facilitate the population transfer14–21. In the limit of large donor-acceptor separation (i.e. >10 Å) inter-chromophore orbital overlaps are negligible, and couplings are dictated by the long-range electrostatic interactions. As the resulting electronic couplings are weaker than the coupling to the environment, the coherently prepared excited state first thermalizes through vibrational cooling on the picosecond time scale22,23. The EET then occurs incoherently on tens of ps-to-ns time scale, with rates that can be estimated by utilizing Förster resonance energy transfer theory. In the case of closely packed heteroaggregates, on the other hand, variations of the inter-chromophore orbital overlaps and mixing with charge-transfer (CT) configurations24,25 occurring during the photoinduced vibrational dynamics give rise to electronic and vibronic couplings which can be as pronounced as their intra-molecular counterparts, accelerating the EET down to the sub-ps regime. In this limit, the coherent motion of the nuclear wavepacket before thermalization (referred to as classical coherence) induces an inherently coherent EET competitive even with possible sub-ps decay channels such as intramolecular internal conversion (IC) to the ground state (GS) mediated by conical intersections (CIs).

NADH, the reduced form of nicotinamide adenine dinucleotide (NAD), is an important coenzyme found in the mitochondria that is part of the electron transport chain leading to the generation of adenosine triphosphate26. It is a dimer consisting of the two chromophores adenine (Ade) and nicotinamide (Nic), which absorb light in the ultraviolet (UV) region, connected through a phosphate bridge. This bichromophoric structure provides a playground rich with photoinduced processes, which are used as spectroscopic markers for cellular metabolism, for example, in label-free multimodal microscopy27,28.

The high conformational flexibility afforded by the phosphate bridge linking the two monomers leads to DNA-like stacked (folded) or unstacked (unfolded) conformers whose ratio depends on solvent polarity, pH, and temperature29–33. While in alcoholic environments, like propylene glycol or methanol, there is a high propensity of unstacked conformations, a considerable percentage of stacked conformers is observed in water34–37. In its folded state, NADH is characterized by a photoinduced ultrafast EET process from Ade (donor) to Nic (acceptor). Since its first observation by Weber in 195734, the EET process, which involves the frontier π-orbitals of the locally excited (LE) states of the two monomers38,39, has been extensively studied. Upon impulsive resonant excitation of the La-state of Ade, immediate (within the ≈100-fs temporal resolution of the ultrafast experiments) emission was observed from the lowest ππ* state of Nic (referred henceforth as Nic*). Previous time-resolved fluorescence and transient absorption (TA) experiments revealed the ultrafast nature of the process, placing the EET timescale in the sub-100-fs range36,37, but with significant uncertainty due to their limited temporal resolution.

In this letter, using the example of NADH, we illuminate the structural, energetic and environmental conditions which favor a sub-100-fs coherent EET mechanism against IC in a closely stacked molecular heteroaggregate. Using UV TA spectroscopy with sub-30-fs temporal resolution, we resolve the EET from Ade to Nic* in water and find that, with a time constant τEET = 54 ± 11 fs, it is much faster than intra-Ade IC to the GS, which occurs with a time constant τIC ~ 160 fs40–46. This observation places the mechanism through which the EET occurs completely outside the applicability of Förster theory still employed in the analysis of NADH photophysics36,37. Through multidimensional nonadiabatic quantum dynamics, we anticipate that this ultrafast EET in folded conformers is channeled through coherent vibrational motion along the same planar modes that facilitate intra-Ade IC. These coherent molecular vibrations contribute to the EET being faster than IC to the GS as, for the latter process, additional out-of-plane structural deformations are needed to bring the La and GS electronic potential energy surfaces to degeneracy. This is a general system-independent observation which implies that EET is a competitive deactivation process in closely stacked aggregates. Furthermore, our simulations reveal a dark Nic→Ade CT state, highly susceptible to structural and environmental disorder. In favorably stacked dimers and in cooperation with the fluctuations of the polar environment, the CT state is stabilized in the vicinity of the La state and acts as an intermediary that further boosts the EET yield.

Results

Transient absorption spectra of adenosine and NADH

Figure 1a, b depict the TA maps of Adenosine (a) and NADH (b) in water upon resonant excitation of the Ade moiety at 4.7 eV with sub-20-fs pulses and broadband (1.9-3.2 eV) probing. The excited state of the Ade monomer exhibits an intense photo-induced absorption (PA) band that spans our entire probe photon energy range, characterized by a sharp peak at >3 eV (Fig. 1a). This PA band decays on the <300-fs timescale as a result of Ade’s La-state IC pathway back to the ground state. This is consistent with the vibronically driven 100–300 fs IC, reported for Ade and its derivatives40–47.Fig. 1 Ultrafast energy transfer in NADH.

Transient absorption (TA) maps of a Adenosine (Ade) and b NADH in piperazine-N,N′-bis(2-ethanesulfonic acid) (PIPES) aqueous buffer upon photoexcitation with sub-20-fs pulses at 4.7 eV. The dashed lines indicate the selected time-traces shown in Fig. 1c, d for detection at 3.1 eV (black dashed line) and 2.5 eV (gray dashed line). c, d Time-traces of Ade (green) and NADH (blue) for probing at c 3.1 eV and d 2.5 eV respectively. The faint gray line in Fig. 1c, d defines the zero line (ΔA = 0). Schematics depicting the ultrafast vibronic processes after Ade resonant excitation which lead to e internal conversion (IC) in Ade and f direct/through-CT coherent EET in NADH. Source data are provided as a Source Data file.

The 4.7 eV UV pulse predominantly photoexcites the Ade moiety in NADH, resulting in an initial PA band similar to the one reported for the monomeric Ade (Fig. 2b). As mentioned earlier, the equilibrium between folded/unfolded conformers in NADH heavily depends on the solvent. In methanol, NADH exists almost entirely in its open form and thus the EET pathway is negligible, resulting in excited state dynamics similar to isolated Ade in water (see Suppl. Fig. 1 and Suppl. Fig. 2 in the Supplementary Information (SI)). However, in water, NADH has an approximate 30/70 population ratio between folded and unfolded configurations37. The significant population of folded structures activates the EET from Ade to Nic, resulting in an additional ultrafast de-activation channel of Ade’s excited state (Fig. 1b). At 3.1 eV probe photon energy, the EET process results in a long-lived signal, due to a PA band from the Nic* state (Fig. 1c). In contrast, for a 2.5 eV probe, only Ade displays a significant PA. Comparison of the dynamics of Adenosine and NADH in water at this probe photon energy (Fig. 1d) reveals a faster decay for NADH, due to the additional deactivation channel for the Laexcited state of Ade through EET to Nic*.

Estimation of experimental EET timescale

Global fitting of the TA data, as shown in Table 1, gives time constants of τ1Ade = 157 ± 4 fs and τ1NADH, water = 122 ± 4 fs for Ade and NADH in water, respectively. While in isolated Ade and NADH solvated in methanol, the corresponding time constants of 157 fs and 167 fs correspond to the IC back to the GS, in NADH solvated in water the EET mechanism opens up an additional channel of excited state deactivation, increasing its rate. Taking into account the 30/70 equilibrium between folded and unfolded forms, the faster lifetime τ1NADH, water = 122 fs extracted by global analysis can be considered as a weighted average of the folded and unfolded ultrafast lifetime components:τ1NADH,water=0.7τ1NADH,unfolded+0.3τ1NADH,folded

Table 1 Time constants, accompanied by their experimental uncertainty, as obtained from the global fitting of the TA data in the case of Adenosine in water and NADH in water

	τ1 (fs)	τ2 (fs)	τ3 (fs)	
Adenosine in water	157 ± 4	894 ± 50	-	
NADH in water	122 ± 4	878 ± 40	long	

Assuming similar excited-state deactivations in unfolded water-solvated NADH and isolated Ade, one can take τNADH, unfolded = τ1Ade, water = 157 fs, which results in a time constant τNADH, folded = 40 ± 6 fs for the folded NADH conformers. This total decay rate obtained for the folded NADH components (ktotal = (40 fs)−1) reflects the sum of the rates for Ade IC (kIC = (157 fs)−1) and ultrafast EET from the Ade La-state to Nic* (kET), resulting in ktotal = kIC + kET, thus giving a time constant of τET, folded = 54 ± 11 fs for the EET process (see a detailed discussion in Suppl. Note 1.2 and 1.3).

Quantum dynamics simulations of the EET in NADH

To gain insight into the ~50-fs Ade →Nic EET mechanism in the folded NADH forms, we have modeled the photoinduced time-evolution of the system by nonadiabatic quantum dynamics of multidimensional wavepackets using the MCTDH method48–52. Due to the large number of normal modes in the bichromophoric molecule, the multilayer formulation of the MCTDH method (ML-MCTDH) was employed53–56. The coupled dynamics of the electronic states of the system through the vibrational degrees of freedom are described by a Linear Vibronic Coupling (LVC) Hamiltonian57,58, parametrized with the energies and gradients computed at the multiconfigurational wavefunction-based XMS-CASPT2 level of theory59. Calculations were performed on different water solvated stacked NADH aggregates, obtained from a 20 ns of Replica Exchange Molecular Dynamics (REMD)60, in a quantum mechanics/molecular mechanics (QM/MM)61 setup allowing to account for the conformational diversity and solvent effects with quantitative accuracy. Details are provided in the Suppl. Note 2 to 7.

The structural flexibility of NADH leads to diverse relative orientations of the two bases in the stacked conformers, whose relative populations in REMD dynamics are detailed in Suppl. Table 2. Figure 2d–f shows the time-evolution of electronic states population in representative structures of the three most populated stacked conformers, capturing 65% of stacked population. The dynamics reveal that EET is an ultrafast process mediated by coherent vibrational motions. The coupling between the electronic states induced by nuclear motions along these vibrations facilitates a near complete Ade →Nic population transfer on a sub-100-fs time scale. As shown in Fig. 2 and additional dynamics reported in Suppl. Note. 8 of the SI, the EET mechanism can be either direct (Fig. 2f) or proceed involving a Nic→Ade CT state as an intermediary (Fig. 2d,e). The propensity for a direct or CT-mediated EET mechanism is dictated by the energetic position of the CT state versus the La. We shall show later how this is affected by the structural and solvent heterogeneity present in the system.Fig. 2 Ultrafast energy transfer dynamics in highest populated clusters.

a–c Conformational representatives of the three largest clusters of Replica-Exchange Molecular Dynamics. The relative orientation of two chromophores in different clusters is shown through arrows colored blue (on Nic) and red (on Ade). Carbon atoms are colored in green, nitrogens in blue, hydrogens in white, oxygen in red and phosphorus in yellow. Details about clustering and the population of clusters are given in Suppl. Note. 2. d–f Quantum Dynamics (using ML-MCTDH) on a representative structure of the three largest clusters. Additional dynamics on other representative structures are shown in Suppl. Fig. 10. Source data are provided as a Source Data file.

Vibronically coherent mechanisms are ubiquitous in ultrafast photophysical processes involving population transfer between electronic states62,63. Impulsively excited molecular vibrations, called tuning modes, modulate the energy gaps between electronic states, whereas vibrations known as coupling modes, together with purely electronic coupling, induce wave function mixing which dictates the rate of the non-adiabatic population transfer between electronic states. The direct transfer in NADH is mediated by high-frequency stretching modes of the initially excited Ade moiety. In particular, the C-C and C-N stretchings, with frequencies of ~1400 cm−1 and ~1600 cm−1, respectively, show the highest tuning and, simultaneously, significant coupling between La and Nic*, irrespective of stacking conformation (Suppl. Fig. 15 in the SI). The oscillatory pattern in the first 100 fs of the population dynamics (Fig. 2f) with a period of ≈10 fs is the result of the tuning modes leading to energy degeneracy (real La/Nic* CI) between La and Nic* every half a period (Suppl. Fig. 12 and Suppl. Fig. 16 in SI).

To gain further insight, we computed the quantum electronic coherences between La and Nic* along the dynamics (Suppl. Note. 15) and observed that they last for at least 70 fs while exhibiting recurrences with a period comparable to the frequencies of the stretching modes (Suppl. Fig. 18). Their eventual decay is caused by the concurrent effects of: a) progessive population depletion in the Lastate accompanied by the accumulation in the Nic* state; b) progessive decrease of the overlap of the wavepackets evolving in the La and Nic* states. The direct EET is facilitated to equal amounts by electronic (Eij0) and vibronic (λij) couplings, which can be verified by switching off their respective contribution (Suppl. Fig. 18). The population transfer becomes slower in both cases, yet, notably, population inversion still occurs within 100 fs. Interestingly, the observation of ultrafast EET in the absence of electronic couplings makes apparent that in closely stacked aggregates even small localized deformations such as stretchings induce vibronic couplings large enough to facilitate a sub-100-fs EET.

We note that intra-monomer IC in Ade is initially driven by the same impusively excited high-frequency modes. However, as the energy gap (with the GS) to overcome is significantly bigger (ca. 5 eV with respect to ca. 1.5 eV), additional deformations, such as ring puckering, are necessary to facilitate IC to the GS42,46,64,65. Since these deformations need activation and are associated with more pronounced displacements, IC is naturally slower, taking a few hundreds of fs. Source data are provided as a Source Data file.

Factors affecting the participation of CT state in EET process

The extent of participation of the dark CT state in the EET process is determined by its energetic position and coupling to the two monomer LE states La and Nic*. The structural and solvent dynamics impart a large amount of heterogeneity to the system, leading to various types of stacked conformers with different relative orientations of the two bases (Fig. 2a–c). While the energetic positions of the La and Nic* are relatively unaffected by the different arrangements of the two chromophores, the CT state is markedly more sensitive. Figure 3 shows how both the stacking and solvent configuration individually affect the energetic position of the CT state. The transfer of an electron from Nic to Ade in the CT state creates an electron-hole dipole in the system (Fig. 3b). The relative orientation of the two bases modulates the distance between the centers of (positive and negative) charge of this electron-hole pair in stacked conformations from 3 to 5 Å. In return, this causes a variation of the CT state energy by ~1 eV in vacuo (Fig. 3c) at the time of photoexcitation. Introducing the solvent, the interaction of the CT dipole with the electric field of water additionally modulates the fluctuation of CT state energy by ~1 eV, as seen from computations in solvent ensemble around two geometries selected from the manifold of stacked conformations (Fig. 3d, e).Fig. 3 Stacking and solvent effects on energies of charge-transfer (CT) and locally excited (LE) states.

a The natural transition orbitals of the two LE and lowest CT state involved in the EET process. b The electron-hole distance77,78 in the Nic→Ade CT state in a representative stacked conformer depicted with a dashed line. The two dots specify the position of the center of charges of hole and electron created in the CT state. c Stacking Effect: Vertical transition energies of the CT state vs the electron-hole distance of the CT state in gas-phase. The black and green bands indicate the energies of Nic* and La state which show minimal fluctuation compared to CT state. The Gaussian profile shows the spread of the values of the CT state with FWHM of ca. 1 eV. Two structures with markedly different CT energies are highlighted with black dots for studying solvent effects. d, e Solvent Effect: The effect of the solvent fluctuations on the energies of the diabatic states for the two structures highlighted in (c), which represent respectively a case in which the CT is near degenerate with the La-state of Ade (d) and a case in which the CT state is significantly higher in energy (e). The two profiles were obtained by running an equilibrated solvent MD around the two frozen NADH solute geometry and computing vertical transition energies at XMS-CASPT2 level. Details in Suppl. Note. 14.

The combined effects of stacking and solvent orientation dictate whether the CT state can actively participate in ultrafast EET process. If a certain stacking places the CT state more than 0.5 eV above the La state in gas-phase, its energy fluctuations due to the interaction with the solvent cannot render the CT accessible from the La state after photoexcitation (Fig. 3e). In such conformations, populating predominantly cluster C2 – with a median interbase distance of ~4.6 Å – the EET occurs via a direct transfer. Instead, for stacking where the CT state is energetically closer to the La (i.e. within less than 0.5 eV) in gas-phase, the solvent fluctuations can induce stabilization of the CT band around and even below the La state (Fig. 3d), eventually enabling its active participation in the EET dynamics. Such is the case for the most populated stacked cluster C0 – with a median interbase distance of ~4.12 Å. The average ratio between La/CT electronic coupling and energy gap |V/Δ| = 0.268 calculated for the conformations in this cluster indicates a strong wave function mixing and, consequently, admixing of CT character in the bright adiabatic state in the Franck-Condon (FC) region.

Figure 2d shows a representative dynamics for a snapshot of cluster C0 in which the CT state is 0.5 eV below the La in the FC region. A virtually immediate increase of the population of the CT state to 50% on a 10 fs timescale is observed, which is accompanied by the instantaneous formation of a quantum coherence which lasts for about 70 fs (column LVC in Suppl. Fig. 19). Tuning along the high-frequency C-C and C-N stretching modes quickly overcomes the La/CT energy gap, so that electronic and vibronic couplings – exhibiting larger magnitudes compared to their counterparts in the direct EET due to the one-electron nature of the process – efficiently promote La→CT population transfer. Notably, if tuning modes are switched off, thus forcing purely electronic dynamics between stationary wavepackets, no population transfer takes place despite the formation of coherences of sizeable magnitude (column “no λii All” in Suppl. Fig. 19). This makes evident that tuning modes are essential to circumvent the La/CT energy gap for the electronic and vibronic couplings to have effect. Furthermore, we selectively switched off either coupling to investigate its impact on the Nic→Ade electron transfer event (columns “no Eij0 All” and “no λij All” in Suppl. Fig. 19). Purely electronic couplings are found to promote La/CT population inversion on a 10-fs timescale. Electronic coherence formation is vibrationally assisted as can be noted by the coherence recurrences with times compatible with the periods of the high-frequency tuning modes, and by the fact that the coherence magnitude moderately increases when switching off electronic couplings.

The population of the CT state is only transient and is transferred to Nic* creating, at the same time, an electronic coherence between CT and Nic* that lasts till 150 fs. Contrary to the high frequency modes driving the direct EET, we find that multiple conformation-dependent low frequency modes delocalized over both moieties facilitate this second step of the CT-mediated EET. These low-frequency modes show the largest vibronic coupling between the CT and Nic* state amongst the normal modes of the system. When the vibronic couplings along these modes are switched off (see Suppl. Fig. 14 in SI for a representative case), the transfer to Nic* from the CT state is significantly slowed down and the population remains trapped in the CT state.

To explore the role of solvent heterogeneity in the ultrafast EET process, we carried out independent wavepacket dynamics on a representative structure from the most populated conformer embedded in 100 different solvent configurations obtained through equilibrated molecular dynamics at room temperature (Fig. 4a). Thereby, for each solvent arrangement the LVC Hamiltonian was reparametrized to take into account the solvent effect on the energetics and electronic coupling (details in Suppl. Note. 14.1). Figure 4b shows a false color plot of the state-specific population evolution upon instantaneous La photoexcitation along the 100 wavepacket dynamics ordered according to the energy gap between CT and La in the FC point. An ultrafast population transfer to Nic* is observed when the CT state is energetically below the La state at the FC geometry. The CT state acts as a doorway populated within 10-20 fs after excitation of the La state and depopulated within 60 fs by transferring the population to the Nic* state. The acceleration provided by the CT state becomes apparent when comparing against dynamics performed in the same solvent ensemble after switching off LE-CT couplings (Suppl. Fig. 17 in SI, also see Suppl. Figs. 10/S11) where the EET process gets substantially slowed down. A similar retardation is observed for those solvent/stacking configurations where the CT state is substantially above the La (Fig. 4b) and, therefore, cannot act as a doorway in the EET mechanism. Thus, the acceleration provided by the CT state can play a paramount role in making the EET mechanism the dominant relaxation path immediately following photoexcitation.Fig. 4 Wavepacket dynamics on most populated conformer in 100 equilibrated solvent ensembles.

a An ensemble of equilibrated solvent configurations is created around a representative structure of most populated cluster of the REMD. The solute NADH is kept frozen during the solvent equilibrium dynamics. b–d Populations dynamics in the three diabatic electronic states - La (left), CT (middle), Nic* (right) – during the first 200 fs after photoexcitation of the La state as a function of the relative energy difference between the CT and Lastates (i.e. ECT-ELa). The plots are generated from individual ML-MCTDH dynamics performed on 100 snapshots from the solvent equilibrium dynamics as indicated in (a), each one with characteristic solvent-modulated electronic energies, electronic and vibronic couplings. Details on the simulation can be found in Suppl. Note. 14. Source data are provided as a Source Data file.

Besides the ultrafast population transfer to Nic*, we observe population trapping in the CT state in our dynamics reported in Fig. 4 when La and CT are energetically close in the FC region. In this case, after an ultrafast population of the CT state, a subsequent transfer to Nic* state is substantially slowed down as the CT and Nic* surfaces do not cross in the course of the vibrational dynamics in the CT state. However, it should be noted that our model lacks the dynamic relaxation of the electronic states induced by the temporal response of the polar solvent. Indeed, a charge-separated state can access lower energies upon relaxation of the molecules of the surrounding polar solvent which can occur on a timescale as fast as sub-50-fs in water66–68. Thus, we expect that the CT state stabilization induced by modeling the response of the polar solvent will further put these dynamics into the realm of ultrafast population transfer to Nic* state.

Discussion

In this work, we unambiguously demonstrate, by means of TA spectroscopy with sub-30-fs temporal resolution, that the ultrafast Ade to Nic EET in water-solvated NADH occurs on a 50-fs timescale. By means of multimode wavepacket quantum dynamics, we reveal that this is an electronically and vibrationally coherent process. Coherent molecular vibrations – high-frequency stretchings – activated upon optical excitation of Ade (the donor) tune the energy gap where electronic and vibronic couplings directly channel the electronic population to Nic (the acceptor). Depending on the spatial separation and relative orientation of the two monomers (which affects the magnitude of the couplings) the direct pathway exhibits timescales ranging from few tens of fs to several ps. Sub-100 fs direct EET pathways are only possible for very close ( < 5 Å) stacking, since this activates inter-chromophore vibronic couplings of the same order of magnitude as intra-chromophore ones.

The EET process can be further accelerated via the mediation of a low-lying CT state coupled to the locally excited states of the two monomers. The relative orientation of the two monomers in closely stacked conformations, along with the solvent thermal fluctuations, affect the energetics of CT state and determine the propensity of the CT-mediated EET pathway. When the CT state is near-degenerate or below the Lastate of Ade, the first step of the CT-mediated EET, i.e. the Ade →Nic hole transfer, occurs quasi instantaneously. Completion of the EET process is observed within 100 fs, through the subsequent Ade→Nic electron transfer, which is dependent on the relaxation timescale of the solvent.

Our findings are of general validity for describing the dependence of energy transfer and charge separation processes on the coupling between electronic and nuclear degrees of freedom, conformational heterogeneity and solvent fluctuations in closely packed hetero-aggregates. The EET depends on the interplay of many factors: nature of the tuning and coupling modes and respectively the amplitude of their motion or the strength of coupling they carry; conformational freedom and spatial distance between chromophores in the aggregate; solvent polarity and thermalizaiton timescale. We identify the parameter window which facilitates ultrafast EET. At an inter-chromophore distance of <5 Å: a) electronic and vibronic couplings are equally relevant and cooperate in promoting efficient EET; b) even localized vibrational modes are capable of generating sufficiently strong orbital mixing and, thus, strong vibronic couplings; c) the short distance between the centers of negative and positive charge of the two chromophores stabilizes the CT state opening an additional EET pathway, favored in polar solvents. Most importantly, despite the conformational heterogeneity of such macro-structures, the EET is an intrinsically coherent process, i.e. governed by coherently oscillating wavepackets on the potential energy surfaces of the involved electronic states which tune the energy gaps and drive the EET unidirectionally from the donor to the acceptor. Such coherent EET cannot be described by Förster theory, which is suited to describe incoherent EET mechanisms.

Classical coherences manifest themselves through quantum beating in the TA spectra69–71, yet, the high frequencies of the dominant modes in NADH prevent them from being registered with current state-of-the-art experimental setups. This makes the theoretical analysis much more relevant. We hope that our findings will be a stimulus to increase the temporal resolution of UV TA spectroscopy and to apply transient spectroscopies which interrogate other spectral ranges in which the CT state may show unambiguous fingerprints. In particular, we propose photoelectron spectroscopy as a more sensitive way to single out CT states, as ionization from the intermediately created negatively charged Ade would give rise to characteristic signatures at lower ionization energies compared to the neutral form. Eventually, this study discloses a learning path on how to underdstand, and thereafter tune, the control knobs that are crucial for the design of systems with tailored EET efficiency.

Methods

Experimental setup

Ultrafast TA experiments72 were performed using a Ti:sapphire laser generating 100 fs pulses at 800 nm wavelength and 1 kHz repetition rate. Deep UV pump pulses tunable in the 4.2-4.8 eV range were generated as the second harmonic of a visible non-collinear optical parametric amplifier and compressed to sub-20-fs duration with a prism pair. Probe pulses covering 1.9-3.2 eV were obtained through white-light continuum generation by focusing a fraction of the fundamental beam in a calcium fluoride plate. Pump and probe polarizations were set at the magic angle (54.7°).

Sample preparation

β-NADH, in the form of reduced disodium salt hydrate, and adenosine were purchased from Sigma-Aldrich and used as received, dissolved in either methanol or 0.1 M piperazine-N,N′-bis(2-ethanesulfonic acid)(PIPES) aqueous buffer at pH 7.0. The samples were flown in a 150-μm-thick laminar liquid jet configuration, resulting in the absorbance of 2 OD at 4.7 eV pump photon energy. The used pump fluence was below 300 μJ/cm2 to minimize the coherent artifact and solvated electron signals.

Estimation of experimental EET timescale

The timescale of the EET process was estimated by employing a sequential kinetic model, taking into account the diverse relaxation pathways undertaken in stacked and unstacked conformers. The timescale of the EET pathway was estimated from the obtained lifetimes of the various evolution-associated spectra (EAS) obtained through global analysis of the TA map. The full details are reported in Suppl. Note. 1.3.

Computational methods

The conformational heterogeneity of solvated NADH was sampled through replica-exchange molecular dynamics performed through AMBER molecular dynamics package as detailed in Suppl. Note. 2. The major conformers were identified through a cluster analysis using of internal distance matrix of the solute atoms. The population fraction and inter-base distances of major clusters are reported in Suppl Table 2. Hybrid QM/MM calculations were carried out with the COBRAMM73,74 program, interfacing Gaussian1675 and openMOLCAS76 QM codes with the AMBER molecular dynamics package. The QM and MM partitioning involves three layers: High, Medium, and Low. The two chromophore bases in NADH solute are treated at QM level (High level), while the rest of the backbone along with the solvent is treated at MM level (details in Suppl. Note. 4). All geometry optimizations were done while allowing the backbone and the hydrogen-bonded solvent molecules to relax to the QM solute (Medium level), while the rest of the MM water molecules are frozen (Low Layer). The ground state geometry minimum was obtained at the Møller−Plesset second order perturbation theory (MP2) level. All computations utilized the Pople 6-31 G* basis set employing the Cholesky decomposition method to speed up the computation of atomic integrals. The vertical energies and gradients for the LVC model were obtained at XMS-CASPT2 level utilizing an active space of 4 electrons in 4 orbitals, comprising the frontier HOMO and LUMO of the two bases. The validity of the minimal active space was established by benchmarking against larger active space computations discussed in Suppl. Note. 5. Structure analysis (NTOs, attachment and detachment densities, electron-hole distance) was performed with the WFA module of Molcas77,78. All CASPT2 computations were done with zero IPEA shift, and an imaginary shift of 0.2.

Transformation of the adiabatic XMS-CASPT2 states to locally-excited and pure-CT states was done by the Fragment Excitation Difference procedure outlined by Hsu et al.79. In this method an excitation matrix Δxmn is built which characterizes the excitation difference between the donor and acceptor. In our case, the diagonalization of this matrix separates out the locally excited states with eigenvalues +/−1 and the CT-state with an eigenvalue of 0. The transformation matrix thus obtained can be used to transform the Hamiltonian to the diabatic basis, to compute the static electronic coupling and the vibronically induced coupling between the diabatic states. Further details outlining the procedure are reported in Suppl. Note. 7.

ML-MCTDH dynamics were performed using a variant of the multilayer algorithm of Heidelberg MCTDH package80 as implemented in Quantics81 program. The system was parametrized with a Linear Vibronic Coupling model using energies, gradients and couplings computed at XMS-CASPT2 level. The wavepacket propagations were done with ML-MCTDH adopting a variable mean field scheme with a Runge-Kutta integrator up to a final time of 200 fs. The dynamics on different cluster representatives were performed with 62 photoactive vibrational modes selected from the total 84 normal modes of the bichromophoric system included in the QM part of QM/MM setup. To select the photoactive modes, the maximum from the absolute value of either the gradients of the electronic states or interstate couplings of every pair of electronic states at Franck-Condon for each mode was chosen. All the 84 modes were sorted in decreasing order based on this value and the first 62 modes were selected. In this way, the discarded modes are the least active upon excitation as they have smallest gradients/interstate couplings at Franck-Condon (less than ~ 0.02 eV) and don’t promote population transfer. The selected 62 modes were partitioned into system (16) and bath modes (46). For the primitive basis-set we adopted Hermite DVR functions. The 16 coordinates containing the largest coupling/gradients were grouped together by pairs of coordinates with similar frequency in the ML tree and described by a larger number of primitive (30) and single-particle functions, while the remainder 46 (the “bath”) were grouped in pairs and described with less primitives (15) and single-particle functions. The 100 dynamics in different solvent ensembles were performed with 62 normal modes partitioned into a main system (8 modes) and bath (54 modes). The system modes were described 30 primitive bases and bath modes 15 primitive basis functions.

Supplementary information

Supplementary information

Peer Review File

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-48871-4.

Acknowledgements

D.A. acknowledges Fundación Ramón Areces and Generalitat Valenciana/European Social Fund (APOSTD/2021/025) for fundings and ICCOM-CNR (Pisa)/ICMol-MolMatTC (Valencia) for hospitality. Support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences and Biosciences Division under award no. DE-SC0022225 (A.N., F.M., S.M., V.K.J., and M.G.) is gratefully acknowledged. AN acknowledge financial support from PNRR MUR project ECS_00000033_ECOSISTER and from the Si-Fi-MYSTERY project, PRIN: PROGETTI DI RICERCA DI RILEVANTE INTERESSE NAZIONALE-Bando 2022, Prot. P2022WSJKS. G.C., F.S. and M.G. acknowledge financial support from the CRESCENDO project, PRIN: PROGETTI DI RICERCA DI RILEVANTE INTERESSE NAZIONALE-Bando 2022, Prot. 2022HL9PRP. G.C. and M.G. acknowledge financial support from the European Union – Next Generation EU – PNRR – M4C2, investimento 1.1 – “Fondo PRIN 2022” – “Understanding the pHotochemistry of sulfur substituted dnA bases by advanced ultrafast spectroscopy for phototherapeutic applications (HAPPY) ID P20224AWLB – CUP D53D23016720001”.

Author contributions

P.K. build the experimental setup, V.P., P.K., L.U. performed the measurements, V.P., P.K. analyzed and interpreted the experimental data. V.K.J, D.A., F.M., and S.U. performed the computational simulations. V.K.J., D.A., F.S., M.G., S.M., and A.N. analyzed and discussed the computational simulations. D.A., S.M., A.N., M.G., and G.C. procured funding for the work. F.S., A.N., M.G., and G.C. supervised the project. V.K.J. and A.N. wrote the manuscript with inputs from all authors.

Peer review

Peer review information

Nature Communications thanks Dongho Kim and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Data availability

The QM/MM optimized structures for the clusters representatives employed in this work, the solvent configurations for the dynamics displayed in Fig. 2 and the parameters used for the Hamiltonian employed in the MLMCTDH wavepacket dynamics have been deposited in the Zenodo database under the accession code 10987710. Source data for Figs. 1–4 are provided in the Source Data file82 (10.6084/m9.figshare.25640532). Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Daniel Aranda Ruiz, Vasilis Petropoulos, Piotr Kabaciński.

Change history

9/10/2024

A Correction to this paper has been published: 10.1038/s41467-024-52280-y
==== Refs
References

1. Ashkenazi G Kosloff R Ratner MA Photoexcited electron transfer:  short-time dynamics and turnover control by dephasing, relaxation, and mixing J. Am. Chem. Soc. 1999 121 3386 3395 10.1021/ja981998p
Ashkenazi, G., Kosloff, R. & Ratner, M. A. Photoexcited electron transfer:  short-time dynamics and turnover control by dephasing, relaxation, and mixing. J. Am. Chem. Soc. 121, 3386–3395 (1999).10.1021/ja981998p
2. Xu D Schulten K Coupling of protein motion to electron transfer in a photosynthetic reaction center: investigating the low temperature behavior in the framework of the spin—boson model Chem. Phys. 1994 182 91 117 10.1016/0301-0104(94)00016-6
Xu, D. & Schulten, K. Coupling of protein motion to electron transfer in a photosynthetic reaction center: investigating the low temperature behavior in the framework of the spin—boson model. Chem. Phys. 182, 91–117 (1994).10.1016/0301-0104(94)00016-6
3. Ando K Sumi H Nonequilibrium oscillatory electron transfer in bacterial photosynthesis J. Phys. Chem. B 1998 102 10991 11000 10.1021/jp982659l
Ando, K. & Sumi, H. Nonequilibrium oscillatory electron transfer in bacterial photosynthesis. J. Phys. Chem. B 102, 10991–11000 (1998).10.1021/jp982659l
4. Olaya-Castro A Scholes GD Energy transfer from Förster–Dexter theory to quantum coherent light-harvesting Int. Rev. Phys. Chem. 2011 30 49 77 10.1080/0144235X.2010.537060
Olaya-Castro, A. & Scholes, G. D. Energy transfer from Förster–Dexter theory to quantum coherent light-harvesting. Int. Rev. Phys. Chem. 30, 49–77 (2011).10.1080/0144235X.2010.537060
5. Simpson WT Peterson DL Coupling strength for resonance force transfer of electronic energy in van der waals solids J. Chem. Phys. 1957 26 588 593 10.1063/1.1743351
Simpson, W. T. & Peterson, D. L. Coupling strength for resonance force transfer of electronic energy in van der waals solids. J. Chem. Phys. 26, 588–593 (1957).10.1063/1.1743351
6. Robinson GW Frosch RP Electronic excitation transfer and relaxation J. Chem. Phys. 1963 38 1187 1203 10.1063/1.1733823
Robinson, G. W. & Frosch, R. P. Electronic excitation transfer and relaxation. J. Chem. Phys. 38, 1187–1203 (1963).10.1063/1.1733823
7. Scholes GD Ghiggino KP Rate expressions for excitation transfer I. Radiationless transition theory perspective J. Chem. Phys. 1994 101 1251 1261 10.1063/1.467817
Scholes, G. D. & Ghiggino, K. P. Rate expressions for excitation transfer I. Radiationless transition theory perspective. J. Chem. Phys. 101, 1251–1261 (1994).10.1063/1.467817
8. Clegg, R. M., Sener, M. & Govindjee. From Förster resonance energy transfer to coherent resonance energy transfer and back. Proc. SPIE 7561, Optical Biopsy VII, 75610C (2010).
9. Mendive-Tapia D Multidimensional quantum mechanical modeling of electron transfer and electronic coherence in plant cryptochromes: the role of initial bath conditions J. Phys. Chem. B 2018 122 126 136 10.1021/acs.jpcb.7b10412 29216421
Mendive-Tapia, D. et al. Multidimensional quantum mechanical modeling of electron transfer and electronic coherence in plant cryptochromes: the role of initial bath conditions. J. Phys. Chem. B 122, 126–136 (2018).29216421 10.1021/acs.jpcb.7b10412
10. Engel GS Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems Nature 2007 446 782 786 10.1038/nature05678 17429397
Engel, G. S. et al. Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature 446, 782–786 (2007).17429397 10.1038/nature05678
11. Collini E Coherently wired light-harvesting in photosynthetic marine algae at ambient temperature Nature 2010 463 644 647 10.1038/nature08811 20130647
Collini, E. et al. Coherently wired light-harvesting in photosynthetic marine algae at ambient temperature. Nature 463, 644–647 (2010).20130647 10.1038/nature08811
12. Panitchayangkoon G Long-lived quantum coherence in photosynthetic complexes at physiological temperature Proc. Natl Acad. Sci. USA 2010 107 12766 12770 10.1073/pnas.1005484107 20615985
Panitchayangkoon, G. et al. Long-lived quantum coherence in photosynthetic complexes at physiological temperature. Proc. Natl Acad. Sci. USA 107, 12766–12770 (2010).20615985 10.1073/pnas.1005484107
13. Duan H-G Nature does not rely on long-lived electronic quantum coherence for photosynthetic energy transfer Proc. Natl Acad. Sci. 2017 114 8493 8498 10.1073/pnas.1702261114 28743751
Duan, H.-G. et al. Nature does not rely on long-lived electronic quantum coherence for photosynthetic energy transfer. Proc. Natl Acad. Sci. 114, 8493–8498 (2017).28743751 10.1073/pnas.1702261114
14. Arsenault EA Yoneda Y Iwai M Niyogi KK Fleming GR The role of mixed vibronic Qy-Qx states in green light absorption of light-harvesting complex II Nat. Commun. 2020 11 6011 10.1038/s41467-020-19800-y 33243997
Arsenault, E. A., Yoneda, Y., Iwai, M., Niyogi, K. K. & Fleming, G. R. The role of mixed vibronic Qy-Qx states in green light absorption of light-harvesting complex II. Nat. Commun. 11, 6011 (2020).33243997 10.1038/s41467-020-19800-y
15. Halpin A Two-dimensional spectroscopy of a molecular dimer unveils the effects of vibronic coupling on exciton coherences Nat. Chem. 2014 6 196 201 10.1038/nchem.1834 24557133
Halpin, A. et al. Two-dimensional spectroscopy of a molecular dimer unveils the effects of vibronic coupling on exciton coherences. Nat. Chem. 6, 196–201 (2014).24557133 10.1038/nchem.1834
16. Sil S Excitation energy transfer and vibronic coherence in intact phycobilisomes Nat. Chem. 2022 14 1286 1294 10.1038/s41557-022-01026-8 36123451
Sil, S. et al. Excitation energy transfer and vibronic coherence in intact phycobilisomes. Nat. Chem. 14, 1286–1294 (2022).36123451 10.1038/s41557-022-01026-8
17. Fuller FD Vibronic coherence in oxygenic photosynthesis Nat. Chem. 2014 6 706 711 10.1038/nchem.2005 25054941
Fuller, F. D. et al. Vibronic coherence in oxygenic photosynthesis. Nat. Chem. 6, 706–711 (2014).25054941 10.1038/nchem.2005
18. Lin C Kim T Schultz JD Young RM Wasielewski MR Accelerating symmetry-breaking charge separation in a perylenediimide trimer through a vibronically coherent dimer intermediate Nat. Chem. 2022 14 786 793 10.1038/s41557-022-00927-y 35469005
Lin, C., Kim, T., Schultz, J. D., Young, R. M. & Wasielewski, M. R. Accelerating symmetry-breaking charge separation in a perylenediimide trimer through a vibronically coherent dimer intermediate. Nat. Chem. 14, 786–793 (2022).35469005 10.1038/s41557-022-00927-y
19. Romero E Quantum coherence in photosynthesis for efficient solar-energy conversion Nat. Phys. 2014 10 676 682 10.1038/nphys3017 26870153
Romero, E. et al. Quantum coherence in photosynthesis for efficient solar-energy conversion. Nat. Phys. 10, 676–682 (2014).26870153 10.1038/nphys3017
20. Perlík V Vibronic coupling explains the ultrafast carotenoid-to-bacteriochlorophyll energy transfer in natural and artificial light harvesters J. Chem. Phys. 2015 142 212434 10.1063/1.4919548 26049454
Perlík, V. et al. Vibronic coupling explains the ultrafast carotenoid-to-bacteriochlorophyll energy transfer in natural and artificial light harvesters. J. Chem. Phys. 142, 212434 (2015).26049454 10.1063/1.4919548
21. Jumper CC van Stokkum IHM Mirkovic T Scholes GD Vibronic wavepackets and energy transfer in cryptophyte light-harvesting complexes J. Phys. Chem. B 2018 122 6328 6340 10.1021/acs.jpcb.8b02629 29847127
Jumper, C. C., van Stokkum, I. H. M., Mirkovic, T. & Scholes, G. D. Vibronic wavepackets and energy transfer in cryptophyte light-harvesting complexes. J. Phys. Chem. B 122, 6328–6340 (2018).29847127 10.1021/acs.jpcb.8b02629
22. Forster Th. energiewanderung und fluoreszenz Naturwissenschaften 1946 33 166 175 10.1007/BF00585226
Forster Th. energiewanderung und fluoreszenz. Naturwissenschaften 33, 166–175 (1946).10.1007/BF00585226
23. Theodor Förster. Energy migration and fluorescence J. Biomed. Opt. 2012 17 011002 10.1117/1.JBO.17.1.011002 22352636
Theodor Förster. Energy migration and fluorescence. J. Biomed. Opt. 17, 011002 (2012).22352636 10.1117/1.JBO.17.1.011002
24. Harcourt RD Scholes GD Ghiggino KP Rate expressions for excitation transfer. II. Electronic considerations of direct and through–configuration exciton resonance interactions J. Chem. Phys. 1994 101 10521 10525 10.1063/1.467869
Harcourt, R. D., Scholes, G. D. & Ghiggino, K. P. Rate expressions for excitation transfer. II. Electronic considerations of direct and through–configuration exciton resonance interactions. J. Chem. Phys. 101, 10521–10525 (1994).10.1063/1.467869
25. Feng X Luzanov AV Krylov AI Fission of entangled spins: an electronic structure perspective J. Phys. Chem. Lett. 2013 4 3845 3852 10.1021/jz402122m
Feng, X., Luzanov, A. V. & Krylov, A. I. Fission of entangled spins: an electronic structure perspective. J. Phys. Chem. Lett. 4, 3845–3852 (2013).10.1021/jz402122m
26. Ghukasyan, V. V. & Heikal, A. A. (Eds.) Natural Biomarkers for Cellular Metabolism: Biology, Techniques, and Applications. (1st ed.). (CRC Press, London, 2014).
27. Ranjit S Lanzanò L Libby AE Gratton E Levi M Advances in fluorescence microscopy techniques to study kidney function Nat. Rev. Nephrol. 2021 17 128 144 10.1038/s41581-020-00337-8 32948857
Ranjit, S., Lanzanò, L., Libby, A. E., Gratton, E. & Levi, M. Advances in fluorescence microscopy techniques to study kidney function. Nat. Rev. Nephrol. 17, 128–144 (2021).32948857 10.1038/s41581-020-00337-8
28. Lakowicz JR Szmacinski H Nowaczyk K Johnson ML Fluorescence lifetime imaging of free and protein-bound NADH Proc. Natl Acad. Sci. 1992 89 1271 1275 10.1073/pnas.89.4.1271 1741380
Lakowicz, J. R., Szmacinski, H., Nowaczyk, K. & Johnson, M. L. Fluorescence lifetime imaging of free and protein-bound NADH. Proc. Natl Acad. Sci. 89, 1271–1275 (1992).1741380 10.1073/pnas.89.4.1271
29. Tropp J Redfield AG Proton magnetic resonance of NADH in water-methanol mixtures. Conformational change and behavior of exchangeable proton resonances as a function of temperature J. Am. Chem. Soc. 1980 102 534 538 10.1021/ja00522a016
Tropp, J. & Redfield, A. G. Proton magnetic resonance of NADH in water-methanol mixtures. Conformational change and behavior of exchangeable proton resonances as a function of temperature. J. Am. Chem. Soc. 102, 534–538 (1980).10.1021/ja00522a016
30. Hull RV Conger PS Hoobler RJ Conformation of NADH studied by fluorescence excitation transfer spectroscopy Biophys. Chem. 2001 90 9 16 10.1016/S0301-4622(00)00239-8 11321678
Hull, R. V., Conger, P. S. & Hoobler, R. J. Conformation of NADH studied by fluorescence excitation transfer spectroscopy. Biophys. Chem. 90, 9–16 (2001).11321678 10.1016/S0301-4622(00)00239-8
31. Li H Ultrafast fluorescence dynamics of NADH in aprotic solvents: Quasi-static self-quenching unmasked J. Photochem. Photobiol. Chem. 2023 436 114384 10.1016/j.jphotochem.2022.114384
Li, H. et al. Ultrafast fluorescence dynamics of NADH in aprotic solvents: Quasi-static self-quenching unmasked. J. Photochem. Photobiol. Chem. 436, 114384 (2023).10.1016/j.jphotochem.2022.114384
32. Cao S Femtosecond fluorescence spectra of NADH in solution: ultrafast solvation dynamics J. Phys. Chem. B 2020 124 771 776 10.1021/acs.jpcb.9b10656 31941277
Cao, S. et al. Femtosecond fluorescence spectra of NADH in solution: ultrafast solvation dynamics. J. Phys. Chem. B 124, 771–776 (2020).31941277 10.1021/acs.jpcb.9b10656
33. Reza MM Primary photophysics of nicotinamide chromophores in their oxidized and reduced forms J. Phys. Chem. B 2023 127 8432 8445 10.1021/acs.jpcb.3c03246 37733881
Reza, M. M. et al. Primary photophysics of nicotinamide chromophores in their oxidized and reduced forms. J. Phys. Chem. B 127, 8432–8445 (2023).37733881 10.1021/acs.jpcb.3c03246
34. Weber G Intramolecular transfer of electronic energy in dihydro diphosphopyridine nucleotide Nature 1957 180 1409 1409 10.1038/1801409a0
Weber, G. Intramolecular transfer of electronic energy in dihydro diphosphopyridine nucleotide. Nature 180, 1409–1409 (1957).10.1038/1801409a0
35. Freed S Neyfakh EA Tumerman LA Influence of solvents on the intramolecular energy transfer in NADH and NADPH Biochim. Biophys. Acta BBA - Bioenerg. 1967 143 432 434 10.1016/0005-2728(67)90099-0
Freed, S., Neyfakh, E. A. & Tumerman, L. A. Influence of solvents on the intramolecular energy transfer in NADH and NADPH. Biochim. Biophys. Acta BBA - Bioenerg. 143, 432–434 (1967).10.1016/0005-2728(67)90099-0
36. Heiner Z Roland T Leonard J Haacke S Groma GI Kinetics of light-induced intramolecular energy transfer in different conformational states of NADH J. Phys. Chem. B 2017 121 8037 8045 10.1021/acs.jpcb.7b04753 28731346
Heiner, Z., Roland, T., Leonard, J., Haacke, S. & Groma, G. I. Kinetics of light-induced intramolecular energy transfer in different conformational states of NADH. J. Phys. Chem. B 121, 8037–8045 (2017).28731346 10.1021/acs.jpcb.7b04753
37. Cadena-Caicedo A Gonzalez-Cano B López-Arteaga R Esturau-Escofet N Peon J Ultrafast fluorescence signals from β-dihydronicotinamide adenine dinucleotide: resonant energy transfer in the folded and unfolded forms J. Phys. Chem. B 2020 124 519 530 10.1021/acs.jpcb.9b10012 31876417
Cadena-Caicedo, A., Gonzalez-Cano, B., López-Arteaga, R., Esturau-Escofet, N. & Peon, J. Ultrafast fluorescence signals from β-dihydronicotinamide adenine dinucleotide: resonant energy transfer in the folded and unfolded forms. J. Phys. Chem. B 124, 519–530 (2020).31876417 10.1021/acs.jpcb.9b10012
38. Kulkarni GV Ray A Patel CC Molecular orbital studies on nicotinic acid, nicotinamide and related compounds J. Mol. Struct. 1978 49 373 382 10.1016/0022-2860(78)87275-5
Kulkarni, G. V., Ray, A. & Patel, C. C. Molecular orbital studies on nicotinic acid, nicotinamide and related compounds. J. Mol. Struct. 49, 373–382 (1978).10.1016/0022-2860(78)87275-5
39. Fülscher MP Serrano-Andrés L Roos BO A theoretical study of the electronic spectra of adenine and guanine J. Am. Chem. Soc. 1997 119 6168 6176 10.1021/ja964426i
Fülscher, M. P., Serrano-Andrés, L. & Roos, B. O. A theoretical study of the electronic spectra of adenine and guanine. J. Am. Chem. Soc. 119, 6168–6176 (1997).10.1021/ja964426i
40. Cohen B Hare PM Kohler B Ultrafast excited-state dynamics of adenine and monomethylated adenines in solution:  implications for the nonradiative decay mechanism J. Am. Chem. Soc. 2003 125 13594 13601 10.1021/ja035628z 14583057
Cohen, B., Hare, P. M. & Kohler, B. Ultrafast excited-state dynamics of adenine and monomethylated adenines in solution:  implications for the nonradiative decay mechanism. J. Am. Chem. Soc. 125, 13594–13601 (2003).14583057 10.1021/ja035628z
41. Conti I Garavelli M Orlandi G Deciphering low energy deactivation channels in adenine J. Am. Chem. Soc. 2009 131 16108 16118 10.1021/ja902311y 19845361
Conti, I., Garavelli, M. & Orlandi, G. Deciphering low energy deactivation channels in adenine. J. Am. Chem. Soc. 131, 16108–16118 (2009).19845361 10.1021/ja902311y
42. Brøndsted Nielsen S Sølling TI Are conical intersections responsible for the ultrafast processes of adenine, protonated adenine, and the corresponding nucleosides? ChemPhysChem 2005 6 1276 1281 10.1002/cphc.200400644 15929162
Brøndsted Nielsen, S. & Sølling, T. I. Are conical intersections responsible for the ultrafast processes of adenine, protonated adenine, and the corresponding nucleosides? ChemPhysChem 6, 1276–1281 (2005).15929162 10.1002/cphc.200400644
43. Evans NL Ullrich S Wavelength dependence of electronic relaxation in isolated adenine using uv femtosecond time-resolved photoelectron spectroscopy J. Phys. Chem. A 2010 114 11225 11230 10.1021/jp1029097 20961159
Evans, N. L. & Ullrich, S. Wavelength dependence of electronic relaxation in isolated adenine using uv femtosecond time-resolved photoelectron spectroscopy. J. Phys. Chem. A 114, 11225–11230 (2010).20961159 10.1021/jp1029097
44. Kang H Jung B Kim SK Mechanism for ultrafast internal conversion of adenine J. Chem. Phys. 2003 118 6717 6719 10.1063/1.1566438
Kang, H., Jung, B. & Kim, S. K. Mechanism for ultrafast internal conversion of adenine. J. Chem. Phys. 118, 6717–6719 (2003).10.1063/1.1566438
45. Gustavsson T Sharonov A Onidas D Markovitsi D Adenine, deoxyadenosine and deoxyadenosine 5′-monophosphate studied by femtosecond fluorescence upconversion spectroscopy Chem. Phys. Lett. 2002 356 49 54 10.1016/S0009-2614(02)00290-7
Gustavsson, T., Sharonov, A., Onidas, D. & Markovitsi, D. Adenine, deoxyadenosine and deoxyadenosine 5′-monophosphate studied by femtosecond fluorescence upconversion spectroscopy. Chem. Phys. Lett. 356, 49–54 (2002).10.1016/S0009-2614(02)00290-7
46. Improta R Santoro F Blancafort L Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases Chem. Rev. 2016 116 3540 3593 10.1021/acs.chemrev.5b00444 26928320
Improta, R., Santoro, F. & Blancafort, L. Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases. Chem. Rev. 116, 3540–3593 (2016).26928320 10.1021/acs.chemrev.5b00444
47. Stange UC Temps F Ultrafast electronic deactivation of UV-excited adenine and its ribo- and deoxyribonucleosides and -nucleotides: A comparative study. Ultrafast Photoinduced Process Polyat. Mol. Struct. Dyn. Spectrosc. Dedic. Wolfgang Domcke Occas. His 70th Birthd. 2018 515 441 451
Stange, U. C. & Temps, F. Ultrafast electronic deactivation of UV-excited adenine and its ribo- and deoxyribonucleosides and -nucleotides: A comparative study. Ultrafast Photoinduced Process. Polyat. Mol. Struct. Dyn. Spectrosc. Dedic. Wolfgang Domcke Occas. His 70th Birthd. 515, 441–451 (2018).
48. Meyer H-D Manthe U Cederbaum LS The multi-configurational time-dependent Hartree approach Chem. Phys. Lett. 1990 165 73 78 10.1016/0009-2614(90)87014-I
Meyer, H.-D., Manthe, U. & Cederbaum, L. S. The multi-configurational time-dependent Hartree approach. Chem. Phys. Lett. 165, 73–78 (1990).10.1016/0009-2614(90)87014-I
49. Beck MH Jäckle A Worth GA Meyer H-D The multiconfiguration time-dependent Hartree (MCTDH) method: a highly efficient algorithm for propagating wavepackets Phys. Rep. 2000 324 1 105 10.1016/S0370-1573(99)00047-2
Beck, M. H., Jäckle, A., Worth, G. A. & Meyer, H.-D. The multiconfiguration time-dependent Hartree (MCTDH) method: a highly efficient algorithm for propagating wavepackets. Phys. Rep. 324, 1–105 (2000).10.1016/S0370-1573(99)00047-2
50. Meyer, H.-D., Gatti, F. & Worth, G. A. Multidimensional Quantum Dynamics: MCTDH Theory and Applications. (Wiley-VCH, Weinheim, 2009).
51. Meyer H-D Worth GA Quantum molecular dynamics: propagating wavepackets and density operators using the multiconfiguration time-dependent Hartree method Theor. Chem. Acc. 2003 109 251 267 10.1007/s00214-003-0439-1
Meyer, H.-D. & Worth, G. A. Quantum molecular dynamics: propagating wavepackets and density operators using the multiconfiguration time-dependent Hartree method. Theor. Chem. Acc. 109, 251–267 (2003).10.1007/s00214-003-0439-1
52. Worth GA Meyer H-D Köppel H Cederbaum LS Burghardt I Using the MCTDH wavepacket propagation method to describe multimode non-adiabatic dynamics Int. Rev. Phys. Chem. 2008 27 569 606 10.1080/01442350802137656
Worth, G. A., Meyer, H.-D., Köppel, H., Cederbaum, L. S. & Burghardt, I. Using the MCTDH wavepacket propagation method to describe multimode non-adiabatic dynamics. Int. Rev. Phys. Chem. 27, 569–606 (2008).10.1080/01442350802137656
53. Wang H Thoss M Multilayer formulation of the multiconfiguration time-dependent Hartree theory J. Chem. Phys. 2003 119 1289 1299 10.1063/1.1580111
Wang, H. & Thoss, M. Multilayer formulation of the multiconfiguration time-dependent Hartree theory. J. Chem. Phys. 119, 1289–1299 (2003).10.1063/1.1580111
54. Manthe U A multilayer multiconfigurational time-dependent Hartree approach for quantum dynamics on general potential energy surfaces J. Chem. Phys. 2008 128 164116 10.1063/1.2902982 18447430
Manthe, U. A multilayer multiconfigurational time-dependent Hartree approach for quantum dynamics on general potential energy surfaces. J. Chem. Phys. 128, 164116 (2008).18447430 10.1063/1.2902982
55. Manthe U Layered discrete variable representations and their application within the multiconfigurational time-dependent Hartree approach J. Chem. Phys. 2009 130 054109 10.1063/1.3069655 19206960
Manthe, U. Layered discrete variable representations and their application within the multiconfigurational time-dependent Hartree approach. J. Chem. Phys. 130, 054109 (2009).19206960 10.1063/1.3069655
56. Vendrell O Meyer H-D Multilayer multiconfiguration time-dependent Hartree method: Implementation and applications to a Henon–Heiles Hamiltonian and to pyrazine J. Chem. Phys. 2011 134 044135 10.1063/1.3535541 21280715
Vendrell, O. & Meyer, H.-D. Multilayer multiconfiguration time-dependent Hartree method: Implementation and applications to a Henon–Heiles Hamiltonian and to pyrazine. J. Chem. Phys. 134, 044135 (2011).21280715 10.1063/1.3535541
57. Köppel, H., Domcke, W. & Cederbaum, L. S. Multimode molecular dynamics beyond the born-oppenheimer approximation. in Adv. Chem. Phys. 57, 59–246 (1984).
58. Köppel, H. Diabatic Representation: Methods For The Construction Of Diabatic Electronic States. in Conical Intersection s 15 175–204 (World Scientific, 2004).
59. Shiozaki T Győrffy W Celani P Werner H-J Communication: extended multi-state complete active space second-order perturbation theory: Energy and nuclear gradients J. Chem. Phys. 2011 135 81106 10.1063/1.3633329
Shiozaki, T., Győrffy, W., Celani, P. & Werner, H.-J. Communication: extended multi-state complete active space second-order perturbation theory: Energy and nuclear gradients. J. Chem. Phys. 135, 81106 (2011).10.1063/1.3633329
60. Sugita Y Okamoto Y Replica-exchange molecular dynamics method for protein folding Chem. Phys. Lett. 1999 314 141 151 10.1016/S0009-2614(99)01123-9
Sugita, Y. & Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 314, 141–151 (1999).10.1016/S0009-2614(99)01123-9
61. Warshel A Levitt M Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme J. Mol. Biol. 1976 103 227 249 10.1016/0022-2836(76)90311-9 985660
Warshel, A. & Levitt, M. Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. J. Mol. Biol. 103, 227–249 (1976).985660 10.1016/0022-2836(76)90311-9
62. Crespo-Hernández CE Cohen B Hare PM Kohler B Ultrafast excited-state dynamics in nucleic acids Chem. Rev. 2004 104 1977 2020 10.1021/cr0206770 15080719
Crespo-Hernández, C. E., Cohen, B., Hare, P. M. & Kohler, B. Ultrafast excited-state dynamics in nucleic acids. Chem. Rev. 104, 1977–2020 (2004).15080719 10.1021/cr0206770
63. Penfold TJ Gindensperger E Daniel C Marian CM Spin-vibronic mechanism for intersystem crossing Chem. Rev. 2018 118 6975 7025 10.1021/acs.chemrev.7b00617 29558159
Penfold, T. J., Gindensperger, E., Daniel, C. & Marian, C. M. Spin-vibronic mechanism for intersystem crossing. Chem. Rev. 118, 6975–7025 (2018).29558159 10.1021/acs.chemrev.7b00617
64. Yamazaki S Kato S Solvent effect on conical intersections in excited-state 9h-adenine:  radiationless decay mechanism in polar solvent J. Am. Chem. Soc. 2007 129 2901 2909 10.1021/ja0669169 17298065
Yamazaki, S. & Kato, S. Solvent effect on conical intersections in excited-state 9h-adenine:  radiationless decay mechanism in polar solvent. J. Am. Chem. Soc. 129, 2901–2909 (2007).17298065 10.1021/ja0669169
65. Tuna D Sobolewski AL Domcke W Mechanisms of ultrafast excited-state deactivation in adenosine J. Phys. Chem. A 2014 118 122 127 10.1021/jp410121h 24320624
Tuna, D., Sobolewski, A. L. & Domcke, W. Mechanisms of ultrafast excited-state deactivation in adenosine. J. Phys. Chem. A 118, 122–127 (2014).24320624 10.1021/jp410121h
66. Jaiswal VK Environment-driven coherent population transfer governs the ultrafast photophysics of tryptophan J. Am. Chem. Soc. 2022 144 12884 12892 10.1021/jacs.2c04565 35796759
Jaiswal, V. K. et al. Environment-driven coherent population transfer governs the ultrafast photophysics of tryptophan. J. Am. Chem. Soc. 144, 12884–12892 (2022).35796759 10.1021/jacs.2c04565
67. Jimenez R Fleming GR Kumar PV Maroncelli M Femtosecond solvation dynamics of water Nature 1994 369 471 473 10.1038/369471a0
Jimenez, R., Fleming, G. R., Kumar, P. V. & Maroncelli, M. Femtosecond solvation dynamics of water. Nature 369, 471–473 (1994).10.1038/369471a0
68. Cerezo J Mixed quantum/classical method for nonadiabatic quantum dynamics in explicit solvent models: the ππ*/nπ* decay of thymine in water as a test case J. Chem. Theory Comput. 2018 14 820 832 10.1021/acs.jctc.7b01015 29207245
Cerezo, J. et al. Mixed quantum/classical method for nonadiabatic quantum dynamics in explicit solvent models: the ππ*/nπ* decay of thymine in water as a test case. J. Chem. Theory Comput. 14, 820–832 (2018).29207245 10.1021/acs.jctc.7b01015
69. Borrego-Varillas R Tracking excited state decay mechanisms of pyrimidine nucleosides in real time Nat. Commun. 2021 12 7285 10.1038/s41467-021-27535-7 34907186
Borrego-Varillas, R. et al. Tracking excited state decay mechanisms of pyrimidine nucleosides in real time. Nat. Commun. 12, 7285 (2021).34907186 10.1038/s41467-021-27535-7
70. Rouxel JR Coupled electronic and nuclear motions during azobenzene photoisomerization monitored by ultrafast electron diffraction J. Chem. Theory Comput. 2022 18 605 613 10.1021/acs.jctc.1c00792 35073085
Rouxel, J. R. et al. Coupled electronic and nuclear motions during azobenzene photoisomerization monitored by ultrafast electron diffraction. J. Chem. Theory Comput. 18, 605–613 (2022).35073085 10.1021/acs.jctc.1c00792
71. Borrego-Varillas R Observation of the sub-100 femtosecond population of a dark state in a thiobase mediating intersystem crossing J. Am. Chem. Soc. 2018 140 16087 16093 10.1021/jacs.8b07057 30380844
Borrego-Varillas, R. et al. Observation of the sub-100 femtosecond population of a dark state in a thiobase mediating intersystem crossing. J. Am. Chem. Soc. 140, 16087–16093 (2018).30380844 10.1021/jacs.8b07057
72. Borrego-Varillas R Ganzer L Cerullo G Manzoni C Ultraviolet transient absorption spectrometer with sub-20-fs time resolution Appl. Sci. 2018 8 989 10.3390/app8060989
Borrego-Varillas, R., Ganzer, L., Cerullo, G. & Manzoni, C. Ultraviolet transient absorption spectrometer with sub-20-fs time resolution. Appl. Sci. 8, 989 (2018).10.3390/app8060989
73. Weingart O COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations J. Mol. Model. 2018 24 271 10.1007/s00894-018-3769-6 30178229
Weingart, O. et al. COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations. J. Mol. Model. 24, 271 (2018).30178229 10.1007/s00894-018-3769-6
74. Avagliano D Bonfanti M Nenov A Garavelli M Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM J. Comput. Chem. 2022 43 1641 1655 10.1002/jcc.26966 35815854
Avagliano, D., Bonfanti, M., Nenov, A. & Garavelli, M. Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM. J. Comput. Chem. 43, 1641–1655 (2022).35815854 10.1002/jcc.26966
75. Frisch, M. J. et al. Gaussian 16. Gaussian, Inc. (2016).
76. Li Manni, G. et al. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J. Chem. Theory Comput. 19, 6933–6991 (2023).
77. Plasser F Krylov AI Dreuw A libwfa: Wavefunction analysis tools for excited and open-shell electronic states WIREs Comput. Mol. Sci. 2022 12 e1595 10.1002/wcms.1595
Plasser, F., Krylov, A. I. & Dreuw, A. libwfa: Wavefunction analysis tools for excited and open-shell electronic states. WIREs Comput. Mol. Sci. 12, e1595 (2022).10.1002/wcms.1595
78. Plasser F Mewes SA Dreuw A González L Detailed wave function analysis for multireference methods: implementation in the molcas program package and applications to tetracene J. Chem. Theory Comput. 2017 13 5343 5353 10.1021/acs.jctc.7b00718 28972759
Plasser, F., Mewes, S. A., Dreuw, A. & González, L. Detailed wave function analysis for multireference methods: implementation in the molcas program package and applications to tetracene. J. Chem. Theory Comput. 13, 5343–5353 (2017).28972759 10.1021/acs.jctc.7b00718
79. Hsu C-P You Z-Q Chen H-C Characterization of the short-range couplings in excitation energy transfer J. Phys. Chem. C. 2008 112 1204 1212 10.1021/jp076512i
Hsu, C.-P., You, Z.-Q. & Chen, H.-C. Characterization of the short-range couplings in excitation energy transfer. J. Phys. Chem. C. 112, 1204–1212 (2008).10.1021/jp076512i
80. Worth, G. A., Beck, M. H., Jäckle, A. & Meyer, H.-D. The MCTDH Package, Version 8.2, (2000). H.-D. Meyer, Version 8.3 (2002), Version 8.4 (2007). O. Vendrell and H.-D. Meyer Version 8.5 (2013).
81. Worth, G. A. et al. The QUANTICS Package, Version 1.1, (University of Birmingham, Birmingham, U.K, 2015).
82. Jaiswal, V. K., et al. Sub-100-fs energy transfer in coenzyme NADH is a coherent process assisted by a charge-transfer state. Nat. Commun., 10.6084/m9.figshare.25640532. (2024).
