Significant impact of deprotonated status on the photoisomerization dynamics of bacteriophytochrome chromophore

Haiyi Huang Chao Xu Kunni Lin Jiawei Peng Feng Long Gu Zhenggang Lan

Citation:  Haiyi Huang, Chao Xu, Kunni Lin, Jiawei Peng, Feng Long Gu, Zhenggang Lan. Significant impact of deprotonated status on the photoisomerization dynamics of bacteriophytochrome chromophore[J]. Chinese Chemical Letters, 2023, 34(7): 107850. doi: 10.1016/j.cclet.2022.107850 shu

Significant impact of deprotonated status on the photoisomerization dynamics of bacteriophytochrome chromophore

English

  • Phytochromes are widely found in plants and microorganisms, which are responsible for fundamental photo-response functions [1-4]. They play crucial roles in a variety of light-regulated biological processes, ranging from seed germination, growth and flowering in plants to pigmentation and phototaxis in bacteria [1, 5, 6]. The phytochromes can be engineered into optical probes for the deep-tissue imaging of living cells, thus their photoinduced properties are considerable significance in photobiological researches [7-9].

    The chromophores in phytochromes, as the photoreaction centers, are mainly classified as three categories, namely biliverdin (BV), phytochromobilin (PΦB) and phycocyanobilin (PCB). Their structures are very similar and all share a methane-bridged open-chain conjugated tetrapyrrole structures with differences on the two terminal rings (Scheme 1). Their basic photochemical reaction mechanism is the reversible photo- and thermal-conversions between two isomers, i.e., the red-absorbing form Pr and the far-red-absorbing form Pfr. The primary step of the Pr-to-Pfr photoisomerization mechanism is governed by the nonadiabatic transitions. Once red-light absorption at ~660 nm, the chromophore in Pr state undergoes an ultrafast Z-to-E photoisomerization accompanied by the internal rotation of pyrrole rings, then follows conformational relaxations and finally leads to the formation of the Pfr isomer via the thermal processes. The reverse Pfr-to-Pr conversion experiences the E-to-Z isomerization upon the photoexcitation with far red light at ~730 nm. The forward and backward isomerization reactions induce the large-sized conformational evolution of the whole photoreceptor protein, giving regulating signals to adjust relevant biological functions. Therefore, their photoinduced reactions have been extensively investigated experimentally and theoretically [10-27].

    Scheme 1

    Scheme 1.  Chromophores, photoreaction centers in phytochromes. (1) R1: =CH-CH2-Cys, R2: -CH=CH2 for BV; (2) R1: -CH-(CH3)Cys, R2: -CH2-CH3 for PCB; (3) R1: -CH-(CH3)Cys, R2: -CH=CH2 for PΦB.

    Although numerous researches were conducted for many years, there are still several fundamental open questions in the photoinduced Pr-to-Pfr processes. The first critical argument is which molecular motion dominates the Z-to-E isomerization in the primary excited-state reaction. Several studies suggested that in protein environments the photoisomization takes place around the C15C16 methine bridge, resulting in the torsion of nearly-free side Ring D along this bond [11, 13, 28, 29]. However, different isomerization mechanisms were also proposed. For instance, Ulijasz et al. discovered that the internal rotation of the Ring A via the twisting motion of the C4C5 double bond triggers the excited-state dynamics of the PCB chromophore in photoreceptor protein [21]. It is also important to clarify the microscopic photoisomerization mechanism of the isolated chromophore at the all-atomic level. Some studies pointed out that the initial excited-state evolution of the isolated chromophore may be governed by the torsional motion at the C10C11 bond [20, 30, 31]. Zhuang et al. confirmed the pivotal role of the C10C11 double bond, and proposed that C9C10 single bond is also an key rotational degree of freedom in the photoisomerization of the isolated PΦB chromophore [20]. The existence of different photoisomerization mechanisms for the isolated chromophore and chromophore in phytochrome protein pocket is possibly due to the geometrical constrains added by the interactions between the chromophore and the surrounding protein. In this sense, the deep understanding of the role of the protein pocket in the modification of the photoinduced dynamics of phytochrome chromophore still requires additional efforts. Furthermore, hydrogen bonding plays a crucial role in many biological processes and accounts for many important physical properties of compounds [32-36]. Several works pointed out that the photochemistry of phytochrome chromophore is also strongly dependent on the hydrogen-bond network, as well as the protonation status of the pyrrole rings and surrounding residues [16, 22, 37-41]. A vibrational spectroscopic study on bathy phytochromes demonstrated that the proton transfer between the tetrapyrrole and a nearby histidine residue directly modules the photobiological functions by the generations of the photoactivated and deactivated states [22]. Osoegawa et al. revealed that the deprotonation of the Ring B/C distorts the methane bridge between the Ring B and the Ring C of the PCB chromophore, and such structural difference further adjusts the green or red light absorption [40]. Modi et al. clarified that the protonation/deprotonation status may regulate the light aborption of the BV chromophore in the protein pocket [24]. Kennis and co-workers proposed that the excited-state hydrogen transfer from the Ring A, B or C directly leads to the nonadiabatic decay to the ground state in the BV chromophore, resulting in the decreasing of the fluorescence signals [37]. Recently, Liu et al. demonstrated that the different orientations of the aspartate residue give the distinctive hydrogen-bond network connectivity that regulates the excited-state reaction channels [16]. As the summary, the protonation/deprotonation of the hydrogens attached to the pyrrole rings plays a critical role in the nonadiabatic dynamics of the phytochrome chromophore. It is interesting to clarify whether the removing of the hydrogen atom connected with the pyrrole nitrogen significantly modifies the nonadiabatic dynamics of the phytochrome or not. Particularly, it is also necessary to identify which key coordinates are responsible for the photoinduced dynamics of different phytochrome derivatives obtained by simple chemical modifications. These questions can be clearly addressed by the simulation of the nonadiabatic dynamics at all-atomic level, in which the active degrees of freedom in the excited-state dynamics are clarified by tracking the real-time molecular motion. In recent years, the on-the-fly ab initio nonadiabatic dynamics simulations were employed extensively to explore the ultrafast nonradiative reactions [42-47]. These simulations directly give the rich information, such as the time-dependent population dynamics, the nonadiabatic decay pathways and the essential molecular structure evolution.

    In present work, we wish to provide the detailed description for the photoisomerization of the isolated BV chromophore by using the on-the-fly nonadiabatic dynamics. The simulation was run with the trajectory surface hopping (TSH) approach [48] and the electronic structure calculations were performed at the semi-empirical OM2/MRCI level (the orthogonalization-corrected OM2 Hamiltonian combined with multireference configuration interaction) [49, 50]. Many studied have verified that the nonadiabatic dynamics at the OM2/MRCI level is a powerful simulation tool in the description of excited-state processes for media-sized and large polyatomic systems with the balance of computational accuracy and efficiency [51-56]. By following the excited-state molecular motions, we found that the protonation/deprotonation states of the pyrrole rings show the strong impact on the excited-state lifetime and reaction pathways. The underline reason is that the hydrogen bond network deeply modules the intramolecular motion, leading to different profiles of excited-state potential energy surfaces. As the consequence, the distinct population dynamics appears and different conical intersections are involved. Such finding confirms the important role of the hydrogen bond in the photoinduced nonadiabatic dynamics of the phytochrome chromophore. In other words, the understanding of the photoisomerization mechanism of phytochrome chromophere must take the hydrogen bond network into account, beside the steric effects caused by its surrounding protein environments. This work provides new insights into the photobiological functions of phytochrome chromophere and suggests the possible ideas to control its photoconversion processes for the further rational engineering of near-infrared fluorescent proteins.

    We focus on the nonadiabatic dynamics of the Pr conformation of the BV chromophore in DrBphP (Deinococcus radiodurans) bacterial phytochrome. The ZsZsZa (C5-Z, syn C10-Z, syn C15-Z, anti) stereoisomer (Fig. S1 in Supporting information) is chosen as our target configuration, which was identified as the dominant isomer by crystallography and NMR spectroscopy [57-59]. We built the simplified BV model to perform the nonadiabatic dynamics simulation, in which some side groups (-CH2COO and -Cys-S) were replaced by hydrogen atoms to reduce the computational cost. As it is rather difficult to assign the protonation status by experimental studies, we consider all possible single deprotonation derivatives. In order to systematically rationalize the influence of phytochrome chromophore protonation state on the photoinduced processes, the four most probable protonation models (Fig. 1 and Fig. S2 in Supporting information) were investigated, namely:

    Figure 1

    Figure 1.  OM2 optimized S0 minimum-energy structures of four most probable protonated BV-chromophore models. Selected bond lengths are given in Angstroms. See the Supporting information for Cartesian coordinates.

    (1) BV_N: all nitrogen atoms protonated;

    (2) BV_C: nitrogen in the Ring C deprotonated and the Rings A, B and D protonated;

    (3) BV_B: nitrogen in the Ring B deprotonated and the Rings A, C and D protonated;

    (4) BV_A: nitrogen in the Ring A deprotonated and the Rings B, C and D protonated.

    We first optimized the structures of the S0 and S1 minima (S0_min and S1_min) for the four protonated/deprotonated BV chromophore models (BV_N, BV_C, BV_B and BV_A), and the key internal-coordinate parameters of these geometries are listed in Tables S1-S4 (Supporting information). The most striking feature here is that the different protonation patterns clearly build distinguishing hydrogen-bond networks, as shown in Fig. 1, Figs. S3 and S4 (Supporting information). Here, the interaction region indicator (IRI) method was employed to investigate the intramolecular hydrogen bond (Fig. S3) [60]. For instance, at the S0_min, only a weak hydrogen bond between OA of the Ring A and HCD atom (bridged the Ring C and the Ring D) is found in BV_N. However, the more complex hydrogen-bonding interactions are observed due to the presence of the multi-hydrogen-bond connectivity in the other three models (BV_C, BV_B and BV_A), as the results of the deprotonation of different pyrrole ring.

    Fig. 2 presents the time-dependent fractional occupations of the S0 and S1 electronic states for the above four models, assuming that all trajectories start from the S1 state. Their lifetimes are far from each other and their corresponding time scales are even differed by the order of magnitude. The BV_N model with the saturated protonation displays the ultrafast excited-state population decay. The S1 population decay to ~50% around 400 fs and the S1 population becomes very low at about 1000 fs. As the contrast, the excited-state population decay of the BV_C model takes place at a much longer time scale and less than 50% of S1 population decays to the S1 state within 2.0 ps. The S1 lifetimes of the latter two models (BV_B and BV_A) are extremely long and only the very minor population decay is observed within the simulation time duration. The remarkably distinctive population decay behaviours among these four BV models imply that their excited-state reaction mechanism may be drastically different. In this scenario, two open questions naturally arise: (1) How to explain the different population dynamics in these BV models and do they display different decay channels? (2) How to correlate their decay mechanism and the protonation status? Next, we try to clarify both questions.

    Figure 2

    Figure 2.  Time-dependent fractional occupations of the S0 and S1 electronic states of the four BV models (BV_N, BV_C, BV_B and BV_A) in the non-adiabatic dynamics staring from the S1 state.

    We begin with the BV_N model to explore the mechanism of the nonadiabatic dynamics. By analysing all S0/S1 hopping geometries and the CI optimization results, we found that two S0/S1 minimum-energy conical intersections (CI1 and CI2 in Figs. 3a and b) govern the nonadiabatic transitions. The former CI1 is the main decay channel, since 93.6% trajectories decay via it, as shown in Fig. 3c. In the beginning of the photoinduced process, both HOMO and LUMO are delocalized on two central rings (Rings B and C) at the Frank-Condon region. When the torsional motion drives the system to access the CI1 and CI2, the HOMO and LUMO tend to become the localized one. As shown in Fig. S5 (Supporting information), when such torsional motion take place at the C10C11 or C9C10 bond, HOMO and LUMO start to locate in two moieties separated by these two bonds. This type of biradical or charge-transfer character widely exists in photoisomerization of double-bond conjugated systems [61].

    Figure 3

    Figure 3.  (a, b) OM2/MRCI optimized two S0/S1 conical intersections of the BV_N model. Cyan, red and gold represent S0_min, CI1 and CI2, respectively. The CI1 is characterized by the rotation of the whole C-D double-rings moiety around the C10C11 bond, while the CI2 is featured by the rotation of the whole A-B double-rings moiety around the C9C10 bond. (c) Distributions of two key dihedral angles at the initial and hopping geometries of the BV_N model. (d) Potential energy curves along the linear interpolated pathway from S0_min to the two conical intersections CI1 and CI2 of the BV_N model.

    The first conical intersection CI1 is characterized by the strong twisting motion of the C10C11 bond, along with the minor torsion of the C15C16 bond. From S0_min to CI1, the dihedral angle of C9C10C11NC between the Ring B and the Ring C changes from −12.4° to −85.4°, and the dihedral angle of C14C15C16ND between the Ring C and the Ring D also changes slightly from 6.4° to 19.0°. For the trajectories via this channel, the rotation of whole C-D double-rings moiety takes place in the nonadiabatic dynamics, while the minor torsion of the terminal Ring D with respect to Ring C also exists.

    The second conical intersection CI2 is characterized by the significant twisting of the C9C10 bond and also the minor torsion of the C5C6 bond. From S0_min to CI2, the dihedral angle of NBC9C10C11 between the Ring B and the Ring C increases from −16.5° to −93.4°, while the dihedral angle of C4C5C6NB between the Ring A and the Ring B shows the minor decreasing from 28.2° to 17.6°. In this decay channel, the obvious rotation of the whole A-B double-ring moiety and the minor torsion of the terminal Ring A exist.

    The reaction pathways from the S0_min to two conical intersections were built by employing the linear interpolated method. Fig. 3d predicts that the CI1 pathway should be more preferable, although both pathways are barrierless. Moreover, the S1_min optimization directly gives the CI1 geometry. These findings clearly explain why the CI1 channel is dominant in the nonadiabatic dynamics of the BV_N model (Fig. 3c).

    After the internal conversion of BV_N, the excessive kinetic energy leads to the high excitation of vibrational motions, and the system can easily overcome barriers on the electronic ground state to form various photoproducts. To explore possible reactions products, we ran the additional 5-ps Born-Oppenheimer molecular dynamics on the electronic ground state after the 2-ps nonadiabatic dynamics. By connecting the nonadiabatic and adiabatic dynamics, a few of representative trajectories are given in Fig. 4. Starting from the S1 state, the trajectories on the excited state move towards two S0/S1 CIs. Since the CI1 channel is dominant, we mainly show the geometry evolution of some representative trajectories decaying via it. At first, the trajectories propagate on the S1 state. Once approaching the CI1 where two electronic states show similar potential energies, the trajectories decay to the ground state. After the internal conversion, the initial BV_N isomer can be reformed though the reversed torsion of the whole C-D double-rings moiety, returning back to the starting BV_N configuration. Next, in some trajectories the vibrational energy flow may furtherly induce the internal rotation of the Ring D along the C15C16 bond, or the mixed Ring A and the Ring D rotations (Figs. S6-S8 in Supporting information). At the end of simulation, the ratio of these three channels is about 14.4%: 32.1%: 19.3%. Other isomers were also found while their probabilities are small (Fig. S9 in Supporting information). Since only 5.88% trajectories passing the CI2, it is not necessary to discuss their final photoproducts. The current finding clearly confirms that the rotation of other rings (such as Ring D) may be highly possible even when such motion is not dominant in the excited-state dynamics. Nevertheless, such multiple-step reaction mechanism provides a possible explanation on why different isomerization patterns were discovered in previous studies [11, 13, 20, 21, 28-31].

    Figure 4

    Figure 4.  Representative trajectories towards to various reaction pathways. The cyan and white structural skeleton represent the snapshot at the current and starting time steps, respectively. Top: the trajectory back to initial BV_N structure; Middle: the trajectory ends with the rotation of the Ring D; Bottom: the trajectory evolution shows the mixed the Ring A and the Ring D rotations.

    Next, let us consider the deprotonated BV_C model, which shows the longer excited state lifetime than the BV_N model (Fig. 2). Different to BV_N, no significant torsional motion at the single and double carbon and carbon bonds between the Ring B and the Ring C was observed for BV_C in the trajectories. This difference can be understood by the presence of the strong hydrogen bond between deprotonated NC and atom HB in BV_C, which locks these torsional motions. As the consequence, the preferred isomerization pathways in BV_N are not observed in the excited-state dynamics of BV_C.

    To clarify which motion is responsible for the nonadiabatic decay of the BV_C model, we collected the hopping geometries and located the key internal coordinates. As shown in Figs. 5a and b, two dihedral angles NCC14C15C16 and HBNBC6C9 are found to play the significant roles. Starting from these hopping geometries, the optimization of the minimum-energy S0/S1 CI geometry also confirms the importance of these two internal coordinates. Therefore, the synergy of the rotation of the Ring D and the out-of-plane motion of the HB atom associated with the Ring B (HBOOP) opens a new pathway in the nonadiabatic dynamics in the deprotonated BV_C model, different to the BV_N model.

    Figure 5

    Figure 5.  (a, b) Distributions of key dihedral angles at the initial and hopping points of the BV_C model. (c) Potential energy curves along the linear interpolated pathway from S0_min to the S1_min and from S1_min to the CI for the BV_C model. (d) OM2/MRCI optimized S0/S1 CI of the BV_C model. Cyan, gold and red represent S0_min, S1_min and CI, respectively.

    The excited state reaction pathway of the BV_C model was given in Fig. 5c. Starting from Frank-Condon region, there is a barrierless pathway to the S1_min. According to the linear-interpolated pathway from the S1_min to the minimum-energy geometry in the S0/S1 CI (Fig. 5d), we noticed that the CI lies much higher than the S1 minimum. Also considering the large size of the current chromophore, the vibrational relaxation on the excited state may lead to that the system tends to stay at the S1 minimum region for very long time before it enters the CI region. This may explain the much slower nonadiabatic decay of the BV_C model with respect to the BV_N model. At the same time, this CI is accessible from the S1_min only after the rather large geometrical adjustment on the excited state. Several motions must be involved, which include the significant HBOOP motion (∠HBNBC6C9 is 125.7° at CI, 159.9° at S1_min and 176.9° at S0_min), the breakdown of the strong hydrogen bond NCHB and the torsion of the NCC14C15C16 dihedral angle by ~10°. As the consequence, it takes a longer time for trajectories to access this CI region, resulting in the slow population decay of the BV_C model.

    Now we pay attention to the deprotonated BV_B and BV_A models, which show significantly longer excited state lifetime (Fig. 2). As the excited-state reaction pathways were clarified for BV_N and BV_C systems, it is necessary to examine whether these channels are possible in the BV_B and BV_A models.

    In BV_B model, the hydrogen-bond network is observed clearly (Fig. S3). Similar to the BV_C case, the NB-HC hydrogen bond limits the torsional motion between the Ring B and the Ring C, preventing the ultrafast decay channel observed in the BV_N model. Meanwhile, the rotational motion of the Ring D is also blocked by the hydrogen bond between atom OA and atom HC-D. Moreover, as the nitrogen atom in the Ring B is deprotonated, there is no HBOOP movement in this system. Therefore, BV_B cannot experience the similar decay channel as BV_C.

    In BV_A model, the hydrogen bond between the Ring A and the Ring B freezes the rotation of C4C5C6NB. However, the twisting motion of the C5C6 bond is clearly involved in the early stage of the nonadiabatic dynamics for the BV_N model (Fig. S10 in Supporting information), although the torsional motions at the bonds between the Ring B and the Ring C are dominant. As such torsional motion is fully locked in the BV_A model (Fig. S11 in Supporting information), the corresponding channel is closed. In addition, the same hydrogen bond also restricts the HBOOP motion of the Ring B, and the presence of the other strong hydrogen bond (S1_min, OAHC-D = 1.629 Å in Fig. S3) prevents the torsional motion of the Ring D. Therefore, the decay pathway in the BV_C model cannot exist for the BV_A model.

    Finally, we give some comments on the role of the hydrogen transfer in the nonadiabatic dynamics of the deprotonated BV_C, BV_B and BV_A models. Although the proton transfer may take place due to the existence of the hydrogen-bond network in the above three model (Table S5 in Supporting information), this motion is not important in the photoisomerization dynamics. Let us take the BV_C model with the highest nonadiabatic transition probability as an example, most hopping structures still remain as the BV_C configurations, and only 6.1% of them become the BV_B configurations. Even in this case, the existence of the hydrogen bond still locks the torsional motion at the single and double CC bonds between the Ring B and the Ring C. For the BV_A and BV_B models, the hydrogen transfer is also not important for the nonadiabatic decay dynamics and thus we do not repeat the discussions here. Overall, the proton transfer has very minor effects on the non-adiabatic transitions processes for the deprotonated models, possibly due to the fact that the overall hydrogen bond network is not destroyed after the proton transfer motion.

    In this work, we employed the on-the-fly surface hopping dynamics simulation at the semiempirical OM2/MRCI level to understand the photoinduced processes for the BV phytochrome chromophore. Several critical geometries, including electronic state minima and conical intersections, were obtained. We found that the deprotonation of different pyrrole nitrogen leads to distinct photophysical and photochemical behaviors, such as the substantial changes in the excited state lifetime and nonadiabatic decay pathways.

    The BV_N model shows the fastest excited-state decay because there are barrierless pathway from Franck-Condon region to the two CIs, and its excited-state process is dominated by the rotation of whole A-B double-rings moiety or whole C-D double-rings moiety. As the contrast, the excited-state population decay of the BV_C model takes place at a much longer time scale than that of BV_N model. The Ring D rotation and HOOP motion in the Ring B play efficient roles in the nonadiabatic process in the BV_C model. Interestingly, the other two models (BV_A and BV_B) show only extremely minor excited-state population decay within the simulation time duration.

    We attribute the striking difference to the formation of different intramolecular hydrogen bond connective patterns, which affect the preferable reaction pathways and accessible conical intersections, and give diverse decay channels from the excited state to the ground state. At the end, we emphasize that in addition to the protein steric effect, the variation in deprotonation state also greatly impacts on the excited-state dynamics of the phytochrome chromophore. This finding not only deepens our understanding of the primary photoinduced reaction of the phytochrome chromophore, but also provides new design ideas for the efficient engineering near-infrared fluorescent proteins based on the phytochrome. Certainly, it is highly interesting to examine whether the protein environments can modify the protonation states of the BV chromophores and the protein pocket can regulate the photoisomerization dynamics. Clearly, the nonadiabatic dynamics mechanism of different phytochrome chromophore (BV, PΦB and PCB) is strongly modified by hydrogen bond networks, as well as their interactions with the surrounding protein and solution environments. Such theoretical simulations may involve the further large scale QM/MM nonadiabatic dynamics simulations. This represents the challenging and on-going topic in the future study.

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    This work is supported by National Natural Science Foundation of China (NSFC, Nos. 21903030, 21933011 and 21873112).

    Theoretical methods, BV chromophore, S1_min of four models, typical trajectories, reaction products of the BV_N model, distributions of selected internal coordinates for all hopping geometries in the BV_N model, time-dependent evolution of the distributions of a few key dihedral angles in the nonadiabatic dynamics, configurations at hop events, benchmark of the vertical energies, reaction pathways in the OM2/MRCI(24, 22) level, cartesian coordinates and additional implementation details are available. Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2022.107850.


    1. [1]

      N.C. Rockwell, Y.S. Su, J.C. Lagarias, Annu. Rev. Plant Biol. 57 (2006) 837–858. doi: 10.1146/annurev.arplant.56.032604.144208

    2. [2]

      Z. Yu, R. Fischer, Nat. Rev. Microbiol. 17 (2019) 25–36. doi: 10.1038/s41579-018-0109-x

    3. [3]

      Y. Wang, G. Veglia, D. Zhong, J. Gao, Sci. Adv. 7 (2021) eabg3815. doi: 10.1126/sciadv.abg3815

    4. [4]

      C. Bongards, W. Gärtner, Acc. Chem. Res. 43 (2010) 485–495. doi: 10.1021/ar800133x

    5. [5]

      T. Kottke, A. Xie, D.S. Larsen, W.D. Hoff, Annu. Rev. Biophys. 47 (2018) 291–313. doi: 10.1146/annurev-biophys-070317-033047

    6. [6]

      W.R. Briggs, J.L. Spudich, Handbook of Photosensory Receptors, Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2005.

    7. [7]

      X. Shu, A. Royant, M.Z. Lin, et al., Science 324 (2009) 804–807. doi: 10.1126/science.1168683

    8. [8]

      K.G. Chernov, T.A. Redchuk, E.S. Omelina, V.V. Verkhusha, Chem. Rev. 117 (2017) 6423–6446. doi: 10.1021/acs.chemrev.6b00700

    9. [9]

      A.A. Kaberniuk, A.A. Shemetov, V.V. Verkhusha, Nat. Method 13 (2016) 591–597. doi: 10.1038/nmeth.3864

    10. [10]

      M.A. Mroginski, S. Adam, G.S. Amoyal, et al., Photochem. Photobiol. 97 (2021) 243–269. doi: 10.1111/php.13372

    11. [11]

      M.A. Mroginski, D.H. Murgida, P. Hildebrandt, Acc. Chem. Res. 40 (2007) 258–266. doi: 10.1021/ar6000523

    12. [12]

      P. Altoe, T. Climent, G.C. De Fusco, et al., J. Phys. Chem. B 113 (2009) 15067–15073. doi: 10.1021/jp904669x

    13. [13]

      X. Yang, Z. Ren, J. Kuk, K. Moffat, Nature 479 (2011) 428–432. doi: 10.1038/nature10506

    14. [14]

      J. Dasgupta, R.R. Frontiera, K.C. Taylor, J.C. Lagarias, R.A. Mathies, Proc. Natl. Acad. Sci. U. S. A. 106 (2009) 1784–1789. doi: 10.1073/pnas.0812056106

    15. [15]

      X. Xu, A. Höppner, C. Wiebeler, et al., Proc. Natl. Acad. Sci. U. S. A. 117 (2020) 2432–2440. doi: 10.1073/pnas.1910208117

    16. [16]

      X.Y. Liu, T.S. Zhang, Q. Fang, et al., Angew. Chem. 133 (2021) 18836–18841. doi: 10.1002/ange.202104853

    17. [17]

      D. Wang, Y. Qin, S. Zhang, et al., J. Phys. Chem. Lett. 10 (2019) 6197–6201. doi: 10.1021/acs.jpclett.9b02446

    18. [18]

      D. Wang, X. Li, L. Wang, X. Yang, D. Zhong, J. Phys. Chem. Lett. 11 (2020) 8819–8824. doi: 10.1021/acs.jpclett.0c02467

    19. [19]

      D. Wang, Y. Qin, M. Zhang, et al., J. Phys. Chem. Lett. 11 (2020) 5913–5919. doi: 10.1021/acs.jpclett.0c01394

    20. [20]

      X. Zhuang, J. Wang, Z. Lan, J. Phys. Chem. B 117 (2013) 15976–15986. doi: 10.1021/jp408799b

    21. [21]

      A.T. Ulijasz, G. Cornilescu, C.C. Cornilescu, et al., Nature 463 (2010) 250–254. doi: 10.1038/nature08671

    22. [22]

      F.V. Escobar, P. Piwowarski, J. Salewski, et al., Nat. Chem. 7 (2015) 423–430. doi: 10.1038/nchem.2225

    23. [23]

      B. Durbeej, O.A. Borg, L.A. Eriksson, Chem. Phys. Lett. 416 (2005) 83–88. doi: 10.1016/j.cplett.2005.09.050

    24. [24]

      V. Modi, S. Donnini, G. Groenhof, D. Morozov, J. Phys. Chem. B 123 (2019) 2325–2334. doi: 10.1021/acs.jpcb.9b01117

    25. [25]

      D. Morozov, V. Modi, V. Mironov, G. Groenhof, J. Phys. Chem. Lett. 13 (2022) 4538–4542. doi: 10.1021/acs.jpclett.2c00899

    26. [26]

      O. Falklöf, B. Durbeej, ChemPhysChem 17 (2016) 954–957. doi: 10.1002/cphc.201501080

    27. [27]

      V. Macaluso, G. Salvadori, L. Cupellini, B. Mennucci, Chem. Sci. 12 (2021) 5555–5565. doi: 10.1039/d1sc00186h

    28. [28]

      J. Hahn, H.M. Strauss, P. Schmieder, J. Am. Chem. Soc. 130 (2008) 11170–11178. doi: 10.1021/ja8031086

    29. [29]

      Y. Yang, M. Linke, T. Von Haimberger, et al., J. Am. Chem. Soc. 134 (2012) 1408–1411. doi: 10.1021/ja209413d

    30. [30]

      A. Strambi, B. Durbeej, Photochem. Photobiol. Sci. 10 (2011) 569–579. doi: 10.1039/c0pp00307g

    31. [31]

      B. Durbeej, Phys. Chem. Chem. Phys. 11 (2009) 1354–1361. doi: 10.1039/b811813b

    32. [32]

      Y. Liu, S. Han, J. Wang, et al., J. Am. Chem. Soc. 143 (2021) 2130–2137. doi: 10.1021/jacs.0c12513

    33. [33]

      H. Xu, W. Guo, J. Wang, et al., J. Am. Chem. Soc. 143 (2021) 14379–14385. doi: 10.1021/jacs.1c07521

    34. [34]

      W. Guo, X. Liu, S. Han, et al., Angew. Chem. 59 (2020) 13879–13884. doi: 10.1002/anie.202004235

    35. [35]

      D. Herschlag, M.M. Pinney, Biochemistry 57 (2018) 3338–3352. doi: 10.1021/acs.biochem.8b00217

    36. [36]

      M.M. Pinney, A. Natarajan, F. Yabukarski, et al., J. Am. Chem. Soc. 140 (2018) 9827–9843. doi: 10.1021/jacs.8b01596

    37. [37]

      K. Toh, E.A. Stojković, I.H. van Stokkum, K. Moffat, J.T. Kennis, Proc. Natl. Acad. Sci. U. S. A. 107 (2010) 9170–9175. doi: 10.1073/pnas.0911535107

    38. [38]

      S. Bhattacharya, M.E. Auldridge, H. Lehtivuori, J.A. Ihalainen, K.T. Forest, J. Bio. Chem. 289 (2014) 32144–32152. doi: 10.1074/jbc.M114.589739

    39. [39]

      J. Zhu, D.M. Shcherbakova, Y. Hontani, V.V. Verkhusha, J. Kennis, Sci. Rep. 5 (2015) 1–8. doi: 10.1007/s11468-014-9769-1

    40. [40]

      S. Osoegawa, R. Miyoshi, K. Watanabe, et al., J. Phys. Chem. B 123 (2019) 3242–3247. doi: 10.1021/acs.jpcb.9b00965

    41. [41]

      O. Falklöf, B. Durbeej, ChemPhotoChem 2 (2018) 453–457. doi: 10.1002/cptc.201800014

    42. [42]

      R. Crespo-Otero, M. Barbatti, Chem. Rev. 118 (2018) 7026–7068. doi: 10.1021/acs.chemrev.7b00577

    43. [43]

      L. Du, Z. Lan, J. Chem. Theor. Comput. 11 (2015) 1360–1374. doi: 10.1021/ct501106d

    44. [44]

      M. Barbatti, G. Granucci, M. Persico, et al., J. Photochem. Photobiol. A: Chem. 190 (2007) 228–240. doi: 10.1016/j.jphotochem.2006.12.008

    45. [45]

      S. Mai, L. González, Angew. Chem. 59 (2020) 16832–16846. doi: 10.1002/anie.201916381

    46. [46]

      G. Cui, W. Thiel, J. Chem. Phys. 141 (2014) 124101. doi: 10.1063/1.4894849

    47. [47]

      B.F. Curchod, T.J. Martínez, Chem. Rev. 118 (2018) 3305–3336. doi: 10.1021/acs.chemrev.7b00423

    48. [48]

      J.C. Tully, J. Chem. Phys. 93 (1990) 1061–1071. doi: 10.1063/1.459170

    49. [49]

      A. Koslowski, M.E. Beck, W. Thiel, J. Comput. Chem. 24 (2003) 714–726. doi: 10.1002/jcc.10210

    50. [50]

      W. Weber, W. Thiel, Theor. Chem. Acc. 103 (2000) 495–506. doi: 10.1007/s002149900083

    51. [51]

      A. Kazaryan, Z. Lan, L.V. Schafer, W. Thiel, M. Filatov, J. Chem. Theor. Comput. 7 (2011) 2189–2199. doi: 10.1021/ct200199w

    52. [52]

      T.S. Zhang, Z.W. Li, Q. Fang, et al., J. Phys. Chem. A 123 (2019) 6144–6151. doi: 10.1021/acs.jpca.9b04372

    53. [53]

      G. Cui, W. Thiel, Angew. Chem. 125 (2013) 451–454. doi: 10.1002/ange.201207628

    54. [54]

      G. Cui, Z. Lan, W. Thiel, J. Am. Chem. Soc. 134 (2012) 1662–1672. doi: 10.1021/ja208496s

    55. [55]

      X. Pang, X. Cui, D. Hu, et al., J. Phys. Chem. A 121 (2017) 1240–1249. doi: 10.1021/acs.jpca.6b12253

    56. [56]

      Z. Lan, E. Fabiano, W. Thiel, J. Phys. Chem. B 113 (2009) 3548–3555. doi: 10.1021/jp809085h

    57. [57]

      J.R. Wagner, J.S. Brunzelle, K.T. Forest, R.D. Vierstra, Nature 438 (2005) 325–331. doi: 10.1038/nature04118

    58. [58]

      X. Yang, E.A. Stojković, J. Kuk, K. Moffat, Proc. Natl. Acad. Sci. U. S. A. 104 (2007) 12571–12576. doi: 10.1073/pnas.0701737104

    59. [59]

      J.R. Wagner, J. Zhang, J.S. Brunzelle, R.D. Vierstra, K.T. Forest, J. Bio. Chem. 282 (2007) 12298–12309. doi: 10.1074/jbc.M611824200

    60. [60]

      T. Lu, Q. Chen, Chem. Methods 1 (2021) 231–239. doi: 10.1002/cmtd.202100007

    61. [61]

      A. Cembran, F. Bernardi, M. Olivucci, M. Garavelli, J. Am. Chem. Soc. 126 (2004) 16018–16037. doi: 10.1021/ja048782+

  • Scheme 1  Chromophores, photoreaction centers in phytochromes. (1) R1: =CH-CH2-Cys, R2: -CH=CH2 for BV; (2) R1: -CH-(CH3)Cys, R2: -CH2-CH3 for PCB; (3) R1: -CH-(CH3)Cys, R2: -CH=CH2 for PΦB.

    Figure 1  OM2 optimized S0 minimum-energy structures of four most probable protonated BV-chromophore models. Selected bond lengths are given in Angstroms. See the Supporting information for Cartesian coordinates.

    Figure 2  Time-dependent fractional occupations of the S0 and S1 electronic states of the four BV models (BV_N, BV_C, BV_B and BV_A) in the non-adiabatic dynamics staring from the S1 state.

    Figure 3  (a, b) OM2/MRCI optimized two S0/S1 conical intersections of the BV_N model. Cyan, red and gold represent S0_min, CI1 and CI2, respectively. The CI1 is characterized by the rotation of the whole C-D double-rings moiety around the C10C11 bond, while the CI2 is featured by the rotation of the whole A-B double-rings moiety around the C9C10 bond. (c) Distributions of two key dihedral angles at the initial and hopping geometries of the BV_N model. (d) Potential energy curves along the linear interpolated pathway from S0_min to the two conical intersections CI1 and CI2 of the BV_N model.

    Figure 4  Representative trajectories towards to various reaction pathways. The cyan and white structural skeleton represent the snapshot at the current and starting time steps, respectively. Top: the trajectory back to initial BV_N structure; Middle: the trajectory ends with the rotation of the Ring D; Bottom: the trajectory evolution shows the mixed the Ring A and the Ring D rotations.

    Figure 5  (a, b) Distributions of key dihedral angles at the initial and hopping points of the BV_C model. (c) Potential energy curves along the linear interpolated pathway from S0_min to the S1_min and from S1_min to the CI for the BV_C model. (d) OM2/MRCI optimized S0/S1 CI of the BV_C model. Cyan, gold and red represent S0_min, S1_min and CI, respectively.

  • 加载中
计量
  • PDF下载量:  6
  • 文章访问数:  1118
  • HTML全文浏览量:  50
文章相关
  • 发布日期:  2023-07-15
  • 收稿日期:  2022-06-08
  • 接受日期:  2022-09-22
  • 修回日期:  2022-06-27
  • 网络出版日期:  2022-09-25
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

/

返回文章